A review of sparsity vs stability in LU updates
|
|
- Charity Cox
- 5 years ago
- Views:
Transcription
1 A review of sparsity vs stability in LU updates Michael Saunders Dept of Management Science and Engineering (MS&E) Systems Optimization Laboratory (SOL) Institute for Computational Mathematics and Engineering (ICME) Stanford University, California Session: 40 Years of Forrest and Tomlin INFORMS 2012 Annual Meeting Phoenix AZ, Oct 14 18, 2012 Michael Saunders Sparsity vs stability in LU updates 1/31
2 Abstract The Forrest-Tomlin update has stood the test of time within many generations of commercial mathematical programming systems. Its ease of implementation leads to high efficiency and evidently acceptable reliability. We review its relation to Reid s version of the Bartels-Golub update as implemented in LA05, LA15, and LUSOL. In particular, we examine the extent to which FT implementations must live dangerously in order to achieve the desired efficiency. Michael Saunders Sparsity vs stability in LU updates 2/31
3 Outline 1 The Bartels-Golub update 2 The Forrest-Tomlin update 3 Sparse Bartels-Golub 4 Two tolerances 5 A numerical experiment 6 Numerical results 7 Conclusion Michael Saunders Sparsity vs stability in LU updates 3/31
4 The Bartels-Golub update Michael Saunders Sparsity vs stability in LU updates 4/31
5 Dense Bartels-Golub p ū Michael Saunders Sparsity vs stability in LU updates 5/31
6 Dense Bartels-Golub X X X X X X Michael Saunders Sparsity vs stability in LU updates 5/31
7 Dense Bartels-Golub X X X X X X Forget the Hessenberg Michael Saunders Sparsity vs stability in LU updates 5/31
8 The Forrest-Tomlin update Michael Saunders Sparsity vs stability in LU updates 6/31
9 Forrest-Tomlin (conceptually) X X Few nonzeros to eliminate U stored by columns Michael Saunders Sparsity vs stability in LU updates 7/31
10 Sparse Bartels-Golub Michael Saunders Sparsity vs stability in LU updates 8/31
11 Sparse Bartels-Golub (Reid, Saunders) X p X X ū l Michael Saunders Sparsity vs stability in LU updates 9/31
12 Sparse Bartels-Golub (Reid, Saunders) X X l LA05, LA15, LUSOL U stored by rows Michael Saunders Sparsity vs stability in LU updates 9/31
13 The connection X X FT: Eliminate X s with diags as pivots BG: Allow row interchanges ( ) 1 Use a product of and permutations µ 1 How big can µ safely be?? Michael Saunders Sparsity vs stability in LU updates 10/31
14 A Tale of Two Tolerances LU factor tol α LU update tol β Control and/or monitor stability Michael Saunders Sparsity vs stability in LU updates 11/31
15 LU factor tol α Threshold partial pivoting controls cond(l) B = LU L ii = 1 L ij α α = 10.0 or usually sparse and stable Michael Saunders Sparsity vs stability in LU updates 12/31
16 LU factor tol α Threshold partial pivoting controls cond(l) B = LU L ii = 1 L ij α α = 10.0 or usually sparse and stable Threshold rook pivoting controls cond(l), cond(u) B = LDU L ii = 1 L ij α and U ii = 1 U ij α 1 < α 2.0 rank-revealing (for basis repair) Michael Saunders Sparsity vs stability in LU updates 12/31
17 Factor and Update involve triangular matrices M = ( ) cond(m) =?? Michael Saunders Sparsity vs stability in LU updates 13/31
18 Factor and Update involve triangular matrices M = ( ) cond(m) 10 4 Michael Saunders Sparsity vs stability in LU updates 13/31
19 Factor and Update involve triangular matrices M = ( ) 1 µ 1 cond(m) µ 2 µ > 1 Michael Saunders Sparsity vs stability in LU updates 13/31
20 FT update involves triangular matrices 1 R = r p+1... r m 1 cond(r) r 2 max Michael Saunders Sparsity vs stability in LU updates 14/31
21 Product-form update involves triangular matrices 1 v 1 1. T = v p. v m... I cond(t ) v 2 max/ v p Michael Saunders Sparsity vs stability in LU updates 15/31
22 LU update tol β B = B(I + (v e p )e T p ) Use β to control or monitor basis updates Product-form B = BT monitor cond(t ) Bartels-Golub L = LM 1 M 2... control cond(m j ) Forrest-Tomlin L = LR monitor cond(r) Michael Saunders Sparsity vs stability in LU updates 16/31
23 Product-form update LU update tol β B = BT T = I + (v e p )ep T monitor cond(t ) Update if vmax/ v 2 p β else refactor How big should β be?? 10 6? Michael Saunders Sparsity vs stability in LU updates 17/31
24 Product-form update LU update tol β B = BT T = I + (v e p )e T p monitor cond(t ) LU updates Update if vmax/ v 2 p β else refactor How big should β be?? 10 6? BG L = LM 1 M 2... control M ij β FT L = LR monitor R ij β? For BG we can enforce β = 10 say For FT how big should β be?? 10 6? Michael Saunders Sparsity vs stability in LU updates 17/31
25 A numerical experiment Michael Saunders Sparsity vs stability in LU updates 18/31
26 Simulate FT via BG Run MINOS on a numerically challenging LP (e.g. pilot87) Set LUSOL s LU update tol β = 10, 10 2, 10 3, 10 4, 10 5, 10 6,... Count number of times something broke Michael Saunders Sparsity vs stability in LU updates 19/31
27 Simulate FT via BG Run MINOS on a numerically challenging LP (e.g. pilot87) Set LUSOL s LU update tol β = 10, 10 2, 10 3, 10 4, 10 5, 10 6,... Count number of times something broke Normally Factorization freq = 100 LU factor tol α = LU update tol β = 10.0 Michael Saunders Sparsity vs stability in LU updates 19/31
28 Simulate FT via BG Run MINOS on a numerically challenging LP (e.g. pilot87) Set LUSOL s LU update tol β = 10, 10 2, 10 3, 10 4, 10 5, 10 6,... Count number of times something broke Normally Factorization freq = 100 LU factor tol α = LU update tol β = 10.0 How would MINOS measure broke? Michael Saunders Sparsity vs stability in LU updates 19/31
29 Simulate FT via BG Run MINOS on a numerically challenging LP (e.g. pilot87) Set LUSOL s LU update tol β = 10, 10 2, 10 3, 10 4, 10 5, 10 6,... Count number of times something broke Normally Factorization freq = 100 LU factor tol α = LU update tol β = 10.0 How would MINOS measure broke? Row check failed? Never happened Every 60 itns, see if Ax b is too big Michael Saunders Sparsity vs stability in LU updates 19/31
30 Simulate FT via BG Run MINOS on a numerically challenging LP (e.g. pilot87) Set LUSOL s LU update tol β = 10, 10 2, 10 3, 10 4, 10 5, 10 6,... Count number of times something broke Normally Factorization freq = 100 LU factor tol α = LU update tol β = 10.0 How would MINOS measure broke? Row check failed? Never happened Every 60 itns, see if Ax b is too big LU became singular? Every update, check if U has small diag Increasingly often Michael Saunders Sparsity vs stability in LU updates 19/31
31 Tests for accepting LU update ( ) I L = L r T 1 Ū = ( ) U1 u δ Tomlin 1975 Our experiment r max u max δ 16 7 ( ) r max β and u 1 δ ɛ 2/3 ( ) Michael Saunders Sparsity vs stability in LU updates 20/31
32 Tests for accepting LU update ( ) I L = L r T 1 Ū = ( ) U1 u δ Tomlin 1975 Our experiment r max u max δ r max β and u 1 δ ɛ 2/3 ( ) Roughly equivalent to r max u max δ β10 10 = 10 11, 10 12,... Michael Saunders Sparsity vs stability in LU updates 20/31
33 Numerical results Michael Saunders Sparsity vs stability in LU updates 21/31
34 MINOS simulation of FT updates LPnetlib problem pilot ɛ = Featol = Opttol = 10 6 LU update tol β = 10, 10 2, 10 3,... singular average beta itns time LUfacs warnings updates 1e e e e e e e e Michael Saunders Sparsity vs stability in LU updates 22/31
35 MINOS simulation of FT updates 120 pilots: LU update tol = 1e+8 Average 31 updates 100 Updates before new LU LU number Michael Saunders Sparsity vs stability in LU updates 23/31
36 MINOS simulation of FT updates LPnetlib problem pilot87 ɛ = Featol = Opttol = 10 6 LU update tol β = 10, 10 2, 10 3,... singular average beta itns time LUfacs warnings updates 1e e e e e e e e Michael Saunders Sparsity vs stability in LU updates 24/31
37 MINOS simulation of FT updates LPnetlib problem mod2pre ɛ = Featol = Opttol = 10 6 LU update tol β = 10, 10 2, 10 3,... singular beta itns time LUfacs warnings 1e e e e * 1e * 1e infeas 1e * 1e infeas Refactorization gave singular LU Rook pivoting invoked, β reduced to 10.0 Michael Saunders Sparsity vs stability in LU updates 25/31
38 Quad SQOPT simulation of FT updates LPnetlib problem pilot ɛ = Featol = Opttol = LU update tol β = 10 4, 10 6, 10 8,... singular average dual beta itns time LUfacs warnings updates infeas 1e E-25 1e E-27 1e E-21 1e E-24 1e E-27 1e E-29 1e E-12 1e E-19 1e * 1e * 1e * *Refac singular, beta reduced to 50.0 Michael Saunders Sparsity vs stability in LU updates 26/31
39 Conclusion Michael Saunders Sparsity vs stability in LU updates 27/31
40 Conclusion Test for accepting an FT update r max u max WHIZARD 10 8 δ CLP, CPLEX, Gurobi, Xpress, Mozek,... r max u max δ? or? Probably safer r max β and u max δ Average updates per LU β = β = β = β = Michael Saunders Sparsity vs stability in LU updates 28/31
41 References (in chronological order) R. H. Bartels (1971). A stabilization of the simplex method, Numerische Mathematik 16, Eliminated subdiag of dense Hessenberg H = ŪP Michael Saunders Sparsity vs stability in LU updates 29/31
42 References (in chronological order) R. H. Bartels (1971). A stabilization of the simplex method, Numerische Mathematik 16, Eliminated subdiag of dense Hessenberg H = ŪP J. J. H. Forrest and J. A. Tomlin (1972). Updating triangular factors of the basis to maintain sparsity in the product form simplex method, Mathematical Programming 2, Eliminated pth row of sparse H (not subdiag) Michael Saunders Sparsity vs stability in LU updates 29/31
43 References (in chronological order) R. H. Bartels (1971). A stabilization of the simplex method, Numerische Mathematik 16, Eliminated subdiag of dense Hessenberg H = ŪP J. J. H. Forrest and J. A. Tomlin (1972). Updating triangular factors of the basis to maintain sparsity in the product form simplex method, Mathematical Programming 2, Eliminated pth row of sparse H (not subdiag) J. A. Tomlin (1975). An accuracy test for updating triangular factors, Math. Prog. Study 4, Test for an unstable update Michael Saunders Sparsity vs stability in LU updates 29/31
44 References (in chronological order) R. H. Bartels (1971). A stabilization of the simplex method, Numerische Mathematik 16, Eliminated subdiag of dense Hessenberg H = ŪP J. J. H. Forrest and J. A. Tomlin (1972). Updating triangular factors of the basis to maintain sparsity in the product form simplex method, Mathematical Programming 2, Eliminated pth row of sparse H (not subdiag) J. A. Tomlin (1975). An accuracy test for updating triangular factors, Math. Prog. Study 4, Test for an unstable update M. A. Saunders (1976). The complexity of LU updating, in Anderssen and Brent (eds.), The Complexity of Computational Problem Solving, Picture of symmetrically permuted P T ŪP with row spike to be eliminated Michael Saunders Sparsity vs stability in LU updates 29/31
45 References (contd) J. K. Reid (1982). A sparsity-exploiting variant of the Bartels-Golub decomposition for linear programming bases, Mathematical Programming 24, U stored by rows to facilitate BG update Michael Saunders Sparsity vs stability in LU updates 30/31
46 References (contd) J. K. Reid (1982). A sparsity-exploiting variant of the Bartels-Golub decomposition for linear programming bases, Mathematical Programming 24, U stored by rows to facilitate BG update P. E. Gill, W. Murray, M. A. Saunders, and M. H. Wright (1987). Maintaining LU factors of a general sparse matrix, Linear Algebra and its Applics. 88/89, Rectangular LU, more general updates, cond(l) controlled throughout Michael Saunders Sparsity vs stability in LU updates 30/31
47 References (contd) J. K. Reid (1982). A sparsity-exploiting variant of the Bartels-Golub decomposition for linear programming bases, Mathematical Programming 24, U stored by rows to facilitate BG update P. E. Gill, W. Murray, M. A. Saunders, and M. H. Wright (1987). Maintaining LU factors of a general sparse matrix, Linear Algebra and its Applics. 88/89, Rectangular LU, more general updates, cond(l) controlled throughout P. E. Gill, W. Murray, and M. A. Saunders (2005). SNOPT: An SQP algorithm for large-scale constrained optimization, SIAM Review 47:1, LU with Threshold Rook Pivoting, same general updates Michael Saunders Sparsity vs stability in LU updates 30/31
48 Forrest-Tomlin Michael Saunders Sparsity vs stability in LU updates 31/31
LUSOL: A basis package for constrained optimization
LUSOL: A basis package for constrained optimization IFORS Triennial Conference on OR/MS Honolulu, HI, July 11 15, 2005 Michael O Sullivan and Michael Saunders Dept of Engineering Science Dept of Management
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME MS&E 38 (CME 338 Large-Scale Numerical Optimization Course description Instructor: Michael Saunders Spring 28 Notes : Review The course teaches
More informationSparse Rank-Revealing LU Factorization
Sparse Rank-Revealing LU Factorization SIAM Conference on Optimization Toronto, May 2002 Michael O Sullivan and Michael Saunders Dept of Engineering Science Dept of Management Sci & Eng University of Auckland
More informationSparse Rank-Revealing LU Factorization
Sparse Rank-Revealing LU Factorization Householder Symposium XV on Numerical Linear Algebra Peebles, Scotland, June 2002 Michael O Sullivan and Michael Saunders Dept of Engineering Science Dept of Management
More informationA convex QP solver based on block-lu updates
Block-LU updates p. 1/24 A convex QP solver based on block-lu updates PP06 SIAM Conference on Parallel Processing for Scientific Computing San Francisco, CA, Feb 22 24, 2006 Hanh Huynh and Michael Saunders
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization Instructor: Michael Saunders Spring 2015 Notes 4: The Primal Simplex Method 1 Linear
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME MS&E 318 (CME 338 Large-Scale Numerical Optimization 1 Origins Instructor: Michael Saunders Spring 2018 Notes 7: LUSOL: a Basis Factorization
More informationImplementation of a KKT-based active-set QP solver
Block-LU updates p. 1/27 Implementation of a KKT-based active-set QP solver ISMP 2006 19th International Symposium on Mathematical Programming Rio de Janeiro, Brazil, July 30 August 4, 2006 Hanh Huynh
More informationSatisfying Flux Balance and Mass-Action Kinetics in a Network of Biochemical Reactions
Satisfying Flux Balance and Mass-Action Kinetics in a Network of Biochemical Reactions Speakers Senior Personnel Collaborators Michael Saunders and Ines Thiele Ronan Fleming, Bernhard Palsson, Yinyu Ye
More informationLinear Algebra Section 2.6 : LU Decomposition Section 2.7 : Permutations and transposes Wednesday, February 13th Math 301 Week #4
Linear Algebra Section. : LU Decomposition Section. : Permutations and transposes Wednesday, February 1th Math 01 Week # 1 The LU Decomposition We learned last time that we can factor a invertible matrix
More informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 13
STAT 309: MATHEMATICAL COMPUTATIONS I FALL 208 LECTURE 3 need for pivoting we saw that under proper circumstances, we can write A LU where 0 0 0 u u 2 u n l 2 0 0 0 u 22 u 2n L l 3 l 32, U 0 0 0 l n l
More informationScientific Computing
Scientific Computing Direct solution methods Martin van Gijzen Delft University of Technology October 3, 2018 1 Program October 3 Matrix norms LU decomposition Basic algorithm Cost Stability Pivoting Pivoting
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6
CME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 6 GENE H GOLUB Issues with Floating-point Arithmetic We conclude our discussion of floating-point arithmetic by highlighting two issues that frequently
More informationSolving Linear and Integer Programs
Solving Linear and Integer Programs Robert E. Bixby ILOG, Inc. and Rice University Ed Rothberg ILOG, Inc. DAY 2 2 INFORMS Practice 2002 1 Dual Simplex Algorithm 3 Some Motivation Dual simplex vs. primal
More informationI-v k e k. (I-e k h kt ) = Stability of Gauss-Huard Elimination for Solving Linear Systems. 1 x 1 x x x x
Technical Report CS-93-08 Department of Computer Systems Faculty of Mathematics and Computer Science University of Amsterdam Stability of Gauss-Huard Elimination for Solving Linear Systems T. J. Dekker
More informationNovel update techniques for the revised simplex method (and their application)
Novel update techniques for the revised simplex method (and their application) Qi Huangfu 1 Julian Hall 2 Others 1 FICO 2 School of Mathematics, University of Edinburgh ERGO 30 November 2016 Overview Background
More information1. Introduction. We are concerned with the linear programming (LP) problem in the standard form. minimize subject to Ax = b, x 0,
A GENERALIZATION OF THE REVISED SIMPLEX ALGORITHM FOR LINEAR PROGRAMMING PING-QI PAN Abstract. Recently the concept of basis, which plays a fundamental role in the simplex methodology, were generalized
More informationPivoting. Reading: GV96 Section 3.4, Stew98 Chapter 3: 1.3
Pivoting Reading: GV96 Section 3.4, Stew98 Chapter 3: 1.3 In the previous discussions we have assumed that the LU factorization of A existed and the various versions could compute it in a stable manner.
More informationChapter 8 Cholesky-based Methods for Sparse Least Squares: The Benefits of Regularization
In L. Adams and J. L. Nazareth eds., Linear and Nonlinear Conjugate Gradient-Related Methods, SIAM, Philadelphia, 92 100 1996. Chapter 8 Cholesky-based Methods for Sparse Least Squares: The Benefits of
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Logistics Notes for 2016-09-14 1. There was a goof in HW 2, problem 1 (now fixed) please re-download if you have already started looking at it. 2. CS colloquium (4:15 in Gates G01) this Thurs is Margaret
More informationMATH 3511 Lecture 1. Solving Linear Systems 1
MATH 3511 Lecture 1 Solving Linear Systems 1 Dmitriy Leykekhman Spring 2012 Goals Review of basic linear algebra Solution of simple linear systems Gaussian elimination D Leykekhman - MATH 3511 Introduction
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization Instructor: Michael Saunders Spring 2017 Notes 10: MINOS Part 1: the Reduced-Gradient
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization Instructor: Michael Saunders Spring 2015 Notes 11: NPSOL and SNOPT SQP Methods 1 Overview
More informationThe Dual Flow Between Linear Algebra and Optimization
The Dual Flow Between Linear Algebra and Optimization Margaret H. Wright Computer Science Department Courant Institute of Mathematical Sciences New York University History of Numerical Linear Algebra Minisymposium
More information1 Multiply Eq. E i by λ 0: (λe i ) (E i ) 2 Multiply Eq. E j by λ and add to Eq. E i : (E i + λe j ) (E i )
Direct Methods for Linear Systems Chapter Direct Methods for Solving Linear Systems Per-Olof Persson persson@berkeleyedu Department of Mathematics University of California, Berkeley Math 18A Numerical
More informationIntroduction to Mathematical Programming
Introduction to Mathematical Programming Ming Zhong Lecture 6 September 12, 2018 Ming Zhong (JHU) AMS Fall 2018 1 / 20 Table of Contents 1 Ming Zhong (JHU) AMS Fall 2018 2 / 20 Solving Linear Systems A
More informationMain matrix factorizations
Main matrix factorizations A P L U P permutation matrix, L lower triangular, U upper triangular Key use: Solve square linear system Ax b. A Q R Q unitary, R upper triangular Key use: Solve square or overdetrmined
More informationNumerical Linear Algebra
Numerical Linear Algebra Direct Methods Philippe B. Laval KSU Fall 2017 Philippe B. Laval (KSU) Linear Systems: Direct Solution Methods Fall 2017 1 / 14 Introduction The solution of linear systems is one
More informationNumerical Linear Algebra
Numerical Linear Algebra Decompositions, numerical aspects Gerard Sleijpen and Martin van Gijzen September 27, 2017 1 Delft University of Technology Program Lecture 2 LU-decomposition Basic algorithm Cost
More informationProgram Lecture 2. Numerical Linear Algebra. Gaussian elimination (2) Gaussian elimination. Decompositions, numerical aspects
Numerical Linear Algebra Decompositions, numerical aspects Program Lecture 2 LU-decomposition Basic algorithm Cost Stability Pivoting Cholesky decomposition Sparse matrices and reorderings Gerard Sleijpen
More informationReview Questions REVIEW QUESTIONS 71
REVIEW QUESTIONS 71 MATLAB, is [42]. For a comprehensive treatment of error analysis and perturbation theory for linear systems and many other problems in linear algebra, see [126, 241]. An overview of
More informationPermutations in the factorization of simplex bases
Permutations in the factorization of simplex bases Ricardo Fukasawa, Laurent Poirrier {rfukasawa,lpoirrier}@uwaterloo.ca December 13, 2016 Abstract The basis matrices corresponding to consecutive iterations
More informationRevised Simplex Method
DM545 Linear and Integer Programming Lecture 7 Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 2 Motivation Complexity of single pivot operation
More informationAN EFFICIENT APPROACH TO UPDATING SIMPLEX MULTIPLIERS IN THE SIMPLEX ALGORITHM
AN EFFICIENT APPROACH TO UPDATING SIMPLEX MULTIPLIERS IN THE SIMPLEX ALGORITHM JIAN-FENG HU AND PING-QI PAN Abstract. The simplex algorithm computes the simplex multipliers by solving a system (or two
More informationWhat s New in Active-Set Methods for Nonlinear Optimization?
What s New in Active-Set Methods for Nonlinear Optimization? Philip E. Gill Advances in Numerical Computation, Manchester University, July 5, 2011 A Workshop in Honor of Sven Hammarling UCSD Center for
More informationLinear Programming The Simplex Algorithm: Part II Chapter 5
1 Linear Programming The Simplex Algorithm: Part II Chapter 5 University of Chicago Booth School of Business Kipp Martin October 17, 2017 Outline List of Files Key Concepts Revised Simplex Revised Simplex
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationLinear Programming: Simplex
Linear Programming: Simplex Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. IMA, August 2016 Stephen Wright (UW-Madison) Linear Programming: Simplex IMA, August 2016
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 12: Gaussian Elimination and LU Factorization Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 10 Gaussian Elimination
More informationSparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices
Sparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices Jaehyun Park June 1 2016 Abstract We consider the problem of writing an arbitrary symmetric matrix as
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME) MS&E 318 (CME 338) Large-Scale Numerical Optimization 1 Origins Instructor: Michael Saunders Spring 2015 Notes 9: Augmented Lagrangian Methods
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 19: Computing the SVD; Sparse Linear Systems Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationSection 5.6. LU and LDU Factorizations
5.6. LU and LDU Factorizations Section 5.6. LU and LDU Factorizations Note. We largely follow Fraleigh and Beauregard s approach to this topic from Linear Algebra, 3rd Edition, Addison-Wesley (995). See
More informationLinear Algebra Linear Algebra : Matrix decompositions Monday, February 11th Math 365 Week #4
Linear Algebra Linear Algebra : Matrix decompositions Monday, February 11th Math Week # 1 Saturday, February 1, 1 Linear algebra Typical linear system of equations : x 1 x +x = x 1 +x +9x = 0 x 1 +x x
More informationMODULE 7. where A is an m n real (or complex) matrix. 2) Let K(t, s) be a function of two variables which is continuous on the square [0, 1] [0, 1].
Topics: Linear operators MODULE 7 We are going to discuss functions = mappings = transformations = operators from one vector space V 1 into another vector space V 2. However, we shall restrict our sights
More informationMatrices and systems of linear equations
Matrices and systems of linear equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T.
More informationLU Factorization. Marco Chiarandini. DM559 Linear and Integer Programming. Department of Mathematics & Computer Science University of Southern Denmark
DM559 Linear and Integer Programming LU Factorization Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark [Based on slides by Lieven Vandenberghe, UCLA] Outline
More informationLecture 12 (Tue, Mar 5) Gaussian elimination and LU factorization (II)
Math 59 Lecture 2 (Tue Mar 5) Gaussian elimination and LU factorization (II) 2 Gaussian elimination - LU factorization For a general n n matrix A the Gaussian elimination produces an LU factorization if
More informationDM545 Linear and Integer Programming. Lecture 7 Revised Simplex Method. Marco Chiarandini
DM545 Linear and Integer Programming Lecture 7 Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 2 Motivation Complexity of single pivot operation
More informationAck: 1. LD Garcia, MTH 199, Sam Houston State University 2. Linear Algebra and Its Applications - Gilbert Strang
Gaussian Elimination CS6015 : Linear Algebra Ack: 1. LD Garcia, MTH 199, Sam Houston State University 2. Linear Algebra and Its Applications - Gilbert Strang The Gaussian Elimination Method The Gaussian
More informationAx=b. Zack 10/4/2013
Ax=b Zack 10/4/2013 Iteration method Ax=b v 1, v 2 x k = V k y k Given (A, b) standard orthonormal base v 1, v 2 x k = V k y k For symmetric A, the method to generate V k is called Lanczos Method Lanczos
More informationLU Factorization. LU factorization is the most common way of solving linear systems! Ax = b LUx = b
AM 205: lecture 7 Last time: LU factorization Today s lecture: Cholesky factorization, timing, QR factorization Reminder: assignment 1 due at 5 PM on Friday September 22 LU Factorization LU factorization
More informationMatrix decompositions
Matrix decompositions How can we solve Ax = b? 1 Linear algebra Typical linear system of equations : x 1 x +x = x 1 +x +9x = 0 x 1 +x x = The variables x 1, x, and x only appear as linear terms (no powers
More informationNumerical Linear Algebra
Chapter 3 Numerical Linear Algebra We review some techniques used to solve Ax = b where A is an n n matrix, and x and b are n 1 vectors (column vectors). We then review eigenvalues and eigenvectors and
More informationNumerical Linear Algebra
Numerical Linear Algebra The two principal problems in linear algebra are: Linear system Given an n n matrix A and an n-vector b, determine x IR n such that A x = b Eigenvalue problem Given an n n matrix
More information9. Numerical linear algebra background
Convex Optimization Boyd & Vandenberghe 9. Numerical linear algebra background matrix structure and algorithm complexity solving linear equations with factored matrices LU, Cholesky, LDL T factorization
More informationRank Revealing QR factorization. F. Guyomarc h, D. Mezher and B. Philippe
Rank Revealing QR factorization F. Guyomarc h, D. Mezher and B. Philippe 1 Outline Introduction Classical Algorithms Full matrices Sparse matrices Rank-Revealing QR Conclusion CSDA 2005, Cyprus 2 Situation
More informationGAUSSIAN ELIMINATION AND LU DECOMPOSITION (SUPPLEMENT FOR MA511)
GAUSSIAN ELIMINATION AND LU DECOMPOSITION (SUPPLEMENT FOR MA511) D. ARAPURA Gaussian elimination is the go to method for all basic linear classes including this one. We go summarize the main ideas. 1.
More information9 Appendix. Determinants and Cramer s formula
LINEAR ALGEBRA: THEORY Version: August 12, 2000 133 9 Appendix Determinants and Cramer s formula Here we the definition of the determinant in the general case and summarize some features Then we show how
More informationMatrix decompositions
Matrix decompositions How can we solve Ax = b? 1 Linear algebra Typical linear system of equations : x 1 x +x = x 1 +x +9x = 0 x 1 +x x = The variables x 1, x, and x only appear as linear terms (no powers
More informationKey words. numerical linear algebra, direct methods, Cholesky factorization, sparse matrices, mathematical software, matrix updates.
ROW MODIFICATIONS OF A SPARSE CHOLESKY FACTORIZATION TIMOTHY A. DAVIS AND WILLIAM W. HAGER Abstract. Given a sparse, symmetric positive definite matrix C and an associated sparse Cholesky factorization
More informationMATH 2030: MATRICES. Example 0.2. Q:Define A 1 =, A. 3 4 A: We wish to find c 1, c 2, and c 3 such that. c 1 + c c
MATH 2030: MATRICES Matrix Algebra As with vectors, we may use the algebra of matrices to simplify calculations. However, matrices have operations that vectors do not possess, and so it will be of interest
More informationCS 323: Numerical Analysis and Computing
CS 323: Numerical Analysis and Computing MIDTERM #1 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 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 informationMatrix Multiplication Chapter IV Special Linear Systems
Matrix Multiplication Chapter IV Special Linear Systems By Gokturk Poyrazoglu The State University of New York at Buffalo BEST Group Winter Lecture Series Outline 1. Diagonal Dominance and Symmetry a.
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationAM 205: lecture 8. Last time: Cholesky factorization, QR factorization Today: how to compute the QR factorization, the Singular Value Decomposition
AM 205: lecture 8 Last time: Cholesky factorization, QR factorization Today: how to compute the QR factorization, the Singular Value Decomposition QR Factorization A matrix A R m n, m n, can be factorized
More informationExperiments with iterative computation of search directions within interior methods for constrained optimization
Experiments with iterative computation of search directions within interior methods for constrained optimization Santiago Akle and Michael Saunders Institute for Computational Mathematics and Engineering
More informationCS412: Lecture #17. Mridul Aanjaneya. March 19, 2015
CS: Lecture #7 Mridul Aanjaneya March 9, 5 Solving linear systems of equations Consider a lower triangular matrix L: l l l L = l 3 l 3 l 33 l n l nn A procedure similar to that for upper triangular systems
More informationIntroduction to Applied Linear Algebra with MATLAB
Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationThe practical revised simplex method (Part 2)
The practical revised simplex method (Part 2) Julian Hall School of Mathematics University of Edinburgh January 25th 2007 The practical revised simplex method Overview (Part 2) Practical implementation
More informationCS227-Scientific Computing. Lecture 4: A Crash Course in Linear Algebra
CS227-Scientific Computing Lecture 4: A Crash Course in Linear Algebra Linear Transformation of Variables A common phenomenon: Two sets of quantities linearly related: y = 3x + x 2 4x 3 y 2 = 2.7x 2 x
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More informationSolution of Linear Equations
Solution of Linear Equations (Com S 477/577 Notes) Yan-Bin Jia Sep 7, 07 We have discussed general methods for solving arbitrary equations, and looked at the special class of polynomial equations A subclass
More informationLinear Systems of n equations for n unknowns
Linear Systems of n equations for n unknowns In many application problems we want to find n unknowns, and we have n linear equations Example: Find x,x,x such that the following three equations hold: x
More informationGaussian Elimination without/with Pivoting and Cholesky Decomposition
Gaussian Elimination without/with Pivoting and Cholesky Decomposition Gaussian Elimination WITHOUT pivoting Notation: For a matrix A R n n we define for k {,,n} the leading principal submatrix a a k A
More informationGaussian Elimination and Back Substitution
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 4 Notes These notes correspond to Sections 31 and 32 in the text Gaussian Elimination and Back Substitution The basic idea behind methods for solving
More informationNumerical Linear Algebra Primer. Ryan Tibshirani Convex Optimization /36-725
Numerical Linear Algebra Primer Ryan Tibshirani Convex Optimization 10-725/36-725 Last time: proximal gradient descent Consider the problem min g(x) + h(x) with g, h convex, g differentiable, and h simple
More information7. LU factorization. factor-solve method. LU factorization. solving Ax = b with A nonsingular. the inverse of a nonsingular matrix
EE507 - Computational Techniques for EE 7. LU factorization Jitkomut Songsiri factor-solve method LU factorization solving Ax = b with A nonsingular the inverse of a nonsingular matrix LU factorization
More informationAMS 209, Fall 2015 Final Project Type A Numerical Linear Algebra: Gaussian Elimination with Pivoting for Solving Linear Systems
AMS 209, Fall 205 Final Project Type A Numerical Linear Algebra: Gaussian Elimination with Pivoting for Solving Linear Systems. Overview We are interested in solving a well-defined linear system given
More informationAlgorithm 853: an Efficient Algorithm for Solving Rank-Deficient Least Squares Problems
Algorithm 853: an Efficient Algorithm for Solving Rank-Deficient Least Squares Problems LESLIE FOSTER and RAJESH KOMMU San Jose State University Existing routines, such as xgelsy or xgelsd in LAPACK, for
More informationNumerical Methods I Solving Square Linear Systems: GEM and LU factorization
Numerical Methods I Solving Square Linear Systems: GEM and LU factorization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 18th,
More informationarxiv:math/ v1 [math.oc] 22 Dec 2000
The simplest examples where the simplex method cycles and conditions where expand fails to prevent cycling 1 arxiv:math/0012242v1 [math.oc] 22 Dec 2000 J.A.J. Hall and K.I.M. McKinnon 2 Dept. of Mathematics
More informationGaussian Elimination for Linear Systems
Gaussian Elimination for Linear Systems Tsung-Ming Huang Department of Mathematics National Taiwan Normal University October 3, 2011 1/56 Outline 1 Elementary matrices 2 LR-factorization 3 Gaussian elimination
More information3 QR factorization revisited
LINEAR ALGEBRA: NUMERICAL METHODS. Version: August 2, 2000 30 3 QR factorization revisited Now we can explain why A = QR factorization is much better when using it to solve Ax = b than the A = LU factorization
More informationProcess Model Formulation and Solution, 3E4
Process Model Formulation and Solution, 3E4 Section B: Linear Algebraic Equations Instructor: Kevin Dunn dunnkg@mcmasterca Department of Chemical Engineering Course notes: Dr Benoît Chachuat 06 October
More informationAlgebraic Equations. 2.0 Introduction. Nonsingular versus Singular Sets of Equations. A set of linear algebraic equations looks like this:
Chapter 2. 2.0 Introduction Solution of Linear Algebraic Equations A set of linear algebraic equations looks like this: a 11 x 1 + a 12 x 2 + a 13 x 3 + +a 1N x N =b 1 a 21 x 1 + a 22 x 2 + a 23 x 3 +
More informationOn the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination
On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination J.M. Peña 1 Introduction Gaussian elimination (GE) with a given pivoting strategy, for nonsingular matrices
More informationMath 22AL Lab #4. 1 Objectives. 2 Header. 0.1 Notes
Math 22AL Lab #4 0.1 Notes Green typewriter text represents comments you must type. Each comment is worth one point. Blue typewriter text represents commands you must type. Each command is worth one point.
More informationMTH 464: Computational Linear Algebra
MTH 464: Computational Linear Algebra Lecture Outlines Exam 2 Material Prof. M. Beauregard Department of Mathematics & Statistics Stephen F. Austin State University February 6, 2018 Linear Algebra (MTH
More informationLinear Equations and Matrix
1/60 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Gaussian Elimination 2/60 Alpha Go Linear algebra begins with a system of linear
More informationNumerical Methods I: Numerical linear algebra
1/3 Numerical Methods I: Numerical linear algebra Georg Stadler Courant Institute, NYU stadler@cimsnyuedu September 1, 017 /3 We study the solution of linear systems of the form Ax = b with A R n n, x,
More informationDense LU factorization and its error analysis
Dense LU factorization and its error analysis Laura Grigori INRIA and LJLL, UPMC February 2016 Plan Basis of floating point arithmetic and stability analysis Notation, results, proofs taken from [N.J.Higham,
More informationSection 3.5 LU Decomposition (Factorization) Key terms. Matrix factorization Forward and back substitution LU-decomposition Storage economization
Section 3.5 LU Decomposition (Factorization) Key terms Matrix factorization Forward and back substitution LU-decomposition Storage economization In matrix analysis as implemented in modern software the
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationLU factorization with Panel Rank Revealing Pivoting and its Communication Avoiding version
1 LU factorization with Panel Rank Revealing Pivoting and its Communication Avoiding version Amal Khabou Advisor: Laura Grigori Université Paris Sud 11, INRIA Saclay France SIAMPP12 February 17, 2012 2
More informationThe following steps will help you to record your work and save and submit it successfully.
MATH 22AL Lab # 4 1 Objectives In this LAB you will explore the following topics using MATLAB. Properties of invertible matrices. Inverse of a Matrix Explore LU Factorization 2 Recording and submitting
More information