Dune. Patrick Leidenberger. 13. May Distributed and Unified Numerics Environment. 13. May 2009
|
|
- Maud Mosley
- 5 years ago
- Views:
Transcription
1 13. May / 25 Dune Distributed and Unified Numerics Environment Patrick Leidenberger 13. May 2009
2 13. May / 25 Table of contents 1 Introduction 2 Programming Techniques 3 Dune Modules 4 Examples 5 Parallel Example 6 Conclusion
3 13. May / 25 What is Dune toolbox for solving PDEs modular for grid-based methods supports FE, FV, FD, DG written in modern C++ slim and common interface to different libraries open source (GPL v2 + runtime exception)
4 13. May / 25 Design Principles Efficiency: low overhead Modularity: maintainability and software reuse Flexibility: separate data structures and algorithms
5 13. May / 25 Developers Group of Peter Bastian, IWR Heidelberg Peter Bastian Markus Blatt Christian Engwer Group of Dietmar Kröner, IAM Freiburg Andreas Dedner Robert Klöfkorn Martin Nolte Group of Ralf Kornhuber, Institut für Mathematik, FU Berlin Oliver Sander Group of Mario Ohlberger, Inst. Num. Angew. Math., Münster Mario Ohlberger
6 13. May / 25 Templates template <typename T> T max(t x, T y ) { i f ( x < y ) return y ; e l s e return x ; } cout << max ( 7, 4 2 ) ; cout << max<double >(7.42, ) ;
7 13. May / 25 Template Metaprogramming template<long B, unsigned long E> s t r u c t pow { s t a t i c long const v a l u e = B*pow<B, E 1 >:: v a l u e ; } ; template<long B> s t r u c t pow<b, 0> { s t a t i c long const v a l u e = 1 ; } ; s t d : : cout << pow<2, 8 >:: v a l u e << s t d : : e n d l ; } s t d : : cout << 256 << s t d : : e n d l ;
8 Dune Modules 13. May / 25
9 13. May / 25 Grid Interfaces AlbertaGrid finite element toolbox ALBERTA ( 1d/2d/3/ simplicial meshes, bisection refinement ALU3dGrid 2d/3d hexa/tetra meshes, adaptive and parallel SGrid structured grid, arbitrary dimension UGGrid 2d/3d unstructured, adaptive and parallel meshes YaspGrid structured grid, parallel and arbitrary dimension periodic boundary conditions...
10 13. May / 25 Grid Traversal Iterate over all elements on the leaf grid: typedef GridType : : Codim <0>:: L e a f I t e r a t o r E l e m e n t L e a f I t e r a t o r ; f o r { } ( E l e m e n t L e a f I t e r a t o r i t=g r i d. template l e a f b e g i n <0>(); i t!= g r i d. template l e a f e n d <0>(); ++i t ) s t d : : cout << element : << i t >type ( ) << s t d : : e n d l ;
11 13. May / 25 Quadrature: Integrate a function f over an element *it Dune : : GeometryType gt = i t >type ( ) ; const Dune : : QuadratureRule <double, dim>& r u l e = Dune : : QuadratureRules <double, dim >:: r u l e ( gt, p ) ; double r e s u l t = 0 ; f o r ( i n t i = 0 ; i < r u l e. s i z e ( ) ; ++i ) { F i e l d V e c t o r <double, dim> g l o b a l P o s i t i o n = i t >geometry ( ). g l o b a l ( r u l e [ i ]. p o s i t i o n ( ) ) ; double f V a l = f ( g l o b a l P o s i t i o n ) ; double weight = r u l e [ i ]. weight ( ) ; double d e t j a c = i t >geometry ( ). i n t e g r a t i o n E l e m e n t ( r u l e [ i ]. p o s i t i o n ( ) ) ; r e s u l t += f v a l * weight * d e t j a c ; }
12 13. May / 25 Iterative Solver Template Library ISTL vector matrix, dense and sparse (BCRS) containers are hierarchic solvers realized via operator concept
13 13. May / 25 Iterative Solver Template Library ISTL vector matrix, dense and sparse (BCRS) containers are hierarchic solvers realized via operator concept Preconditioners Jacobi method successive overrelaxation (SOR) symmetric SOR incomplete LU decomposition (ILU) algebraic multigrid method additive overlapping Schwarz
14 13. May / 25 Iterative Solver Template Library ISTL vector matrix, dense and sparse (BCRS) containers are hierarchic solvers realized via operator concept Preconditioners Jacobi method successive overrelaxation (SOR) symmetric SOR incomplete LU decomposition (ILU) algebraic multigrid method additive overlapping Schwarz Solvers apply preconditioner multiple time preconditioned gradient method preconditioned conjugate gradient method preconditioned biconjugate gradient stabilized method
15 13. May / 25 Laplace Equation with Dune PDE Lab ϕ = 0 fast implementation few code lines define boundary condition define value on boundary (define FE operator) write main program
16 13. May / 25 Laplace Equation with Dune PDE Lab ϕ = 0 fast implementation few code lines define boundary condition define value on boundary (define FE operator) write main program
17 13. May / 25 Laplace Equation: Implement Boundary Condition 1 template<typename GV> c l a s s B : p u b l i c Dune : : PDELab : : BoundaryGridFunctionBase <Dune : : PDELab : : B o u n d a r y G r i d F u n c t i o n T r a i t s <GV, int, 1, Dune : : F i e l d V e c t o r <int,1 >, B<GV> >
18 13. May / 25 Laplace Equation: Implement Boundary Condition 2 template<typename I > i n l i n e void e v a l u a t e ( const I& ig, const typename T r a i t s : : DomainType& x, typename T r a i t s : : RangeType& y ) const { Dune : : F i e l d V e c t o r <typename GV : : G r i d : : ctype,gv : : dimension > xg = i g. geometry ( ). g l o b a l ( x ) ; } y = 0 ; // no D i r i c h l e t // top : D i r i c h l e t i f ( xg [ 1 ] >= (1.0 1 e 6)) { y = 1 ; } // bottom : D i r i c h l e t i f ( xg [ 1 ] <= 1e 6) { y = 1 ; }
19 13. May / 25 Laplace Equation: Value on Boundary y = x [ 1 ] ; i f ( x [ 0 ] > 0. 5 ) { y *= 2 ; }
20 13. May / 25 Laplace Equation Main Program: Define Grid Dune : : F i e l d V e c t o r <double,2> L ( 1. 0 ) ; Dune : : F i e l d V e c t o r <int,2> N( 1 ) ; Dune : : F i e l d V e c t o r <bool,2> P( f a l s e ) ; Dune : : YaspGrid <2> g r i d ( L, N, P, 0 ) ; g r i d. g l o b a l R e f i n e ( 5 ) ; // get view typedef Dune : : YaspGrid <2>:: L e a f G r i d V i ew GV; const GV& gv=g r i d. l e a f V i e w ( ) ;
21 13. May / 25 Laplace Equation Main Program: Execution // make f i n i t e element map typedef GV : : G r i d : : c t y p e DF ; typedef Q1LocalFiniteElementMap <DF, double> FEM; FEM fem ; l a p l a c e d i r i c h l e t <GV,FEM, Q1Constraints > ( gv, fem, 3, l a p l a c e y a s p Q 1 2 d ) ;
22 Laplace Equation: Result on 2d Simplex Mesh 13. May / 25
23 Laplace Equation: Result on 2d Simplex Mesh 13. May / 25
24 Laplace Equation: Result on 2d Quadrilateral Mesh 13. May / 25
25 Laplace Equation: Result on 2d Quadrilateral Mesh 13. May / 25
26 Laplace Equation: Result on 3d Quadrilateral Mesh 13. May / 25
27 13. May / 25 2d Frequency Domain FE solve: ( ) 1 µ E z (x, y) ( ω 2 ɛ iωσ ) E z (x, y) = iωj ext,z higher order scalar base functions arbitrary meshes boundary conditions: PEC, ABC, waveguide port
28 13. May / 25 3d Time Domain FE solve: 1 µ E + σ E t 2 + ɛ 0 ɛ r t 2 E = t J ext vector base functions arbitrary meshes dispersion (Debye and Drude model) parallel (distributed memory)
29 3d Time Domain FE Work by B. Oswald and P. Leidenberger 13. May / 25
30 3d Time Domain FE Work by B. Oswald and P. Leidenberger 13. May / 25
31 Convective flow 13. May / 25
32 Convective flow Density-driven flow (by Peter Bastian) cell-centered finite volume scheme YaspGrid, cells, 384 processors 9000 timesteps, 3 days of computation time 13. May 2009 leidenberger@ifh.ee.ethz.ch 24 / 25
33 13. May / 25 Conclusion modern and efficient C++ framework different mesh backends suitable for large simulations suitable for teaching good developer community
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationGENERIC IMPLEMENTATION OF FINITE ELEMENT METHODS IN THE DISTRIBUTED AND UNIFIED NUMERICS ENVIRONMENT (DUNE)
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 2, P AGES 294 315 GENERIC IMPLEMENTATION OF FINITE ELEMENT METHODS IN THE DISTRIBUTED AND UNIFIED NUMERICS ENVIRONMENT (DUNE) Peter Bastian, Felix Heimann and Sven
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 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 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 informationNumerik 2 Motivation
Numerik 2 Motivation P. Bastian Universität Heidelberg Interdisziplinäres Zentrum für Wissenschaftliches Rechnen Im Neuenheimer Feld 368, D-69120 Heidelberg email: Peter.Bastian@iwr.uni-heidelberg.de April
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 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 informationMultipole-Based Preconditioners for Sparse Linear Systems.
Multipole-Based Preconditioners for Sparse Linear Systems. Ananth Grama Purdue University. Supported by the National Science Foundation. Overview Summary of Contributions Generalized Stokes Problem Solenoidal
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 informationITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS
ITERATIVE METHODS FOR SPARSE LINEAR SYSTEMS YOUSEF SAAD University of Minnesota PWS PUBLISHING COMPANY I(T)P An International Thomson Publishing Company BOSTON ALBANY BONN CINCINNATI DETROIT LONDON MADRID
More informationParallel Discontinuous Galerkin Method
Parallel Discontinuous Galerkin Method Yin Ki, NG The Chinese University of Hong Kong Aug 5, 2015 Mentors: Dr. Ohannes Karakashian, Dr. Kwai Wong Overview Project Goal Implement parallelization on Discontinuous
More informationFlow an open source research tool for reservoir simulation
Flow an open source research tool for reservoir simulation Robert Klöfkorn IRIS & The National IOR Centre of Norway April 27, 2016 Task 6 Reservoir Simulation Tools Pål (PostDoc UiS) Trine (PostDoc IRIS)
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
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 informationIMPLEMENTATION OF A PARALLEL AMG SOLVER
IMPLEMENTATION OF A PARALLEL AMG SOLVER Tony Saad May 2005 http://tsaad.utsi.edu - tsaad@utsi.edu PLAN INTRODUCTION 2 min. MULTIGRID METHODS.. 3 min. PARALLEL IMPLEMENTATION PARTITIONING. 1 min. RENUMBERING...
More informationAdvanced numerical methods for nonlinear advectiondiffusion-reaction. Peter Frolkovič, University of Heidelberg
Advanced numerical methods for nonlinear advectiondiffusion-reaction equations Peter Frolkovič, University of Heidelberg Content Motivation and background R 3 T Numerical modelling advection advection
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 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 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 informationAPPLIED NUMERICAL LINEAR ALGEBRA
APPLIED NUMERICAL LINEAR ALGEBRA James W. Demmel University of California Berkeley, California Society for Industrial and Applied Mathematics Philadelphia Contents Preface 1 Introduction 1 1.1 Basic Notation
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 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 informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
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 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 informationFinite Elements (Numerical Solution of Partial Differential Equations)
Finite Elements (Numerical Solution of Partial Differential Equations) P. Bastian Universität Heidelberg Interdisziplinäres Zentrum für Wissenschaftliches Rechnen Im Neuenheimer Feld 205, D-69120 Heidelberg
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 Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More 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 informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
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 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 informationTask 6: Reservoir simulation tools.
Task 6: Reservoir simulation tools Randi.Valestrand@iris.no Theme 2: Tasks and Task Leaders Task6 Robert Köfkorn Senior Researcher IRIS Robert.Kloefkorn@iris.no DEVELOPMENT & IMPROVED METHOLOGY FOR FULL
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 informationUsing DUNE-ACFEM for Non-smooth Minimization of Bounded Variation Functions
Archive of Numerical Software vol. 5, no. 1, 2017, pages 3 19 DOI: 10.11588/ans.2017.1.27475 Using DUNE-ACFEM for Non-smooth Minimization of Bounded Variation Functions Martin Alkämper 1 and Andreas Langer
More information- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline
- Part 4 - Multicore and Manycore Technology: Chances and Challenges Vincent Heuveline 1 Numerical Simulation of Tropical Cyclones Goal oriented adaptivity for tropical cyclones ~10⁴km ~1500km ~100km 2
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
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 informationLecture 17: Iterative Methods and Sparse Linear Algebra
Lecture 17: Iterative Methods and Sparse Linear Algebra David Bindel 25 Mar 2014 Logistics HW 3 extended to Wednesday after break HW 4 should come out Monday after break Still need project description
More informationDirect and Incomplete Cholesky Factorizations with Static Supernodes
Direct and Incomplete Cholesky Factorizations with Static Supernodes AMSC 661 Term Project Report Yuancheng Luo 2010-05-14 Introduction Incomplete factorizations of sparse symmetric positive definite (SSPD)
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
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 informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationAn Empirical Comparison of Graph Laplacian Solvers
An Empirical Comparison of Graph Laplacian Solvers Kevin Deweese 1 Erik Boman 2 John Gilbert 1 1 Department of Computer Science University of California, Santa Barbara 2 Scalable Algorithms Department
More 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 informationA Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow
A Space-Time Multigrid Solver Methodology for the Optimal Control of Time-Dependent Fluid Flow Michael Köster, Michael Hinze, Stefan Turek Michael Köster Institute for Applied Mathematics TU Dortmund Trier,
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 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 informationUniversität Dortmund UCHPC. Performance. Computing for Finite Element Simulations
technische universität dortmund Universität Dortmund fakultät für mathematik LS III (IAM) UCHPC UnConventional High Performance Computing for Finite Element Simulations S. Turek, Chr. Becker, S. Buijssen,
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 informationConjugate Gradients: Idea
Overview Steepest Descent often takes steps in the same direction as earlier steps Wouldn t it be better every time we take a step to get it exactly right the first time? Again, in general we choose a
More informationNUMERICAL METHODS FOR ENGINEERING APPLICATION
NUMERICAL METHODS FOR ENGINEERING APPLICATION Second Edition JOEL H. FERZIGER A Wiley-Interscience Publication JOHN WILEY & SONS, INC. New York / Chichester / Weinheim / Brisbane / Singapore / Toronto
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 informationA Numerical Study of Some Parallel Algebraic Preconditioners
A Numerical Study of Some Parallel Algebraic Preconditioners Xing Cai Simula Research Laboratory & University of Oslo PO Box 1080, Blindern, 0316 Oslo, Norway xingca@simulano Masha Sosonkina University
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 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 informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
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 informationBreaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points
Breaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points Michael Griebel Christian Rieger Peter Zaspel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationMPI parallel implementation of CBF preconditioning for 3D elasticity problems 1
Mathematics and Computers in Simulation 50 (1999) 247±254 MPI parallel implementation of CBF preconditioning for 3D elasticity problems 1 Ivan Lirkov *, Svetozar Margenov Central Laboratory for Parallel
More informationIndex ": :", 50 "? :", 14 <, 14 =, 8, 14 >, 14, 55., 14., 49 ;, 6 %, 9, 37, 176 &, 58, 150 :, 52, 79 *, 8, 22, 150 "++", 10
Index , 14, 55., 14., 49 ;, 6 %, 9, 37, 176 &, 58, 150 :, 52, 79 *, 8, 22, 150 "++", 10 "+=", 10 in complex, 68 in dynamic vector, 672 in matrix, 670 in point, 63 in sparse matrix, 692
More informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
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 informationJ.I. Aliaga 1 M. Bollhöfer 2 A.F. Martín 1 E.S. Quintana-Ortí 1. March, 2009
Parallel Preconditioning of Linear Systems based on ILUPACK for Multithreaded Architectures J.I. Aliaga M. Bollhöfer 2 A.F. Martín E.S. Quintana-Ortí Deparment of Computer Science and Engineering, Univ.
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 informationA DG-based incompressible Navier-Stokes solver
A DG-based incompressible Navier-Stokes solver Higer-order DG metods and finite element software for modern arcitectures - Bat 2016 Marian Piatkowski 1 Peter Bastian 1 1 Interdisciplinary Center for Scientific
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 informationAn Implementation of Hybrid Discontinuous Galerkin Methods in DUNE
Aachen Institute for Advanced Study in Computational Engineering Science Preprint: AICES-2010/12-1 15/December/2010 An Implementation of Hybrid Discontinuous Galerkin Methods in DUNE C. Waluga and H. Egger
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
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 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 informationMultilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota
Multilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota SIAM CSE Boston - March 1, 2013 First: Joint work with Ruipeng Li Work
More informationPreconditioners for the incompressible Navier Stokes equations
Preconditioners for the incompressible Navier Stokes equations C. Vuik M. ur Rehman A. Segal Delft Institute of Applied Mathematics, TU Delft, The Netherlands SIAM Conference on Computational Science and
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 informationIFISS: A computational laboratory for investigating incompressible flow problems. Elman, Howard C. and Ramage, Alison and Silvester, David J.
IFISS: A computational laboratory for investigating incompressible flow problems Elman, Howard C. and Ramage, Alison and Silvester, David J. 22 MIMS EPrint: 22.8 Manchester Institute for Mathematical Sciences
More informationThe Deflation Accelerated Schwarz Method for CFD
The Deflation Accelerated Schwarz Method for CFD J. Verkaik 1, C. Vuik 2,, B.D. Paarhuis 1, and A. Twerda 1 1 TNO Science and Industry, Stieltjesweg 1, P.O. Box 155, 2600 AD Delft, The Netherlands 2 Delft
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 informationAtomistic effect of delta doping layer in a 50 nm InP HEMT
J Comput Electron (26) 5:131 135 DOI 1.7/s1825-6-8832-3 Atomistic effect of delta doping layer in a 5 nm InP HEMT N. Seoane A. J. García-Loureiro K. Kalna A. Asenov C Science + Business Media, LLC 26 Abstract
More information1.0. T (ev) 0.5. z [m] y [m] x [m]
W7-X edge modelling with the 3D SOL fluid code BoRiS M. Borchardt, J. Riemann, R. Schneider, X. Bonnin Max-Planck-Institut für Plasmaphysik, Teilinstitut Greifswald EURATOM Association, D 17491 Greifswald,
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 informationThe FEniCS Project Philosophy, current status and future plans
The FEniCS Project Philosophy, current status and future plans Anders Logg logg@simula.no Simula Research Laboratory FEniCS 06 in Delft, November 8-9 2006 Outline Philosophy Recent updates Future plans
More informationA robust multilevel approximate inverse preconditioner for symmetric positive definite matrices
DICEA DEPARTMENT OF CIVIL, ENVIRONMENTAL AND ARCHITECTURAL ENGINEERING PhD SCHOOL CIVIL AND ENVIRONMENTAL ENGINEERING SCIENCES XXX CYCLE A robust multilevel approximate inverse preconditioner for symmetric
More informationSolving PDEs with CUDA Jonathan Cohen
Solving PDEs with CUDA Jonathan Cohen jocohen@nvidia.com NVIDIA Research PDEs (Partial Differential Equations) Big topic Some common strategies Focus on one type of PDE in this talk Poisson Equation Linear
More informationReview: From problem to parallel algorithm
Review: From problem to parallel algorithm Mathematical formulations of interesting problems abound Poisson s equation Sources: Electrostatics, gravity, fluid flow, image processing (!) Numerical solution:
More informationIterative Solution methods
p. 1/28 TDB NLA Parallel Algorithms for Scientific Computing Iterative Solution methods p. 2/28 TDB NLA Parallel Algorithms for Scientific Computing Basic Iterative Solution methods The ideas to use iterative
More informationIncomplete Block LU Preconditioners on Slightly Overlapping. E. de Sturler. Delft University of Technology. Abstract
Incomplete Block LU Preconditioners on Slightly Overlapping Subdomains for a Massively Parallel Computer E. de Sturler Faculty of Technical Mathematics and Informatics Delft University of Technology Mekelweg
More informationOn the performance of the Algebraic Optimized Schwarz Methods with applications
On the performance of the Algebraic Optimized Schwarz Methods with applications Report 12-1-18 Department of Mathematics Temple University October 212 This report is available at http://www.math.temple.edu/~szyld
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
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 informationIterative methods for Linear System of Equations. Joint Advanced Student School (JASS-2009)
Iterative methods for Linear System of Equations Joint Advanced Student School (JASS-2009) Course #2: Numerical Simulation - from Models to Software Introduction In numerical simulation, Partial Differential
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More 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 informationLarge Scale Simulation of Flow and Transport in Porous Media
Large Scale Simulation of Flow and Transport in Porous Media Olaf Ippisch Institut de Mathematiques École Polytechnique de Clausthal Séminaire de la Modélisation 13 Février 2018 O. Ippisch (IfM, TU Clausthal)
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationANALYSIS OF AUGMENTED LAGRANGIAN-BASED PRECONDITIONERS FOR THE STEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
ANALYSIS OF AUGMENTED LAGRANGIAN-BASED PRECONDITIONERS FOR THE STEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS MICHELE BENZI AND ZHEN WANG Abstract. We analyze a class of modified augmented Lagrangian-based
More information