Computational Linear Algebra
|
|
- Junior Tyler
- 5 years ago
- Views:
Transcription
1 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
2 Part 4: Iterative Methods PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 2
3 overview definitions splitting methods projection and KRYLOV subspace methods multigrid methods PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 3
4 some considerations let s first consider conventional iterative or relaxation methods, hence let Au = f (4.3.1) denote a linear system derived from a finite-difference discretisation of a one-dimensional second-order boundary value problem u (x) = f(x), 0 < x < 1 with u(0) = u(1) = 0 A = PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 4
5 some considerations (cont d) let u denote the exact solution of (4.3.1) and v an approximation to the exact solution, then the error (or algebraic error) is given by e = u v the error is a vector and its magnitude measured by the maximum norm or the EUCLEDIAN norm defined as e := max e i and e 2 := i = 1...n unfortunately, the error is just as inaccessible as the exact solution itself however, a computable measure of how well v approximates u is the residual, given by r = f Av PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 5
6 some considerations (cont d) by uniqueness of the solution, r = 0 if and only if e = 0 from (4.3.1), we can derive the relation between error and residual Ae = r (4.3.2) which is called the residual equation and plays a vital role in multigrid methods by solving (4.3.2) for e we can compute a new approximation using the definition of the error u = v + e for further considerations, it is sufficient to work with the homogeneous linear system Au = 0 using arbitrary initial guesses in order to start a JACOBI relaxation method as the exact solution is known (u = 0), the error in some approximation v is simply v PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 6
7 some considerations (cont d) the one-dimensional model problem with f = 0 appears as u 0 u 1 u 2 u n 1 u n h u i 1 + 2u i u i+1 = 0, 1 i n 1 u 0 = u n = 0, with grid points x i = ih, where h = 1/n to obtain some valuable insight, we apply various iterations to this system with an initial guess consisting of the vectors (or FOURIER modes) m i = sin, 0 i n, 1 k n 1, where integer k is called wave number (or frequency) indicating the number of half sine waves that constitute m on the problem domain PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 7
8 some considerations (cont d) let m k denote the entire vector m with wave number k, shown below are initial guesses m 1, m 3, and m 6 k = 1 k = 3 k = 6 small values of k belong to long, smooth waves while large values of k correspond to highly oscillatory waves PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 8
9 some considerations (cont d) now applying a JACOBI relaxation method with ω = 2/3 to the model problem on a grid with n = 64 points for initial guesses m 1, m 3, and m 6 the iteration is applied 100 times 1,00 k = 1 0,75 error 0,50 k = 3 0,25 0, k = 6 recall that the error is just m iterations PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 9
10 some considerations (cont d) plotting log of e for the JACOBI relaxation method with ω = 2/3 shows a linear decrease in the log of the error norm, indicating that the error decreases geometrically with each iteration 1,00 k = 1 k = 3 log error 0,10 k = 6 0, iterations PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 10
11 some considerations (cont d) in general, most initial guesses would not consist of a single mode hence, a more realistic situation is given by an overlay of two modes: a low-frequency wave (k = 2) and a high frequency wave (k = 16) due to m i = on a grid with n = 64 points 1,00 0,75 error 0,50 0,25 (m 2 + m 16 )/2 0, iterations PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 11
12 some considerations (cont d) live demo playing with coins source: wn.de PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 12
13 analytical approach discussed iteration method can be represented in the (simplified) form v 1 = Rv 0 + g (4.3.3) where R denotes the iteration matrix furthermore, the exact solution u is a fixed-point of (4.3.3), that means u = Ru + g (4.3.4) subtracting (4.3.3) from (4.3.4) yields e 1 = Re 0 repeating this argument, it follows that after m relaxation steps the error in m-th approximation is given by e m = R m e 0 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 13
14 analytical approach (cont d) if we choose a particular vector norm (and its associated matrix norm), it is possible to bound the error after m iterations by e m R m e 0 as for relaxation methods we assume ρ(r) < 1, the error is forced to zero as the iteration proceeds considering the weighted JACOBI iteration applied to the one-dimensional model problem, we have R ω = (1 ω)i + ω R J with R ω = I PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 14
15 analytical approach (cont d) written in this form follows that eigenvalues of R ω and A are related by σ (R ω ) = 1 ω 2 σ (A), with σ (A) ={4 sin 2 }, 1 k n 1 let w k,j denote the j-th component of k-th eigenvector w k, then with eigenvectors of A given by w k,j = sin, 1 k n 1, 0 j n, we see that eigenvectors of A are simply the FOURIER modes discussed earlier in this section PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 15
16 analytical approach (cont d) with these results, we find that eigenvalues of R ω are λ k (R ω ) = 1 2ω sin 2, 1 k n 1, while eigenvectors of R ω are the same as eigenvectors of A it is important to note that only for 0 <ω 1 follows λ k (Rω) <1 and the weighted JACOBI method converges let e 0 be the error in an initial guess used in the weighted JACOBI method, then it is possible to expand e 0 using the eigenvectors of A in the form e 0 = c k w k, where coefficients c k error denote a weighting factor for each mode in the PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 16
17 analytical approach (cont d) we have seen that after m iterations, the error is given by e m = Rm ω e 0 using the eigenvector expansion for e 0, we have R ω w k = λ k (R ω )w k m e m = R ω e 0 = c k Rm m ω w k = c k λ k (R ω )w k this expansion for e m shows that after m iterations, the k-th mode of the m initial error has been reduced by a factor of λ k (R ω ) it should further be noted that the weighted JACOBI method does not mix modes, i.e. amplitudes of single modes can change, but they cannot be converted into different modes we would like to the find the optimal value of ω with 0 < ω 1 that makes λ k (R ω ) as small as possible for all 1 k n 1 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 17
18 analytical approach (cont d) eigenvalues λ k = 1 2ω sin 2 of iteration matrix R ω for ω = ⅓, ½, ⅔, 1 λ k (R ω ) 1 ω = 1/3 ω = 1/2 n/2 n k ω = 2/3 1 ω = 1 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 18
19 analytical approach (cont d) for all values of ω satisfying 0 < ω 1 applies λ 1 = 1 2ω sin 2 = 1 2ω sin 2 1 this implies that λ 1, the eigenvalue associated with the smoothest mode, will always be close to 1, hence no value of ω will reduce the smooth components of the error effectively even worse, the smaller the grid spacing h, the closer λ 1 is to 1, thus deteriorating convergence as no value of ω damps the smooth components satisfactorily, we are looking for values that provide the best damping of the oscillatory components (those with n/2 k n 1) solving condition λ n/2 (R ω ) = λ n (R ω ) leads to the optimal value of ω = ⅔ PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 19
20 analytical approach (cont d) again, we consider the weighted JACOBI method applied to the onedimensional model problem Au = 0 on a grid with n = 64 points for initial guesses consisting of single modes m k (1 k 63) only, the figures show the number of iterations required to reduce the norm of the initial error by a factor of 100 for weighting factors of ω = 1 and ω = ⅔ iterations ω=1 iterations ω=⅔ wave number k wave number k PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 20
21 analytical approach (cont d) another perspective of this convergence behaviour is provided below plotted are approximations for different initial guesses consisting of m 3, m 16, and (m 2 + m 16 )/2 shown after one and after 10 relaxation steps initial guess m 3 initial guess m 16 initial guess (m 2 + m 16 )/2 after 10 relax. steps after 10 relax. steps after 10 relax. steps PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 21
22 elements of multigrid how do smooth components look like on coarser grids? consider some fine (Ω h ) and coarse (Ω 2h ) grid with double grid spacing given some smooth wave with k = 4 on Ω h with n = 12 points Ω 2h representation with n = 6 points via direct projection wave becomes oscillatory on Ω 2h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 22
23 elements of multigrid (cont d) consequence: passing from the fine to the coarse grid, a mode becomes more oscillatory this is true provided that 1 k < n/2; the k = n/2 mode on Ω h becomes the zero vector on Ω 2h idea: when relaxation begins to stall, signalling the predominance of smooth error modes, move to a coarser grid as smooth error modes appear oscillatory there this leads to the following strategy, a procedure which is the basis of what is called the correction scheme relax on Au = f on Ω h to obtain an approximation v h compute residual r = f Av h relax on Ae = r on Ω 2h to obtain an approximation to the error e 2h correct v h v h + e 2h on Ω h with error estimate e 2h obtained on Ω 2h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 23
24 elements of multigrid (cont d) question: how to transfer residual r h from Ω h to Ω 2h (called restriction) and how to transfer the error estimate e 2h back from Ω 2h to Ω h (called interpolation or prolongation)? for further discussions, we consider only cases in which the coarse grid has twice the grid spacing of the next finest one ( universal practice, as there is no advantage in using grid spacings with ratios other than 2) the simplest prolongation method is quite effective and, thus, used for most multigrid purposes here the linear prolongation operator takes coarse-grid vectors and produces fine-grid vectors according to v 2h = v h, where, 0 j n/2 1 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 24
25 elements of multigrid (cont d) Ω h Ω 2h for n = 7 the linear prolongation has the form PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 25
26 elements of multigrid (cont d) quality of the prolongation (I) assume, the real error is a smooth vector on the fine grid Ω h and the coarse-grid approximation to the error obtained on Ω 2h is exact at the coarse grid points hence, when prolongated to the fine grid, the interpolant is also smooth if the exact error on Ω h (indicated by and ) is smooth, an interpolant of the coarse-grid error e 2h (solid line connecting points) should give a good representation of the exact error PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 26
27 elements of multigrid (cont d) quality of the prolongation (II) by contrast, if the real error is oscillatory on the fine grid Ω h, even a very good coarse-grid approximation may produce an interpolant that is not very accurate thus, interpolation is most effective when the error is smooth (in contrast to relaxation, which is most effective when the error is oscillatory) if the exact error on Ω h (indicated by and ) is oscillatory, an interpolant of the coarse-grid error e 2h (solid line connecting points) may give a poor representation of the exact error PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 27
28 elements of multigrid (cont d) second class of intergrid transfer operations restriction operator involves moving vectors from finer to coarser grids the most obvious restriction operator is injection, which is defined by v h = v 2h, where Ω h Ω 2h trivial restriction called injection PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 28
29 elements of multigrid (cont d) for n = 7, the restriction (injection) has the form an alternative restriction operator, called full weighting, is defined by, 1 j n/2 1 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 29
30 elements of multigrid (cont d) Ω h Ω 2h for n = 7, the restriction (full weighting) has the form PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 30
31 variational properties one reason for choosing full weighting as restriction operator is the fact, c, called variational property when expressing the coarse grid problem, A 2h v 2h = f 2h, we assume that A 2h is the Ω 2h version of the original operator A h, i.e. the result of discretising the problem on Ω 2h question: how to generate A 2h...? let us assume that the error e h = u h v h of some approximation v h lies entirely in the range of interpolation, denoted as this means that for some vector v 2h Ω 2h, e h = v 2h therefore, the residual equation on Ω h may be written A h e h = A h v 2h = r h (4.3.5) PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 31
32 variational properties (cont d) in (4.3.5), A h acts on a vector that lies entirely in the range of interpolation but how does A h act on let v 2h be some arbitrary vector on Ω 2h v 2h Ω 2h Ω h v 2h A h v 2h Ω h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 32
33 variational properties (cont d) we may conclude, odd rows of A h in (4.3.5) are zero and even rows correspond to the coarse-grid points of Ω 2h therefore, we can find a coarse-grid form of the residual equation by dropping the odd rows of (4.3.5) formally done by applying the restriction operator to both sides of (4.3.5) hence the residual equations becomes A h v 2h = r h this gives a plausible definition for the coarse-grid operator: A 2h = A h the terms of A 2h may be computed explicitly when applying A h term by term to the j-th unit vector on Ω 2h this establishes that the j-th column of A 2h and, by symmetry, also the j-th row of A 2h are given by 1 (2h) 2 ( 1 2 1) PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 33
34 variational properties (cont d) j 1 j j ½ 1 ½ 0 1 2h 2 1 4h h 2 2h 2 1 2h 2 1 4h 2 calculation of j-th row of A 2h = A h hence, we would get the same result if the original problem were simply discretised on Ω 2h using second-order finite differences therefore, by this definition, A 2h really is the Ω 2h version of A h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 34
35 two-grid correction scheme now using well-defined ways to transfer vectors between grids with previous correction scheme we define a two-grid method; parameters ν 1, ν 2 controlling number of relaxation steps in practice often 1, 2, or 3 relax ν 1 times on A h v h = f h on Ω h with initial guess v h compute residual r h = f h A h v h restrict residual r h to coarse grid by r 2h = r h solve A 2h e 2h = r 2h on Ω 2h prolongate coarse-grid error e 2h to fine grid by e h = e 2h correct fine-grid approximation v h v h + e h relax ν 2 times on A h v h = f h on Ω h with corrected approximation v h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 35
36 two-grid correction scheme (cont d) example: we consider the weighted JACOBI method with ω = ⅔ applied to our one-dimensional model problem Au = 0 on a grid with n = 64 points as initial guess (m 16 + m 40 )/2 is used, consisting of modes k = 16 and k = 40 the following two-grid correction scheme is applied relax three times on A h v h = 0 on Ω h with initial guess (m 16 + m 40 )/2 compute r 2h = r h relax three times on A 2h e 2h = r 2h on Ω 2h with initial guess e 2h = 0 correct fine-grid approximation v h v h + e 2h relax three times on A h v h = 0 on Ω h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 36
37 two-grid correction scheme (cont d) initial guess (m 16 + m 40 )/2 after one relaxation step after three relaxation steps after coarse-grid correction after one full 2-grid cycle after two full 2-grid cycles PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 37
38 V-cycle scheme question: best way to solve the coarse-grid problem A 2h e 2h = r 2h notice: coarse-grid problem not much different from original problem therefore, we can apply two-grid scheme to residual equation on Ω 2h, which means relaxing there (on Ω 2h ) and then moving forward to Ω 4h for the correction step recursive applications leads to successively coarser grids until a direct solution of the residual equation is possible simplified notation: for m > 1 grids with grid spacing l {2h, 4h,..., 2 m 1 h} RHS vector of residual equation on Ω l is also called f l instead of r l approximations on Ω l are also called v l instead of e l initial guesses on first visit to Ω l are chosen as v l = 0 PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 38
39 V-cycle scheme (cont d) idea of recursive algorithm with ν 1 pre- and ν 2 post-smoothing steps Ω h Ω 2h Ω 4h relax ν 1 times A h u h = f h f 2h = r h relax ν 1 times A 2h v 2h = f 2h f 4h = r 2h relax ν 1 times A 4h v 4h = f 4h relax ν 2 times A 4h e 4h = f 4h relax ν 2 times A h u h = f h relax ν 2 times A 2h e 2h = f 2h e 2h e 2h + u h u h + e 4h e 2h f 8h = r 4h e 4h e 4h + e 8h Ω 8h solve e 8h = (A 8h ) 1 f 8h PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 39
40 V-cycle scheme (cont d) Ω h Ω 2h v l MG V (v l, f l ) 1. relax ν 1 times on A l v l = f l with initial guess v l 2. if Ω l = coarsest grid, then go to step 4 else f 2l (f l A l v l ) v 2l 0 v 2l MG V (v 2l, f 2l ) Ω 4h Ω 8h 3. correct v l v l + v 2l 4. relax ν 2 times on A l v l = f l PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 40
41 μ-cycle scheme Ω h Ω 2h Ω 4h v l MG μ (v l, f l ) Ω 8h W-cycle (μ = 2) 1. relax ν 1 times on A l v l = f l with initial guess v l 2. if Ω l = coarsest grid, then go to step 4 else f 2l (f l A l v l ) v 2l 0 v 2l MG μ (v 2l, f 2l ) μ times 3. correct v l v l + v 2l In practice, only μ = 1 (which gives the V-cycle) and μ = 2 (which gives the W-cycle) are used. 4. relax ν 2 times on A l v l = f l PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 41
42 full multigrid V-cycle Ω h Ω 2h Ω 4h Ω 8h full multigrid with ν 0 = 1 v l FMG( f l ) 1. if Ω l = coarsest grid, set v l 0 and go to step 3 else f 2l (f l ) v 2l FMG( f 2l ) 2. correct v l v 2l 3. v l MG V (v l, f l ) ν 0 times Here, the idea is to use coarse grids in order to obtain better initial guesses, a strategy called nested iteration. PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 42
43 complexity storage: consider d-dimensional grid with n d grid points on Ω h on each coarser level, number of grid points decreases by a factor of (½) d i.e. Ω 2h = (½) d Ω h, Ω 4h = (½) d Ω 2h = (¼) d Ω h,... hence, to store all grids Ω h, Ω 2h, Ω 4h,... there is memory for grid points necessary using the sum of the geometric series as an upper bound gives a good estimation of memory requirements, namely storage < n d 1 2 d PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 43
44 computational cost it is convenient to measure the computational cost in terms of work units (WU), which is the cost performing one relaxation sweep on the finest grid Ω h neglecting the cost of restriction and prolongation operations, which typically amounts to 10 20% of the cost of the entire cycle consider a V-cycle with one relaxation sweep on each level (i.e. ν 1 = ν 2 = 1) each level is visited twice and on each coarser level work units decrease by a factor of (½) d adding these costs and again using the geometric series for an upper bound gives the V-cycle computation cost as PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 44
45 example let Au = 0 denote a linear system derived from a finite-difference discretisation of a one-dimensional second-order boundary value problem u (x) = 0, 0 < x < 1 u(0) = 0 u(1) = 1 with grid spacing h = 2 14, thus we have n = 1/h = 2 14 = 16,384 grid points matrix A is SPD and diagonally dominant, hence the JACOBI method will converge for arbitrary start vector v 0 as initial guess we choose v 0 = v h = 0 for a multigrid solver (V-cycle) PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 45
46 example (cont d) convergence plot of (standard) JACOBI method PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 46
47 example (cont d) convergence plot of multigrid method (V-cycle, 11 levels, ν 1 = ν 2 = 3) PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 47
48 example (cont d) comparison JACOBI, multigrid, CG PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 48
49 overview definitions splitting methods projection and KRYLOV subspace methods multigrid methods PD Dr. Ralf-Peter Mundani Computational Linear Algebra Winter Term 2018/19 49
Kasetsart 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 informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 1: Introduction to Multigrid 1 12/02/09 MG02.prz Error Smoothing 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 Initial Solution=-Error 0 10 20 30 40 50 60 70 80 90 100 DCT:
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 informationIntroduction to Multigrid Method
Introduction to Multigrid Metod Presented by: Bogojeska Jasmina /08/005 JASS, 005, St. Petersburg 1 Te ultimate upsot of MLAT Te amount of computational work sould be proportional to te amount of real
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 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 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 information6. Multigrid & Krylov Methods. June 1, 2010
June 1, 2010 Scientific Computing II, Tobias Weinzierl page 1 of 27 Outline of This Session A recapitulation of iterative schemes Lots of advertisement Multigrid Ingredients Multigrid Analysis Scientific
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 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 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 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 2: Direct Methods PD Dr.
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 informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationGeometric Multigrid Methods
Geometric Multigrid Methods Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University IMA Tutorial: Fast Solution Techniques November 28, 2010 Ideas
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationAlgebraic 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 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 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 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 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 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 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 informationMULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW
MULTIGRID METHODS FOR NONLINEAR PROBLEMS: AN OVERVIEW VAN EMDEN HENSON CENTER FOR APPLIED SCIENTIFIC COMPUTING LAWRENCE LIVERMORE NATIONAL LABORATORY Abstract Since their early application to elliptic
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 informationAlgebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov
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 information6. Iterative Methods for Linear Systems. The stepwise approach to the solution...
6 Iterative Methods for Linear Systems The stepwise approach to the solution Miriam Mehl: 6 Iterative Methods for Linear Systems The stepwise approach to the solution, January 18, 2013 1 61 Large Sparse
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 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 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 informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
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 informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 1: Introduction to Multigrid 2000 Eric de Sturler 1 12/02/09 MG01.prz Basic Iterative Methods (1) Nonlinear equation: f(x) = 0 Rewrite as x = F(x), and iterate x i+1
More informationGeometric Multigrid Methods for the Helmholtz equations
Geometric Multigrid Methods for the Helmholtz equations Ira Livshits Ball State University RICAM, Linz, 4, 6 November 20 Ira Livshits (BSU) - November 20, Linz / 83 Multigrid Methods Aim: To understand
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 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 informationPreconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators
Preconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University
More informationMultigrid Methods for Discretized PDE Problems
Towards Metods for Discretized PDE Problems Institute for Applied Matematics University of Heidelberg Feb 1-5, 2010 Towards Outline A model problem Solution of very large linear systems Iterative Metods
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 informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationHW4, Math 228A. Multigrid Solver. Fall 2010
HW4, Math 228A. Multigrid Solver Date due 11/30/2010 UC Davis, California Fall 2010 Nasser M. Abbasi Fall 2010 Compiled on January 20, 2019 at 4:13am [public] Contents 1 Problem 1 3 1.1 Restriction and
More informationThe Fundamentals and Advantages of Multi-grid Techniques
Joseph Kovac 18.086 Final Project Spring 2005 Prof. Gilbert Strang The Fundamentals and Advantages of Multi-grid Techniques Introduction The finite difference method represents a highly straightforward
More informationA Generalized Eigensolver Based on Smoothed Aggregation (GES-SA) for Initializing Smoothed Aggregation Multigrid (SA)
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2007; 07: 6 [Version: 2002/09/8 v.02] A Generalized Eigensolver Based on Smoothed Aggregation (GES-SA) for Initializing Smoothed Aggregation
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 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 informationEFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC ANALYSIS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK EFFICIENT MULTIGRID BASED SOLVERS FOR ISOGEOMETRIC
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 informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationDomain decomposition on different levels of the Jacobi-Davidson method
hapter 5 Domain decomposition on different levels of the Jacobi-Davidson method Abstract Most computational work of Jacobi-Davidson [46], an iterative method suitable for computing solutions of large dimensional
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationSparse Linear Systems. Iterative Methods for Sparse Linear Systems. Motivation for Studying Sparse Linear Systems. Partial Differential Equations
Sparse Linear Systems Iterative Methods for Sparse Linear Systems Matrix Computations and Applications, Lecture C11 Fredrik Bengzon, Robert Söderlund We consider the problem of solving the linear system
More informationChapter 5. Methods for Solving Elliptic Equations
Chapter 5. Methods for Solving Elliptic Equations References: Tannehill et al Section 4.3. Fulton et al (1986 MWR). Recommended reading: Chapter 7, Numerical Methods for Engineering Application. J. H.
More informationNew Multigrid Solver Advances in TOPS
New Multigrid Solver Advances in TOPS R D Falgout 1, J Brannick 2, M Brezina 2, T Manteuffel 2 and S McCormick 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, P.O.
More informationITERATIVE METHODS BASED ON KRYLOV SUBSPACES
ITERATIVE METHODS BASED ON KRYLOV SUBSPACES LONG CHEN We shall present iterative methods for solving linear algebraic equation Au = b based on Krylov subspaces We derive conjugate gradient (CG) method
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 6: Some Other Stuff PD Dr.
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
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 informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
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 informationA SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS
A SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS MICHAEL HOLST AND FAISAL SAIED Abstract. We consider multigrid and domain decomposition methods for the numerical
More informationComparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes
Comparison of V-cycle Multigrid Method for Cell-centered Finite Difference on Triangular Meshes Do Y. Kwak, 1 JunS.Lee 1 Department of Mathematics, KAIST, Taejon 305-701, Korea Department of Mathematics,
More information9. Iterative Methods for Large Linear Systems
EE507 - Computational Techniques for EE Jitkomut Songsiri 9. Iterative Methods for Large Linear Systems introduction splitting method Jacobi method Gauss-Seidel method successive overrelaxation (SOR) 9-1
More informationK.S. Kang. The multigrid method for an elliptic problem on a rectangular domain with an internal conductiong structure and an inner empty space
K.S. Kang The multigrid method for an elliptic problem on a rectangular domain with an internal conductiong structure and an inner empty space IPP 5/128 September, 2011 The multigrid method for an elliptic
More informationA Line search Multigrid Method for Large-Scale Nonlinear Optimization
A Line search Multigrid Method for Large-Scale Nonlinear Optimization Zaiwen Wen Donald Goldfarb Department of Industrial Engineering and Operations Research Columbia University 2008 Siam Conference on
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
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 informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
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 informationCLASSICAL ITERATIVE METHODS
CLASSICAL ITERATIVE METHODS LONG CHEN In this notes we discuss classic iterative methods on solving the linear operator equation (1) Au = f, posed on a finite dimensional Hilbert space V = R N equipped
More informationEECS 275 Matrix Computation
EECS 275 Matrix Computation Ming-Hsuan Yang Electrical Engineering and Computer Science University of California at Merced Merced, CA 95344 http://faculty.ucmerced.edu/mhyang Lecture 20 1 / 20 Overview
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 informationMultigrid finite element methods on semi-structured triangular grids
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real, -5 septiembre 009 (pp. 8) Multigrid finite element methods on semi-structured triangular grids F.J.
More informationIntroduction to Iterative Solvers of Linear Systems
Introduction to Iterative Solvers of Linear Systems SFB Training Event January 2012 Prof. Dr. Andreas Frommer Typeset by Lukas Krämer, Simon-Wolfgang Mages and Rudolf Rödl 1 Classes of Matrices and their
More informationAn improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media
An improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media Henri Calandra, Serge Gratton, Xavier Pinel and Xavier Vasseur Technical Report TR/PA/12/2
More informationSparse Matrix Techniques for MCAO
Sparse Matrix Techniques for MCAO Luc Gilles lgilles@mtu.edu Michigan Technological University, ECE Department Brent Ellerbroek bellerbroek@gemini.edu Gemini Observatory Curt Vogel vogel@math.montana.edu
More informationME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Solution to Final Examination Date: May 4, 2011 1:00 3:00 PM Instructor: J. Murthy Open Book, Open Notes Total: 100 points Use the finite volume
More informationMULTI-LEVEL TECHNIQUES FOR THE SOLUTION OF THE KINETIC EQUATIONS IN CONDENSING FLOWS SIMON GLAZENBORG
MULTI-LEVEL TECHNIQUES FOR THE SOLUTION OF THE KINETIC EQUATIONS IN CONDENSING FLOWS SIMON GLAZENBORG CONTENTS Introduction Theory Test case: Nucleation pulse Conclusions & recommendations 2 WHAT IS CONDENSATION
More informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
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 informationMultigrid Solution of the Debye-Hückel Equation
Rochester Institute of Technology RIT Scholar Works Theses Thesis/Dissertation Collections 5-17-2013 Multigrid Solution of the Debye-Hückel Equation Benjamin Quanah Parker Follow this and additional works
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 informationIntroduction to Multigrid Methods Sampling Theory and Elements of Multigrid
Introduction to Multigrid Methods Sampling Theory and Elements of Multigrid Gustaf Söderlind Numerical Analysis, Lund University Contents V3.15 What is a multigrid method? Sampling theory Digital filters
More information1 Extrapolation: A Hint of Things to Come
Notes for 2017-03-24 1 Extrapolation: A Hint of Things to Come Stationary iterations are simple. Methods like Jacobi or Gauss-Seidel are easy to program, and it s (relatively) easy to analyze their convergence.
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
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 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 informationScientific Computing II
Technische Universität München ST 008 Institut für Informatik Dr. Miriam Mehl Scientific Computing II Final Exam, July, 008 Iterative Solvers (3 pts + 4 extra pts, 60 min) a) Steepest Descent and Conjugate
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
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 informationOn nonlinear adaptivity with heterogeneity
On nonlinear adaptivity with heterogeneity Jed Brown jed@jedbrown.org (CU Boulder) Collaborators: Mark Adams (LBL), Matt Knepley (UChicago), Dave May (ETH), Laetitia Le Pourhiet (UPMC), Ravi Samtaney (KAUST)
More informationOn the Development of Implicit Solvers for Time-Dependent Systems
School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,
More informationarxiv: v1 [math.oc] 17 Dec 2018
ACCELERATING MULTIGRID OPTIMIZATION VIA SESOP TAO HONG, IRAD YAVNEH AND MICHAEL ZIBULEVSKY arxiv:82.06896v [math.oc] 7 Dec 208 Abstract. A merger of two optimization frameworks is introduced: SEquential
More informationMultigrid Techniques for Nonlinear Eigenvalue Problems; Solutions of a. Sorin Costiner and Shlomo Ta'asan
Multigrid Techniques for Nonlinear Eigenvalue Problems; Solutions of a Nonlinear Schrodinger Eigenvalue Problem in 2D and 3D Sorin Costiner and Shlomo Ta'asan Department of Applied Mathematics and Computer
More informationAn Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering
An Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering Bruno Turcksin & Jean C. Ragusa & Jim E. Morel Texas A&M University, Dept. of Nuclear Engineering Bruno
More informationResearch Article Evaluation of the Capability of the Multigrid Method in Speeding Up the Convergence of Iterative Methods
International Scholarly Research Network ISRN Computational Mathematics Volume 212, Article ID 172687, 5 pages doi:1.542/212/172687 Research Article Evaluation of the Capability of the Multigrid Method
More information