A Continuation Approach to a Quadratic Matrix Equation
|
|
- Ethan Chambers
- 6 years ago
- Views:
Transcription
1 A Continuation Approach to a Quadratic Matrix Equation Nils Wagner nwagner@mecha.uni-stuttgart.de Institut A für Mechanik, Universität Stuttgart GAMM Workshop Applied and Numerical Linear Algebra September 22 24, 2005 Technische Universität Dresden Quadratic Matrix... Predictor-Corrector... Algorithms for... Generalized Sylvester... Bartels-Stewart type... Examples Applications Summary and Outlook Acknowledgments References Page 1 of 18
2 Contents Quadratic Matrix... Predictor-Corrector... Algorithms for... Generalized Sylvester... Bartels-Stewart type... Examples Applications Summary and Outlook Acknowledgments References Page 2 of 18
3 1. Quadratic Matrix Equation 1.1. Definition F (X, t) = A X 2 + t B X + C = O, A, B, C C n n (1) ( λ 2 A + t λ B + C ) x = 0 (2) The quadratic eigenvalue problem (2) is closely related to the matrix equation (1). λ 2 A + t λ B + C = (t B + A X + λ A) (X λ I) = F (λ I n, t) (3) When (1) has a solution X, the 2 n eigenvalues of F (λ I n, t) can be found by finding the eigenvalues of the matrix X and the matrix pencil (t B + A X, A). This solvent approach has been explored by Davis (1983), Dennis et al. (1976/1978) and more recently by Guo (2004) and Kim & Higham (2000, 2001). Page 3 of 18
4 1.2. Newton s method The algorithm proposed is Newton s method applied to the matrix function F (X, t). Having chosen an initial guess X 0 at t = t 0 = 0, we produce a sequence of iterates X (k) which will hopefully converge to S such that F (S, t) = O. Successive X (k) are generated by the formula where X (k+1) = X (k) + H (k), k = 0, 1, 2,... (4) H (k) = F (X, t) 1 F (X (k), t), k = 0, 1, 2,... (5) 1.3. Frechét Derivative Examining the form of the Frechét derivative F (X, t) we find F (X + H, t) = F (X, t) + A X H + A H X + t B H +A H }{{} 2 (6) F X (H) Page 4 of 18
5 2. Predictor-Corrector scheme Here, we assume that a solution X 0 does exist for t 0 = 0. Predictor A X C = O (7) Euler step A X 0 H + A H X 0 + t 0 B H = B X 0 (8) X (0) = X 0 + h H (9) Corrector Each H (k+1) is the solution of the system [ A X (k) H (k+1) + H (k+1) X (k)] + t B H (k+1) = F (X (k), t) (10) Until convergence Stopping Criteria X (k+1) = X (k) + H (k+1), k = 0, 1,... (11) F (X (k) ) F < ε (12) Page 5 of 18
6 3. Algorithms for Sylvester Equations An overview of different algorithms is given by Antoulas (2005) The Kronecker product method Bartels-Stewart Algorithm Complex integration method Characteristic polynomial methods Invariant subspace method Sign function method Infinite sum method Square root method... Page 6 of 18
7 4. Generalized Sylvester Equation The brute force way to approach the solution of (8,10) is to rewrite it as a linear matrix vector system in standard form. The system becomes [ ] I (t B + A X 0 ) + (X 0 ) T A vec H = (I B) vec X 0, (13) X (0) = X 0 + h H (14) [ ( I t B + A X (k)) + (X (k)) ] T A vec H (k+1) = vec F (X (k), t) (15) X (k+1) = X (k) + H (k+1) (16) Unfortunately, the coefficient matrix has dimension n 2 n 2, making this approach impractical except for small systems. A second alternative is to cast (8,10) in a form for which effective algorithms already exists. If we premultiply both sides of (8,10) by A 1 we obtain [ t0 A 1 B + X 0 ] H + H X0 = A 1 B X 0 (17) and [ t A 1 B + X (k)] H (k+1) + H (k+1) X (k) = A 1 F (X (k), t) (18) respectively. This method is is not generally satisfactory, if A is singular or illconditioned with respect to inversion. A two equation form could also be used (Kagström and Westin, 1989). Page 7 of 18
8 5. Bartels-Stewart type algorithm We study the linear matrix equation A X X B = C. (19) where A C n n, B C n n, and C C n n. This equation is now often called a Sylvester equation. This equation has a unique solution if and only if A and B have no eigenvalues in common. First B is reduced to Schur form. Schur form B V = V D, V H V = I (20) A Z Z D = R, R = C V, Z = X V, (21) k 1 (A d kk I) z k = r k + d ik z i (22) i=1 X = Z V H (23) Page 8 of 18
9 6. Examples 6.1. Example 1 F (X) has no solvents at all. Pereira (2003) F (X) = X 2 + A 1 X + A 2 (24) A 1 = 98/ /25 112/25 4/5 24/5 4/5 22/25 /38/25 182/ Example 2, A 2 = 89/25 294/25 316/25 7/5 42/5 8/5 46/25 59/25 251/25 Modelling oscillations of an airplane wing - Higham and Kim (2001) F (X) = A X 2 + B X + C (26) (25) A = C = , B = , (27) (28) Page 9 of 18
10 6.3. Example 3 Overdamped eigenvalue problem Higham and Kim (2000) F (X) = A X 2 + B X + C (29) A = I n, B = C = , (30) (31) Page 10 of 18
11 7. Applications F (X, t) = M X 2 + t D X + K = O, M, D, K R n n, M = M T (32) Cholesky Factorization Transformed matrices M = L L T (33) X = L T X L T, (34) D = L 1 D L T, (35) K = L 1 K L T (36) [ t 0 D + X0 ] H + H X0 = D X 0, (37) [ t D + X (k)] H(k+1) + H(k+1) X(k) X (0) 0 = X 0 + h H (38) = F ( X (k), t), (39) X (k+1) = X (k) + H (k+1) (40) Page 11 of 18
12 7.1. Fluid-Conveying Pipes 7.2. Hamilton s principle δ t 1 t 0 t 1 ( L mf v 2 ) w L dt t 0 m f v (ẇ L + v w L) δw L dt = 0, (41) Page 12 of 18 L L T = T p + T f = 1 m p ẇ 2 dx + 1 m f {v 2 + [ẇ + v w ] 2} dx, (42)
13 L L U = U p + U f = 1 EI w 2 dx (m p + m f ) g (L x) w 2 dx. (43) 0 0 EI w + { m f v 2 (m f + m p ) g(l x) } w + 2 m f vẇ + (m p + m f )gw + (m p + m f )ẅ = 0. (44) Nondimensional Quantities ξ = x L, η = w L, τ = EI m f + m p t L 2, u = mf v L, (45) EI γ = m p + m f L 3 m f g, β =, (46) EI m p + m f Page 13 of 18 η + u 2 η + γ (ξ 1) η + 2 βu η + γη + η = 0, (47) 7.3. Fourier-Galerkin Scheme N η(ξ, τ) η N (ξ, τ) = Ψ i (ξ) q i (τ), (48) i=1
14 Eigenfunctions of the cantilever beam Φ i (ξ) = cosh λ i ξ cos λ i ξ σ i (sinh λ i ξ sin λ i ξ), σ i := sinh λ i sin λ i cosh λ i + cos λ i. (49) 7.4. Second-Order Differential Equation I N q + C(u) q + K(u) q = 0 (50) C(u) = 2 β u B, (51) K(u) = A + ( u 2 γ ) G + γ (D + B). (52) F (X, u) = X β u B X + γ (B + D G) + A + u 2 G (53) Page 14 of 18
15 8. Summary and Outlook Quadratic eigenvalue problems are typically solved using the companion matrix technique in which the nonlinear eigenvalue problem is linearized. The obvious disadvantage of such an approach is that the size of the problem essentially doubles. This contribution introduces a method that tackles the quadratic eigenvalue problem directly using an approach suggested by Lancaster and Higham, to name a few. It is well known that the quadratic eigenvalue problem is closely related to a (parameterized) quadratic matrix equation. When solutions, so-called solvents, of this nonlinear matrix equation exist the eigenvalues of the associated quadratic eigenvalue problem can be found by the eigenvalues of the solvent and a matrix pencil. The parameterized quadratic matrix equation is then solved by a continuation approach with a standard Euler-Newton predictor-corrector scheme, where (generalized) Sylvester equations occur. The Sylvester equations can be solved using the Bartels-Stewart algorithm or the Hessenberg-Schur method. These methods require certain relatively expensive initial factorizations (Schur or Hessenberg) of the coefficient matrices. In this context we should make use of iterative methods (see e.g. Benner (2004), Ding (2005),...). Page 15 of 18
16 9. Acknowledgments Thanks to Vera Thümmler for helpful discussions on Bartels-Stewart like algorithms and continuation. Page 16 of 18
17 10. References Bartels and Stewart, Solution of the matrix equation AX+XB=C, Communications of the ACM, (1972) Kagström and Westin, Generalized Schur methods with condition estimators for solving the generalized Sylvester equation, IEEE Transactions on Automatic Control, (1989) Gardiner, Laub, Amato and Moler, Solution of the Sylvester matrix equation AXB T + CXD T = E, ACM Transactions of Mathematical Software, (1992) Dennis, Traub and Weber, The algebraic theory of matrix polynomials, SIAM Journal on Numerical Analysis, (1976) Dennis, Traub and Weber, Algorithms for solvents of matrix polynomials, SIAM Journal on Numerical Analysis, (1978) Davis, An algorithm to compute solvents of the matrix equation AX 2 +BX + C = 0, ACM Transactions on Mathematical Software, (1983) Kratz and Stickel, Numerical solution of matrix polynomial equations by Newton s method, IMA Journal of Numerical Analysis (1987) Bras and de Lima, A spectral approach to polynomial matrices solvents, Appl. Math. Lett (1996) Higham and Kim, Numerical analysis of a quadratic matrix equation, IMA Journal of Numerical Analysis, (2000) Page 17 of 18
18 Higham and Kim, Solving a quadratic matrix equation by Newton s method with exact line searches, SIAM J. Matrix Anal. Appl (2001) Pereira, On Solvents of matrix polynomials, Applied Numerical Mathematics, (2003) Guo and Lancaster, Algorithms for hyperbolic quadratic eigenvalue problems, Mathematics of Computation, (2005) Benner, Peter ; Quintana-Orti, Enrique ; Quintana-Orti, Gregorio, Solving Linear Matrix Equations via Rational Iterative Schemes Technische Universität Chemnitz, SFB 393 (Germany), SFB393-Preprint 8 (2004) Ding and Chen, Gradient based iterative algorithms for solving a class of matrix equations, IEEE Transactions on Automatic Control, (2005) Antoulas, Approximation of large-scale dynamical systems, SIAM (2005) Page 18 of 18
Improved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationPositive Solution to a Generalized Lyapunov Equation via a Coupled Fixed Point Theorem in a Metric Space Endowed With a Partial Order
Filomat 29:8 2015, 1831 1837 DOI 10.2298/FIL1508831B Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Positive Solution to a Generalized
More informationOn a quadratic matrix equation associated with an M-matrix
Article Submitted to IMA Journal of Numerical Analysis On a quadratic matrix equation associated with an M-matrix Chun-Hua Guo Department of Mathematics and Statistics, University of Regina, Regina, SK
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationAlgorithms for Solving the Polynomial Eigenvalue Problem
Algorithms for Solving the Polynomial Eigenvalue Problem Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with D. Steven Mackey
More information1 Sylvester equations
1 Sylvester equations Notes for 2016-11-02 The Sylvester equation (or the special case of the Lyapunov equation) is a matrix equation of the form AX + XB = C where A R m m, B R n n, B R m n, are known,
More informationFactorized Solution of Sylvester Equations with Applications in Control
Factorized Solution of Sylvester Equations with Applications in Control Peter Benner Abstract Sylvester equations play a central role in many areas of applied mathematics and in particular in systems and
More informationAn error estimate for matrix equations
Applied Numerical Mathematics 50 (2004) 395 407 www.elsevier.com/locate/apnum An error estimate for matrix equations Yang Cao, Linda Petzold 1 Department of Computer Science, University of California,
More informationApproximate Low Rank Solution of Generalized Lyapunov Matrix Equations via Proper Orthogonal Decomposition
Applied Mathematical Sciences, Vol. 4, 2010, no. 1, 21-30 Approximate Low Rank Solution of Generalized Lyapunov Matrix Equations via Proper Orthogonal Decomposition Amer Kaabi Department of Basic Science
More informationKrylov Subspace Type Methods for Solving Projected Generalized Continuous-Time Lyapunov Equations
Krylov Subspace Type Methods for Solving Proected Generalized Continuous-Time Lyapunov Equations YUIAN ZHOU YIQIN LIN Hunan University of Science and Engineering Institute of Computational Mathematics
More informationContour integral solutions of Sylvester-type matrix equations
Contour integral solutions of Sylvester-type matrix equations Harald K. Wimmer Mathematisches Institut, Universität Würzburg, 97074 Würzburg, Germany Abstract The linear matrix equations AXB CXD = E, AX
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationChebfun2: Exploring Constant Coefficient PDEs on Rectangles
Chebfun2: Exploring Constant Coefficient PDEs on Rectangles Alex Townsend Oxford University with Sheehan Olver, Sydney University FUN13 workshop, 11th of April 2013 Alex Townsend Constant coefficient PDEs
More informationAn Exact Line Search Method for Solving Generalized Continuous-Time Algebraic Riccati Equations
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 43, NO., JANUARY 998 0 An Exact Line Search Method for Solving Generalized Continuous-Time Algebraic Riccati Equations Peter Benner and Ralph Byers Abstract
More informationParallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation
Parallel Model Reduction of Large Linear Descriptor Systems via Balanced Truncation Peter Benner 1, Enrique S. Quintana-Ortí 2, Gregorio Quintana-Ortí 2 1 Fakultät für Mathematik Technische Universität
More informationOn GPU Acceleration of Common Solvers for (Quasi-) Triangular Generalized Lyapunov Equations
Max Planck Institute Magdeburg Preprints Martin Köhler Jens Saak On GPU Acceleration of Common Solvers for (Quasi-) Triangular Generalized Lyapunov Equations MAX PLANCK INSTITUT FÜR DYNAMIK KOMPLEXER TECHNISCHER
More informationAvailable online at ScienceDirect. Procedia Engineering 100 (2015 ) 56 63
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 100 (2015 ) 56 63 25th DAAAM International Symposium on Intelligent Manufacturing and Automation, DAAAM 2014 Definite Quadratic
More informationSolving projected generalized Lyapunov equations using SLICOT
Solving projected generalized Lyapunov equations using SLICOT Tatjana Styel Abstract We discuss the numerical solution of projected generalized Lyapunov equations. Such equations arise in many control
More informationAn Arnoldi method with structured starting vectors for the delay eigenvalue problem
An Arnoldi method with structured starting vectors for the delay eigenvalue problem Elias Jarlebring, Karl Meerbergen, Wim Michiels Department of Computer Science, K.U. Leuven, Celestijnenlaan 200 A, 3001
More informationAPPLIED NUMERICAL LINEAR ALGEBRA
APPLIED NUMERICAL LINEAR ALGEBRA James W. Demmel University of California Berkeley, California Society for Industrial and Applied Mathematics Philadelphia Contents Preface 1 Introduction 1 1.1 Basic Notation
More informationTheorem Let X be a symmetric solution of DR(X) = 0. Let X b be a symmetric approximation to Xand set V := X? X. b If R b = R + B T XB b and I + V B R
Initializing Newton's method for discrete-time algebraic Riccati equations using the buttery SZ algorithm Heike Fabender Zentrum fur Technomathematik Fachbereich - Mathematik und Informatik Universitat
More informationOn the Computations of Eigenvalues of the Fourth-order Sturm Liouville Problems
Int. J. Open Problems Compt. Math., Vol. 4, No. 3, September 2011 ISSN 1998-6262; Copyright c ICSRS Publication, 2011 www.i-csrs.org On the Computations of Eigenvalues of the Fourth-order Sturm Liouville
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationCS760, S. Qiao Part 3 Page 1. Kernighan and Pike [5] present the following general guidelines for testing software.
CS760, S. Qiao Part 3 Page 1 Testing 1 General Guidelines Kernighan and Pike [5] present the following general guidelines for testing software. Know exactly what you are testing and what results you expect.
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon's method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 2. First Results and Algorithms Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Friday, 7 June, 2013 1:30 pm to 4:30 pm PAPER 4 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question
More informationResearch Matters. February 25, The Nonlinear Eigenvalue Problem. Nick Higham. Part III. Director of Research School of Mathematics
Research Matters February 25, 2009 The Nonlinear Eigenvalue Problem Nick Higham Part III Director of Research School of Mathematics Françoise Tisseur School of Mathematics The University of Manchester
More informationBalanced Truncation Model Reduction of Large and Sparse Generalized Linear Systems
Balanced Truncation Model Reduction of Large and Sparse Generalized Linear Systems Jos M. Badía 1, Peter Benner 2, Rafael Mayo 1, Enrique S. Quintana-Ortí 1, Gregorio Quintana-Ortí 1, A. Remón 1 1 Depto.
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationLow Rank Solution of Data-Sparse Sylvester Equations
Low Ran Solution of Data-Sparse Sylvester Equations U. Baur e-mail: baur@math.tu-berlin.de, Phone: +49(0)3031479177, Fax: +49(0)3031479706, Technische Universität Berlin, Institut für Mathemati, Straße
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationOrder reduction numerical methods for the algebraic Riccati equation. V. Simoncini
Order reduction numerical methods for the algebraic Riccati equation V. Simoncini Dipartimento di Matematica Alma Mater Studiorum - Università di Bologna valeria.simoncini@unibo.it 1 The problem Find X
More informationNUMERICAL METHODS USING MATLAB
NUMERICAL METHODS USING MATLAB Dr John Penny George Lindfield Department of Mechanical Engineering, Aston University ELLIS HORWOOD NEW YORK LONDON TORONTO SYDNEY TOKYO SINGAPORE Preface 1 An introduction
More informationEfficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers
Efficient Implementation of Large Scale Lyapunov and Riccati Equation Solvers Jens Saak joint work with Peter Benner (MiIT) Professur Mathematik in Industrie und Technik (MiIT) Fakultät für Mathematik
More informationModel Reduction for Dynamical Systems
Otto-von-Guericke Universität Magdeburg Faculty of Mathematics Summer term 2015 Model Reduction for Dynamical Systems Lecture 8 Peter Benner Lihong Feng Max Planck Institute for Dynamics of Complex Technical
More information2 Computing complex square roots of a real matrix
On computing complex square roots of real matrices Zhongyun Liu a,, Yulin Zhang b, Jorge Santos c and Rui Ralha b a School of Math., Changsha University of Science & Technology, Hunan, 410076, China b
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More informationAn Accelerated Jacobi-gradient Based Iterative Algorithm for Solving Sylvester Matrix Equations
Filomat 31:8 (2017), 2381 2390 DOI 10.2298/FIL1708381T Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat An Accelerated Jacobi-gradient
More informationOn matrix equations X ± A X 2 A = I
Linear Algebra and its Applications 326 21 27 44 www.elsevier.com/locate/laa On matrix equations X ± A X 2 A = I I.G. Ivanov,V.I.Hasanov,B.V.Minchev Faculty of Mathematics and Informatics, Shoumen University,
More informationOn the computation of the Jordan canonical form of regular matrix polynomials
On the computation of the Jordan canonical form of regular matrix polynomials G Kalogeropoulos, P Psarrakos 2 and N Karcanias 3 Dedicated to Professor Peter Lancaster on the occasion of his 75th birthday
More informationS.F. Xu (Department of Mathematics, Peking University, Beijing)
Journal of Computational Mathematics, Vol.14, No.1, 1996, 23 31. A SMALLEST SINGULAR VALUE METHOD FOR SOLVING INVERSE EIGENVALUE PROBLEMS 1) S.F. Xu (Department of Mathematics, Peking University, Beijing)
More informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationGAUSS-SIDEL AND SUCCESSIVE OVER RELAXATION ITERATIVE METHODS FOR SOLVING SYSTEM OF FUZZY SYLVESTER EQUATIONS
GAUSS-SIDEL AND SUCCESSIVE OVER RELAXATION ITERATIVE METHODS FOR SOLVING SYSTEM OF FUZZY SYLVESTER EQUATIONS AZIM RIVAZ 1 AND FATEMEH SALARY POUR SHARIF ABAD 2 1,2 DEPARTMENT OF MATHEMATICS, SHAHID BAHONAR
More informationChemnitz Scientific Computing Preprints
Peter Benner Heike Faßbender On the numerical solution of large-scale sparse discrete-time Riccati equations CSC/09-11 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationTwo improved classes of Broyden s methods for solving nonlinear systems of equations
Available online at www.isr-publications.com/jmcs J. Math. Computer Sci., 17 (2017), 22 31 Research Article Journal Homepage: www.tjmcs.com - www.isr-publications.com/jmcs Two improved classes of Broyden
More informationBalancing-Related Model Reduction for Large-Scale Systems
Balancing-Related Model Reduction for Large-Scale Systems Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz D-09107 Chemnitz benner@mathematik.tu-chemnitz.de
More informationSolvability of Linear Matrix Equations in a Symmetric Matrix Variable
Solvability of Linear Matrix Equations in a Symmetric Matrix Variable Maurcio C. de Oliveira J. William Helton Abstract We study the solvability of generalized linear matrix equations of the Lyapunov type
More informationConvergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations
Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations Wen-Wei Lin Shu-Fang u Abstract In this paper a structure-preserving transformation of a symplectic pencil
More informationMATLAB TOOLS FOR SOLVING PERIODIC EIGENVALUE PROBLEMS 1
MATLAB TOOLS FOR SOLVING PERIODIC EIGENVALUE PROBLEMS 1 Robert Granat Bo Kågström Daniel Kressner,2 Department of Computing Science and HPC2N, Umeå University, SE-90187 Umeå, Sweden. {granat,bokg,kressner}@cs.umu.se
More informationStructure preserving Krylov-subspace methods for Lyapunov equations
Structure preserving Krylov-subspace methods for Lyapunov equations Matthias Bollhöfer, André Eppler Institute Computational Mathematics TU Braunschweig MoRePas Workshop, Münster September 17, 2009 System
More informationConsistency and efficient solution of the Sylvester equation for *-congruence
Electronic Journal of Linear Algebra Volume 22 Volume 22 (2011) Article 57 2011 Consistency and efficient solution of the Sylvester equation for *-congruence Fernando de Teran Froilan M. Dopico Follow
More informationNonlinear palindromic eigenvalue problems and their numerical solution
Nonlinear palindromic eigenvalue problems and their numerical solution TU Berlin DFG Research Center Institut für Mathematik MATHEON IN MEMORIAM RALPH BYERS Polynomial eigenvalue problems k P(λ) x = (
More informationParametrische Modellreduktion mit dünnen Gittern
Parametrische Modellreduktion mit dünnen Gittern (Parametric model reduction with sparse grids) Ulrike Baur Peter Benner Mathematik in Industrie und Technik, Fakultät für Mathematik Technische Universität
More informationPONDI CHERRY UNI VERSI TY
B.Sc. ALLIED MATHEMATICS SYLLABUS 2009-2010 onwards PONDI CHERRY UNI VERSI TY PUDUCHERRY 605 014 B.Sc. ALLIED MATHEMATICS Syllabus for Allied Mathematics for B.Sc. Physics Main/Chemistry Main/Electronics
More informationNonlinear Eigenvalue Problems: An Introduction
Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction
More informationON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS
ON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS ALEKSANDRA KOSTIĆ AND HEINRICH VOSS Key words. eigenvalue, variational characterization, principle, Sylvester s law of inertia AMS subject
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationAutomatica, 33(9): , September 1997.
A Parallel Algorithm for Principal nth Roots of Matrices C. K. Koc and M. _ Inceoglu Abstract An iterative algorithm for computing the principal nth root of a positive denite matrix is presented. The algorithm
More informationIndex. for generalized eigenvalue problem, butterfly form, 211
Index ad hoc shifts, 165 aggressive early deflation, 205 207 algebraic multiplicity, 35 algebraic Riccati equation, 100 Arnoldi process, 372 block, 418 Hamiltonian skew symmetric, 420 implicitly restarted,
More informationApproximation algorithms for nonnegative polynomial optimization problems over unit spheres
Front. Math. China 2017, 12(6): 1409 1426 https://doi.org/10.1007/s11464-017-0644-1 Approximation algorithms for nonnegative polynomial optimization problems over unit spheres Xinzhen ZHANG 1, Guanglu
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationS N. hochdimensionaler Lyapunov- und Sylvestergleichungen. Peter Benner. Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz
Ansätze zur numerischen Lösung hochdimensionaler Lyapunov- und Sylvestergleichungen Peter Benner Mathematik in Industrie und Technik Fakultät für Mathematik TU Chemnitz S N SIMULATION www.tu-chemnitz.de/~benner
More informationStructured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries
Structured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries Fakultät für Mathematik TU Chemnitz, Germany Peter Benner benner@mathematik.tu-chemnitz.de joint work with Heike Faßbender
More informationGeneralized Shifted Inverse Iterations on Grassmann Manifolds 1
Proceedings of the Sixteenth International Symposium on Mathematical Networks and Systems (MTNS 2004), Leuven, Belgium Generalized Shifted Inverse Iterations on Grassmann Manifolds 1 J. Jordan α, P.-A.
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationDirect Methods for Matrix Sylvester and Lyapunov Equations
Direct Methods for Matrix Sylvester and Lyapunov Equations D. C. Sorensen and Y. Zhou Dept. of Computational and Applied Mathematics Rice University Houston, Texas, 77005-89. USA. e-mail: {sorensen,ykzhou}@caam.rice.edu
More informationAN ITERATIVE METHOD TO SOLVE SYMMETRIC POSITIVE DEFINITE MATRIX EQUATIONS
AN ITERATIVE METHOD TO SOLVE SYMMETRIC POSITIVE DEFINITE MATRIX EQUATIONS DAVOD KHOJASTEH SALKUYEH and FATEMEH PANJEH ALI BEIK Communicated by the former editorial board Let A : R m n R m n be a symmetric
More informationParameter-Dependent Eigencomputations and MEMS Applications
Parameter-Dependent Eigencomputations and MEMS Applications David Bindel UC Berkeley, CS Division Parameter-Dependent Eigencomputations and MEMS Applications p.1/36 Outline Some applications Mathematical
More informationThe matrix sign function
The matrix sign function 1 Re x > 0, sign(x) = 1 Re x < 0, undefined Re x = 0. Suppose the Jordan form of A is reblocked as A = [ V 1 V 2 ] [ J 1 J 2 ] [ V 1 V 2 ] 1, where J 1 contains all eigenvalues
More informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
More informationA BOUNDARY VALUE PROBLEM FOR LINEAR PDAEs
Int. J. Appl. Math. Comput. Sci. 2002 Vol.12 No.4 487 491 A BOUNDARY VALUE PROBLEM FOR LINEAR PDAEs WIESŁAW MARSZAŁEK ZDZISŁAW TRZASKA DeVry College of Technology 630 US Highway One North Brunswick N8902
More informationSOLUTION OF ALGEBRAIC RICCATI EQUATIONS ARISING IN CONTROL OF PARTIAL DIFFERENTIAL EQUATIONS
SOLUTION OF ALGEBRAIC RICCATI EQUATIONS ARISING IN CONTROL OF PARTIAL DIFFERENTIAL EQUATIONS Kirsten Morris Department of Applied Mathematics University of Waterloo Waterloo, CANADA N2L 3G1 kmorris@uwaterloo.ca
More informationA Smoothing Newton Method for Solving Absolute Value Equations
A Smoothing Newton Method for Solving Absolute Value Equations Xiaoqin Jiang Department of public basic, Wuhan Yangtze Business University, Wuhan 430065, P.R. China 392875220@qq.com Abstract: In this paper,
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface
More information1 Quasi-definite matrix
1 Quasi-definite matrix The matrix H is a quasi-definite matrix, if there exists a permutation matrix P such that H qd P T H11 H HP = 1 H1, 1) H where H 11 and H + H1H 11 H 1 are positive definite. This
More informationBOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 4, P AGES 655 664 BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION Guang-Da Hu and Qiao Zhu This paper is concerned with bounds of eigenvalues of a complex
More informationTotal least squares. Gérard MEURANT. October, 2008
Total least squares Gérard MEURANT October, 2008 1 Introduction to total least squares 2 Approximation of the TLS secular equation 3 Numerical experiments Introduction to total least squares In least squares
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationEquations in Quadratic Form
Equations in Quadratic Form MATH 101 College Algebra J. Robert Buchanan Department of Mathematics Summer 2012 Objectives In this lesson we will learn to: make substitutions that allow equations to be written
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More informationAlexander Ostrowski
Ostrowski p. 1/3 Alexander Ostrowski 1893 1986 Walter Gautschi wxg@cs.purdue.edu Purdue University Ostrowski p. 2/3 Collected Mathematical Papers Volume 1 Determinants Linear Algebra Algebraic Equations
More informationResidual iterative schemes for largescale linear systems
Universidad Central de Venezuela Facultad de Ciencias Escuela de Computación Lecturas en Ciencias de la Computación ISSN 1316-6239 Residual iterative schemes for largescale linear systems William La Cruz
More informationMatrix Algorithms. Volume II: Eigensystems. G. W. Stewart H1HJ1L. University of Maryland College Park, Maryland
Matrix Algorithms Volume II: Eigensystems G. W. Stewart University of Maryland College Park, Maryland H1HJ1L Society for Industrial and Applied Mathematics Philadelphia CONTENTS Algorithms Preface xv xvii
More information1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl)
DDM Preprint Editors: editor1, editor2, editor3, editor4 c DDM.org 1. Fast Solvers and Schwarz Preconditioners for Spectral Nédélec Elements for a Model Problem in H(curl) Bernhard Hientzsch 1 1. Introduction.
More informationAn Algorithm for. Nick Higham. Françoise Tisseur. Director of Research School of Mathematics.
An Algorithm for the Research Complete Matters Solution of Quadratic February Eigenvalue 25, 2009 Problems Nick Higham Françoise Tisseur Director of Research School of Mathematics The School University
More informationPolynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations
Polynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations Elias Jarlebring a, Michiel E. Hochstenbach b,1, a Technische Universität Braunschweig,
More informationFIXED POINT ITERATION
FIXED POINT ITERATION The idea of the fixed point iteration methods is to first reformulate a equation to an equivalent fixed point problem: f (x) = 0 x = g(x) and then to use the iteration: with an initial
More informationRecursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q
Applied Mathematical Sciences, Vol. 2, 2008, no. 38, 1855-1872 Recursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q Nicholas Assimakis Department of Electronics, Technological
More informationArnoldi Methods in SLEPc
Scalable Library for Eigenvalue Problem Computations SLEPc Technical Report STR-4 Available at http://slepc.upv.es Arnoldi Methods in SLEPc V. Hernández J. E. Román A. Tomás V. Vidal Last update: October,
More informationCondition Numbers and Backward Error of a Matrix Polynomial Equation Arising in Stochastic Models
J Sci Comput (2018) 76:759 776 https://doi.org/10.1007/s10915-018-0641-x Condition Numbers and Backward Error of a Matrix Polynomial Equation Arising in Stochastic Models Jie Meng 1 Sang-Hyup Seo 1 Hyun-Min
More information1.1: The bisection method. September 2017
(1/11) 1.1: The bisection method Solving nonlinear equations MA385/530 Numerical Analysis September 2017 3 2 f(x)= x 2 2 x axis 1 0 1 x [0] =a x [2] =1 x [3] =1.5 x [1] =b 2 0.5 0 0.5 1 1.5 2 2.5 1 Solving
More informationMICHIEL E. HOCHSTENBACH
VARIATIONS ON HARMONIC RAYLEIGH RITZ FOR STANDARD AND GENERALIZED EIGENPROBLEMS MICHIEL E. HOCHSTENBACH Abstract. We present several variations on the harmonic Rayleigh Ritz method. First, we introduce
More informationFinal Year M.Sc., Degree Examinations
QP CODE 569 Page No Final Year MSc, Degree Examinations September / October 5 (Directorate of Distance Education) MATHEMATICS Paper PM 5: DPB 5: COMPLEX ANALYSIS Time: 3hrs] [Max Marks: 7/8 Instructions
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationDownloaded 05/27/14 to Redistribution subject to SIAM license or copyright; see
SIAM J. MATRIX ANAL. APPL. Vol. 31, No. 2, pp. 272 288 c 2009 Society for Industrial and Applied Mathematics INEXACT KLEINMAN NEWTON METHOD FOR RICCATI EQUATIONS F. FEITZINGER, T. HYLLA, AND E. W. SACHS
More information