Multigrid Method ZHONG-CI SHI. Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China. Joint work: Xuejun Xu
|
|
- Emory Brown
- 6 years ago
- Views:
Transcription
1 Multigrid Method ZHONG-CI SHI Institute of Computational Mathematics Chinese Academy of Sciences, Beijing, China Joint work: Xuejun Xu
2 Outline Introduction Standard Cascadic Multigrid method(cmg) Economical Cascadic Multigrid method(ecmg) Summary Appendix
3 1. Introduction Multigrid techniques give algorithms that solve sparse linear systems Ax = b of N unknowns with O(N) computational complexity for large classes of problems. The systems arise from approximation of partial differential equations. The condition number is usually very large when the mesh size becomes small.
4 Background: Idea (Finite difference methods), 60.s: [Fedorenko, 1964], [Bachvalov, 1966] Applications, 70.s: [Brandt, 1973], etc. Theoretical Analysis, 80.s: [Bank and Dupont,1981], [Braess and Hackbusch,1983], [Bramble, Pasciak and Xu, 1987], etc. Features: For many iterative methods one can prove that the iterates {x i } satisfy: x i+1 x δ x i x for δ (0, 1) in certain norm. Jacobi, Gauss-Seidel: δ = 1 O(h 2 ). MG: δ < 1 is independent of the mesh size h.
5 Observation Two solution techniques Iterative method: Jacobi, Gauss-Seidel Easy to perform (matrix vector), but slow convergence with operations O(N 2 ) Reason: the higher frequencies of the residual can be supressed quickly, but lower frequencies are reduced very slowly. After few iterations, the residual becomes smooth, but not small. Direct method: Gauss elimination Good for small-size problem Rather difficult to perform for large-size problem, the accuracy is not high (round error accumulation)
6 MGM: Combination of the above two techniques Using different techniques on different levels of grid 1. Fine grid (smoothing) Iterative method to smooth the residual 2. Coarse grid (correction) Direct method to solve the residual equation 3. Add the correction to the initial residual a more precise residual 4. Back to 1. Total operations: O(N)
7 Two types of MGM: 1. W cycle: two corrections per cycle (1980 ) Difficult to carry out, easy to analyze 2. V cycle: one correction per cycle (1990 ) Easy to carry out, difficult to analyze New scheme Cascadic scheme No correction at all, only iterations on each level. Very easy to carry out, difficult to analyze
8 Model problem: The second order elliptic problem 2 i,j=1 The differential operator is uniformly elliptic: x i (a ij u x j ) = f in Ω, u = 0 on Ω, cξ t ξ d a ij ξ i ξ j Cξ t ξ (x, y) Ω, ξ R d. i,j=1 The variational form: find u H 1 0(Ω) such that a(u, v) = (f, v) v H 1 0(Ω), where the bilinear form is a(u, v) = d Ω i,j=1 a ij u x j v x i dx u, v H 1 (Ω). First Prev Next Last Go Back Full Screen Close Quit
9 Let Γ l (l 0) denote triangulation of Ω, which is generated uniformly from the initial triangulation Γ 0. The mesh size is h l = h 0 2 l. For instance, in 2D case, Γ l is obtained by linking the midpoints of three edges of triangle or the midpoints of counteredges of rectangle on Γ l 1. Let V l L 2 (Ω) denote the finite element space on Γ l. On each level, we define a l (u, v) = K Γ l d K i,j=1 a ij u x j v x i dx and v i,l = K Γ l v 2 i,k i = 0, 1, here v 0 = v 0,l.
10 It is easy to see that the following properties are satisfied for V l : (1) a l (u, v) C u 1,l v 1,l u, v V l. (2) a l (v, v) C v 2 1,l v V l. We assume that v 1,l is a norm over V l. The finite element approximation on V l : Find u l V l such that a l (u l, v) = (f, v) v V l. It is known that the above equation has a unique solution u l V l. Define the operator A l : V l V l as follows (A l u, v) = a l (u, v) u, v V l. Then we get a linear system A l u l = f l, here f l V l and (f l, v) = (f, v) v V l.
11 2. Cascadic Multigrid Method Bornemann and Deuflhard (1996) Shaidurov (1996) Shi and Xu (1998) Braess and Dahmen (1999) Stevensen (2002) Advantages: Simplicity No coarse grid corrections One way multigrid
12 Cascadic Multigrid Algorithm (1) Set u 0 0 = u 0 ˆ=u 0 and let u 0 l = I l u l 1. (2) For l = 1,..., L, do iterations: (3) Set u l ˆ=um l l. u m l l = C m l l u 0 l. Choose the iterative operator C l : V l V l on the level l and assume that there exists a linear operator T l : V l V l such that u l C m l l u 0 l = T m l l (u l u 0 l ). C l : Gauss-Seidel, Jacobi, CG.
13 Full Multigrid Algorithm (1) Set u 0 0 = u 0 ˆ=u 0 and let u 0 l = I l u l 1. (2) For l = 1,..., L, do iterations: u m l l = MG(l, u 0 l, r). (3) Set u l ˆ=um l l, where MG(l, u 0 l, r) denotes to do r-th multigrid iterations. Difference: No coarse grid corrections in CMG.
14 Optimal Cascadic Multigrid Method: If we obtain both the Accuracy u L u L L u u L L which means that the iterative error is comparable to the approximation error of the finite element method and the multigrid Complexity amount of work = O(n L ), n L = dimv L.
15 General Framework: (Shi and Xu, 1998) Three assumptions: Define the intergrid transfer operator I l : V l 1 V l which is assumed to satisfy the condition: ( H1) (1). v I l v t 1,l Ch l v t,l 1 v V l 1, (2). u l I l u l 1 t 1,l Ch 2 l f 0, where u l is the finite element solution on V l.
16 Choose the iterative operator C l : V l V l on the level l and assume that there exists a linear operator T l : V l V l such that u l C m l l u 0 l = T m l l (u l u 0 l ) ( H2) (1). T m l l v l C h 1 l m γ v t 1,l v V l, l (2). T m l l v l v l v V l, where m l is the number of iteration steps on the level l, and γ is a positive number depending on the given iteration.
17 Introduce a projection operator P l : V l 1 + V l V l defined by a l (P l u, v) = a l (u, v) v V l. From the definition, it is easily seen that P l v l v l 1 v V l 1. For the operator P l, we assume that (H3) v P l v t 1,l Ch l v t,l 1 v V l 1. Note that the project operator P l is used only in the convergence analysis, it is not needed in real computation.
18 Note that the mesh size on the level l is h l = h L 2 L l. Let m l (0 l L) be the smallest integer satisfying m l β L l m L (1) for some fixed β > 1, where m L is the number of iterations on the finest level L. Under the assumptions (H1), (H2) and (H3), if m l, the number of iterations on level l is given by (1), then the accuracy of the cascadic multigrid is u L u L L 1 C 1 ( 2 h L β ) m γ γ L f 0 for β > 2 1 γ, CL h L m γ f 0 for β = 2 1 γ. L
19 Complexity of CMG: The computational cost of the cascadic multigrid is proportional to L C 1 m l n l 1 β m L n L for β < 2 d, 2 l=0 d CLm L n L for β = 2 d, where d is the dimension of the domain Ω.
20 Main Results: Suppose (H1), (H2) and (H3) hold. (1) If γ = 1 2, d = 3, then the cascadic multigrid is optimal. (2) If γ = 1, d = 2, 3, then the cascadic multigrid is optimal. (3) If γ = 1 2, d = 2, and the number of iterations on the level L is then the error in the energy norm is m L = [m L 2 ], and the complexity is u L u L L C h L m 1 2 f 0, L m l n l cm n L (1 + logn L ) 3. l=0 It means that the cascadic multigrid is nearly optimal in this case.
21 Applications: 1. Lagrange conforming elements In this case, V l 1 V l. Therefore, both the transfer operator I l and the projection operator P l are the identity I, so that (H1)-1 and (H3) are trivial. By finite element error estimates, so (H1)-2 also holds. u l I l u l 1 0 = u l u l 1 0 u l u 0 + u u l 1 0 Ch 2 l f 0, (H2) has been proved by Bank and Dupont (1981)
22 2. Nonconforming elements Let V l be the P 1 nonconforming finite element space. Define the intergrid transfer operator I l : V l 1 V l as follows: let m be the midpoint of a side of triangles in Γ l. (1) If m lies in the interior of a triangle in Γ l 1, then (I l v)(m) := v(m). (2) If m lies on the common edge of two adjacent triangles T 1 and T 2 in Γ l 1, then (I l v)(m) := 1 2 (v T 1 (m) + v T2 (m)). (H1) is valid for the P1 nonconforming element. Using a duality argument, we can show that (H3) is valid for the P1 nonconforming element. (H2) is obvious. Other nonconforming elements: the Wilson element, the Carey element.
23 3. The plate bending problem Let Ω be a convex polygon in R 2, the variational form of the plate bending problem is to find u H 2 0(Ω) such that a(u, v) = Ω a(u, v) = (f, v) v H 2 0(Ω), { u v + (1 σ)(2 2 u 2 v 2 u 2 v x 1 x 2 x 1 x 2 x 2 1 x 2 2 (f, v) = fvdx, 2 u x 2 2 where f L 2 (Ω), and 0 < σ < 1 2 is the Possion ratio. In this case V = H2 0(Ω). Ω 2 v )}dx, x 2 1
24 Morley element The Morley finite element space is nonconforming and nonnested. We define the intergrid transfer operator as follows: (1) I l v = 0 at the vertex of Γ l along Ω and I lv n = 0 at the midpoint of each edge of Γ l along Ω. (2) Let p be a vertex of Γ l inside Ω. If p is also a vertex of Γ l 1, then (I l v)(p) = v(p). If p is the midpoint of a common edge of two triangles K 1 and K 2 in Γ l 1, then (I l v)(p) = 1 2 (v K 1 (p) + v K2 (p)).
25 (3) Let m be a midpoint of the edge e of Γ l inside Ω and n denote the unit normal of e. If m is inside of a triangle in Γ l 1, then I l v v (m) = n n (m). If m is on the common edge of two triangles K 1 and K 2 in Γ l 1, then We can prove (Shi and Xu 1998) I l v n (m) = 1 2 [ v K 1 n (m) + v K 2 n (m)]. (H1)-(H3) hold for the Morley element.
26 3. Economical Cascadic Multigrid Method Advantages: Less operations on the each level. Computational costs can be greatly reduced.
27 Example: CMG: m l = [m L β L l ] l = 1,, L, ECMG: a new criterion for choosing smoothing steps on each level. When d = 2, that is: (i) If l > L 0, then (ii) If l L 0, then m l = [m L β L l ]. m l = [m 1 2 (L (2 ε 0 )l)h 2 l ], here L 0 is a positive integer which will be determined later, and 0 < ε 0 1 is a fixed positive number, m L = m 0 (L L 0 ) 2. Note that in the standard cascadic algorithm m L = m 0 L 2.
28 Consider the second order elliptic problem and use Gauss-Seidel iteration as a smoother. For simplicity, we only consider 2D case. Then L 0 = L/2 = 4, ε 0 = 1 2, h 0 = 1 4, and m 0 = 2, m = 1. The following table shows the iteration steps on each level for the two methods. Level CMG ECMG We will prove that the new cascadic multigrid algorithm is still optimal in both accuracy and complexity as the standard CMG.
29 Basic smoothers: For the basic iterative smoothers, we have where T m l l v l ɛ(m l, h l ) v 0 v V l, C h 1 l, for m ɛ(m l, h l ) = m 1 l < κ l, 2 l C2 m l 1 κ l, for m l κ l, where κ l denotes the condition number of the matrix A l.
30 New criteria for choosing iteration parameters Determine the largest positive integer L 0, which satisfies the following inequality β L L 0 m L κ L0. For simplicity, we assume that κ l = h 2 l, then by a simple manipulation, L 0 Llogβ + logm L + 2logh 0. logβ + 2log2 Define the level parameter L 0 as the largest positive integer, which satisfies the following inequality: L 0 min{ Llogβ + logm L + 2logh 0 dl, }. (2) logβ + 2log2 2 + d
31 New criteria: 1. When d = 2, (i) If l > L 0, then (ii) If l L 0, then m l = [m L β L l ]. m l = [m 1 2 (L (2 ε 0 )l)κ l ]. In practical implementation, because κ l h 2 l, the above terms can be replaced by: m l = [m 1 2 (L (2 ε 0 )l)h 2 l ].
32 2. When d = 3, (i) If l > L 0, then m l = [m L β L l ]. (ii) If l L 0, there are two cases: (a) If (2 ε 0 )L 0 L, then m l = [m 1 2 (L (2 ε 0 )l)κ l ]. (b) If(2 ε 0 )L 0 > L, then there exists a largest positive integer L < L 0 such that (2 ε 0 )L L. For all l L, we choose m l as follows: m l = [m 1 2 (L (2 ε 0 )l)κ l ], for all L < l L 0, we choose m l as follows: m l = [m 1 2 κ l ].
33 Economical Cascadic Multigrid (1) Set u 0 0 = u 0 ˆ=u 0 and let u 0 l = I l u l 1. (2) For l = 1,..., L, do u m l l (3) Set u l ˆ=um l l, where the m l is determined by the new criteria. = C m l l u 0 l.
34 Main Results: Suppose (H1) and (H3) hold. (1) If d = 3, then the cascadic multigrid is optimal, i.e.,the error in the energy norm is u L u L L Ch L f 0 ( h 0 m 1 2 L + h 0 m 1 2 ) and the complexity is L l=0 m l n l C(h 2 0 m β 2 d m L )n L.
35 (2) If d = 2, and the number of iterations on the level L is m L = [m 0 (L L 0 ) 2 ], then the error in the energy norm is and the complexity of computation is L l=0 u L u L L Ch L f 0 ( h 0 m h 0 m ) m l n l C(m h m 0 (L L 0 ) 3 )n L. It means that the cascadic multigrid is nearly optimal in this case.
36 CG method as a smoother: Main results: Suppose (H1) and (H2) hold. For d = 2, 3, the cascadic multigrid with CG smoother is optimal, i.e., the error in the energy norm is and the complexity of computation is u L u L L Ch L f 0 ( h 0 m L + h 0 m ) L l=0 m l n l C(h β 2 d )n L. It means that the cascadic multigrid is optimal in this case.
37 Applications: Conforming and Nonconforming elements for the second order problem. Nonconforming plate elements, Morley element, Adini element.
38 Numerical experiments: We use ECMG and standard CMG algorithms to solve the following problem: { u = f, in Ω = ( 1, 1) ( 1, 1), u = 0, on Ω, where f is chosen such that the exact solution of the problem is u(x, y) = (1 x 2 )(1 y 2 ).
39 ECMG and CMG with GS smoother (P1 conforming element) CMG ECMG Mesh Energy error CPU Energy error CPU e-03 38(s) e-03 20(s) e (s) e (s) e (s) e (s)
40 ECMG and CMG with CG smoother (P1 conforming element) CMG ECMG Mesh Energy error CPU Energy error CPU e-03 14(s) e-03 10(s) e-03 66(s) e-03 42(s) e (s) e (s)
41 ECMG and CMG with GS smoother (P1 nonconforming element) CMG ECMG Mesh Energy error CPU Energy error CPU e-02 40(s) e-02 19(s) e (s) e-03 98(s) e (s) e (s)
42 ECMG and CMG with CG smoother (P1 nonconforming element) CMG ECMG Mesh Energy error CPU Energy error CPU e-02 18(s) e-03 14(s) e-03 72(s) e-03 48(s) e (s) e (s)
43 4. Summary Simplicity No coarse grid corrections, one way mutigrid ECMG: less operations on the each level. ECMG: computational costs can be greatly reduced.
44 5. Appendix Efficiency Computer vs Algorithm Computer from 10 3 (Kilo) (T era), increase 10 9 (Giga) Algorithm (solution of linear system) 1950 (Gauss elimination) 2000 (multigrid method) Convergence rate: from O(N 3 ) O(N), increase O(N 2 ) 3D problem: N = (10 2 ) 3 (Mega) Operations: , reduce (T era)!!
45 THANK YOU!
Geometric 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 informationA MULTILEVEL SUCCESSIVE ITERATION METHOD FOR NONLINEAR ELLIPTIC PROBLEMS
MATHEMATICS OF COMPUTATION Volume 73, Number 246, Pages 525 539 S 0025-5718(03)01566-7 Article electronically published on July 14, 2003 A MULTILEVEL SUCCESSIVE ITERATION METHOD FOR NONLINEAR ELLIPTIC
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 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 informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Xuejun Xu, Huangxin Chen, Ronald H.W. Hoppe Local Multilevel Methods for Adaptive Nonconforming Finite Element Methods Preprint Nr. 21/2009
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 informationLOCAL MULTILEVEL METHODS FOR ADAPTIVE NONCONFORMING FINITE ELEMENT METHODS
LOCAL MULTILEVEL METHODS FOR ADAPTIVE NONCONFORMING FINITE ELEMENT METHODS XUEUN XU, HUANGXIN CHEN, AND RONALD H.W. HOPPE Abstract. In this paper, we propose a local multilevel product algorithm and its
More informationTHE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE FOURTH ORDER PROBLEM *
Journal of Computational Mathematics Vol.30, No.5, 2012, 449 460. http://www.global-sci.org/jcm doi:10.4208/jcm.1203-m3855 THE BEST L 2 NORM ERROR ESTIMATE OF LOWER ORDER FINITE ELEMENT METHODS FOR THE
More informationJournal of Computational and Applied Mathematics. Multigrid method for solving convection-diffusion problems with dominant convection
Journal of Computational and Applied Mathematics 226 (2009) 77 83 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More informationAn Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions
An Iterative Substructuring Method for Mortar Nonconforming Discretization of a Fourth-Order Elliptic Problem in two dimensions Leszek Marcinkowski Department of Mathematics, Warsaw University, Banacha
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationMultigrid Methods for Maxwell s Equations
Multigrid Methods for Maxwell s Equations Jintao Cui Institute for Mathematics and Its Applications University of Minnesota Outline Nonconforming Finite Element Methods for a Two Dimensional Curl-Curl
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 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 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 for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
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 informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
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 informationAn additive average Schwarz method for the plate bending problem
J. Numer. Math., Vol. 10, No. 2, pp. 109 125 (2002) c VSP 2002 Prepared using jnm.sty [Version: 02.02.2002 v1.2] An additive average Schwarz method for the plate bending problem X. Feng and T. Rahman Abstract
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 informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationNumerische Mathematik
Numer. Math. (1999) 82: 179 191 Numerische Mathematik c Springer-Verlag 1999 Electronic Edition A cascadic multigrid algorithm for the Stokes equations Dietrich Braess 1, Wolfgang Dahmen 2 1 Fakultät für
More informationIt is known that Morley element is not C 0 element and it is divergent for Poisson equation (see [6]). When Morley element is applied to solve problem
Modied Morley Element Method for a ourth Order Elliptic Singular Perturbation Problem Λ Wang Ming LMAM, School of Mathematical Science, Peking University Jinchao u School of Mathematical Science, Peking
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
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 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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationFinite Element Multigrid Framework for Mimetic Finite Difference Discretizations
Finite Element Multigrid Framework for Mimetic Finite ifference iscretizations Xiaozhe Hu Tufts University Polytopal Element Methods in Mathematics and Engineering, October 26-28, 2015 Joint work with:
More informationVARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR.
VARIATIONAL AND NON-VARIATIONAL MULTIGRID ALGORITHMS FOR THE LAPLACE-BELTRAMI OPERATOR. ANDREA BONITO AND JOSEPH E. PASCIAK Abstract. We design and analyze variational and non-variational multigrid algorithms
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
More informationMultigrid solvers for equations arising in implicit MHD simulations
Multigrid solvers for equations arising in implicit MHD simulations smoothing Finest Grid Mark F. Adams Department of Applied Physics & Applied Mathematics Columbia University Ravi Samtaney PPPL Achi Brandt
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 informationRobust multigrid methods for nonsmooth coecient elliptic linear systems
Journal of Computational and Applied Mathematics 123 (2000) 323 352 www.elsevier.nl/locate/cam Robust multigrid methods for nonsmooth coecient elliptic linear systems Tony F. Chan a; ; 1, W.L. Wan b; 2
More informationOPTIMAL MULTILEVEL METHODS FOR GRADED BISECTION GRIDS
OPTIMAL MULTILEVEL METHODS FOR GRADED BISECTION GRIDS LONG CHEN, RICARDO H. NOCHETTO, AND JINCHAO XU Abstract. We design and analyze optimal additive and multiplicative multilevel methods for solving H
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 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 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 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 informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
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 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 informationEfficient smoothers for all-at-once multigrid methods for Poisson and Stokes control problems
Efficient smoothers for all-at-once multigrid methods for Poisson and Stoes control problems Stefan Taacs stefan.taacs@numa.uni-linz.ac.at, WWW home page: http://www.numa.uni-linz.ac.at/~stefant/j3362/
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 PRECONDITIONING FOR THE BIHARMONIC DIRICHLET PROBLEM M. R. HANISCH
MULTIGRID PRECONDITIONING FOR THE BIHARMONIC DIRICHLET PROBLEM M. R. HANISCH Abstract. A multigrid preconditioning scheme for solving the Ciarlet-Raviart mixed method equations for the biharmonic Dirichlet
More informationETNA Kent State University
g ' Electronic Transactions on umerical Analysis. Volume 17 pp. 112-132 2004. Copyright 2004. SS 10-913. COVERGECE OF V-CYCE AD F-CYCE MTGRD METHODS FOR THE BHARMOC PROBEM SG THE MOREY EEMET JE ZHAO Abstract.
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 informationInstitut für Mathematik
U n i v e r s i t ä t A u g s b u r g Institut für Mathematik Xuejun Xu, Huangxin Chen, Ronald H.W. Hoppe Optimality of Local Multilevel Methods on Adaptively Refined Meshes for Elliptic Boundary Value
More informationFinite Element Methods for Fourth Order Variational Inequalities
Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2013 Finite Element Methods for Fourth Order Variational Inequalities Yi Zhang Louisiana State University and Agricultural
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationHigh order, finite volume method, flux conservation, finite element method
FLUX-CONSERVING FINITE ELEMENT METHODS SHANGYOU ZHANG, ZHIMIN ZHANG, AND QINGSONG ZOU Abstract. We analyze the flux conservation property of the finite element method. It is shown that the finite element
More informationSolving Negative Order Equations by the Multigrid Method Via Variable Substitution
J Sci Comput (214) 59:371 385 DOI 1.17/s1915-13-9762-4 Solving Negative Order Equations by the Multigrid Method Via Variable Substitution George C. Hsiao Liwei Xu Shangyou Zhang Received: 26 December 212
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 informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
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 informationTREATMENTS OF DISCONTINUITY AND BUBBLE FUNCTIONS IN THE MULTIGRID METHOD
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 1055 1072 S 0025-5718(97)00853-3 TREATMENTS OF DISCONTINUITY AND BUBBLE FUNCTIONS IN THE MULTIGRID METHOD SHANGYOU ZHANG AND ZHIMIN ZHANG
More informationOversampling for the partition of unity parallel finite element algorithm
International conference of Computational Mathematics and its application Urumqi, China, July 25-27, 203 Oversampling for the partition of unity parallel finite element algorithm HAIBIAO ZHENG School of
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 informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationFEniCS Course. Lecture 0: Introduction to FEM. Contributors Anders Logg, Kent-Andre Mardal
FEniCS Course Lecture 0: Introduction to FEM Contributors Anders Logg, Kent-Andre Mardal 1 / 46 What is FEM? The finite element method is a framework and a recipe for discretization of mathematical problems
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More information3D Space Charge Routines: The Software Package MOEVE and FFT Compared
3D Space Charge Routines: The Software Package MOEVE and FFT Compared Gisela Pöplau DESY, Hamburg, December 4, 2007 Overview Algorithms for 3D space charge calculations Properties of FFT and iterative
More informationLecture Notes: African Institute of Mathematics Senegal, January Topic Title: A short introduction to numerical methods for elliptic PDEs
Lecture Notes: African Institute of Mathematics Senegal, January 26 opic itle: A short introduction to numerical methods for elliptic PDEs Authors and Lecturers: Gerard Awanou (University of Illinois-Chicago)
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationA DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT
A DELTA-REGULARIZATION FINITE ELEMENT METHOD FOR A DOUBLE CURL PROBLEM WITH DIVERGENCE-FREE CONSTRAINT HUOYUAN DUAN, SHA LI, ROGER C. E. TAN, AND WEIYING ZHENG Abstract. To deal with the divergence-free
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 informationarxiv: v2 [math.na] 29 Nov 2016
Noname manuscript No. (will be inserted by the editor) Multigrid algorithms for hp-version Interior Penalty Discontinuous Galerkin methods on polygonal and polyhedral meshes P. F. Antonietti P. Houston
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Numerical Methods for Partial Differential Equations Finite Difference Methods
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 informationChapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction
Chapter 6 A posteriori error estimates for finite element approximations 6.1 Introduction The a posteriori error estimation of finite element approximations of elliptic boundary value problems has reached
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 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 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 informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
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 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 informationOn Surface Meshes Induced by Level Set Functions
On Surface Meshes Induced by Level Set Functions Maxim A. Olshanskii, Arnold Reusken, and Xianmin Xu Bericht Nr. 347 Oktober 01 Key words: surface finite elements, level set function, surface triangulation,
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
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 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 informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationGoal-oriented error control of the iterative solution of finite element equations
J. Numer. Math., Vol. 0, No. 0, pp. 1 31 (2009) DOI 10.1515/ JNUM.2009.000 c de Gruyter 2009 Prepared using jnm.sty [Version: 20.02.2007 v1.3] Goal-oriented error control of the iterative solution of finite
More informationLocal Inexact Newton Multilevel FEM for Nonlinear Elliptic Problems
Konrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D-10711 Berlin-Wilmersdorf Peter Deuflhard Martin Weiser Local Inexact Newton Multilevel FEM for Nonlinear Elliptic Problems Preprint
More informationNumerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element
Numerical analysis and a-posteriori error control for a new nonconforming quadrilateral linear finite element M. Grajewski, J. Hron, S. Turek Matthias.Grajewski@math.uni-dortmund.de, jaroslav.hron@math.uni-dortmund.de,
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 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 informationA Mixed Nonconforming Finite Element for Linear Elasticity
A Mixed Nonconforming Finite Element for Linear Elasticity Zhiqiang Cai, 1 Xiu Ye 2 1 Department of Mathematics, Purdue University, West Lafayette, Indiana 47907-1395 2 Department of Mathematics and Statistics,
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More information