Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
|
|
- Lynn Fox
- 6 years ago
- Views:
Transcription
1 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.1/19
2 Overview Two Real World Problems: 1. Model of the blood circulation in a human body 2. Google TM Matrix - determining PageRank TM Introduction to Markov chains Problem setting and solution approaches Introduction to algebraic multigrid (AMG) Some numerical tests Summary and outlook Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.2/19
3 Real World Problem: Biotech Blood Circulation in a Human Body Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.3/19
4 Real World Problem: Biotech Blood Circulation in a Human Body Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.3/19
5 Real World Problem: Biotech Blood Circulation in a Human Body (I T T )x = 0 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.3/19
6 Example of a blood flow matrix Entry t ij represents the probability that a blood particle passes from organ i to organ j within one time step BloodFlow Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.4/19
7 Google TM Matrix Example: S1 S2 S3 The link structure of the Internet is transformed into a column stochastic matrix as follows: S4 We search the web from a starting page following the rules: 1. With probability 0.85 follow a link on the page 2. With probability 0.15 choose a random site Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.5/19
8 Google TM Matrix Example: S1 S2 S3 The link structure of the Internet is transformed into a column stochastic matrix as follows: T T = S Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.5/19
9 Google TM Matrix Example: S1 S2 S3 The link structure of the Internet is transformed into a column stochastic matrix as follows: S4 The eigenvector x to the eigenvalue 1 is the PageRank TM. T T x = x Here, the solution is x = (0.26, 0.29, 0.12, 0.33) T. S4 is the site with the highest PageRank. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.5/19
10 Introduction to Markov Chains [Markov Chain] (Named after Andrei Markov) A process is called finite homogeneous Markov chain (MC) if s 1,..,s n are the states of the process P[s i s j ] =: t ij is time independent Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.6/19
11 Introduction to Markov Chains [Markov Chain] (Named after Andrei Markov) A process is called finite homogeneous Markov chain (MC) if s 1,..,s n are the states of the process P[s i s j ] =: t ij is time independent The matrix T = (t ij ) i,j=1,..,n satisfies t ij 0, n t ij = 1, i,j = 1,..n j=1 and is called (state) transition matrix. In many applications T is singular. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.6/19
12 Probability Distribution The probability distribution vector x k = (x k i ) is defined by x k i = probability that the system is in state s i after k steps Properties: x k i 0 and n i=1 xk i = 1 for each k. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.7/19
13 Probability Distribution The probability distribution vector x k = (x k i ) is defined by x k i = probability that the system is in state s i after k steps Properties: x k i 0 and n i=1 xk i = 1 for each k. Q.: What is the probability distribution at time k? Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.7/19
14 Probability Distribution The probability distribution vector x k = (x k i ) is defined by x k i = probability that the system is in state s i after k steps Properties: x k i 0 and n i=1 xk i = 1 for each k. Q.: What is the probability distribution at time k? A: Let M = (T,x 0 ). Then for each k > 0 (x k ) T = (x k 1 ) T T = (x 0 ) T T k. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.7/19
15 Is there a stable state? A distribution vector x is said to be stationary if x T T = x T. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
16 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
17 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Iterate (x k ) T = (x k 1 ) T T = (x 0 ) T T k If x := lim k (x 0 ) T T k exists, then it is stationary. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
18 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Iterate (x k ) T = (x k 1 ) T T = (x 0 ) T T k If x := lim k (x 0 ) T T k exists, then it is stationary. Problem: Does lim k (x 0 ) T T k always exist? Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
19 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Iterate (x k ) T = (x k 1 ) T T = (x 0 ) T T k If x := lim k (x 0 ) T T k exists, then it is stationary. Problem: Does lim k (x 0 ) T T k always exist? Answer:NO! Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
20 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Iterate (x k ) T = (x k 1 ) T T = (x 0 ) T T k If x := lim (x 0 ) T T k exists, then it is stationary. k Problem: Does lim ] [ ] 1 0 Example: T = (x 0 ) T T k always exist? Answer:NO! [ k [ ] , T =, T =. But 1 0 ] [ ] = [ ] there is a stable distribution [ 1 2 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
21 Is there a stable state? A distribution vector x is said to be stationary if How to compute it? x T T = x T. Iterate (x k ) T = (x k 1 ) T T = (x 0 ) T T k If x := lim k (x 0 ) T T k exists, then it is stationary. Problem: Does lim k (x 0 ) T T k always exist? Answer:NO! Transform fixed point equation into a linear system: Ax := (I T T )x = 0 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.8/19
22 How do we solve (I T T )x = 0? Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
23 How do we solve (I T T )x = 0? 1. Solvability of a singular system: Ax = b is solvable if b im(a). The system is then called consistent. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
24 How do we solve (I T T )x = 0? 1. Solvability of a singular system: Ax = b is solvable if b im(a). The system is then called consistent. our system is consistent! Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
25 How do we solve (I T T )x = 0? 1. Solvability of a singular system: Ax = b is solvable if b im(a). The system is then called consistent. our system is consistent! 2. Possible approaches to solve Ax = (I T T )x = 0: Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
26 How do we solve (I T T )x = 0? 1. Solvability of a singular system: Ax = b is solvable if b im(a). The system is then called consistent. our system is consistent! 2. Possible approaches to solve Ax = (I T T )x = 0: Direct methods (e.g. Gauss, LU) Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
27 How do we solve (I T T )x = 0? 1. Solvability of a singular system: Ax = b is solvable if b im(a). The system is then called consistent. our system is consistent! 2. Possible approaches to solve Ax = (I T T )x = 0: Direct methods (e.g. Gauss, LU) Iterative methods Splitting methods (e.g. Jacobi, SOR) Krylov subspace methods (e.g. CG, GMRES) AMG Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.9/19
28 Direct solvers Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.10/19
29 Direct solvers LU works despite singularity if A is an M-Matrix Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.10/19
30 Direct solvers LU works despite singularity if A is an M-Matrix A of the form A = si B, s > 0, B 0 where s ρ(b) is called M-Matrix. If s > ρ(b) then A is nonsingular M-Matrix s = ρ(b) then A is singular M-Matrix Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.10/19
31 Direct solvers LU works despite singularity if A is an M-Matrix A of the form A = si B, s > 0, B 0 where s ρ(b) is called M-Matrix. If s > ρ(b) then A is nonsingular M-Matrix s = ρ(b) then A is singular M-Matrix A = (I T T ) is a singular M-Matrix as ρ(t) = 1 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.10/19
32 Direct solvers LU works despite singularity if A is an M-Matrix A of the form A = si B, s > 0, B 0 where s ρ(b) is called M-Matrix. If s > ρ(b) then A is nonsingular M-Matrix s = ρ(b) then A is singular M-Matrix A = (I T T ) is a singular M-Matrix as ρ(t) = 1 BUT: high computational effort - only useful for small n Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.10/19
33 Krylov subspace methods e.g GMRES Let A R n n, r R n. Then K m := {r,ar,a 2 r,..,a m 1 r} is the m-dimensional Krylov space. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.11/19
34 Krylov subspace methods e.g GMRES Let A R n n, r R n. Then K m := {r,ar,a 2 r,..,a m 1 r} is the m-dimensional Krylov space. Basic Idea of GMRES for solving Ax = b: Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.11/19
35 Krylov subspace methods e.g GMRES Let A R n n, r R n. Then K m := {r,ar,a 2 r,..,a m 1 r} is the m-dimensional Krylov space. Basic Idea of GMRES for solving Ax = b: Let r = b Ax 0 = Ax 0 be the initial residual Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.11/19
36 Krylov subspace methods e.g GMRES Let A R n n, r R n. Then K m := {r,ar,a 2 r,..,a m 1 r} is the m-dimensional Krylov space. Basic Idea of GMRES for solving Ax = b: Let r = b Ax 0 = Ax 0 be the initial residual Choose x such that x = x 0 + V y where columns of V orthogonally span K m so that b Ax 0 2 = r 2 is minimized. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.11/19
37 AMG as a preconditioner in GMRES The concept of a preconditioner M an easily invertible matrix such that M 1 A 1 solve à x = b, where à = M 1 A, and b = M 1 b instead of the original system. Goal: make M 1 A I. Idea: take as M 1 the solution from AMG iteration. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.12/19
38 Introduction to AMG Define coarse grids: a sequence of smaller and smaller systems of equations A m u m A 1 = b m, where = A The interpolation operator P and the restriction operator R are defined as follows: Let A 1 = A R n n. Then R : R n R N P : R N R n. It is A 2 = RA 1 P and A 2 R N N. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.13/19
39 AMG structure The AMG consists of two main parts: Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
40 AMG structure The AMG consists of two main parts: setup phase solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
41 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C define grid transfer P,R define coarse grid operators A H = RA h P solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
42 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C Undirected graph of A define grid transfer P,R Coarse grid nodes Fine grid nodes define coarse grid operators A H = RA h P solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
43 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C define grid transfer P,R define coarse grid operators A H = RA h P solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
44 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C define grid transfer P,R Restriction: Interpolation: define coarse grid operators A H = RA h P solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
45 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C define grid transfer P,R define coarse grid operators A H = RA h P solution phase Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
46 AMG structure The AMG consists of two main parts: setup phase choose coarser grid F,C define grid transfer P,R define coarse grid operators A H = RA h P solution phase : V-cycle smoothing R smoothing R P coarse grid correction smoothing P smoothing Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.14/19
47 Recursive application Essentiality for setting up a multi-level cycle: Properties of the finest grid operator should carry over to the coarsest grid. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.15/19
48 Recursive application Essentiality for setting up a multi-level cycle: Properties of the finest grid operator should carry over to the coarsest grid. J.W.Ruge and K.Stueben ([4]) show this for symmetric, weakly diagonally dominant M-Matrices, Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.15/19
49 Recursive application Essentiality for setting up a multi-level cycle: Properties of the finest grid operator should carry over to the coarsest grid. J.W.Ruge and K.Stueben ([4]) show this for symmetric, weakly diagonally dominant M-Matrices, In our case the main properties are the row sum zero property and the singular M-matrix property. Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.15/19
50 Numerical tests 1.Blood circulation: GMRES AMG+GMRES AMG n nnz TIME ITS grid TIME ITS TIME , ,18 0,11 2 1, , ,44 0,06 2 1, , ,9 0,17 2 2, , ,02 0,17 2 3, , ,26 0,58 2 4,84 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.16/19
51 Numerical tests 2.Blood circulation - nearly decoupled: ǫ ǫ ǫ... ǫ ǫ ǫ ǫ ǫ is taken to be 10 6 in the following computation. ǫ Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.16/19
52 Numerical tests 2.Blood circulation - nearly decoupled: GMRES AMG+GMRES AMG n TIME ITS grid TIME ITS TIME 432 1, ,67 0,29 7 2, , ,73 0,46 7 4, , ,78 1,03 9 6, , ,40 6, , >317 > ,87 64, ,35 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.16/19
53 Numerical tests 3.Google TM Matrix: A 500 nodes example starting at This matrix represents the link structure. An entry t ij represents a link from page j to page i. The right eigenvector corresponding to the eigenvalue 1 is the PageRank nz = 5783 Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.16/19
54 Numerical tests 3.Google TM Matrix: A 500 nodes example starting at PageRank In Out URL Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.16/19
55 Summary and Outlook AMG+GMRES was tested on some Bloodflow matrices Nearly decoupled Markov chains (Bloodflow) Google TM Matrix Generally, AMG preconditioner accelerates the convergence of GMRES for larger matrices For smaller matrices and for random Markov chains, the preconditioner does not accelerate convergence To-do s: Convergence estimate for singular M-matrices Application to other areas Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.17/19
56 Bibliography [1] S. Bächle. Numerische Berechnung der stationären Wahrscheinlichkeiten von Markov Ketten mit dem GMRES Verfahren. Master s thesis, Universität Karlsruhe (TH), [2] A. Berman and R.J. Plemmons. Nonnegative Matrices in the Mathematical Sciences. Academic Press, New York, [3] I-wen Kuo. A note on factorizations of singular m-matrices. Linear Algebra and its Applications, 197/198: , [4] J.W. Ruge and K. Stueben. Algebraic multigrid. Multigrid Methods, pages , [5] Richard S. Varga and Da-Yong Cai. On the LU factorization of M-matrices. Numerische Mathematik, 38: , Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.18/19
57 List of Slides 1. Overview 2. Real World Problem: Biotech 3. Example of a blood flow matrix 4. Google TM Matrix 5. Introduction to Markov Chains 6. Probability Distribution 7. Is there a stable state? 8. How do we solve (I T T )x = 0? 9. Direct solvers 10. Krylov subspace methods e.g GMRES 11. AMG as a preconditioner in GMRES 12. Introduction to AMG 13. AMG structure 14. Recursive application 15. Numerical tests 16. Summary and Outlook 17. Bibliography 18. List of Slides Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes p.19/19
Algebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationLecture 18 Classical Iterative Methods
Lecture 18 Classical Iterative Methods MIT 18.335J / 6.337J Introduction to Numerical Methods Per-Olof Persson November 14, 2006 1 Iterative Methods for Linear Systems Direct methods for solving Ax = b,
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 informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationAggregation Algorithms for K-cycle Aggregation Multigrid for Markov Chains
Aggregation Algorithms for K-cycle Aggregation Multigrid for Markov Chains by Manda Winlaw A research paper presented to the University of Waterloo in partial fulfillment of the requirements for the degree
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationBootstrap AMG. Kailai Xu. July 12, Stanford University
Bootstrap AMG Kailai Xu Stanford University July 12, 2017 AMG Components A general AMG algorithm consists of the following components. A hierarchy of levels. A smoother. A prolongation. A restriction.
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationMarkov Chains and Web Ranking: a Multilevel Adaptive Aggregation Method
Markov Chains and Web Ranking: a Multilevel Adaptive Aggregation Method Hans De Sterck Department of Applied Mathematics, University of Waterloo Quoc Nguyen; Steve McCormick, John Ruge, Tom Manteuffel
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationRestricted Additive Schwarz Methods for Markov Chains
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 23 [Version: 2002/09/18 v1.02] Restricted Additive Schwarz Methods for Markov Chains Michele Benzi Verena Kuhlemann Department
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 informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationThe Conjugate Gradient Method
The Conjugate Gradient Method Classical Iterations We have a problem, We assume that the matrix comes from a discretization of a PDE. The best and most popular model problem is, The matrix will be as large
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationc 2010 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 32, No. 1, pp. 40 61 c 2010 Society for Industrial and Applied Mathematics SMOOTHED AGGREGATION MULTIGRID FOR MARKOV CHAINS H. DE STERCK, T. A. MANTEUFFEL, S. F. MCCORMICK, K.
More informationAggregation-based Adaptive Algebraic Multigrid for Sparse Linear Systems. Eran Treister
Aggregation-based Adaptive Algebraic Multigrid for Sparse Linear Systems Eran Treister Aggregation-based Adaptive Algebraic Multigrid for Sparse Linear Systems Research Thesis In Partial Fulfillment of
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationA hybrid reordered Arnoldi method to accelerate PageRank computations
A hybrid reordered Arnoldi method to accelerate PageRank computations Danielle Parker Final Presentation Background Modeling the Web The Web The Graph (A) Ranks of Web pages v = v 1... Dominant Eigenvector
More informationA NEW EFFECTIVE PRECONDITIONED METHOD FOR L-MATRICES
Journal of Mathematical Sciences: Advances and Applications Volume, Number 2, 2008, Pages 3-322 A NEW EFFECTIVE PRECONDITIONED METHOD FOR L-MATRICES Department of Mathematics Taiyuan Normal University
More informationIterative methods for Linear System
Iterative methods for Linear System JASS 2009 Student: Rishi Patil Advisor: Prof. Thomas Huckle Outline Basics: Matrices and their properties Eigenvalues, Condition Number Iterative Methods Direct and
More informationMULTILEVEL ADAPTIVE AGGREGATION FOR MARKOV CHAINS, WITH APPLICATION TO WEB RANKING
MULTILEVEL ADAPTIVE AGGREGATION FOR MARKOV CHAINS, WITH APPLICATION TO WEB RANKING H. DE STERCK, THOMAS A. MANTEUFFEL, STEPHEN F. MCCORMICK, QUOC NGUYEN, AND JOHN RUGE Abstract. A multilevel adaptive aggregation
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 informationScientific Computing WS 2018/2019. Lecture 9. Jürgen Fuhrmann Lecture 9 Slide 1
Scientific Computing WS 2018/2019 Lecture 9 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 9 Slide 1 Lecture 9 Slide 2 Simple iteration with preconditioning Idea: Aû = b iterative scheme û = û
More information9.1 Preconditioned Krylov Subspace Methods
Chapter 9 PRECONDITIONING 9.1 Preconditioned Krylov Subspace Methods 9.2 Preconditioned Conjugate Gradient 9.3 Preconditioned Generalized Minimal Residual 9.4 Relaxation Method Preconditioners 9.5 Incomplete
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccs Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary c 2008,2010
More informationSolving Sparse Linear Systems: Iterative methods
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 2: Iterative Methods Dianne P. O Leary
More informationUsing an Auction Algorithm in AMG based on Maximum Weighted Matching in Matrix Graphs
Using an Auction Algorithm in AMG based on Maximum Weighted Matching in Matrix Graphs Pasqua D Ambra Institute for Applied Computing (IAC) National Research Council of Italy (CNR) pasqua.dambra@cnr.it
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 informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationConditioning of the Entries in the Stationary Vector of a Google-Type Matrix. Steve Kirkland University of Regina
Conditioning of the Entries in the Stationary Vector of a Google-Type Matrix Steve Kirkland University of Regina June 5, 2006 Motivation: Google s PageRank algorithm finds the stationary vector of a stochastic
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationIterative methods for Linear System of Equations. Joint Advanced Student School (JASS-2009)
Iterative methods for Linear System of Equations Joint Advanced Student School (JASS-2009) Course #2: Numerical Simulation - from Models to Software Introduction In numerical simulation, Partial Differential
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2017/18 Part 3: Iterative Methods PD
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 informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More informationAccepted Manuscript. Multigrid Algorithm from Cyclic Reduction for Markovian Queueing Networks. Shu-Ling Yang, Jian-Feng Cai, Hai-Wei Sun
Accepted Manuscript Multigrid Algorithm from Cyclic Reduction for Markovian Queueing Networks Shu-Ling Yang, Jian-Feng Cai, Hai-Wei Sun PII: S0096-3003()0050-0 DOI: 006/jamc20008 Reference: AMC 5730 To
More informationJordan Journal of Mathematics and Statistics (JJMS) 5(3), 2012, pp A NEW ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS OF EQUATIONS
Jordan Journal of Mathematics and Statistics JJMS) 53), 2012, pp.169-184 A NEW ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS OF EQUATIONS ADEL H. AL-RABTAH Abstract. The Jacobi and Gauss-Seidel iterative
More informationKrylov Subspace Methods to Calculate PageRank
Krylov Subspace Methods to Calculate PageRank B. Vadala-Roth REU Final Presentation August 1st, 2013 How does Google Rank Web Pages? The Web The Graph (A) Ranks of Web pages v = v 1... Dominant Eigenvector
More informationJACOBI S ITERATION METHOD
ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes
More informationCAAM 454/554: Stationary Iterative Methods
CAAM 454/554: Stationary Iterative Methods Yin Zhang (draft) CAAM, Rice University, Houston, TX 77005 2007, Revised 2010 Abstract Stationary iterative methods for solving systems of linear equations are
More informationeigenvalues, markov matrices, and the power method
eigenvalues, markov matrices, and the power method Slides by Olson. Some taken loosely from Jeff Jauregui, Some from Semeraro L. Olson Department of Computer Science University of Illinois at Urbana-Champaign
More informationMultigrid Methods for Linear Systems with Stochastic Entries Arising in Lattice QCD. Andreas Frommer
Methods for Linear Systems with Stochastic Entries Arising in Lattice QCD Andreas Frommer Collaborators The Dirac operator James Brannick, Penn State University Björn Leder, Humboldt Universität Berlin
More informationComputational Economics and Finance
Computational Economics and Finance Part II: Linear Equations Spring 2016 Outline Back Substitution, LU and other decomposi- Direct methods: tions Error analysis and condition numbers Iterative methods:
More information0.1 Naive formulation of PageRank
PageRank is a ranking system designed to find the best pages on the web. A webpage is considered good if it is endorsed (i.e. linked to) by other good webpages. The more webpages link to it, and the more
More informationUniversity of Illinois at Urbana-Champaign. Multigrid (MG) methods are used to approximate solutions to elliptic partial differential
Title: Multigrid Methods Name: Luke Olson 1 Affil./Addr.: Department of Computer Science University of Illinois at Urbana-Champaign Urbana, IL 61801 email: lukeo@illinois.edu url: http://www.cs.uiuc.edu/homes/lukeo/
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationConjugate Gradients: Idea
Overview Steepest Descent often takes steps in the same direction as earlier steps Wouldn t it be better every time we take a step to get it exactly right the first time? Again, in general we choose a
More informationLab 8: Measuring Graph Centrality - PageRank. Monday, November 5 CompSci 531, Fall 2018
Lab 8: Measuring Graph Centrality - PageRank Monday, November 5 CompSci 531, Fall 2018 Outline Measuring Graph Centrality: Motivation Random Walks, Markov Chains, and Stationarity Distributions Google
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 68 (2014) 1151 1160 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A GPU
More informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationEQUIVALENCE OF CONDITIONS FOR CONVERGENCE OF ITERATIVE METHODS FOR SINGULAR LINEAR SYSTEMS
EQUIVALENCE OF CONDITIONS FOR CONVERGENCE OF ITERATIVE METHODS FOR SINGULAR LINEAR SYSTEMS DANIEL B. SZYLD Department of Mathematics Temple University Philadelphia, Pennsylvania 19122-2585 USA (szyld@euclid.math.temple.edu)
More informationCourse Notes: Week 1
Course Notes: Week 1 Math 270C: Applied Numerical Linear Algebra 1 Lecture 1: Introduction (3/28/11) We will focus on iterative methods for solving linear systems of equations (and some discussion of eigenvalues
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 21: Sensitivity of Eigenvalues and Eigenvectors; Conjugate Gradient Method Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis
More informationMath 1553, Introduction to Linear Algebra
Learning goals articulate what students are expected to be able to do in a course that can be measured. This course has course-level learning goals that pertain to the entire course, and section-level
More informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
More informationA Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations
A Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations Jin Yun Yuan Plamen Y. Yalamov Abstract A method is presented to make a given matrix strictly diagonally dominant
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationAdaptive Multigrid for QCD. Lattice University of Regensburg
Lattice 2007 University of Regensburg Michael Clark Boston University with J. Brannick, R. Brower, J. Osborn and C. Rebbi -1- Lattice 2007, University of Regensburg Talk Outline Introduction to Multigrid
More informationIntroduction to Scientific Computing
(Lecture 5: Linear system of equations / Matrix Splitting) Bojana Rosić, Thilo Moshagen Institute of Scientific Computing Motivation Let us resolve the problem scheme by using Kirchhoff s laws: the algebraic
More informationSolving Ax = b, an overview. Program
Numerical Linear Algebra Improving iterative solvers: preconditioning, deflation, numerical software and parallelisation Gerard Sleijpen and Martin van Gijzen November 29, 27 Solving Ax = b, an overview
More informationEigenvalues and Eigenvectors
November 3, 2016 1 Definition () The (complex) number λ is called an eigenvalue of the n n matrix A provided there exists a nonzero (complex) vector v such that Av = λv, in which case the vector v is called
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 information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
More informationComputers and Mathematics with Applications. Convergence analysis of the preconditioned Gauss Seidel method for H-matrices
Computers Mathematics with Applications 56 (2008) 2048 2053 Contents lists available at ScienceDirect Computers Mathematics with Applications journal homepage: wwwelseviercom/locate/camwa Convergence analysis
More informationA Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator
A Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator Artur Strebel Bergische Universität Wuppertal August 3, 2016 Joint Work This project is joint work with: Gunnar
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
More information7.3 The Jacobi and Gauss-Siedel Iterative Techniques. Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP.
7.3 The Jacobi and Gauss-Siedel Iterative Techniques Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP. 7.3 The Jacobi and Gauss-Siedel Iterative Techniques
More informationLecture 17: Iterative Methods and Sparse Linear Algebra
Lecture 17: Iterative Methods and Sparse Linear Algebra David Bindel 25 Mar 2014 Logistics HW 3 extended to Wednesday after break HW 4 should come out Monday after break Still need project description
More informationIntroduction to Scientific Computing II Multigrid
Introduction to Scientific Computing II Multigrid Miriam Mehl Slide 5: Relaxation Methods Properties convergence depends on method clear, see exercises and 3), frequency of the error remember eigenvectors
More informationMATH36001 Perron Frobenius Theory 2015
MATH361 Perron Frobenius Theory 215 In addition to saying something useful, the Perron Frobenius theory is elegant. It is a testament to the fact that beautiful mathematics eventually tends to be useful,
More informationSplitting Iteration Methods for Positive Definite Linear Systems
Splitting Iteration Methods for Positive Definite Linear Systems Zhong-Zhi Bai a State Key Lab. of Sci./Engrg. Computing Inst. of Comput. Math. & Sci./Engrg. Computing Academy of Mathematics and System
More informationGoogle Page Rank Project Linear Algebra Summer 2012
Google Page Rank Project Linear Algebra Summer 2012 How does an internet search engine, like Google, work? In this project you will discover how the Page Rank algorithm works to give the most relevant
More informationEigenvalue Problems Computation and Applications
Eigenvalue ProblemsComputation and Applications p. 1/36 Eigenvalue Problems Computation and Applications Che-Rung Lee cherung@gmail.com National Tsing Hua University Eigenvalue ProblemsComputation and
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationModified Gauss Seidel type methods and Jacobi type methods for Z-matrices
Linear Algebra and its Applications 7 (2) 227 24 www.elsevier.com/locate/laa Modified Gauss Seidel type methods and Jacobi type methods for Z-matrices Wen Li a,, Weiwei Sun b a Department of Mathematics,
More informationLab 1: Iterative Methods for Solving Linear Systems
Lab 1: Iterative Methods for Solving Linear Systems January 22, 2017 Introduction Many real world applications require the solution to very large and sparse linear systems where direct methods such as
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More information1. Let m 1 and n 1 be two natural numbers such that m > n. Which of the following is/are true?
. Let m and n be two natural numbers such that m > n. Which of the following is/are true? (i) A linear system of m equations in n variables is always consistent. (ii) A linear system of n equations in
More informationMonte Carlo Method for Finding the Solution of Dirichlet Partial Differential Equations
Applied Mathematical Sciences, Vol. 1, 2007, no. 10, 453-462 Monte Carlo Method for Finding the Solution of Dirichlet Partial Differential Equations Behrouz Fathi Vajargah Department of Mathematics Guilan
More informationIterative Solution methods
p. 1/28 TDB NLA Parallel Algorithms for Scientific Computing Iterative Solution methods p. 2/28 TDB NLA Parallel Algorithms for Scientific Computing Basic Iterative Solution methods The ideas to use iterative
More informationADAPTIVE ALGEBRAIC MULTIGRID
ADAPTIVE ALGEBRAIC MULTIGRID M. BREZINA, R. FALGOUT, S. MACLACHLAN, T. MANTEUFFEL, S. MCCORMICK, AND J. RUGE Abstract. Efficient numerical simulation of physical processes is constrained by our ability
More informationOnline Social Networks and Media. Link Analysis and Web Search
Online Social Networks and Media Link Analysis and Web Search How to Organize the Web First try: Human curated Web directories Yahoo, DMOZ, LookSmart How to organize the web Second try: Web Search Information
More informationPreface to Second Edition... vii. Preface to First Edition...
Contents Preface to Second Edition..................................... vii Preface to First Edition....................................... ix Part I Linear Algebra 1 Basic Vector/Matrix Structure and
More informationFEM and sparse linear system solving
FEM & sparse linear system solving, Lecture 9, Nov 19, 2017 1/36 Lecture 9, Nov 17, 2017: Krylov space methods http://people.inf.ethz.ch/arbenz/fem17 Peter Arbenz Computer Science Department, ETH Zürich
More informationGoogle PageRank. Francesco Ricci Faculty of Computer Science Free University of Bozen-Bolzano
Google PageRank Francesco Ricci Faculty of Computer Science Free University of Bozen-Bolzano fricci@unibz.it 1 Content p Linear Algebra p Matrices p Eigenvalues and eigenvectors p Markov chains p Google
More informationAn Empirical Comparison of Graph Laplacian Solvers
An Empirical Comparison of Graph Laplacian Solvers Kevin Deweese 1 Erik Boman 2 John Gilbert 1 1 Department of Computer Science University of California, Santa Barbara 2 Scalable Algorithms Department
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More information6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities
6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities 1 Outline Outline Dynamical systems. Linear and Non-linear. Convergence. Linear algebra and Lyapunov functions. Markov
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationINTRODUCTION TO MULTIGRID METHODS
INTRODUCTION TO MULTIGRID METHODS LONG CHEN 1. ALGEBRAIC EQUATION OF TWO POINT BOUNDARY VALUE PROBLEM We consider the discretization of Poisson equation in one dimension: (1) u = f, x (0, 1) u(0) = u(1)
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More information