Roberto Castelli Basque Center for Applied Mathematics
|
|
- Octavia Shields
- 6 years ago
- Views:
Transcription
1 A rigorous computational method to enclose the eigendecomposition of interval matrices Roberto Castelli Basque Center for Applied Mathematics Joint work with: Jean-Philippe Lessard Basque / Hungarian Workshop on Numerical Methods for Large Systems
2 Why Interval entries? Uncertainty of measurement - inaccuracy of data collection To consider variability of natural phenomena Qualitative indicators with statistical ranges Errors in computations...
3 Computation of bundles of periodic orbits ẏ = g(y) γ(t + T )=γ(t) γ 3T/4 (0) γ T/2 (0) Γ Γ = {γ(t) :t [0,T)} γ 0 (0) γ θ (t) =γ(t + θ) γ T/4 (0)
4 Computation of bundles of periodic orbits Consider γ(t + T )=γ(t) Γ = {γ(t) :t [0,T)} γ θ (t) =γ(t + θ) Φ θ (t) ẏ = g(y) γ 0 (0) Linearize around the periodic orbit: A θ (t) = g(γ θ (t)) Y = A θ (t)y Y (0) = I γ 3T/4 (0) γ T/4 (0) y = A θ (t)y : T- periodic matrix function γ T/2 (0) a fundamental matrix solution of the non-autonomous linear system Definition: The monodromy matrix is defined to be Φ θ (T ). Definition: The eigenvalues of Φ θ (T ) are called the Floquet multipliers. dim W s (Γ) = number of Floquet multipliers with modulus less than one. dim W u (Γ) = number of Floquet multipliers with modulus greater than one. Γ
5 Computation of bundles of periodic orbits
6 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. Ingredients
7 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions Ingredients
8 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions Ingredients 2. Spectral methods, choice of suitable Banach spaces
9 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions Ingredients 2. Spectral methods, choice of suitable Banach spaces 3. Finite dim. projection, bounds on the truncation error terms (Analytical estimates)
10 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions Ingredients 2. Spectral methods, choice of suitable Banach spaces 3. Finite dim. projection, bounds on the truncation error terms (Analytical estimates) 4. Theorems: Banach Fixed Point Theorem, Conley Index...
11 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions Ingredients 2. Spectral methods, choice of suitable Banach spaces 3. Finite dim. projection, bounds on the truncation error terms (Analytical estimates) 4. Theorems: Banach Fixed Point Theorem, Conley Index Numerical analysis (Newton s method, Fast Fourier transform)
12 Rigorous Computations The goal of rigorous computations is to construct algorithms that provide an approximate solution to the problem f(x)=0 together with precise and possibly efficient bounds within which the exact solution is guaranteed to exist in the mathematically rigorous sense. 1. Smoothness of the solutions 2. Spectral methods, choice of suitable Banach spaces 3. Finite dim. projection, bounds on the truncation error terms (Analytical estimates) 4. Theorems: Banach Fixed Point Theorem, Conley Index Numerical analysis (Newton s method, Fast Fourier transform) 6. Interval Arithmetic Ingredients
13 Banach Fixed Point Argument Topological Arguments Rigorous integration Some Techniques - Sarah Day, Jean-Philippe Lessard, and Konstantin Mischaikow. Validated continuation for equilibria of PDEs. SIAM J. Numer. Anal., 45(4): , Marcio Gameiro, Jean-Philippe Lessard, and Konstantin Mischaikow. Validated continuation over large parameter ranges for equilibria of PDEs. Math. Comput. Simulation, 79(4): , Roberto Castelli and Jean-Philippe Lessard. Rigorous numerics in Floquet theory: computing stable and unstable bundles of periodic orbits. Submitted., Piotr Zgliczynski and Konstantin Mischaikow. Rigorous numerics for partial differential equations: the Kuramoto-Sivashinsky equation. Found. Comput. Math., 1(3): , Gianni Arioli, Piotr Zgliczynski, Symbolic dynamics for the Henon-Heiles Hamiltonian on the critical level. J. Diff. Equations 171, , Piotr Zgliczynski. C1 -Lohner algorithm. Found. Comput. Math., 2(4): , Tomasz Kapela and Carles Simo. Computer assisted proofs for nonsymmetric planar choreographies and for stability of the Eight. Nonlinearity, 20(5):1241{1255, 2007.
14 Rigorous Computations: Banach Fixed Point Theorem Let (B, B ) be a Banach space, that is a complete normed vector space. Assume the existence of a contraction mapping T on B, that is a mapping such that 1. T (B) B; 2. there exists 0 < κ < 1 such that, for every x, y B, T (x) T (y) B κx y B. Then there exists a unique x B such that T (x )=x.
15 Rigorous Computations: Banach Fixed Point Theorem Let (B, B ) be a Banach space, that is a complete normed vector space. Assume the existence of a contraction mapping T on B, that is a mapping such that 1. T (B) B; 2. there exists 0 < κ < 1 such that, for every x, y B, T (x) T (y) B κx y B. Then there exists a unique x B such that T (x )=x. can be verified by strict inequalities
16 Rigorous Computations: Banach Fixed Point Theorem Let (B, B ) be a Banach space, that is a complete normed vector space. Assume the existence of a contraction mapping T on B, that is a mapping such that 1. T (B) B; 2. there exists 0 < κ < 1 such that, for every x, y B, T (x) T (y) B κx y B. Then there exists a unique x B such that T (x )=x. can be verified by strict inequalities can be verified using interval arithmetic!
17 Rigorous computations: Formulation in Banach Space
18 Rigorous computations: Formulation in Banach Space F (u, ν) =0 (E.g : Differential Equations) Knowledge about regularity x Ωs = f(x, Spectral methods ν) =0 x ν : modes : parameter { (x k ) k sup ωk x s k < k }
19 Rigorous computations: Formulation in Banach Space F (u, ν) =0 (E.g : Differential Equations) Knowledge about regularity x Ωs = f(x, Spectral methods ν) =0 x ν : modes : parameter { (x k ) k sup ωk x s k < k } x Consider such that f (m) ( x, ν 0 ) 0. Finite dim. reduction f(x, ν) =0 T ν (x) =x Newton-like operator at x
20 Rigorous computations: Formulation in Banach Space F (u, ν) =0 (E.g : Differential Equations) Knowledge about regularity x Ωs = f(x, Spectral methods ν) =0 x ν : modes : parameter { (x k ) k sup ωk x s k < k } x Consider such that f (m) ( x, ν 0 ) 0. Finite dim. reduction f(x, ν) =0 T ν (x) =x 8) Newton-like operator at x T ν :Ω s Ω s T ν (x) =x Jf(x, ν) J D x f( x, ν 0 ) 1 The chances of contracting a small set B around x depends on the magnitude of the eigenvalues of. J
21 Q: How to find a ball B x (r) such that T : B x (r) B x (r) is a contraction? ν ˆx x B x (r) Ω s
22 Q: How to find a ball B x (r) such that T : B x (r) B x (r) is a contraction? ν B x (r) = x + B(r) ˆx x B x (r) Ω s
23 Q: How to find a ball B x (r) such that T : B x (r) B x (r) is a contraction? ν B x (r) = x + B(r) Ball of radius r centered at 0 in the space Ω s ˆx x B x (r) Ω s
24 Q: How to find a ball B x (r) such that T : B x (r) B x (r) is a contraction? ν B x (r) = x + B(r) Ball of radius r centered at 0 in the space Ω s ˆx x B x (r) Ω s A: Radii polynomials Suppose the bounds Y, Z(r) satisfy D x T ν ( x + b)c T ν ( x) x Y Z(r) k k sup b,c B(r) p k (r) =Y + Z(r) r w s k
25 Q: How to find a ball B x (r) such that T : B x (r) B x (r) is a contraction? ν B x (r) = x + B(r) Ball of radius r centered at 0 in the space Ω s ˆx x B x (r) Ω s A: Radii polynomials Suppose the bounds Y, Z(r) satisfy D x T ν ( x + b)c T ν ( x) x Y Z(r) k k Lemma: If there exists r>0such that p k (r) < 0 for all k, then there is a unique ˆx B x (r) s.t. f(ˆx, ν) =0. proof. Banach fixed point theorem. sup b,c B(r) p k (r) =Y + Z(r) r w s k
26 Rigorous Computations An Example : Compute solution of x 2 2=0
27 Rigorous Computations An Example : Compute solution of x 2 2=0
28 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Choose the Banach Space Compute num. solution Define the fixed point operator x Define Define radii polynomial Find r sup b,c B(r) T ( x) x Y k D x T ( x + b)c Z(r) k
29 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Choose the Banach Space Compute num. solution Define the fixed point operator x Define Define radii polynomial Find r sup b,c B(r) T ( x) x Y k D x T ( x + b)c Z(r) k
30 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
31 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
32 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) T ( x) x DT( x + b)c r +0.7 r 2 Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
33 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) T ( x) x DT( x + b)c r +0.7 r 2 Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
34 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) T ( x) x DT( x + b)c r +0.7 r 2 p(r) =0.7 r r Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
35 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) T ( x) x DT( x + b)c r +0.7 r 2 p(r) =0.7 r r Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k
36 Rigorous Computations An Example : Compute solution of x 2 2=0 Ω = R x =1.41 Df( x) T (x) =x 0.35(x 2 2) T ( x) x DT( x + b)c r +0.7 r 2 p(r) =0.7 r r I =[ , ], r I,p(r) < 0 Choose the Banach Space Compute num. solution Define the fixed point operator Define Define radii polynomial Find r sup b,c B(r) x T ( x) x Y k D x T ( x + b)c Z(r) k ±
37 Eigendecomposition of interval matrices Notation: bold case interval normal case scalar A A A i,j =[a i,j a i,j ], Â is the center of A
38 Eigendecomposition of interval matrices Notation: bold case interval normal case scalar A A A i,j =[a i,j a i,j ], Â is the center of A GOAL: enclosing eigendecomposition Given A IC n n construct a set of balls B i C C n, i =1...n such that for any A A and any i =1...n there exists x i =(λ i,v i ) B i s.t Av i = λ i v i, where x i is unique up to a scaling factor of v i.
39 Some literature Tetsuro Yamamoto. Error bounds for computed eigenvalues and eigenvectors. Numer. Math., 34 (2): , Siegfried M. Rump and Jens-Peter M. Zemke. On eigenvector bounds. BIT, 43(4): ,2003. Milan Hladik, David Daney, and Elias Tsigaridas. Bounds on real eigenvalues and singular values of interval matrices. SIAM J. Matrix Anal. Appl., 31(4): , 2009/10. G. Mayer. Result verification for eigenvectors and eigenvalues. In Topics in validated computations (Oldenburg, 1993), volume 5 of Stud. Comput. Math., pages North- Holland, Amsterdam, G. Alefeld and H. Spreuer. Iterative improvement of componentwise error bounds for invariant subspaces belonging to a double or nearly double eigenvalue. Computing, 36(4): , Siegfried M. Rump. Computational error bounds for multiple or nearly multiple eigenvalues. Linear Algebra Appl., 324(1-3): , Special issue on linear algebra in self-validating methods.
40 Rigorous computation of eigendecomposition: scalar case Define the problem f(x)=0
41 Rigorous computation of eigendecomposition: scalar case Define the problem f(x)=0 Suppose ( λ, v) is an approximate eigenpair of A, i.e.a v λ v f(x) =A x =(λ,v 1,v 2,...,v k 1,v k+1,...,v n ) v 1. v ḳ. v n f(x) :C n C n λ v 1. v ḳ. v n, v k = max{ v i } f(x) =0 (λ,v) is an eigenpair of A.
42 Rigorous computation of eigendecomposition: scalar case Define the problem f(x)=0 Suppose ( λ, v) is an approximate eigenpair of A, i.e.a v λ v f(x) =A x =(λ,v 1,v 2,...,v k 1,v k+1,...,v n ) v 1. v ḳ. v n f(x) :C n C n λ v 1. v ḳ. v n, v k = max{ v i } f(x) =0 (λ,v) is an eigenpair of A. Note that the roots of f(x) are isolated
43 Rigorous computation of eigendecomposition: scalar case Numerical solution : x =( λ, v 1, v 2,..., v k 1, v k+1,..., v n ) Banach space: Ω = {x C n : x Ω < }, x Ω = max{ x i }, i Operator T(x): T (x) =x Rf(x) v 1. R Df( x) 1, Df( x) = v ḳ (A I n )ˆk.. v n
44 Rigorous computation of eigendecomposition: scalar case Definition of the bounds Y, Z(r) ( T(x)=x-Rf(x))
45 Rigorous computation of eigendecomposition: scalar case Definition of the bounds Y, Z(r) Y T ( x) x Y := Rf( x) ( T(x)=x-Rf(x))
46 Rigorous computation of eigendecomposition: scalar case Definition of the bounds Y, Z(r) Y T ( x) x Y := Rf( x) Z(r) sup b,c B(r) DT( x + b)c ( T(x)=x-Rf(x)) DT( x + b)c (I n R Df( x))c + R[(Df( x) Df( x + b))c] r I n R Df( x) 1 n + r 2 2 R ˆ1 n ( b, c r) Z(r) :=rz 0 + r 2 Z 1
47 Rigorous computation of eigendecomposition: scalar case Definition of the bounds Y, Z(r) Y T ( x) x Y := Rf( x) Z(r) sup b,c B(r) DT( x + b)c ( T(x)=x-Rf(x)) DT( x + b)c (I n R Df( x))c + R[(Df( x) Df( x + b))c] r I n R Df( x) 1 n + r 2 2 R ˆ1 n ( b, c r) Z(r) :=rz 0 + r 2 Z 1 p k (r) =(Y + Z(r)) k r.
48 Rigorous computation of eigendecomposition: scalar case Definition of the bounds Y, Z(r) Y T ( x) x Y := Rf( x) Z(r) sup b,c B(r) DT( x + b)c ( T(x)=x-Rf(x)) DT( x + b)c (I n R Df( x))c + R[(Df( x) Df( x + b))c] r I n R Df( x) 1 n + r 2 2 R ˆ1 n ( b, c r) Z(r) :=rz 0 + r 2 Z 1 p k (r) =(Y + Z(r)) k r. Suppose that for any r I, p k (r) < 0 for any k. Then, for any r I, thereexistsagenuine eigenpair x =(λ,v) of A in the ball B x (r) around ( λ, v)
49 Rigorous computation of eigendecomposition: extension to interval matrix Consider interval matrix A and interval arithmetics
50 Rigorous computation of eigendecomposition: extension to interval matrix Consider interval matrix A and interval arithmetics A IC n n, f : C n IC n, Df IC n n
51 Rigorous computation of eigendecomposition: extension to interval matrix Consider interval matrix A and interval arithmetics A IC n n, f : C n IC n, Df IC n n Define x =( λ, v) a numerical eigenpair of  Let R be a numerical inverse of ˆDf and define T = x Rf(x)
52 Rigorous computation of eigendecomposition: extension to interval matrix Consider interval matrix A and interval arithmetics A IC n n, f : C n IC n, Df IC n n Define x =( λ, v) a numerical eigenpair of  Let R be a numerical inverse of ˆDf and define T = x Rf(x) p k (r) < 0 for any A A there exists a unique (λ,v) so that λ λ <r, v j v j <r, v k = v k and Av = λv.
53 The Matlab Code 1. function r=enc_eig(a,v,l) 2. n= size(v,1); 3. iv=intval(v); 4. [m,iv]=max(abs(v)); 5. O=ones(n,1); 6. df=intval(a)-intval(l)*speye(n); 7. F=df*iV; 8. df(:,iv)=-iv; 9. R=inv(mid(df)); 10. Y=sup(abs(R*F)); 11. Z0=sup(abs(eye(n)-R*df))*O; 12. O(IV)=0; 13. Z1=2*abs(R)*O; 14. rd=[z1,z0-ones(n,1),y]; % coefficients of radii polynomials [r^2,r,const] 15. r=evaluate_intersection_int(rd); % evaluate where all the polynomials are negative 16. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 17. function neg_inter=evaluate_intersection_int(rd) 18. i2=intval(2); 19. rd=intval(rd); 20. neg_inter=nan; 21. Rm=(-rd(:,2)-sqrt(rd(:,2).^i2-intval(4)*rd(:,1).*rd(:,3)))./(i2*rd(:,1)); 22. Rp=-Rm-rd(:,2)./rd(:,1); 23. ISCOMP=find(isreal(Rp)+isreal(Rm)<2, 1); 24. if isempty(iscomp)==0 25. 'error: a polynomial is positive everywhere' 26. return 27. end 28. intersection=[max(rm.sup),min(rp.inf)]; 29. if intersection(1)<intersection(2) 30. neg_inter=intersection(1); 31. end 32. end 33. end 34. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 35. Definition of radii polynomilas Compute intersections of intervals
54 Some computational results Enclosure of one eigenpair for large matrices D = diag( N/2,...,N/2), X = rand(size(d)), A = X D X 1 [q, w] =eig(a) x =(q(:,k),w(k, k)) N r time (s) e e e e
55 Some computational results Enclosure of the Eigendecomposition for large intervals D = diag(0,e ik 2π 5 ), X = rand(size(d)) + i rand(size(d)) Â = X D X 1, A = Â + 1rad
56 Some computational results Comparison of the computational time verifyeig - verification method, Brouwer s fix p. verifynlss - Krawczyk operator for the computation of zeros of nonlinear functions  is a matrix with eigenvalue 0 and N equispaced on the unit circle. On the left rad = 10 15
57 Some computational results Comparison of the computational time Our approach is satisfactory from the point of view of accuracy of the result and the algorithm is faster then verifyeig, especially for small N. This achievement has been possible thanks to the analytical estimates, in the form of radii polynomials, that allow minimizing the number of computations done with intervals.
58 Thank you Köszönöm - Eskerrik asko for your attention
Rigorous verification of saddle-node bifurcations in ODEs
Rigorous verification of saddle-node bifurcations in ODEs Jean-Philippe Lessard Abstract In this paper, we introduce a general method for the rigorous verification of saddlenode bifurcations in ordinary
More informationComputation of homoclinic and heteroclinic orbits for flows
Computation of homoclinic and heteroclinic orbits for flows Jean-Philippe Lessard BCAM BCAM Mini-symposium on Computational Math February 1st, 2011 Rigorous Computations Connecting Orbits Compute a set
More informationDynamic Systems and Applications 11 (2002) TOWARDS A RIGOROUS NUMERICAL STUDY OF THE KOT SCHAFFER MODEL
Dynamic Systems and Applications (2002) 87-97 TOWARDS A RIGOROUS NUMERICAL STUDY OF THE KOT SCHAFFER MODEL SARAH DAY Georgia Institute of Technology, Department of Mathematics, Atlanta, GA 30332 090, USA.
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationExistence of Secondary Bifurcations or Isolas for PDEs
Existence of Secondary Bifurcations or Isolas for PDEs Marcio Gameiro Jean-Philippe Lessard Abstract In this paper, we introduce a method to conclude about the existence of secondary bifurcations or isolas
More informationRigorous numerics for nonlinear operators with tridiagonal dominant linear parts
Rigorous numerics for nonlinear operators with tridiagonal dominant linear parts Maxime Breden Laurent Desvillettes Jean-Philippe Lessard Abstract We propose a method to compute solutions of infinite dimensional
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationIn particular, if A is a square matrix and λ is one of its eigenvalues, then we can find a non-zero column vector X with
Appendix: Matrix Estimates and the Perron-Frobenius Theorem. This Appendix will first present some well known estimates. For any m n matrix A = [a ij ] over the real or complex numbers, it will be convenient
More informationValidated numerics for equilibria of analytic vector fields: invariant manifolds and connecting orbits
Validated numerics for equilibria of analytic vector fields: invariant manifolds and connecting orbits J.D. Mireles James, May 24, 2017 Abstract The goal of these notes is to illustrate the use of validated
More informationRapidly and slowly oscillating periodic solutions of a delayed van der Pol oscillator
Rapidly and slowly oscillating periodic solutions of a delayed van der Pol oscillator Gabor Kiss Jean-Philippe Lessard Abstract In this paper, we introduce a method to prove existence of several rapidly
More informationSpectral Theorem for Self-adjoint Linear Operators
Notes for the undergraduate lecture by David Adams. (These are the notes I would write if I was teaching a course on this topic. I have included more material than I will cover in the 45 minute lecture;
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationRigorous Numerics for Global Dynamics: a study of the Swift-Hohenberg equation
Rigorous Numerics for Global Dynamics: a study of the Swift-Hohenberg equation Sarah Day Yasuaki Hiraoka Konstantin Mischaikow Toshiyuki Ogawa February 24, 2004 Abstract This paper presents a rigorous
More informationValidated numerics for continuation and bifurcation of connecting orbits of maps
Validated numerics for continuation and bifurcation of connecting orbits of maps Ronald Adams 1 and JD Mireles James 2 1 Embry-Riddle Aeronautical University, Department of Mathematics 2 Florida Atlantic
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationValidated computation tool for the Perron-Frobenius eigenvalue using graph theory
Validated computation tool for the Perron-Frobenius eigenvalue using graph theory Kaori Nagatou Kyushu University, JAPAN [ Joint work with Yutaka Ishii, Kyushu University ] AGA Seminar in INI, March 13,
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 informationTHE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR
THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR WEN LI AND MICHAEL K. NG Abstract. In this paper, we study the perturbation bound for the spectral radius of an m th - order n-dimensional
More informationLinear Algebra - Part II
Linear Algebra - Part II Projection, Eigendecomposition, SVD (Adapted from Sargur Srihari s slides) Brief Review from Part 1 Symmetric Matrix: A = A T Orthogonal Matrix: A T A = AA T = I and A 1 = A T
More informationUsing the Karush-Kuhn-Tucker Conditions to Analyze the Convergence Rate of Preconditioned Eigenvalue Solvers
Using the Karush-Kuhn-Tucker Conditions to Analyze the Convergence Rate of Preconditioned Eigenvalue Solvers Merico Argentati University of Colorado Denver Joint work with Andrew V. Knyazev, Klaus Neymeyr
More informationRigorous Numerics in Dynamics
Rigorous Numerics in Dynamics Jan Bouwe van den Berg and Jean-Philippe Lessard Motivation Nonlinear dynamics shape the world around us, from the harmonious movements of celestial bodies, via the swirling
More informationTOWARDS AUTOMATED CHAOS VERIFICATION
TOWARDS AUTOMATED CHAOS VERIFICATION SARAH DAY CDSNS, Georgia Institute of Technology, Atlanta, GA 30332 OLIVER JUNGE Institute for Mathematics, University of Paderborn, 33100 Paderborn, Germany KONSTANTIN
More informationOPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY
published in IMA Journal of Numerical Analysis (IMAJNA), Vol. 23, 1-9, 23. OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY SIEGFRIED M. RUMP Abstract. In this note we give lower
More informationNORMS ON SPACE OF MATRICES
NORMS ON SPACE OF MATRICES. Operator Norms on Space of linear maps Let A be an n n real matrix and x 0 be a vector in R n. We would like to use the Picard iteration method to solve for the following system
More informationOn intermediate value theorem in ordered Banach spaces for noncompact and discontinuous mappings
Int. J. Nonlinear Anal. Appl. 7 (2016) No. 1, 295-300 ISSN: 2008-6822 (electronic) http://dx.doi.org/10.22075/ijnaa.2015.341 On intermediate value theorem in ordered Banach spaces for noncompact and discontinuous
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationBounds on the radius of nonsingularity
Towards bounds on the radius of nonsingularity, Milan Hladík Department of Applied Mathematics, Charles University, Prague and Institute of Computer Science of the Academy of Czech Republic SCAN 2014,
More informationAsymptotic Stability by Linearization
Dynamical Systems Prof. J. Rauch Asymptotic Stability by Linearization Summary. Sufficient and nearly sharp sufficient conditions for asymptotic stability of equiiibria of differential equations, fixed
More informationON FRACTAL DIMENSION OF INVARIANT SETS
ON FRACTAL DIMENSION OF INVARIANT SETS R. MIRZAIE We give an upper bound for the box dimension of an invariant set of a differentiable function f : U M. Here U is an open subset of a Riemannian manifold
More informationThe University of Texas at Austin Department of Electrical and Computer Engineering. EE381V: Large Scale Learning Spring 2013.
The University of Texas at Austin Department of Electrical and Computer Engineering EE381V: Large Scale Learning Spring 2013 Assignment Two Caramanis/Sanghavi Due: Tuesday, Feb. 19, 2013. Computational
More informationPersonal notes on renormalization
Personal notes on renormalization Pau Rabassa Sans April 17, 2009 1 Dynamics of the logistic map This will be a fast (and selective) review of the dynamics of the logistic map. Let us consider the logistic
More informationVALIDATED COMPUTATION OF HETEROCLINIC SETS
VALIDATED COMPUTATION OF HETEROCLINIC SETS MACIEJ J. CAPIŃSKI AND J.D. MIRELES JAMES Abstract. In this work we develop a method for computing mathematically rigorous enclosures of some one dimensional
More informationThroughout these notes we assume V, W are finite dimensional inner product spaces over C.
Math 342 - Linear Algebra II Notes Throughout these notes we assume V, W are finite dimensional inner product spaces over C 1 Upper Triangular Representation Proposition: Let T L(V ) There exists an orthonormal
More informationFor δa E, this motivates the definition of the Bauer-Skeel condition number ([2], [3], [14], [15])
LAA 278, pp.2-32, 998 STRUCTURED PERTURBATIONS AND SYMMETRIC MATRICES SIEGFRIED M. RUMP Abstract. For a given n by n matrix the ratio between the componentwise distance to the nearest singular matrix and
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.2: Fundamentals 2 / 31 Eigenvalues and Eigenvectors Eigenvalues and eigenvectors of
More informationCOMPUTATIONAL ERROR BOUNDS FOR MULTIPLE OR NEARLY MULTIPLE EIGENVALUES
LAA 324, Issue 1-3, pp. 209-226, 2001 COMPUTATIONAL ERROR BOUNDS FOR MULTIPLE OR NEARLY MULTIPLE EIGENVALUES SIEGFRIED M. RUMP Abstract. In this paper bounds for clusters of eigenvalues of non-selfadjoint
More informationarxiv: v1 [math.na] 28 Jun 2013
A New Operator and Method for Solving Interval Linear Equations Milan Hladík July 1, 2013 arxiv:13066739v1 [mathna] 28 Jun 2013 Abstract We deal with interval linear systems of equations We present a new
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationNonlinear Programming Algorithms Handout
Nonlinear Programming Algorithms Handout Michael C. Ferris Computer Sciences Department University of Wisconsin Madison, Wisconsin 5376 September 9 1 Eigenvalues The eigenvalues of a matrix A C n n are
More informationStat 159/259: Linear Algebra Notes
Stat 159/259: Linear Algebra Notes Jarrod Millman November 16, 2015 Abstract These notes assume you ve taken a semester of undergraduate linear algebra. In particular, I assume you are familiar with the
More informationComputer-assisted proofs for radially symmetric solutions of PDEs
Computer-assisted proofs for radially symmetric solutions of PDEs István Balázs Jan Bouwe van den Berg Julien Courtois János Dudás Jean-Philippe Lessard Anett Vörös-Kiss JF Williams Xi Yuan Yin Abstract
More informationAffine iterations on nonnegative vectors
Affine iterations on nonnegative vectors V. Blondel L. Ninove P. Van Dooren CESAME Université catholique de Louvain Av. G. Lemaître 4 B-348 Louvain-la-Neuve Belgium Introduction In this paper we consider
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
More informationFUNCTIONAL ANALYSIS LECTURE NOTES: COMPACT SETS AND FINITE-DIMENSIONAL SPACES. 1. Compact Sets
FUNCTIONAL ANALYSIS LECTURE NOTES: COMPACT SETS AND FINITE-DIMENSIONAL SPACES CHRISTOPHER HEIL 1. Compact Sets Definition 1.1 (Compact and Totally Bounded Sets). Let X be a metric space, and let E X be
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
More informationNumerical Analysis Lecture Notes
Numerical Analysis Lecture Notes Peter J Olver 8 Numerical Computation of Eigenvalues In this part, we discuss some practical methods for computing eigenvalues and eigenvectors of matrices Needless to
More informationA Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators
A Perron-type theorem on the principal eigenvalue of nonsymmetric elliptic operators Lei Ni And I cherish more than anything else the Analogies, my most trustworthy masters. They know all the secrets of
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationarxiv: v2 [math.ds] 28 Jun 2016
VALIDATED COMPUTATION OF HETEROCLINIC SETS MACIEJ J. CAPIŃSKI AND J.D. MIRELES JAMES arxiv:1602.02973v2 [math.ds] 28 Jun 2016 Abstract. In this work we develop a method for computing mathematically rigorous
More informationIntroduction to rigorous numerics in dynamics: general functional analytic setup and an example that forces chaos
Introduction to rigorous numerics in dynamics: general functional analytic setup and an example that forces chaos Jan Bouwe van den Berg October 22, 2017 Abstract In this paper some basic concepts of rigorous
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationThe Eigenvalue Problem: Perturbation Theory
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 13 Notes These notes correspond to Sections 7.2 and 8.1 in the text. The Eigenvalue Problem: Perturbation Theory The Unsymmetric Eigenvalue Problem Just
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationParameterization of invariant manifolds for periodic orbits(i): efficient numerics via the Floquet normal form
Parameterization of invariant manifolds for periodic orbits(i): efficient numerics via the Floquet normal form Roberto Castelli Jean-Philippe Lessard J.D. Mireles James June 28, 2014 Abstract We present
More informationChapter 4. Inverse Function Theorem. 4.1 The Inverse Function Theorem
Chapter 4 Inverse Function Theorem d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d dd d d d d This chapter
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationVERIFIED ERROR BOUNDS FOR MULTIPLE ROOTS OF SYSTEMS OF NONLINEAR EQUATIONS
The original publication is available at www.springerlink.com VERIFIED ERROR BOUNDS FOR MULTIPLE ROOTS OF SYSTEMS OF NONLINEAR EQUATIONS SIEGFRIED M. RUMP AND STEF GRAILLAT Abstract. It is well known that
More informationGeneralized Numerical Radius Inequalities for Operator Matrices
International Mathematical Forum, Vol. 6, 011, no. 48, 379-385 Generalized Numerical Radius Inequalities for Operator Matrices Wathiq Bani-Domi Department of Mathematics Yarmouk University, Irbed, Jordan
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More information08a. Operators on Hilbert spaces. 1. Boundedness, continuity, operator norms
(February 24, 2017) 08a. Operators on Hilbert spaces Paul Garrett garrett@math.umn.edu http://www.math.umn.edu/ garrett/ [This document is http://www.math.umn.edu/ garrett/m/real/notes 2016-17/08a-ops
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationLinear Algebra Practice Problems
Math 7, Professor Ramras Linear Algebra Practice Problems () Consider the following system of linear equations in the variables x, y, and z, in which the constants a and b are real numbers. x y + z = a
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationPerron-Frobenius theorem for nonnegative multilinear forms and extensions
Perron-Frobenius theorem for nonnegative multilinear forms and extensions Shmuel Friedland Univ. Illinois at Chicago ensors 18 December, 2010, Hong-Kong Overview 1 Perron-Frobenius theorem for irreducible
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 informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
More informationOn rigorous integration of piece-wise linear systems
On rigorous integration of piece-wise linear systems Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology 06.2011, Trieste, Italy On rigorous integration of PWL
More informationLecture 1: Review of linear algebra
Lecture 1: Review of linear algebra Linear functions and linearization Inverse matrix, least-squares and least-norm solutions Subspaces, basis, and dimension Change of basis and similarity transformations
More informationA New Proof of Anti-Maximum Principle Via A Bifurcation Approach
A New Proof of Anti-Maximum Principle Via A Bifurcation Approach Junping Shi Abstract We use a bifurcation approach to prove an abstract version of anti-maximum principle. The proof is different from previous
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationComputing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices
Computing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices Alper Demir David Long Jaijeet Roychowdhury Bell Laboratories Murray Hill New Jersey USA Abstract The main effort in oscillator
More informationOn the Local Quadratic Convergence of the Primal-Dual Augmented Lagrangian Method
Optimization Methods and Software Vol. 00, No. 00, Month 200x, 1 11 On the Local Quadratic Convergence of the Primal-Dual Augmented Lagrangian Method ROMAN A. POLYAK Department of SEOR and Mathematical
More informationScientific Computing: Dense Linear Systems
Scientific Computing: Dense Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 February 9th, 2012 A. Donev (Courant Institute)
More informationAN ITERATIVE METHOD FOR COMPUTING ZEROS OF OPERATORS SATISFYING AUTONOMOUS DIFFERENTIAL EQUATIONS
Electronic Journal: Southwest Journal of Pure Applied Mathematics Internet: http://rattler.cameron.edu/swjpam.html ISBN 1083-0464 Issue 1 July 2004, pp. 48 53 Submitted: August 21, 2003. Published: July
More informationFinding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients. Takayuki Yamaguchi
Hokkaido Mathematical Journal Vol. 44 (2015) p. 277 312 Finding numerically Newhouse sinks near a homoclinic tangency and investigation of their chaotic transients Takayuki Yamaguchi (Received March 13,
More informationAn estimation of the spectral radius of a product of block matrices
Linear Algebra and its Applications 379 (2004) 267 275 wwwelseviercom/locate/laa An estimation of the spectral radius of a product of block matrices Mei-Qin Chen a Xiezhang Li b a Department of Mathematics
More informationLebesgue Measure on R n
CHAPTER 2 Lebesgue Measure on R n Our goal is to construct a notion of the volume, or Lebesgue measure, of rather general subsets of R n that reduces to the usual volume of elementary geometrical sets
More informationNUMERICAL ENCLOSURE OF SOLUTIONS FOR TWO DIMENSIONAL DRIVEN CAVITY PROBLEMS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) O. Nevanlinna and R. Rannacher (assoc. eds.) Jyväskylä,
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationMATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,
MATH 205 HOMEWORK #3 OFFICIAL SOLUTION Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. a F = R, V = R 3, b F = R or C, V = F 2, T = T = 9 4 4 8 3 4 16 8 7 0 1
More informationDynamical Systems and Ergodic Theory PhD Exam Spring Topics: Topological Dynamics Definitions and basic results about the following for maps and
Dynamical Systems and Ergodic Theory PhD Exam Spring 2012 Introduction: This is the first year for the Dynamical Systems and Ergodic Theory PhD exam. As such the material on the exam will conform fairly
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationMathematica Bohemica
Mathematica Bohemica Cristian Bereanu; Jean Mawhin Existence and multiplicity results for nonlinear second order difference equations with Dirichlet boundary conditions Mathematica Bohemica, Vol. 131 (2006),
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More informationContinuation of cycle-to-cycle connections in 3D ODEs
HET p. 1/2 Continuation of cycle-to-cycle connections in 3D ODEs Yuri A. Kuznetsov joint work with E.J. Doedel, B.W. Kooi, and G.A.K. van Voorn HET p. 2/2 Contents Previous works Truncated BVP s with projection
More informationComputational math: Assignment 1
Computational math: Assignment 1 Thanks Ting Gao for her Latex file 11 Let B be a 4 4 matrix to which we apply the following operations: 1double column 1, halve row 3, 3add row 3 to row 1, 4interchange
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More information5 Compact linear operators
5 Compact linear operators One of the most important results of Linear Algebra is that for every selfadjoint linear map A on a finite-dimensional space, there exists a basis consisting of eigenvectors.
More informationMATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam. Topics
MATH 213 Linear Algebra and Ordinary Differential Equations Spring 2015 Study Sheet for Final Exam This study sheet will not be allowed during the test. Books and notes will not be allowed during the test.
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationLinear Algebra Methods for Data Mining
Linear Algebra Methods for Data Mining Saara Hyvönen, Saara.Hyvonen@cs.helsinki.fi Spring 2007 1. Basic Linear Algebra Linear Algebra Methods for Data Mining, Spring 2007, University of Helsinki Example
More informationLinear Algebra for Machine Learning. Sargur N. Srihari
Linear Algebra for Machine Learning Sargur N. srihari@cedar.buffalo.edu 1 Overview Linear Algebra is based on continuous math rather than discrete math Computer scientists have little experience with it
More informationOn the Stability of the Best Reply Map for Noncooperative Differential Games
On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University
More informationat time t, in dimension d. The index i varies in a countable set I. We call configuration the family, denoted generically by Φ: U (x i (t) x j (t))
Notations In this chapter we investigate infinite systems of interacting particles subject to Newtonian dynamics Each particle is characterized by its position an velocity x i t, v i t R d R d at time
More information