Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
|
|
- Reynard Hicks
- 6 years ago
- Views:
Transcription
1 Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt 10th Workshop of the INRIA-Illinois-ANL Joint Laboratory November 26th, 2013
2 A short introduction to DDM Consider the linear system: Au = f R n. Ω 2 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
3 A short introduction to DDM Consider the linear system: Au = f R n. Given a decomposition of 1; n, (N 1, N 2 ), define: the restriction operator R i from 1; n into N i, R T i as the extension by 0 from N i into 1; n. Ω 1 Ω 2 2 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
4 A short introduction to DDM Consider the linear system: Au = f R n. Given a decomposition of 1; n, (N 1, N 2 ), define: the restriction operator R i from 1; n into N i, Ri T as the extension by 0 from N i into 1; n. Then solve concurrently: u m+1 1 = u m 1 + A 1 11 R 1 (f Au m ) u m+1 2 = u m 2 + A 1 22 R 2 (f Au m ) where u i = R i u and A ij := R i AR T j. Ω 1 Ω 2 2 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
5 A short introduction II We have effectively divided, but we have yet to conquer. Duplicated unknowns coupled via a partition of unity: I = N i=1 R T i D i R i / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
6 A short introduction II We have effectively divided, but we have yet to conquer. Duplicated unknowns coupled via a partition of unity: N I = Ri T D i R i. i= Then, u m+1 = N i=1 R T i D i u m+1 i. 3 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
7 A short introduction II We have effectively divided, but we have yet to conquer. Duplicated unknowns coupled via a partition of unity: N I = Ri T D i R i. i= Then, u m+1 = N i=1 R T i D i u m+1 i. M 1 = N Ri T i=1 D i A 1 ii R i. 3 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
8 Contributions and goals Based on algebraic results with the p. of u., we propose: 1 a reformulation of the global matrix-vector product eliminating the need of a global ordering, 2 a construction of a so-called coarse operator to enhance a simple preconditioner. We are interested in the solution of various SPD systems, independently of: the discretization order, the contrast in the coefficients, the number of subdomains. 4 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
9 Using the overlap to its fullest extent DDM methods are seldom used as standalone solvers. Krylov methods and overlapping Schwarz methods Au = efficient global matrix-vector product 5 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
10 Using the overlap to its fullest extent DDM methods are seldom used as standalone solvers. Krylov methods and overlapping Schwarz methods Au = N j=1 AR T j D j R j u 5 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
11 Using the overlap to its fullest extent DDM methods are seldom used as standalone solvers. Krylov methods and overlapping Schwarz methods N R i Au = R i ARj T j=1 D j R j u 5 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
12 Using the overlap to its fullest extent DDM methods are seldom used as standalone solvers. Krylov methods and overlapping Schwarz methods N R i Au = R i ARj T j=1 D j R j u = j O i A ij D j R j u O i are the neighbors of Ω i, O i = O i {i}. 5 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
13 Using the overlap to its fullest extent DDM methods are seldom used as standalone solvers. Krylov methods and overlapping Schwarz methods N R i Au = R i ARj T D j R j u = A ij D j R j u j=1 j O i = R i Rj T A jj D j R j u. local unknowns on Ω j j O i no need for the global matrix, only local to neighbors mappings. explicit point-to-point communications via R i R T j. reuse of the operators from the preconditioner, A ii. O i are the neighbors of Ω i, O i = O i {i}. 5 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
14 Limitations of one-level methods One-level methods don t require exchange of global information. This hampers numerical scalability of such preconditioners. 6 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
15 Two-level preconditioners I A common technique in the field of DDM, MG, deflation: introduce an auxiliary coarse problem. Let Z be a rectangular matrix. Define E := Z T AZ. Z has O(N) columns, hence E is much smaller than A. 7 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
16 Two-level preconditioners I A common technique in the field of DDM, MG, deflation: introduce an auxiliary coarse problem. Let Z be a rectangular matrix. Define E := Z T AZ. Z has O(N) columns, hence E is much smaller than A. Enrich the original preconditioner, e.g. additively c.f. (Tang et al. 2009). P 1 = M 1 + ZE 1 Z T, 7 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
17 Two-level preconditioners II The construction of Z and the assembly of E are challenging. Let each domain compute concurrently ν i vectors { Λ ij } νi j=1. Define local dense rectangular matrices: W i = [ D i Λ i1 D i Λ i2 D i Λ iνi ]. Then, define the global deflation matrix as: Z = [ R T 1 W 1 R T 2 W 2 R T N W N ] 8 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
18 Generalized eigenvalue problems For theoretical justification of Z, see (Spillane et al. 2011). Solved by ARPACK concurently: where A N i Λ j = λ j D i R T i,0r i,0 A N i D i Λ j A N i is the local unassembled matrix, ( R i,0 is the restriction from Ω i to Ω i j O i Ω j ). 9 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
19 Workflow during one coarse operator correction How to compute ZE 1 Z T u R n? 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
20 Workflow during one coarse operator correction How to compute Z T u R n? Z T u = = m n n operations & MPI_Gather 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
21 Workflow during one coarse operator correction How to compute E 1 Z T u R n? Z T u = = m n n (Z T AZ) 1 Z T u = \ = operations & MPI_Gather + linear solve 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
22 Workflow during one coarse operator correction How to compute ZE 1 Z T u R n? Z T u = = m n = n (Z T AZ) 1 Z T u = \ = = Z(Z T AZ) 1 Z T u operations & MPI_Gather + linear solve + MPI_Scatter & operations 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
23 Workflow during one coarse operator correction How to compute ZE 1 Z T u R n? Z T u = = m n = n (Z T AZ) 1 Z T u = \ = = Z(Z T AZ) 1 Z T u operations & MPI_Gather + linear solve + MPI_Scatter & operations Communication pattern = global reduction at the coarse level. 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
24 Workflow during one coarse operator correction How to compute ZE 1 Z T u R n? Z T u = = m n = n = (Z T AZ) 1 Z T u = \ = = Z(Z T AZ) 1 Z T u operations & MPI_Gather + linear solve + MPI_Scatter & operations Communication pattern = global reduction at the coarse level. 10 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
25 Distribution of the coarse operator How can one solve E 1 z = c R m? Some constraints: 1 E cannot be centralized on a single MPI process, 2 E cannot be distributed on all MPI processes, 3 the solution must be computed fast and reliably. 11 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
26 Distribution of the coarse operator How can one solve E 1 z = c R m? Some constraints: 1 E cannot be centralized on a single MPI process, 2 E cannot be distributed on all MPI processes, 3 the solution must be computed fast and reliably. = use a direct solver with a distributed matrix on few master processes (number chosen at runtime). 11 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
27 Assembly for Schwarz methods Recalling E = ZAZ T, it can be proven that the block (i, j) E ij = Wi T A ij W j = Wi T R i Rj T A jj W j 1 compute locally T i = A ii W i (csrmm), 2 send to each neighbor, S j = R j R T i T i, 3 receive from each neighbor U j = R i R T j T j, 4 compute locally E i,i = Wi T T i (gemm), 5 compute locally E i,j = Wi T U j (gemm). Note: steps 2 and 3 overlap with step 4, if j O i, R i R T j = / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
28 Example of heterogeneous coefficients (κ u) = f + BC κ(x, y) 13 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
29 2D geometry (E 1, E 2 ) = (200, 0.01) GPa (ν 1, ν 2 ) = (0.25, 0.45) σ = f + BC 14 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
30 3D geometry 15 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
31 Machine used for scaling runs Curie Thin Nodes 5,040 compute nodes. 2 eight-core Intel Sandy Bridge@2.7 GHz per node. IB QDR full fat tree. 1.7 PFLOPs peak performance. 16 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
32 Strong scaling (linear elasticity) 1 subdomain/mpi process, 2 OpenMP threads/mpi process. Runtime (seconds) Linear speedup 3D (P 2 FE) 2D (P 3 FE) #processes
33 Strong scaling (linear elasticity) 1 subdomain/mpi process, 2 OpenMP threads/mpi process. 2.1B d.o.f. in 2D (P 3 FE) 300M d.o.f. in 3D (P 2 FE) 100% 100% Ratio 80% 60% 40% 20% 80% 60% 40% 20% 0% #processes #processes Factorization Deflation vectors Coarse operator Krylov method 0%
34 Weak scaling (scalar diffusion equation) 1 subdomain/mpi process, 2 OpenMP threads/mpi process. Efficiency relative to 256 MPI processes 100 % 80 % 60 % 40 % 20 % 0 % D (P 2 FE) 2D (P 4 FE) #d.o.f. (in millions) #processes
35 Weak scaling (scalar diffusion equation) Time (seconds) 1 subdomain/mpi process, 2 OpenMP threads/mpi process M d.o.f. sbdmn in 2D (P 4 FE) k d.o.f. sbdmn in 3D (P 2 FE) #processes #processes Factorization Deflation vectors Coarse operator Krylov method 0
36 Distributed global matrix Local to global mapping = distribution of the global matrix à la PETSc (split row-wise). Comparing performance of setup and solution phases between our solver against purely algebraic (+ near null space) solvers: GASM one-level domain decomposition method (ANL), Hypre BoomerAMG algebraic multigrid (LLNL), GAMG algebraic multrigrid (ANL/LBL). 20 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
37 Relative residual error Solution of a linear system I Homogeneous 3D Poisson equation discretized by P 1 FE solved on 2,048 MPI processes, 111M d.o.f Schwarz GenEO PETSc GASM Hypre BoomerAMG PETSc GAMG #iterations Time (seconds) " % " " GASM BoomerAMG GAMG Schwarz GenEO Setup Solution 21 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
38 Relative residual error Solution of a linear system II Heterogeneous 3D linear elasticity equation discretized by P 2 FE solved on 2,048 MPI processes, 127M d.o.f. " 100 Setup Solution Schwarz GenEO PETSc GASM Hypre BoomerAMG PETSc GAMG (fail) #iterations Time (seconds) 50 0 % % % GASM BoomerAMG GAMG Schwarz GenEO 22 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
39 Final words Limitations: scaling of the coarse operator in 3D beyond 10k subdomains, deflation vectors need elementary matrices to be computed. Summary: scalable framework for building two-level preconditioners for both Schwarz or substructuring methods (FETI-1), easily interfacable (FEM, FVM) without a global ordering. Outlooks: adaptive (re)construction/recycling of the coarse operator, nonlinear and saddle point problems. 23 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
40 Final words Limitations: scaling of the coarse operator in 3D beyond 10k subdomains, deflation vectors need elementary matrices to be computed. Summary: scalable framework for building two-level preconditioners for both Schwarz or substructuring methods (FETI-1), easily interfacable (FEM, FVM) without a global ordering. Outlooks: adaptive (re)construction/recycling of the coarse operator, nonlinear and saddle point problems. Thank you! 23 / 23 Pierre Jolivet Scalable Domain Decomposition Preconditioners
41 References Backup slides Spillane, N., V. Dolean, P. Hauret, F. Nataf, C. Pechstein, and R. Scheichl (2011). A robust two-level domain decomposition preconditioner for systems of PDEs. In: Comptes Rendus Mathematique , pp Tang, J., R. Nabben, C. Vuik, and Y. Erlangga (2009). Comparison of two-level preconditioners derived from deflation, domain decomposition and multigrid methods. In: Journal of Scientific Computing 39.3, pp / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
42 References Backup slides Solvers parameters Schwarz GenEO: ν i = 20, overlap = 1 (geometric). PETSc GASM: overlap = 10 (algebraic). Hyper BoomerAMG: HMIS coarsening, extended classical interpolation, no CF-relaxation, 2 levels of aggressive coarsening. PETSc GAMG: 1 smoothing step, -mg_levels_ksp_type richardson -mg_levels_pc_type sor. OpenMPI bindings for hybrid runs: --bind-to-socket --bycore. 2 / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
43 References Backup slides Distribution of the coarse operator Uniform distribution Non-uniform distribution Distribution of E when built with N = 16 using 4 masters. On the right, the number of values per master is roughly the same if the values below the diagonal are dropped (symmetric coarse operator). 3 / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
44 References Backup slides Timings for assembling the coarse operator N P dim(e) O i (average) Memory cost of E 1 Time 3D MB 2.78 s MB 3.35 s MB 93 MB 4.42 s s MB 138 MB 6.91 s 5.68 s MB 172 MB s 8.04 s MB 241 MB s s N P dim(e) O i (average) Memory cost of E 1 Time 2D MB 9.39 s MB 9.96 s MB 57 MB 9.92 s s MB 83 MB s 6.20 s MB 73 MB s 5.10 s MB 118 MB s 6.96 s 4 / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
45 References Backup slides Strong scaling (linear elasticity) 3D 2D N Factorization Deflation Solution #it. Total #d.o.f s s s s s s s s s s 9.73 s s s s 6.05 s s s s s s s s s s s s 8.64 s s s s 4.79 s s / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
46 References Backup slides Weak scaling (scalar diffusion equation) 3D 2D N Factorization Deflation Solution #it. Total #d.o.f s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s s / 6 (+ 23) Pierre Jolivet Scalable Domain Decomposition Preconditioners
Multilevel spectral coarse space methods in FreeFem++ on parallel architectures
Multilevel spectral coarse space methods in FreeFem++ on parallel architectures Pierre Jolivet Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann DD 21, Rennes. June 29, 2012 In collaboration with
More informationAvancées récentes dans le domaine des solveurs linéaires Pierre Jolivet Journées inaugurales de la machine de calcul June 10, 2015
spcl.inf.ethz.ch @spcl eth Avancées récentes dans le domaine des solveurs linéaires Pierre Jolivet Journées inaugurales de la machine de calcul June 10, 2015 Introduction Domain decomposition methods The
More informationHPDDM une bibliothèque haute performance unifiée pour les méthodes de décomposition de domaine
HPDDM une bibliothèque haute performance unifiée pour les méthodes de décomposition de domaine P. Jolivet and Frédéric Nataf Laboratory J.L. Lions, Univ. Paris VI, Equipe Alpines INRIA-LJLL et CNRS joint
More informationToward black-box adaptive domain decomposition methods
Toward black-box adaptive domain decomposition methods Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI joint work with Victorita Dolean (Univ. Nice Sophia-Antipolis)
More informationUne méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis
Une méthode parallèle hybride à deux niveaux interfacée dans un logiciel d éléments finis Frédéric Nataf Laboratory J.L. Lions (LJLL), CNRS, Alpines Inria and Univ. Paris VI Victorita Dolean (Univ. Nice
More informationRecent advances in HPC with FreeFem++
Recent advances in HPC with FreeFem++ Pierre Jolivet Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann Fourth workshop on FreeFem++ December 6, 2012 With F. Hecht, F. Nataf, C. Prud homme. Outline
More informationOverlapping Domain Decomposition Methods with FreeFem++
Overlapping Domain Decomposition Methods with FreeFem++ Pierre Jolivet, Frédéric Hecht, Frédéric Nataf, Christophe Prud Homme To cite this version: Pierre Jolivet, Frédéric Hecht, Frédéric Nataf, Christophe
More informationOverlapping domain decomposition methods withfreefem++
Overlapping domain decomposition methods withfreefem++ Pierre Jolivet 1,3, Frédéric Hecht 1, Frédéric Nataf 1, and Christophe Prud homme 2 1 Introduction Developping an efficient and versatile framework
More informationMultipréconditionnement adaptatif pour les méthodes de décomposition de domaine. Nicole Spillane (CNRS, CMAP, École Polytechnique)
Multipréconditionnement adaptatif pour les méthodes de décomposition de domaine Nicole Spillane (CNRS, CMAP, École Polytechnique) C. Bovet (ONERA), P. Gosselet (ENS Cachan), A. Parret Fréaud (SafranTech),
More informationSome examples in domain decomposition
Some examples in domain decomposition Frédéric Nataf nataf@ann.jussieu.fr, www.ann.jussieu.fr/ nataf Laboratoire J.L. Lions, CNRS UMR7598. Université Pierre et Marie Curie, France Joint work with V. Dolean
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
More informationEfficient domain decomposition methods for the time-harmonic Maxwell equations
Efficient domain decomposition methods for the time-harmonic Maxwell equations Marcella Bonazzoli 1, Victorita Dolean 2, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier 4 1 Inria Saclay (Defi
More informationContents. Preface... xi. Introduction...
Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism
More informationANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow
ANR Project DEDALES Algebraic and Geometric Domain Decomposition for Subsurface Flow Michel Kern Inria Paris Rocquencourt Maison de la Simulation C2S@Exa Days, Inria Paris Centre, Novembre 2016 M. Kern
More informationAdaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms
Adaptive Coarse Spaces and Multiple Search Directions: Tools for Robust Domain Decomposition Algorithms Nicole Spillane Center for Mathematical Modelling at the Universidad de Chile in Santiago. July 9th,
More informationScalable domain decomposition preconditioners for heterogeneous elliptic problems 1
Scientific Programming 22 (2014) 157 171 157 DOI 10.3233/SPR-140381 IOS Press Scalable domain decomposition preconditioners for heterogeneous elliptic problems 1 Pierre Jolivet a,b,c,, Frédéric Hecht b,c,
More informationOn the choice of abstract projection vectors for second level preconditioners
On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik
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 informationDomain Decomposition solvers (FETI)
Domain Decomposition solvers (FETI) a random walk in history and some current trends Daniel J. Rixen Technische Universität München Institute of Applied Mechanics www.amm.mw.tum.de rixen@tum.de 8-10 October
More informationA High-Performance Parallel Hybrid Method for Large Sparse Linear Systems
Outline A High-Performance Parallel Hybrid Method for Large Sparse Linear Systems Azzam Haidar CERFACS, Toulouse joint work with Luc Giraud (N7-IRIT, France) and Layne Watson (Virginia Polytechnic Institute,
More informationAlgebraic Coarse Spaces for Overlapping Schwarz Preconditioners
Algebraic Coarse Spaces for Overlapping Schwarz Preconditioners 17 th International Conference on Domain Decomposition Methods St. Wolfgang/Strobl, Austria July 3-7, 2006 Clark R. Dohrmann Sandia National
More informationJOHANNES KEPLER UNIVERSITY LINZ. Abstract Robust Coarse Spaces for Systems of PDEs via Generalized Eigenproblems in the Overlaps
JOHANNES KEPLER UNIVERSITY LINZ Institute of Computational Mathematics Abstract Robust Coarse Spaces for Systems of PDEs via Generalized Eigenproblems in the Overlaps Nicole Spillane Frédérik Nataf Laboratoire
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 informationFrom Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D
From Direct to Iterative Substructuring: some Parallel Experiences in 2 and 3D Luc Giraud N7-IRIT, Toulouse MUMPS Day October 24, 2006, ENS-INRIA, Lyon, France Outline 1 General Framework 2 The direct
More informationAn Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation
An Introduction to Domain Decomposition Methods: algorithms, theory and parallel implementation Victorita Dolean, Pierre Jolivet, Frédéric ataf To cite this version: Victorita Dolean, Pierre Jolivet, Frédéric
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 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 informationAuxiliary space multigrid method for elliptic problems with highly varying coefficients
Auxiliary space multigrid method for elliptic problems with highly varying coefficients Johannes Kraus 1 and Maria Lymbery 2 1 Introduction The robust preconditioning of linear systems of algebraic equations
More informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
More informationProgress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling
Progress in Parallel Implicit Methods For Tokamak Edge Plasma Modeling Michael McCourt 1,2,Lois Curfman McInnes 1 Hong Zhang 1,Ben Dudson 3,Sean Farley 1,4 Tom Rognlien 5, Maxim Umansky 5 Argonne National
More informationConvergence Behavior of a Two-Level Optimized Schwarz Preconditioner
Convergence Behavior of a Two-Level Optimized Schwarz Preconditioner Olivier Dubois 1 and Martin J. Gander 2 1 IMA, University of Minnesota, 207 Church St. SE, Minneapolis, MN 55455 dubois@ima.umn.edu
More informationComparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials
Comput Mech (2012) 50:321 333 DOI 10.1007/s00466-011-0661-y ORIGINAL PAPER Comparison of the deflated preconditioned conjugate gradient method and algebraic multigrid for composite materials T. B. Jönsthövel
More informationMaster Thesis Literature Study Presentation
Master Thesis Literature Study Presentation Delft University of Technology The Faculty of Electrical Engineering, Mathematics and Computer Science January 29, 2010 Plaxis Introduction Plaxis Finite Element
More informationArchitecture Solutions for DDM Numerical Environment
Architecture Solutions for DDM Numerical Environment Yana Gurieva 1 and Valery Ilin 1,2 1 Institute of Computational Mathematics and Mathematical Geophysics SB RAS, Novosibirsk, Russia 2 Novosibirsk State
More informationElectronic Transactions on Numerical Analysis Volume 49, 2018
Electronic Transactions on Numerical Analysis Volume 49, 2018 Contents 1 Adaptive FETI-DP and BDDC methods with a generalized transformation of basis for heterogeneous problems. Axel Klawonn, Martin Kühn,
More informationAdaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers
Adaptive Coarse Space Selection in BDDC and FETI-DP Iterative Substructuring Methods: Towards Fast and Robust Solvers Jan Mandel University of Colorado at Denver Bedřich Sousedík Czech Technical University
More informationC. Vuik 1 R. Nabben 2 J.M. Tang 1
Deflation acceleration of block ILU preconditioned Krylov methods C. Vuik 1 R. Nabben 2 J.M. Tang 1 1 Delft University of Technology J.M. Burgerscentrum 2 Technical University Berlin Institut für Mathematik
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
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 informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
More informationShort title: Total FETI. Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ Ostrava, Czech Republic
Short title: Total FETI Corresponding author: Zdenek Dostal, VŠB-Technical University of Ostrava, 17 listopadu 15, CZ-70833 Ostrava, Czech Republic mail: zdenek.dostal@vsb.cz fax +420 596 919 597 phone
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 informationTwo new enriched multiscale coarse spaces for the Additive Average Schwarz method
346 Two new enriched multiscale coarse spaces for the Additive Average Schwarz method Leszek Marcinkowski 1 and Talal Rahman 2 1 Introduction We propose additive Schwarz methods with spectrally enriched
More informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationOVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY TR
OVERLAPPING SCHWARZ ALGORITHMS FOR ALMOST INCOMPRESSIBLE LINEAR ELASTICITY MINGCHAO CAI, LUCA F. PAVARINO, AND OLOF B. WIDLUND TR2014-969 Abstract. Low order finite element discretizations of the linear
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 informationParallel Preconditioning Methods for Ill-conditioned Problems
Parallel Preconditioning Methods for Ill-conditioned Problems Kengo Nakajima Information Technology Center, The University of Tokyo 2014 Conference on Advanced Topics and Auto Tuning in High Performance
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationKey words. Parallel iterative solvers, saddle-point linear systems, preconditioners, timeharmonic
PARALLEL NUMERICAL SOLUTION OF THE TIME-HARMONIC MAXWELL EQUATIONS IN MIXED FORM DAN LI, CHEN GREIF, AND DOMINIK SCHÖTZAU Numer. Linear Algebra Appl., Vol. 19, pp. 525 539, 2012 Abstract. We develop a
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 informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
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 informationDeflated Krylov Iterations in Domain Decomposition Methods
305 Deflated Krylov Iterations in Domain Decomposition Methods Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 1 Introduction The goal of this research is an investigation of some advanced versions of
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 68 (2014) 1151 1160 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa A GPU
More informationExtending the theory for domain decomposition algorithms to less regular subdomains
Extending the theory for domain decomposition algorithms to less regular subdomains Olof Widlund Courant Institute of Mathematical Sciences New York University http://www.cs.nyu.edu/cs/faculty/widlund/
More informationIterative Helmholtz Solvers
Iterative Helmholtz Solvers Scalability of deflation-based Helmholtz solvers Delft University of Technology Vandana Dwarka March 28th, 2018 Vandana Dwarka (TU Delft) 15th Copper Mountain Conference 2018
More informationA FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem
A FETI-DP Method for Mortar Finite Element Discretization of a Fourth Order Problem Leszek Marcinkowski 1 and Nina Dokeva 2 1 Department of Mathematics, Warsaw University, Banacha 2, 02 097 Warszawa, Poland,
More informationSome Domain Decomposition Methods for Discontinuous Coefficients
Some Domain Decomposition Methods for Discontinuous Coefficients Marcus Sarkis WPI RICAM-Linz, October 31, 2011 Marcus Sarkis (WPI) Domain Decomposition Methods RICAM-2011 1 / 38 Outline Discretizations
More informationSchwarz-type methods and their application in geomechanics
Schwarz-type methods and their application in geomechanics R. Blaheta, O. Jakl, K. Krečmer, J. Starý Institute of Geonics AS CR, Ostrava, Czech Republic E-mail: stary@ugn.cas.cz PDEMAMIP, September 7-11,
More informationA Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation
A Domain Decomposition Based Jacobi-Davidson Algorithm for Quantum Dot Simulation Tao Zhao 1, Feng-Nan Hwang 2 and Xiao-Chuan Cai 3 Abstract In this paper, we develop an overlapping domain decomposition
More informationTHEORETICAL AND NUMERICAL COMPARISON OF PROJECTION METHODS DERIVED FROM DEFLATION, DOMAIN DECOMPOSITION AND MULTIGRID METHODS
THEORETICAL AND NUMERICAL COMPARISON OF PROJECTION METHODS DERIVED FROM DEFLATION, DOMAIN DECOMPOSITION AND MULTIGRID METHODS J.M. TANG, R. NABBEN, C. VUIK, AND Y.A. ERLANGGA Abstract. For various applications,
More informationFakultät für Mathematik und Informatik
Fakultät für Mathematik und Informatik Preprint 2015-17 Axel Klawonn, Martin Lanser and Oliver Rheinbach A Highly Scalable Implementation of Inexact Nonlinear FETI-DP without Sparse Direct Solvers ISSN
More informationParallel sparse linear solvers and applications in CFD
Parallel sparse linear solvers and applications in CFD Jocelyne Erhel Joint work with Désiré Nuentsa Wakam () and Baptiste Poirriez () SAGE team, Inria Rennes, France journée Calcul Intensif Distribué
More informationAggregation-based algebraic multigrid
Aggregation-based algebraic multigrid from theory to fast solvers Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire CEMRACS, Marseille, July 18, 2012 Supported by the Belgian FNRS
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
More informationA User Friendly Toolbox for Parallel PDE-Solvers
A User Friendly Toolbox for Parallel PDE-Solvers Gundolf Haase Institut for Mathematics and Scientific Computing Karl-Franzens University of Graz Manfred Liebmann Mathematics in Sciences Max-Planck-Institute
More informationSolving Ax = b, an overview. Program
Numerical Linear Algebra Improving iterative solvers: preconditioning, deflation, numerical software and parallelisation Gerard Sleijpen and Martin van Gijzen November 29, 27 Solving Ax = b, an overview
More informationSolving Large Nonlinear Sparse Systems
Solving Large Nonlinear Sparse Systems Fred W. Wubs and Jonas Thies Computational Mechanics & Numerical Mathematics University of Groningen, the Netherlands f.w.wubs@rug.nl Centre for Interdisciplinary
More 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 informationAlgebraic Adaptive Multipreconditioning applied to Restricted Additive Schwarz
Algebraic Adaptive Multipreconditioning applied to Restricted Additive Schwarz Nicole Spillane To cite this version: Nicole Spillane. Algebraic Adaptive Multipreconditioning applied to Restricted Additive
More informationAN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS. Gérard MEURANT CEA
Marrakech Jan 2003 AN INTRODUCTION TO DOMAIN DECOMPOSITION METHODS Gérard MEURANT CEA Introduction Domain decomposition is a divide and conquer technique Natural framework to introduce parallelism in the
More informationVarious ways to use a second level preconditioner
Various ways to use a second level preconditioner C. Vuik 1, J.M. Tang 1, R. Nabben 2, and Y. Erlangga 3 1 Delft University of Technology Delft Institute of Applied Mathematics 2 Technische Universität
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationDomain decomposition, hybrid methods, coarse space corrections
Domain decomposition, hybrid methods, coarse space corrections Victorita Dolean, Pierre Jolivet & Frédéric Nataf UNICE, ENSEEIHT & UPMC and CNRS CEMRACS 2016 18 22 July 2016 Outline 1 Introduction 2 Schwarz
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY, USA. Dan_Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More information18. Balancing Neumann-Neumann for (In)Compressible Linear Elasticity and (Generalized) Stokes Parallel Implementation
Fourteenth nternational Conference on Domain Decomposition Methods Editors: smael Herrera, David E Keyes, Olof B Widlund, Robert Yates c 23 DDMorg 18 Balancing Neumann-Neumann for (n)compressible Linear
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationUsing an Auction Algorithm in AMG based on Maximum Weighted Matching in Matrix Graphs
Using an Auction Algorithm in AMG based on Maximum Weighted Matching in Matrix Graphs Pasqua D Ambra Institute for Applied Computing (IAC) National Research Council of Italy (CNR) pasqua.dambra@cnr.it
More informationParallel scalability of a FETI DP mortar method for problems with discontinuous coefficients
Parallel scalability of a FETI DP mortar method for problems with discontinuous coefficients Nina Dokeva and Wlodek Proskurowski University of Southern California, Department of Mathematics Los Angeles,
More informationA decade of fast and robust Helmholtz solvers
A decade of fast and robust Helmholtz solvers Werkgemeenschap Scientific Computing Spring meeting Kees Vuik May 11th, 212 1 Delft University of Technology Contents Introduction Preconditioning (22-28)
More informationEfficient multigrid solvers for mixed finite element discretisations in NWP models
1/20 Efficient multigrid solvers for mixed finite element discretisations in NWP models Colin Cotter, David Ham, Lawrence Mitchell, Eike Hermann Müller *, Robert Scheichl * * University of Bath, Imperial
More informationMATRICES. a m,1 a m,n A =
MATRICES Matrices are rectangular arrays of real or complex numbers With them, we define arithmetic operations that are generalizations of those for real and complex numbers The general form a matrix of
More informationAn efficient multigrid solver based on aggregation
An efficient multigrid solver based on aggregation Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Graz, July 4, 2012 Co-worker: Artem Napov Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationMultispace and Multilevel BDDC. Jan Mandel University of Colorado at Denver and Health Sciences Center
Multispace and Multilevel BDDC Jan Mandel University of Colorado at Denver and Health Sciences Center Based on joint work with Bedřich Sousedík, UCDHSC and Czech Technical University, and Clark R. Dohrmann,
More informationAMG for a Peta-scale Navier Stokes Code
AMG for a Peta-scale Navier Stokes Code James Lottes Argonne National Laboratory October 18, 2007 The Challenge Develop an AMG iterative method to solve Poisson 2 u = f discretized on highly irregular
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 informationAn Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne
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 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 informationIsogEometric Tearing and Interconnecting
IsogEometric Tearing and Interconnecting Christoph Hofer and Ludwig Mitter Johannes Kepler University, Linz 26.01.2017 Doctoral Program Computational Mathematics Numerical Analysis and Symbolic Computation
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 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 informationPreprint ISSN
Fakultät für Mathematik und Informatik Preprint 2015-14 Allison H. Baker, Axel Klawonn, Tzanio Kolev, Martin Lanser, Oliver Rheinbach, and Ulrike Meier Yang Scalability of Classical Algebraic Multigrid
More informationBath Institute For Complex Systems
BICS Bath Institute for Complex Systems Energy Minimizing Coarse Spaces for Two-level Schwarz Methods for Multiscale PDEs J. Van lent, R. Scheichl and I.G. Graham Bath Institute For Complex Systems Preprint
More informationLecture 8: Fast Linear Solvers (Part 7)
Lecture 8: Fast Linear Solvers (Part 7) 1 Modified Gram-Schmidt Process with Reorthogonalization Test Reorthogonalization If Av k 2 + δ v k+1 2 = Av k 2 to working precision. δ = 10 3 2 Householder Arnoldi
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 information