On a root-finding approach to the polynomial eigenvalue problem
|
|
- Marybeth Cobb
- 6 years ago
- Views:
Transcription
1 On a root-finding approach to the polynomial eigenvalue problem Dipartimento di Matematica, Università di Pisa bini Joint work with Vanni Noferini and Meisam Sharify Limoges, May, 10, 2012
2 The problem Given m m matrices A i, i = 0,..., n compute the eigenvalues λ such that A(λ)v = 0, v 0, where A(λ) = n λ i A i, det A(λ) 0 i=0 Standard matrix approach: linearization (A λb)w = 0, for suitable mn mn matrices A, B Polynomial approach: solving the scalar equation a(λ) = 0, where a(λ) = det(a(λ))
3 Structured matrix polynomials Skew-symmetric coefficients of even size: eigenvalues with even multiplicities since a(x) = q(x) 2 roots appear in pairs (λ, λ) Palindromic polynomials: A i = A n i T-palindromic polynomials: A i = A T n i roots appear in pairs (λ, 1/λ) including eigenvalues at infinity Hamiltonian polynomials: A 2i is Hamiltonian, A 2i+1 is skew-hamiltonian roots appear in pairs (λ, λ) General property of the roots: roots appear in pairs (λ, f (λ)), where f (z) is such that f (f (z)) = z Problem: to design a polynomial eigensolver which exploits the structure, and delivers approximations in pairs (λ, f (λ))
4 The Ehrlich-Aberth iteration Given a scalar polynomial a(x) of degree N, given approximations x (0) 1,..., x (0) N to the roots λ 1,..., λ N, the Ehrlich-Aberth method generates the sequence x (ν+1) i = x (ν) i N (x (ν) i ) 1 N (x (ν) i ) N j=1, j i x (ν) i 1 x (ν) j, i = 1,..., N where N (x) = a(x)/a (x) is the Newton correction This iteration is just Newton s iteration applied to the rational functions (implicit deflation) g i (x) = a(x) n j=1, j i (x x j)
5 Some features of the E-A iteration Local convergence to simple roots is cubic Local convergence to simple roots is more than cubic in the Gauss-Seidel fashion Local convergence to multiple roots is linear Global convergence is not proved: practically, convergence always holds if initial approximations are equally placed along a circle Cost per iteration: O(N 2 ) + N Cost of computing N (x) Number of iterations for arriving at numerical convergence is practically bounded under a suitable choice of starting approximations
6 Problems related to the E-A implementation computing N (x) choosing initial approximations design an effective stop condition design a posteriori error bounds design effective strategies for accelerating convergence in case of clusters or multiple roots These issues were treated in the package MPSolve for approximating polynomial zeros up to any number of (guaranteed) digits [Bini, Fiorentino, Numer. Algo. 2000] We try to do the same for (structured) matrix polynomials
7 Exploiting the structure: Newton s correction For a general m m matrix polynomial A(x) of degree n, set a(x) = det A(x). Then N (x) = a(x) a (x) = 1 trace(a 1 (x)a (x)) Computational cost: O(m 3 + m 2 n) ops
8 Exploiting the structure: Newton s correction Structured case (palindromic and T-palindromic polynomials): Assume that N = mn is even, set z = z(x) = x + 1/x then q(z) = x(z) mn/2 a(x(z)) is a polynomial of degree mn/2, where x(z) = (z z 2 4)/2 or x(z) = (z + z 2 4)/2 is any of the two branches of the inverse of the function z(x).
9 Exploiting the structure: Newton s correction Structured case (palindromic and T-palindromic polynomials): Assume that N = mn is even, set z = z(x) = x + 1/x then q(z) = x(z) mn/2 a(x(z)) is a polynomial of degree mn/2, where x(z) = (z z 2 4)/2 or x(z) = (z + z 2 4)/2 is any of the two branches of the inverse of the function z(x). Moreover, q(z) q (z) = 1 1/x 2 g(x) mn/(2x), g(x) = trace(a(x) 1 A (x)) This way, the computation of the roots of the palyndromic polynomial a(x) of degree mn is reduced to computing the roots of the polynomial q(z) of degree mn/2
10 Exploiting the structure: Newton s correction Structured case (palindromic and T-palindromic polynomials): Assume mn odd In this case 1 is root. Moreover the matrix polynomial Â(x) = (x + 1)A(x) is such that (n + 1)m is even and 1 has multiplicity at least m + 1. Aberth iteration can be applied to Â(x) with m + 1 components of the initial vector equal to 1.
11 Exploiting the structure: Newton s correction Structured case: Hamiltonian polynomials (mn is even) Set z = x 2. For x(z) = z or x(z) = z, q(z) = a(x(z)) is a polynomial of degree mn/2. Moreover, q(z) q (z) = 2x g(x) 1, g(x) = trace(a(x) 1 A (x)) x If mn is odd, then q(z) = a(x(z))/x(z) is a polynomial and q(z) q (z) = 2x 3 g(x) 1, g(x) = trace(a(x) 1 A (x)) Once the roots z 1,..., z mn/2 have been computed, the roots of a(x) are given by the pairs ( z i, z i ), i = 1,..., mn/2
12 Exploiting the structure: Newton s correction General case: roots in pairs (x, f (x)), f (x) = (αx + β)/(γx + α), α 2 + βγ 0 Set z(x) = xf (x) = αx2 +βx γx α and denote x(z) one of the two branches of the inverse of z(x). If there are no eigenvalues with odd multiplicity, then is a polynomial q(z) = a(x(z)) (γx(z) α) mn/2
13 Computational cost The overall cost of E-A iteration applied to an m m matrix polynomial of degree n is O(m 4 n + m 3 n 2 ) ops The overall cost of the technique based on linearization is O(m 3 n 3 ) ops From the complexity point of view the E-A approach is more convenient if n > m The cost of E-A, i.e., the computation of trace(a(x) 1 A (x)), can be reduced to O(m 2 n 2 ) ops, by using linearization. This leads to the overall cost of O(m 3 n 3 ) ops The drawback is that linearization may increase the condition number of eigenvalues.
14 Choosing initial approximations The choice of initial approximations is crucial to arrive at numerical convergence with a small number of iterations This is particularly important for polynomials having very unbalanced roots The package MPSolve for computing polynomial roots to any guaranteed precision, relies on a very effective tool: the Newton polygon construction We synthesize the main ideas of this tool and then extend it to the case of matrix polynomials.
15 Choosing initial approximations The Rouché theorem: If p(x) and q(x) are polynomials such that p(x) > q(x) for x = r then p(x) and p(x) + q(x) have the same number of roots in the open disk of center 0 and radius r Given a(x) = n i=0 a ix i apply Rouché theorem with p(x) = a k x k and q(x) = n i=0, i k x i a i. It follows that if p(x) = a k x k > n i=0,i k a i x i q(x), x = r then a(x) = p(x) + q(x) has k roots of modulus less than r.
16 Choosing initial approximations Properties: The equation has a k r k = n i=0,i k a i r i one positive solution t 0 if k = 0 one positive solution s n if k = n either no positive solutions or two positive solutions s t otherwise Pellet s Theorem, 1881 If h 1,..., h p are the values of k for which the above equation has either one or two positive solutions s hi t hi then the open annulus with radii s hi, t hi contains no roots of any polynomial equimodular with p(x) the closed annulus with radii t hi 1, s hi contains h i h i 1 roots
17 Choosing initial approximations
18 Choosing initial approximations: strategy of MPSolve Computing the positive roots of n polynomials is expensive. An alternative solution relies on the following properties: Any positive solution r of the blue equation, if it exists, is such that u k := max i<k a i a k 1 k i < r < min i>k a i a k 1 k i =: v k Let k i, i = 1,... q be such that u ki < v ki. Then v ki = u ki+1, moreover any interval [s hj, t hj ] does not contain any u ki The indices k i are the abscissas of the Newton polygon: upper convex hull of the set {(i, log a i ), i = 0,..., n} The values u ki are the exponential of the slopes of the edges of the Newton polygon. They are called tropical roots
19 Choosing initial approximations: strategy of MPSolve
20 Choosing initial approximations: strategy of MPSolve
21 Starting approximations: an example p(x) = x 10 + (1000x + 1) 3 Coefficient vector: (1, 3000, , , 0, 0, 0, 0, 0, 0, 1) There are: 3 roots of modulus close to 1/1000, 7 roots of modulus close to The Newton polygon has three vertices besides 0 and n. The values of u k are , 0.001, 0.003, initial approximations are placed in a circle of radius initial approximations are placed inside a disk of radius 3/1000
22 The Sharify improvement What happens inside a blue annulus? If a vertex of the polynomial is sufficiently sharp then we may arrive at a strong localization of the roots Denote r i = u ki the value of the radius corresponding to the vertex k i m i = k i+1 k i its multiplicity δ i = r i /r i+1 the ratio of two consecutive radii Theorem (Sharify) If δ i, δ i 1 < 1/9 then the polynomial a(x) has m i roots in the annulus with radii r i /3 and 3r i.
23 Extension to matrix polynomials Let A, B Hermitian matrices, define A B if A B is positive definite, and A B if A B is positive semidefinite Let F (x), G(x) be matrix polynomials Theorem (Rouché theorem for matrix polynomials) If F (x) F (x) G(x) G(x) for x = r then F (x) and F (x) + G(x) have the same number of roots in the disc of center 0 and radius r. Corollary If (A k A k)r 2k ( n i=0, i k A ix i ) ( n i=0, i k A ix i ), for x = r then det A(x) has kn roots inside the disk of center 0 and radius r. The above condition requires that det A k 0 and is implied by n r k A 1 k A i r i i=0, i k
24 Generalization Theorem 1. If 0 < k < n is such that det A k 0 then the equation r k = A 1 k A i r i (1) i=0, i k has either no real positive solution or two real positive solutions s k t k. 2. In the latter case, the polynomial a(x) = det A(x) has no roots in the open annulus of radii s k, t k, while it has mk roots in the disk of radius s k. 3. If k = 0 and det A 0 0 then (1) has only one real positive root t 0, moreover, the polynomial a(x) has no root in the open disk of radius t If k = n and det A n 0, then (1) has only one real positive solution s n and the polynomial a(x) has no root outside the disk of radius s n.
25 Corollary (Generalized Pellet Thm.) Let h 0 < h 1 <... < h q be the values of k such that det A k 0 and there exist positive real solution(s) s hi t hi of (1). Then 1. t hi 1 s hi, i = 1,..., q; 2. there are m(h i h i 1 ) roots of a(x) in the annulus of radii t hi 1, s hi ; 3. there are no roots of the polynomial a(x) in the open annulus of radii s hi, t hi, where i = 0, 1,..., q and we assume that s 0 = 0, t n =.
26 Theorem If k is such that (1) has two real positive solutions s k t k, then u k s k t k v k where u k := max i<k A 1 k A i 1/(k i), v k := min i>k A 1 k A i 1/(k i) If, for k = 0 (1) has a solution t 0 then t 0 v 0 := min i>0 A 1 0 A i 1/i If for k = n (1) has a solution s n then s n u n := max i<n A 1 n A i 1/(n i) Moreover, v ki u ki+1, i = 1,..., p 1.
27 A special class Consider the class Q m,n of matrix polynomials A(x) = n i=0 x i A i, where A i = α i Q i, Q i Q i = I, α i 0 The condition given in the blue equation (1) turns into α k r k = n i=0,i k α i r i and leads to the machinery valid in the scalar case for the location of the roots of a polynomial by means of the Newton polygon. This is the best result that we can obtain for polynomials in the class Q m,n. In fact, for Q i = I we obtain exactly n copies of the scalar polynomial p(x) = n i=0 α ix i for which the inclusion result is optimal.
28 Theorem Let S = {h i, k S. Then, i = 0,..., p} be such that u k < v k if and only if 1. {k 0,..., k q } {h 0,..., h p }, moreover {h 0,..., h p } [s ki, t ki ] = ; 2. v ki u ki+1 ; 3. if A(x) Q m,n, then v ki = u ki+1 and v ki coincide with the vertices of the Newton polygon of the polynomial w(x) = n i=0 A i x i = n i=0 α ix i.
29 The Sharify bound for matrix polynomials Let A(x) = n i=0 x i A i Q m,n, i.e., A i C m m, A i = α i Q i, Q i Q i = I. Let a(x) = det A(x). Define r i the radii (tropical roots) corresponding to the i-th vertex of the Newton polygon of the polynomial n x i A i 2 i=0 and m i their multiplicities, define δ i = r i /r i 1 Theorem For A(x) Q m,n, if δ i, δ i 1 < 1/12.12 then the polynomial a(x) det A(x) has mm i roots in the annulus with radii r i /4.37 and 4.37r i.
30 Some numerical experiments: Counting # of iterations Orthogonal random coefficients, degree n = 13, scaling factors: α = [1, , , 10 9, 0,... 0, 10 4, 0, 0, 1] Newton polygon Unit circle m simul it aver it simul it aver it
31 Some numerical experiments: Counting # of iterations Non-orthogonal random coefficients, degree n = 13, scaling factors: α = [1, , , 10 9, 0,... 0, 10 4, 0, 0, 1] Newton polygon Unit circle m simul it aver it simul it aver it
32 Some numerical experiments: Counting # of iterations General random coefficients, scaling factors: coefficient i, factor: (rand < 1/4) Number of average iterations with m = 10 and n = 30 Newton s Polygon 7.2 Unit circle 111
33 Some numerical experiments: Approximation error Ehrlich-Aberth: Polyeig: + Problem: Sommerfeld problem
34 Some numerical experiments: Approximation error Ehrlich-Aberth: Polyeig: + Problem: Plasma drift
35 Some numerical experiments: Approximation error Ehrlich-Aberth: Polyeig: + QuadEig (Hammarling-Munro-Tisseur): Problem: hospital
36 Some numerical experiments: Approximation error Ehrlich-Aberth: Polyeig: + QuadEig (Hammarling-Munro-Tisseur): Problem: power plant
37 Some numerical experiments: Approximation error Structured Ehrlich-Aberth: Ehrlich-Aberth: Polyeig (QZ): + Structured URV (Schroeder): Problem: butterfly
38 Some numerical experiments: Approximation error Structured Ehrlich-Aberth: Ehrlich-Aberth: Polyeig (QZ): + Structured URV (Schroeder): Problem: wiresaw1
39 Some numerical experiments: Approximation error Structured Ehrlich-Aberth: Ehrlich-Aberth: Polyeig (QZ): + Problem with symmetry (λ, λ+1 λ 1 )
40 Some references O. Aberth. Iteration methods for finding all zeros of a polynomial simultaneously. Math. Comput, D. A. Bini and G. Fiorentino. Design, analysis, and implementation of a multiprecision polynomial rootfinder. Numer. Algorithms, D.A. Bini, V. Noferini. Solving polynomial eigenvalue problems by means of the Ehrlich-Aberth method. TR L. Gemignani and V. Noferini. The Ehrlich-Aberth method for palindromic matrix polynomials represented in the Dickson basis. Linear Algebra Appl H. Guggenheimer, Initial Approximations in Durand-Kerner s Root Finding Method, BIT, 26 (1986). S. Hammarling, C. Munro and F. Tisseur. An Algorithm for the Com- plete Solution of Quadratic Eigenvalue Problems. MIMS EPrint , avaliable online. N. J. Higham, D. S. Mackey and F. Tisseur. The conditioning of linearizations of matrix polynomials. SIAM J. Matrix Anal. Appl., V. Noferini. The application of the Ehrlich-Aberth method to structured polynomial eigenvalue problems. Proc. Appl. Math. Mech., V. Noferini. Polynomial Eigenproblems: A Root-Finding Approach, Ph.D. Thesis, 2012, Dipartimento di Matematica Università di Pisa, Italy. A. Ostrowski, On a Theorem by J.L. Walsh Concerning the Moduli of Roots of Algebraic Equations, Bull. A.M.S., 47 (1941)
41 Some references A.E. Pellet, Sur un mode de séparation des racines des équations et la formule de Lagrange, Bulletin des Sciences Mathématiques, (2), vol 5 (1881), pp M. Sharify. Scaling Algorithms and Tropical Methods in Numerical Ma- trix Analysis: Application to the Optimal Assignment Problem and to the Accurate Computation of Eigenvalues. Ph.D. Thesis, Ecóle Polytechnique Paris Tech, 2011 C. Schröder. Palindromic and Even Eigenvalue Problems Analysis and Numerical Methods. PhD thesis, T.U. Berlin, P.P. Vaidyanathan, S.K. Mitra, A unified structural interpretation of some well-known stability test procedures for linear systems. Proceedings of the IEEE, vol 75, no. 4, J.L. Walsh, On Pellet s theorem concerning the roots of a polynomial, Annals of Mathematics, vol 26 (1924), pp
Tropical roots as approximations to eigenvalues of matrix polynomials. Noferini, Vanni and Sharify, Meisam and Tisseur, Francoise
Tropical roots as approximations to eigenvalues of matrix polynomials Noferini, Vanni and Sharify, Meisam and Tisseur, Francoise 2014 MIMS EPrint: 2014.16 Manchester Institute for Mathematical Sciences
More informationin Numerical Linear Algebra
Exploiting ResearchTropical MattersAlgebra in Numerical Linear Algebra February 25, 2009 Nick Françoise Higham Tisseur Director Schoolof ofresearch Mathematics The School University of Mathematics of Manchester
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Mathematics ( ) Contents lists available at SciVerse ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam Solving
More informationNick Higham. Director of Research School of Mathematics
Exploiting Research Tropical Matters Algebra in Numerical February 25, Linear 2009 Algebra Nick Higham Françoise Tisseur Director of Research School of Mathematics The School University of Mathematics
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 eigenvalue solver based on tropically scaled Lagrange linearization. Van Barel, Marc and Tisseur, Francoise. MIMS EPrint: 2016.
Polynomial eigenvalue solver based on tropically scaled Lagrange linearization Van Barel, Marc and Tisseur, Francoise 2016 MIMS EPrint: 201661 Manchester Institute for Mathematical Sciences School of Mathematics
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 informationA ROOTFINDING ALGORITHM FOR POLYNOMIALS AND SECULAR EQUATIONS
UNIVERSITÀ DEGLI STUDI DI PISA Facoltà di Scienze Matematiche, Fisiche e Naturali Corso di Laurea Magistrale in Matematica Tesi di Laurea A ROOTFINDING ALGORITHM FOR POLYNOMIALS AND SECULAR EQUATIONS Relatore:
More informationSolving the Polynomial Eigenvalue Problem by Linearization
Solving the Polynomial Eigenvalue Problem by Linearization Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with Ren-Cang Li,
More informationarxiv: v2 [math.na] 15 Feb 2013
GENERALIZATION AND VARIATIONS OF PELLET S THEOREM FOR MATRIX POLYNOMIALS arxiv:1210.0172v2 [math.na] 15 Feb 2013 A. Melman Department of Applied Mathematics School of Engineering, Santa Clara University
More informationarxiv: v1 [math.na] 21 Jan 2015
NEW SIMULTANEOUS METHODS FOR FINDING ALL ZEROS OF A POLYNOMIAL JUN-SEOP SONG A, A COLLEGE OF MEDICINE, YONSEI UNIVERSITY, SEOUL, REPUBLIC OF KOREA arxiv:1501.05033v1 [math.na] 21 Jan 2015 Abstract. The
More informationOn the reduction of matrix polynomials to Hessenberg form
Electronic Journal of Linear Algebra Volume 3 Volume 3: (26) Article 24 26 On the reduction of matrix polynomials to Hessenberg form Thomas R. Cameron Washington State University, tcameron@math.wsu.edu
More informationTHE EHRLICH-ABERTH METHOD FOR THE NONSYMMETRIC TRIDIAGONAL EIGENVALUE PROBLEM
THE EHRLICH-ABERTH METHOD FOR THE NONSYMMETRIC TRIDIAGONAL EIGENVALUE PROBLEM DARIO A. BINI, LUCA GEMIGNANI, AND FRANÇOISE TISSEUR Abstract. An algorithm based on the Ehrlich-Aberth iteration is presented
More informationStructured Matrix Methods for Polynomial Root-Finding
Structured Matrix Methods for Polynomial Root-Finding Luca Gemignani Department of Mathematics University of Pisa Largo B. Pontecorvo 5, 5617 Pisa, Italy gemignan@dm.unipi.it ABSTRACT In this paper we
More informationCertified Complex Numerical Root Finding
Intro Subdivision Numerical Comparison Goals Max Planck Institute for Informatics Dep. 1: Algorithms and Complexity Seminar on Computational Geometry and Geometric Computing May 6, 2010 Intro Subdivision
More informationScaling, Sensitivity and Stability in Numerical Solution of the Quadratic Eigenvalue Problem
Scaling, Sensitivity and Stability in Numerical Solution of the Quadratic Eigenvalue Problem Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/
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 informationUMUC MATH-107 Final Exam Information
UMUC MATH-07 Final Exam Information What should you know for the final exam? Here are some highlights of textbook material you should study in preparation for the final exam. Review this material from
More informationKU Leuven Department of Computer Science
Backward error of polynomial eigenvalue problems solved by linearization of Lagrange interpolants Piers W. Lawrence Robert M. Corless Report TW 655, September 214 KU Leuven Department of Computer Science
More informationc 1995 Society for Industrial and Applied Mathematics Vol. 37, No. 1, pp , March
SIAM REVIEW. c 1995 Society for Industrial and Applied Mathematics Vol. 37, No. 1, pp. 93 97, March 1995 008 A UNIFIED PROOF FOR THE CONVERGENCE OF JACOBI AND GAUSS-SEIDEL METHODS * ROBERTO BAGNARA Abstract.
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More information1) The line has a slope of ) The line passes through (2, 11) and. 6) r(x) = x + 4. From memory match each equation with its graph.
Review Test 2 Math 1314 Name Write an equation of the line satisfying the given conditions. Write the answer in standard form. 1) The line has a slope of - 2 7 and contains the point (3, 1). Use the point-slope
More informationMath /Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined
Math 400-001/Foundations of Algebra/Fall 2017 Numbers at the Foundations: Real Numbers In calculus, the derivative of a function f(x) is defined using limits. As a particular case, the derivative of f(x)
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis Volume 38, pp 75-30, 011 Copyright 011, ISSN 1068-9613 ETNA PERTURBATION ANALYSIS FOR COMPLEX SYMMETRIC, SKEW SYMMETRIC, EVEN AND ODD MATRIX POLYNOMIALS SK
More informationRecent Advances in the Numerical Solution of Quadratic Eigenvalue Problems
Recent Advances in the Numerical Solution of Quadratic Eigenvalue Problems Françoise Tisseur School of Mathematics The University of Manchester ftisseur@ma.man.ac.uk http://www.ma.man.ac.uk/~ftisseur/
More informationAn Algorithm for the Complete Solution of Quadratic Eigenvalue Problems. Hammarling, Sven and Munro, Christopher J. and Tisseur, Francoise
An Algorithm for the Complete Solution of Quadratic Eigenvalue Problems Hammarling, Sven and Munro, Christopher J. and Tisseur, Francoise 2011 MIMS EPrint: 2011.86 Manchester Institute for Mathematical
More informationTropical Polynomials
1 Tropical Arithmetic Tropical Polynomials Los Angeles Math Circle, May 15, 2016 Bryant Mathews, Azusa Pacific University In tropical arithmetic, we define new addition and multiplication operations on
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationGraeffe s, Chebyshev-like, and Cardinal s Processes for Splitting a Polynomial into Factors*
JOURNAL OF COMPLEXITY 12, 492 511 (1996) ARTICLE NO 0030 Graeffe s, Chebyshev-like, and Cardinal s Processes for Splitting a Polynomial into Factors* DARIO BINI Dipartimento di Matematica, Università di
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 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 informationAlgebro-geometric aspects of Heine-Stieltjes theory
Dedicated to Heinrich Eduard Heine and his 4 years old riddle Algebro-geometric aspects of Heine-Stieltjes theory Boris Shapiro, Stockholm University, shapiro@math.su.se February, 9 Introduction and main
More informationOptimal Preconditioners for Interval Gauss Seidel Methods
Optimal Preconditioners for Interval Gauss Seidel Methods R B Kearfott and Xiaofa Shi Introduction Consider the following nonlinear system f 1 x 1, x 2,, x n ) F X) f n x 1, x 2,, x n ) where bounds x
More informationImproved 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 informationthat determines x up to a complex scalar of modulus 1, in the real case ±1. Another condition to normalize x is by requesting that
Chapter 3 Newton methods 3. Linear and nonlinear eigenvalue problems When solving linear eigenvalue problems we want to find values λ C such that λi A is singular. Here A F n n is a given real or complex
More informationSection 0.2 & 0.3 Worksheet. Types of Functions
MATH 1142 NAME Section 0.2 & 0.3 Worksheet Types of Functions Now that we have discussed what functions are and some of their characteristics, we will explore different types of functions. Section 0.2
More informationNeville s Method. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Neville s Method
Neville s Method MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have learned how to approximate a function using Lagrange polynomials and how to estimate
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
More informationEigenvector error bound and perturbation for nonlinear eigenvalue problems
Eigenvector error bound and perturbation for nonlinear eigenvalue problems Yuji Nakatsukasa School of Mathematics University of Tokyo Joint work with Françoise Tisseur Workshop on Nonlinear Eigenvalue
More informationThe absolute value (modulus) of a number
The absolute value (modulus) of a number Given a real number x, its absolute value or modulus is dened as x if x is positive x = 0 if x = 0 x if x is negative For example, 6 = 6, 10 = ( 10) = 10. The absolute
More informationFactoring block Fiedler Companion Matrices
Chapter Factoring block Fiedler Companion Matrices Gianna M. Del Corso, Federico Poloni, Leonardo Robol, and Raf Vandebril Abstract When Fiedler published his A note on Companion matrices in 00 in Linear
More informationHomework 8 Solutions to Selected Problems
Homework 8 Solutions to Selected Problems June 7, 01 1 Chapter 17, Problem Let f(x D[x] and suppose f(x is reducible in D[x]. That is, there exist polynomials g(x and h(x in D[x] such that g(x and h(x
More informationThe quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying
I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem QEP is to find scalars λ and nonzero vectors u satisfying where Qλx = 0, 1.1 Qλ = λ 2 M + λd + K, M, D and K are given
More informationAre resultant methods numerically unstable for multidimensional rootfinding? Alex Townsend MIT
Are resultant methods numerically unstable for multidimensional rootfinding? University of Oxford, 16th June 2015 NA internal seminar Alex Townsend MIT Joint work with Vanni Noferini Betteridge s law of
More information1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
Math120 - Precalculus. Final Review. Fall, 2011 Prepared by Dr. P. Babaali 1 Algebra 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
More informationLinear algebra issues in Interior Point methods for bound-constrained least-squares problems
Linear algebra issues in Interior Point methods for bound-constrained least-squares problems Stefania Bellavia Dipartimento di Energetica S. Stecco Università degli Studi di Firenze Joint work with Jacek
More informationLecture 10 - Eigenvalues problem
Lecture 10 - Eigenvalues problem Department of Computer Science University of Houston February 28, 2008 1 Lecture 10 - Eigenvalues problem Introduction Eigenvalue problems form an important class of problems
More informationAdvanced Algebra II 1 st Semester Exam Review
dname Advanced Algebra II 1 st Semester Exam Review Chapter 1A: Number Sets & Solving Equations Name the sets of numbers to which each number belongs. 1. 34 2. 525 3. 0.875 4. Solve each equation. Check
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
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 informationSQUARE-ROOT FAMILIES FOR THE SIMULTANEOUS APPROXIMATION OF POLYNOMIAL MULTIPLE ZEROS
Novi Sad J. Math. Vol. 35, No. 1, 2005, 59-70 SQUARE-ROOT FAMILIES FOR THE SIMULTANEOUS APPROXIMATION OF POLYNOMIAL MULTIPLE ZEROS Lidija Z. Rančić 1, Miodrag S. Petković 2 Abstract. One-parameter families
More informationUp to this point, our main theoretical tools for finding eigenvalues without using det{a λi} = 0 have been the trace and determinant formulas
Finding Eigenvalues Up to this point, our main theoretical tools for finding eigenvalues without using det{a λi} = 0 have been the trace and determinant formulas plus the facts that det{a} = λ λ λ n, Tr{A}
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 informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationA PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
IJMMS 25:6 2001) 397 409 PII. S0161171201002290 http://ijmms.hindawi.com Hindawi Publishing Corp. A PROJECTED HESSIAN GAUSS-NEWTON ALGORITHM FOR SOLVING SYSTEMS OF NONLINEAR EQUATIONS AND INEQUALITIES
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 Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationMathematical Olympiad Training Polynomials
Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,
More information1. The profit from selling local ballet tickets depends on the ticket price. Using past receipts, we find that the profit can be modeled
Warm-up 1. The profit from selling local ballet tickets depends on the ticket price. Using past receipts, we find that the profit can be modeled by the function p x = 15x 2 + 600x + 60, where x is the
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #2 Instructions: This is an open notes exam, i.e., you are allowed to consult any textbook, your class notes, homeworks, or any of the handouts from us.
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationarxiv: v1 [cs.it] 12 Jun 2016
New Permutation Trinomials From Niho Exponents over Finite Fields with Even Characteristic arxiv:606.03768v [cs.it] 2 Jun 206 Nian Li and Tor Helleseth Abstract In this paper, a class of permutation trinomials
More informationLecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 15 Review of Matrix Theory III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Matrix An m n matrix is a rectangular or square array of
More information1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents.
Math120 - Precalculus. Final Review Prepared by Dr. P. Babaali 1 Algebra 1. Use the properties of exponents to simplify the following expression, writing your answer with only positive exponents. (a) 5
More informationA numerical method for polynomial eigenvalue problems using contour integral
A numerical method for polynomial eigenvalue problems using contour integral Junko Asakura a Tetsuya Sakurai b Hiroto Tadano b Tsutomu Ikegami c Kinji Kimura d a Graduate School of Systems and Information
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationComplex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i
Complex Numbers: Definition: A complex number is a number of the form: z = a + bi where a, b are real numbers and i is a symbol with the property: i 2 = 1 Sometimes we like to think of i = 1 We can treat
More informationNotes on generating functions in automata theory
Notes on generating functions in automata theory Benjamin Steinberg December 5, 2009 Contents Introduction: Calculus can count 2 Formal power series 5 3 Rational power series 9 3. Rational power series
More informationOn Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials. University of Bucharest CADE 2007.
On Computation of Positive Roots of Polynomials and Applications to Orthogonal Polynomials Doru Ştefănescu University of Bucharest CADE 2007 21 February 2007 Contents Approximation of the Real Roots Bounds
More informationNumerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity
ANZIAM J. 46 (E) ppc46 C438, 005 C46 Numerical techniques for linear and nonlinear eigenvalue problems in the theory of elasticity Aliki D. Muradova (Received 9 November 004, revised April 005) Abstract
More informationQuadratic Matrix Polynomials
Research Triangularization Matters of Quadratic Matrix Polynomials February 25, 2009 Nick Françoise Higham Tisseur Director School of of Research Mathematics The University of Manchester School of Mathematics
More informationAlgebra II Vocabulary Alphabetical Listing. Absolute Maximum: The highest point over the entire domain of a function or relation.
Algebra II Vocabulary Alphabetical Listing Absolute Maximum: The highest point over the entire domain of a function or relation. Absolute Minimum: The lowest point over the entire domain of a function
More informationMATHEMATICAL ENGINEERING TECHNICAL REPORTS. Eigenvector Error Bound and Perturbation for Polynomial and Rational Eigenvalue Problems
MATHEMATICAL ENGINEERING TECHNICAL REPORTS Eigenvector Error Bound and Perturbation for Polynomial and Rational Eigenvalue Problems Yuji NAKATSUKASA and Françoise TISSEUR METR 2016 04 April 2016 DEPARTMENT
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationApproximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method
Approximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method Antti Koskela KTH Royal Institute of Technology, Lindstedtvägen 25, 10044 Stockholm,
More informationHow to Detect Definite Hermitian Pairs
How to Detect Definite Hermitian Pairs Françoise Tisseur School of Mathematics The University of Manchester ftisseur@ma.man.ac.uk http://www.ma.man.ac.uk/~ftisseur/ Joint work with Chun-Hua Guo and Nick
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationAn LMI description for the cone of Lorentz-positive maps II
An LMI description for the cone of Lorentz-positive maps II Roland Hildebrand October 6, 2008 Abstract Let L n be the n-dimensional second order cone. A linear map from R m to R n is called positive if
More informationTest 2 Review Math 1111 College Algebra
Test 2 Review Math 1111 College Algebra 1. Begin by graphing the standard quadratic function f(x) = x 2. Then use transformations of this graph to graph the given function. g(x) = x 2 + 2 *a. b. c. d.
More informationFirst, we review some important facts on the location of eigenvalues of matrices.
BLOCK NORMAL MATRICES AND GERSHGORIN-TYPE DISCS JAKUB KIERZKOWSKI AND ALICJA SMOKTUNOWICZ Abstract The block analogues of the theorems on inclusion regions for the eigenvalues of normal matrices are given
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 informationPolynomial Jacobi Davidson Method for Large/Sparse Eigenvalue Problems
Polynomial Jacobi Davidson Method for Large/Sparse Eigenvalue Problems Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan April 28, 2011 T.M. Huang (Taiwan Normal Univ.)
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 informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT -09 Computational and Sensitivity Aspects of Eigenvalue-Based Methods for the Large-Scale Trust-Region Subproblem Marielba Rojas, Bjørn H. Fotland, and Trond Steihaug
More informationInverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with J-(Skew) Centrosymmetry
Inverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with -(Skew) Centrosymmetry Zhong-Yun Liu 1 You-Cai Duan 1 Yun-Feng Lai 1 Yu-Lin Zhang 1 School of Math., Changsha University
More informationTHE SUM OF DIGITS OF n AND n 2
THE SUM OF DIGITS OF n AND n 2 KEVIN G. HARE, SHANTA LAISHRAM, AND THOMAS STOLL Abstract. Let s q (n) denote the sum of the digits in the q-ary expansion of an integer n. In 2005, Melfi examined the structure
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
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 informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 10-12 Large-Scale Eigenvalue Problems in Trust-Region Calculations Marielba Rojas, Bjørn H. Fotland, and Trond Steihaug ISSN 1389-6520 Reports of the Department of
More informationMatrix Theory. A.Holst, V.Ufnarovski
Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different
More informationCenterville High School Curriculum Mapping Algebra II 1 st Nine Weeks
Centerville High School Curriculum Mapping Algebra II 1 st Nine Weeks Chapter/ Lesson Common Core Standard(s) 1-1 SMP1 1. How do you use a number line to graph and order real numbers? 2. How do you identify
More informationChebyshev coordinates and Salem numbers
Chebyshev coordinates and Salem numbers S.Capparelli and A. Del Fra arxiv:181.11869v1 [math.co] 31 Dec 018 January 1, 019 Abstract By expressing polynomials in the basis of Chebyshev polynomials, certain
More informationPolynomials. Math 4800/6080 Project Course
Polynomials. Math 4800/6080 Project Course 2. Algebraic Curves. Everyone knows what a curve is, until he has studied enough mathematics to become confused through the countless number of possible exceptions.
More informationCURRICULUM MAP. Course/Subject: Honors Math I Grade: 10 Teacher: Davis. Month: September (19 instructional days)
Month: September (19 instructional days) Numbers, Number Systems and Number Relationships Standard 2.1.11.A: Use operations (e.g., opposite, reciprocal, absolute value, raising to a power, finding roots,
More informationOperators with numerical range in a closed halfplane
Operators with numerical range in a closed halfplane Wai-Shun Cheung 1 Department of Mathematics, University of Hong Kong, Hong Kong, P. R. China. wshun@graduate.hku.hk Chi-Kwong Li 2 Department of Mathematics,
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationA CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION
STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, September 2003 A CLASS OF EVEN DEGREE SPLINES OBTAINED THROUGH A MINIMUM CONDITION GH. MICULA, E. SANTI, AND M. G. CIMORONI Dedicated to
More information