MAA 2007) Dedicated to the Centenary of Mark Krein Odessa, Ukraine, April 9–14, 2007

Total Page:16

File Type:pdf, Size:1020Kb

MAA 2007) Dedicated to the Centenary of Mark Krein Odessa, Ukraine, April 9–14, 2007 International Conference MODERN ANALYSIS AND APPLICATIONS (MAA 2007) dedicated to the centenary of Mark Krein Odessa, Ukraine, April 9–14, 2007 Book of abstracts Kyiv - 2007 International Conference MODERN ANALYSIS AND APPLICATIONS dedicated to the centenary of Mark Krein ORGANIZERS: Ministry of Education and Science of Ukraine • National Academy of Sciences of Ukraine • Municipality of Odessa • Institute of Mathematics, National Academy of Sciences of Ukraine • Odessa National I.I. Mechnikov University • Odessa Nationa Maritime University • South Ukrainian State Pedagogical University “K.D. Ushynsky” • Institute of Mathematics, Economics and Mechanics of Odessa National • I.I. Mechnikov University Department of Mathematical Physics of Odessa National I.I. Mechnikov • University Department of Theoretical Physics of Odessa National I.I. Mechnikov University • SPONSORS: Swedish Institute • European Science Foundation • Municipality of Odessa • Manufactured Goods Market Limited • JV Dipolos Limited • Imexbank (Joint Stock commercial bank) • Pivdennyi Bank • Port-Franco Bank • c 2007 Institute of Mathematics, National Acad. Sci. of Ukraine ABSTRACTS OF TALKS Rustam Abdullaev, Bahrom Hikmetov AN ANALOGUE OF ARENS ALGEBRA CONSTRUCTED BY A STATE ON VON NEUMANN ALGEBRA Let M be a von Neumann algebra with faithful normal finite trace µ. For any faithful normal state ϕ on M exists positive operator h affiliated to M such that ϕ(x)= µ(hx) for every x M + ∈ (see[1]). On M we introduce two norms 1 x µ = (µ( x p)) p , k kp | | 1 1 1 x ϕ = (µ( h 2p xh 2p )p) p . k kp | | µ The completion of M by the norm p denote by Lp(M,µ), the completion of M by the norm ϕ k · k p denote by Lp(M, ϕ). It is easy to see, that Lp(M,µ) is a partical case of Lp(M, ϕ). For k · k 1 regular normal states (i.e when h− is local measurable ) space Lp(M, ϕ) can be described by local measurable operators affiliated to M (see[2]). Next we consider only regular states. The set Lp(M, ϕ)-is called Arens space. Consider on it the topology t, generated by p [1; ) ∈ ∞ ϕ ω system norm T p p [1; ) (see[3]). Denote it by (L (M, ϕ),t). {k · k } ∈ ∞ ω Theorem 1. The space (L (M, ϕ),tϕ) is complete metricable local convex space with dual space, which is isometrical isomorphic to Lq(M, ϕ). q (1; ] ∈S∞ Theorem 2. The space Lω(M,µ) consides with the space Lω(M, ϕ) if and only if 1 h Lq(M,µ) and h− Lq(M, ϕ). ∈ q (1; ] ∈ q (1; ] ∈S∞ ∈S∞ Theorem 3. On noncommutative algebra von Neumann exists a normal state ψ, such that Lω(M, ψ) = Lω(M,ν) for any faithful normal finite trace ν. 6 [1] Pedersen G.K., Takesaki M. The Radon Nikodym theorem for von Neumann algebras, Acta Math, (1973), v.130, no. 1-2, 53-87 [2] Trunov H.B K teorii normalnix vesov na algebrax fon Neymana.,Izv.Vuzov. Matem, (1982),no. 9, 61-70 [3] Abdullaev R.Z. Izomorfizmi algebr Arensa., Sibirskiy jurnal industrialnoy matematiki, (1988), v.1, no. 2, 3-13 Vadym Adamyan, Boris Pavlov KREIN FORMULA IN TRANSPORT PROBLEMS FOR QUANTUM NETWORKS Quantum switches, spin filters and other types of quantum nano- electronic devices are quantum networks. Mathematical design of quantum networks with prescribed transport properties is a challenging problem of modern physics. Contemporary mathematics is already able to contribute to the solution of the simplest type of transport problems- to the one-body scattering problem for networks constructed on a surface of a semiconductor of quantum dots and straight quantum wire connecting them to each other or to infinity. Construction of the network with prescribed 4 transport properties - with given scattering matrix - is a subject of the inverse scattering problem for quantum networks, which is still beyond the limits of modern mathematics. Optimization of design of the network based on straightforward computing is expensive and non-efficient because of large number of geometrical and physical parameters defining the properties of the network. One can expect that substitution of the network by the corresponding quantum graph, possibly with resonances nodes supplied with an inner structure, looks like an attractive alternative to the straightforward computing. This model is solvable, but it requires appropriate fitting. We suggest a procedure of fitting of this model based on semi-analytic approach, involving direct computing for eigenfunctions of discrete spectrum on quantum dots and subsequent rational approximation of the Dirichlet-to-Neumann maps and asymptotic calculation of the relevant scattering matrix. Attachment of the semi-infinite quantum wires to quantum dots defines the transformation of bound states in quantum dots into resonances. This is a typical difficulty for perturbation problems on the continuous spectrum. We overcome this difficulty via introduction of an inter- mediate Hamiltonian which is obtained by Glazman splitting of the Schr¨odinger operator on the quantum network into orthogonal sum of two operators with complementary branches of the absolutely continuous spectra, constituted by the open and closed branches of the continuous spectrum, for given Fermi level. Based on Dirichlet-to-Neumann map of the intermediate Hamil- tonian, we eliminate calculation of coefficients in front of exponentially decreasing modes in the scattered waves, reducing the problem for scattered waves to solution of a finite linear system. The rational approximation of the Dirichlet-to-Neumann map of the Intermediate Hamiltonian is interpreted as Weyl-Titchmarsh function of the corresponding solvable model on the relevant quantum graph which is automatically fitted such that the model scattering matrix serves an approximation of the scattering matrix of the full Hamiltonian on an essential spectral interval near the Fermi level. On the other hand, the scattering matrix of the model is represented by the modified Krein formula, which gives an approximation of the full scattering matrix. The solvable model may serve an intermediate step in the analytic perturbation procedure on the continuous spectrum, because spectral characteristics of the full Hamiltonian are connected to the corresponding spectral characteristics of the solvable model on the essential spectral interval by the convergent standard analytic perturbation procedure. Mikhail Agranovich SPECTRAL PROBLEMS FOR STRONGLY ELLIPTIC SYSTEMS IN LIPSCHITZ DOMAINS We will discuss the regularity of solutions to boundary value problems for strongly elliptic partial differential systems in Lipschitz domains. The results are applied to regularity properties of eigenfunctions of spectral boundary value problems. Some generalizations of the author’s results published in Functional Analysis and Its Applications, 40, 2006, No. 4, will be explained. Sergey Aizikovich BILATERAL ASYMPTOTIC SOLUTION OF ONE CLASS OF DUAL INTEGRAL EQUATIONS OF THE STATIC CONTACT PROBLEMS FOR INHOMOGENEOUS WITH DEPTH FOUNDATION 5 Sergio Albeverio, Rostyslav O. Hryniv, Yaroslav V. Mykytyuk THE KREIN EQUATION AND INVERSE SPECTRAL PROBLEMS In early 1950-ies, M.Krein suggested a method of solving the inverse spectral problems for Sturm– Liouville operators using special integral equations. We show how this method can further be developed to provide solutions to inverse spectral problems for a wide class of singular Sturm– Liouville operators and of Dirac operators with summable potentials. Daniel Alpay, I. Gohberg CONTINUOUS AND DISCRETE SYSTEMS WITH RATIONAL SPECTRAL DATA We first review the theory of canonical differential expressions in the rational case, and in particular discuss the class of potentials introduced by the authors in [1]. Then, we define and study the discrete analogue of canonical differential expressions. We focus on the rational case. Two kinds of discrete systems are to be distinguished: one–sided and two–sided. In both cases the analogue of the potential is a sequence of numbers in the open unit disk (Schur coefficients). We define the characteristic spectral functions of the discrete systems and provide exact realization formulas for them when the Schur coefficients are of a special form called strictly pseudo–exponential and introduced in [2]. The matrix-valued case, to appear in [3] will be briefly outlined. [1] D. Alpay and I. Gohberg. Inverse spectral problem for differential operators with rational scattering matrix functions. Journal of differential equations, 118:1–19, 1995. [2] D. Alpay and I. Gohberg. Connections between the Carath´eodory-Toeplitz and the Nehari extension problems: the discrete scalar case. Integral Equations Operator Theory, 37(2):125– 142, 2000. [3] D. Alpay and I. Gohberg. Discrete systems and their characteristic spectral functions. Mediterranean Journal of Mathematics, 2007 (in press). Arthur Amirshadyan BOUNDARY INDEFINITE NEVANLINNA-PICK INTERPOLATION PROBLEM The Problem ∂IPκ. Given are: real points zj and symmetric matrices Wj, Dj (j=1,..,m). Find the matrix function F (λ) N (dCn), satisfying the conditions ∈ κ lim F (λ)= W (j = 1, .., m) (1) ∧ j λ zj −→ lim F ′(λ) D (j = 1, .., m). (2) ∧ j λ zj ≤ −→ It is assumed that there exist nontangential limits in (1), (2). 6 The condition for the solvablity of this problem will be stated in terms of the Pick matrix m P = (Pjk)j,k=1, of the form [7], [1]: Wj W − k , ifj = k; zj z Pjk = − k 6 ( Dj, ifj = k. In what follows, we assume that det P = 0. 6 The matrices V = (In, .., In), W = (W1, .., Wm), Z = diag(z1In, .., zmIn) are called the data of the Problem ∂IP . Let Φ := (Z λ)P + V W . Now we define the the solution matrix κ λ − ∗ 2 W 1 1 Ω(λ)=(Ω ) =: I + (Z λ)− P− ( V ∗,W ∗). ij i,j=1 2n V − − Theorem 1. Suppose that det P = 0, sq P = κ0 κ, detΦzj = 0 (j=1,..,m). Then the 6 − ≤ 6 formula 1 F (λ)=(Ω (λ)ψ(λ) Ω (λ)φ(λ))(Ω (λ)ψ(λ) Ω (λ)φ(λ))− (3) 12 − 11 22 − 21 establishes a one-to-one correspondence between the solution set of Problem ∂IPκ and the set of Nκ κ0 – pairs φ(λ), ψ(λ) , for which − { } det(Ω (λ)ψ(λ) Ω (λ)φ(λ)) 0 (4) 22 − 21 6≡ 1 1 lim (λ z )φ(λ)(V Φ− V ∗φ(λ)+ ψ(λ))− 0 (j = 1, ..., m).
Recommended publications
  • Tilburg University on the Matrix
    Tilburg University On the Matrix (I + X)-1 Engwerda, J.C. Publication date: 2005 Link to publication in Tilburg University Research Portal Citation for published version (APA): Engwerda, J. C. (2005). On the Matrix (I + X)-1. (CentER Discussion Paper; Vol. 2005-120). Macroeconomics. General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. Download date: 25. sep. 2021 No. 2005–120 -1 ON THE MATRIX INEQUALITY ( I + x ) I . By Jacob Engwerda November 2005 ISSN 0924-7815 On the matrix inequality (I + X)−1 ≤ I. Jacob Engwerda Tilburg University Dept. of Econometrics and O.R. P.O. Box: 90153, 5000 LE Tilburg, The Netherlands e-mail: [email protected] November, 2005 1 Abstract: In this note we consider the question under which conditions all entries of the matrix I −(I +X)−1 are nonnegative in case matrix X is a real positive definite matrix.
    [Show full text]
  • Bart, Hempfling, Kaashoek. (Eds.) Israel Gohberg and Friends.. on The
    Israel Gohberg and Friends On the Occasion of his 80th Birthday Harm Bart Thomas Hempfling Marinus A. Kaashoek Editors Birkhäuser Basel · Boston · Berlin Editors: Harm Bart Marinus A. Kaashoek Econometrisch Instituut Department of Mathematics, FEW Erasmus Universiteit Rotterdam Vrije Universiteit Postbus 1738 De Boelelaan 1081A 3000 DR Rotterdam 1081 HV Amsterdam The Netherlands The Netherlands e-mail: [email protected] e-mail: [email protected] Thomas Hempfling Editorial Department Mathematics Birkhäuser Publishing Ltd. P.O. Box 133 4010 Basel Switzerland e-mail: thomas.hempfl[email protected] Library of Congress Control Number: 2008927170 Bibliographic information published by Die Deutsche Bibliothek Die Deutsche Bibliothek lists this publication in the Deutsche Nationalbibliografie; detailed bibliographic data is available in the Internet at <http://dnb.ddb.de>. ISBN 978-3-7643-8733-4 Birkhäuser Verlag, Basel – Boston – Berlin This work is subject to copyright. All rights are reserved, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, re-use of illustrations, recitation, broadcasting, reproduction on microfilms or in other ways, and storage in data banks. For any kind of use permission of the copyright owner must be obtained. © 2008 Birkhäuser Verlag AG Basel · Boston · Berlin P.O. Box 133, CH-4010 Basel, Switzerland Part of Springer Science+Business Media Printed on acid-free paper produced of chlorine-free pulp. TCF ∞ Printed in Germany ISBN 978-3-7643-8733-4 e-ISBN 978-3-7643-8734-1 9 8 7 6 5 4 3 2 1 www.birkhauser.ch Contents Preface.......................................................................ix CongratulationsfromthePublisher...........................................xii PartI.MathematicalandPhilosophical-MathematicalTales...................1 I.
    [Show full text]
  • Arxiv:1909.13402V1 [Math.CA] 30 Sep 2019 Routh-Hurwitz Array [14], Argument Principle [23] and So On
    ON GENERALIZATION OF CLASSICAL HURWITZ STABILITY CRITERIA FOR MATRIX POLYNOMIALS XUZHOU ZHAN AND ALEXANDER DYACHENKO Abstract. In this paper, we associate a class of Hurwitz matrix polynomi- als with Stieltjes positive definite matrix sequences. This connection leads to an extension of two classical criteria of Hurwitz stability for real polynomials to matrix polynomials: tests for Hurwitz stability via positive definiteness of block-Hankel matrices built from matricial Markov parameters and via matricial Stieltjes continued fractions. We obtain further conditions for Hurwitz stability in terms of block-Hankel minors and quasiminors, which may be viewed as a weak version of the total positivity criterion. Keywords: Hurwitz stability, matrix polynomials, total positivity, Markov parameters, Hankel matrices, Stieltjes positive definite sequences, quasiminors 1. Introduction Consider a high-order differential system (n) (n−1) A0y (t) + A1y (t) + ··· + Any(t) = u(t); where A0;:::;An are complex matrices, y(t) is the output vector and u(t) denotes the control input vector. The asymptotic stability of such a system is determined by the Hurwitz stability of its characteristic matrix polynomial n n−1 F (z) = A0z + A1z + ··· + An; or to say, by that all roots of det F (z) lie in the open left half-plane <z < 0. Many algebraic techniques are developed for testing the Hurwitz stability of matrix polynomials, which allow to avoid computing the determinant and zeros: LMI approach [20, 21, 27, 28], the Anderson-Jury Bezoutian [29, 30], matrix Cauchy indices [6], lossless positive real property [4], block Hurwitz matrix [25], extended arXiv:1909.13402v1 [math.CA] 30 Sep 2019 Routh-Hurwitz array [14], argument principle [23] and so on.
    [Show full text]
  • Implications for Understanding the Role of Affect in Mathematical Thinking
    Mathematicians and music: Implications for understanding the role of affect in mathematical thinking Rena E. Gelb Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy under the Executive Committee of the Graduate School of Arts and Sciences COLUMBIA UNIVERSITY 2021 © 2021 Rena E. Gelb All Rights Reserved Abstract Mathematicians and music: Implications for understanding the role of affect in mathematical thinking Rena E. Gelb The study examines the role of music in the lives and work of 20th century mathematicians within the framework of understanding the contribution of affect to mathematical thinking. The current study focuses on understanding affect and mathematical identity in the contexts of the personal, familial, communal and artistic domains, with a particular focus on musical communities. The study draws on published and archival documents and uses a multiple case study approach in analyzing six mathematicians. The study applies the constant comparative method to identify common themes across cases. The study finds that the ways the subjects are involved in music is personal, familial, communal and social, connecting them to communities of other mathematicians. The results further show that the subjects connect their involvement in music with their mathematical practices through 1) characterizing the mathematician as an artist and mathematics as an art, in particular the art of music; 2) prioritizing aesthetic criteria in their practices of mathematics; and 3) comparing themselves and other mathematicians to musicians. The results show that there is a close connection between subjects’ mathematical and musical identities. I identify eight affective elements that mathematicians display in their work in mathematics, and propose an organization of these affective elements around a view of mathematics as an art, with a particular focus on the art of music.
    [Show full text]
  • Fundamental Theorems in Mathematics
    SOME FUNDAMENTAL THEOREMS IN MATHEMATICS OLIVER KNILL Abstract. An expository hitchhikers guide to some theorems in mathematics. Criteria for the current list of 243 theorems are whether the result can be formulated elegantly, whether it is beautiful or useful and whether it could serve as a guide [6] without leading to panic. The order is not a ranking but ordered along a time-line when things were writ- ten down. Since [556] stated “a mathematical theorem only becomes beautiful if presented as a crown jewel within a context" we try sometimes to give some context. Of course, any such list of theorems is a matter of personal preferences, taste and limitations. The num- ber of theorems is arbitrary, the initial obvious goal was 42 but that number got eventually surpassed as it is hard to stop, once started. As a compensation, there are 42 “tweetable" theorems with included proofs. More comments on the choice of the theorems is included in an epilogue. For literature on general mathematics, see [193, 189, 29, 235, 254, 619, 412, 138], for history [217, 625, 376, 73, 46, 208, 379, 365, 690, 113, 618, 79, 259, 341], for popular, beautiful or elegant things [12, 529, 201, 182, 17, 672, 673, 44, 204, 190, 245, 446, 616, 303, 201, 2, 127, 146, 128, 502, 261, 172]. For comprehensive overviews in large parts of math- ematics, [74, 165, 166, 51, 593] or predictions on developments [47]. For reflections about mathematics in general [145, 455, 45, 306, 439, 99, 561]. Encyclopedic source examples are [188, 705, 670, 102, 192, 152, 221, 191, 111, 635].
    [Show full text]
  • A Note on Nonnegative Diagonally Dominant Matrices Geir Dahl
    UNIVERSITY OF OSLO Department of Informatics A note on nonnegative diagonally dominant matrices Geir Dahl Report 269, ISBN 82-7368-211-0 April 1999 A note on nonnegative diagonally dominant matrices ∗ Geir Dahl April 1999 ∗ e make some observations concerning the set C of real nonnegative, W n diagonally dominant matrices of order . This set is a symmetric and n convex cone and we determine its extreme rays. From this we derive ∗ dierent results, e.g., that the rank and the kernel of each matrix A ∈Cn is , and may b e found explicitly. y a certain supp ort graph of determined b A ∗ ver, the set of doubly sto chastic matrices in C is studied. Moreo n Keywords: Diagonal ly dominant matrices, convex cones, graphs and ma- trices. 1 An observation e recall that a real matrix of order is called diagonal ly dominant if W P A n | |≥ | | for . If all these inequalities are strict, is ai,i j=6 i ai,j i =1,...,n A strictly diagonal ly dominant. These matrices arise in many applications as e.g., discretization of partial dierential equations [14] and cubic spline interp ola- [10], and a typical problem is to solve a linear system where tion Ax = b strictly diagonally dominant, see also [13]. Strict diagonal dominance A is is a criterion which is easy to check for nonsingularity, and this is imp ortant for the estimation of eigenvalues confer Ger²chgorin disks, see e.g. [7]. For more ab out diagonally dominant matrices, see [7] or [13]. A matrix is called nonnegative positive if all its elements are nonnegative p ositive.
    [Show full text]
  • Conditioning Analysis of Positive Definite Matrices by Approximate Factorizations
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector Journal of Computational and Applied Mathematics 26 (1989) 257-269 257 North-Holland Conditioning analysis of positive definite matrices by approximate factorizations Robert BEAUWENS Service de M&rologie NucGaire, UniversitP Libre de Bruxelles, 50 au. F.D. Roosevelt, 1050 Brussels, Belgium Renaud WILMET The Johns Hopkins University, Baltimore, MD 21218, U.S.A Received 19 February 1988 Revised 7 December 1988 Abstract: The conditioning analysis of positive definite matrices by approximate LU factorizations is usually reduced to that of Stieltjes matrices (or even to more specific classes of matrices) by means of perturbation arguments like spectral equivalence. We show in the present work that a wider class, which we call “almost Stieltjes” matrices, can be used as reference class and that it has decisive advantages for the conditioning analysis of finite element approxima- tions of large multidimensional steady-state diffusion problems. Keywords: Iterative analysis, preconditioning, approximate factorizations, solution of finite element approximations to partial differential equations. 1. Introduction The conditioning analysis of positive definite matrices by approximate factorizations has been the subject of a geometrical approach by Axelsson and his coworkers [1,2,7,8] and of an algebraic approach by one of the authors [3-51. In both approaches, the main role is played by Stieltjes matrices to which other cases are reduced by means of spectral equivalence or by a similar perturbation argument. In applications to discrete partial differential equations, such arguments lead to asymptotic estimates which may be quite poor in nonasymptotic regions, a typical instance in large engineering multidimensional steady-state diffusion applications, as illustrated below.
    [Show full text]
  • Similarities on Graphs: Kernels Versus Proximity Measures1
    Similarities on Graphs: Kernels versus Proximity Measures1 Konstantin Avrachenkova, Pavel Chebotarevb,∗, Dmytro Rubanova aInria Sophia Antipolis, France bTrapeznikov Institute of Control Sciences of the Russian Academy of Sciences, Russia Abstract We analytically study proximity and distance properties of various kernels and similarity measures on graphs. This helps to understand the mathemat- ical nature of such measures and can potentially be useful for recommending the adoption of specific similarity measures in data analysis. 1. Introduction Until the 1960s, mathematicians studied only one distance for graph ver- tices, the shortest path distance [6]. In 1967, Gerald Sharpe proposed the electric distance [42]; then it was rediscovered several times. For some col- lection of graph distances, we refer to [19, Chapter 15]. Distances2 are treated as dissimilarity measures. In contrast, similarity measures are maximized when every distance is equal to zero, i.e., when two arguments of the function coincide. At the same time, there is a close relation between distances and certain classes of similarity measures. One of such classes consists of functions defined with the help of kernels on graphs, i.e., positive semidefinite matrices with indices corresponding to the nodes. Every kernel is the Gram matrix of some set of vectors in a Euclidean space, and Schoenberg’s theorem [39, 40] shows how it can be transformed into the matrix of Euclidean distances between these vectors. arXiv:1802.06284v1 [math.CO] 17 Feb 2018 Another class consists of proximity measures characterized by the trian- gle inequality for proximities. Every proximity measure κ(x, y) generates a ∗Corresponding author Email addresses: [email protected] (Konstantin Avrachenkov), [email protected] (Pavel Chebotarev), [email protected] (Dmytro Rubanov) 1This is the authors’ version of a work that was submitted to Elsevier.
    [Show full text]
  • Facts from Linear Algebra
    Appendix A Facts from Linear Algebra Abstract We introduce the notation of vector and matrices (cf. Section A.1), and recall the solvability of linear systems (cf. Section A.2). Section A.3 introduces the spectrum σ(A), matrix polynomials P (A) and their spectra, the spectral radius ρ(A), and its properties. Block structures are introduced in Section A.4. Subjects of Section A.5 are orthogonal and orthonormal vectors, orthogonalisation, the QR method, and orthogonal projections. Section A.6 is devoted to the Schur normal form (§A.6.1) and the Jordan normal form (§A.6.2). Diagonalisability is discussed in §A.6.3. Finally, in §A.6.4, the singular value decomposition is explained. A.1 Notation for Vectors and Matrices We recall that the field K denotes either R or C. Given a finite index set I, the linear I space of all vectors x =(xi)i∈I with xi ∈ K is denoted by K . The corresponding square matrices form the space KI×I . KI×J with another index set J describes rectangular matrices mapping KJ into KI . The linear subspace of a vector space V spanned by the vectors {xα ∈V : α ∈ I} is denoted and defined by α α span{x : α ∈ I} := aαx : aα ∈ K . α∈I I×I T Let A =(aαβ)α,β∈I ∈ K . Then A =(aβα)α,β∈I denotes the transposed H matrix, while A =(aβα)α,β∈I is the adjoint (or Hermitian transposed) matrix. T H Note that A = A holds if K = R . Since (x1,x2,...) indicates a row vector, T (x1,x2,...) is used for a column vector.
    [Show full text]
  • Inverse M-Matrix, a New Characterization
    Linear Algebra and its Applications 595 (2020) 182–191 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Inverse M-matrix, a new characterization Claude Dellacherie a, Servet Martínez b, Jaime San Martín b,∗ a Laboratoire Raphaël Salem, UMR 6085, Université de Rouen, Site du Madrillet, 76801 Saint Étienne du Rouvray Cedex, France b CMM-DIM, Universidad de Chile, UMI-CNRS 2807, Casilla 170-3 Correo 3, Santiago, Chile a r t i c l e i n f o a b s t r a c t Article history: In this article we present a new characterization of inverse M- Received 18 June 2019 matrices, inverse row diagonally dominant M-matrices and Accepted 20 February 2020 inverse row and column diagonally dominant M-matrices, Available online 22 February 2020 based on the positivity of certain inner products. Submitted by P. Semrl © 2020 Elsevier Inc. All rights reserved. MSC: 15B35 15B51 60J10 Keywords: M-matrix Inverse M-matrix Potentials Complete Maximum Principle Markov chains Contents 1. Introduction and main results.......................................... 183 2. Proof of Theorem 1.1 and Theorem 1.3 .................................... 186 3. Some complements.................................................. 189 * Corresponding author. E-mail addresses: [email protected] (C. Dellacherie), [email protected] (S. Martínez), [email protected] (J. San Martín). https://doi.org/10.1016/j.laa.2020.02.024 0024-3795/© 2020 Elsevier Inc. All rights reserved. C. Dellacherie et al. / Linear Algebra and its Applications 595 (2020) 182–191 183 Acknowledgement....................................................... 190 References............................................................ 190 1. Introduction and main results In this short note, we give a new characterization of inverses M-matrices, inverses of row diagonally dominant M-matrices and inverses of row and column diagonally domi- nant M-matrices.
    [Show full text]
  • Localization in Matrix Computations: Theory and Applications
    Localization in Matrix Computations: Theory and Applications Michele Benzi Department of Mathematics and Computer Science, Emory University, Atlanta, GA 30322, USA. Email: [email protected] Summary. Many important problems in mathematics and physics lead to (non- sparse) functions, vectors, or matrices in which the fraction of nonnegligible entries is vanishingly small compared the total number of entries as the size of the system tends to infinity. In other words, the nonnegligible entries tend to be localized, or concentrated, around a small region within the computational domain, with rapid decay away from this region (uniformly as the system size grows). When present, localization opens up the possibility of developing fast approximation algorithms, the complexity of which scales linearly in the size of the problem. While localization already plays an important role in various areas of quantum physics and chemistry, it has received until recently relatively little attention by researchers in numerical linear algebra. In this chapter we survey localization phenomena arising in various fields, and we provide unified theoretical explanations for such phenomena using general results on the decay behavior of matrix functions. We also discuss compu- tational implications for a range of applications. 1 Introduction In numerical linear algebra, it is common to distinguish between sparse and dense matrix computations. An n ˆ n sparse matrix A is one in which the number of nonzero entries is much smaller than n2 for n large. It is generally understood that a matrix is dense if it is not sparse.1 These are not, of course, formal definitions. A more precise definition of a sparse n ˆ n matrix, used by some authors, requires that the number of nonzeros in A is Opnq as n Ñ 8.
    [Show full text]
  • F(X, Y) C[(A, B) X
    SIAM J. ScI. STAT. COMPUT. 1983 Society for Industrial and Applied Mathematics Vol. 4, No. 2, June 1983 0196-5204/83/0402--0008 $1.25/0 A HYBRID ASYMPTOTIC-FINITE ELEMENT METHOD FOR STIFF TWO-POINT BOUNDARY VALUE PROBLEMS* R. C. Y. CHINe AND R. KRASNY' Abstract. An accurate and efficient numerical method has been developed for a nonlinear stiff second order two-point boundary value problem. The scheme combines asymptotic methods with the usual solution techniques for two-point boundary value problems. A new modification of Newton's method or quasilinear- ization is used to reduce the nonlinear problem to a sequence of linear problems. The resultant linear problem is solved by patching local solutions at the knots or equivalently by projecting onto an affine subset constructed from asymptotic expansions. In this way, boundary layers are naturally incorporated into the approximation. An adaptive mesh is employed to achieve an error of O(1/N2)+O(/e). Here, N is the number of intervals and e << is the singular perturbation parameter. Numerical computations are presented. Key words, projection method, quasilinearization, asymptotic expansions, boundary layers, singular perturbations, two-point boundary value problems, y-elliptic splines 1. Introduction. We treat the following stiff boundary value problem: ey" =f(x, y), a <x <b, f(x, y) C[(a, b) x R], f(x, y)>0, O<e << 1, with boundary conditions: (1.2) aoy(a)-aly'(a)=a, (1.3) boy(b + bly'(b fl. Under the assumptions that aoa >0, bobl >0 and la01/lb01 >0 Keller [26] and Bernfield and Lakshmikantham [8] have shown that the boundary value problem (BVP) has a unique solution.
    [Show full text]