Parallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks
|
|
- Katrina Parks
- 5 years ago
- Views:
Transcription
1 Parallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks Gopal Yalla Björn Engquist University of Texas at Austin Institute for Computational Engineering and Sciences April 17, 2018 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
2 Outline 1 Parallel in Time Algorithms 2 Parareal Algorithm 3 Phase-accurate Coarse Solver 4 Numerical Examples 5 Future Work Yalla & Engquist Parallel in Time Algorithms April 17, / 20
3 Parallel in Time Algorithms 1 Required for massively parallel simulations of systems governed by time-dependent dynamical systems, e.g., molecular dynamics. Challenge arises due to causality in time. 1 Martin J Gander. 50 years of time parallel time integration. In: Multiple Shooting and Time Domain Decomposition Methods. Springer, 2015, pp Yalla & Engquist Parallel in Time Algorithms April 17, / 20
4 Parallel in Time Algorithms 1 1 Multiple Shooting Methods 2 Domain Decomposition & Waveform Relaxation Methods 3 Multigrid Based Methods 1 Martin J Gander. 50 years of time parallel time integration. In: Multiple Shooting and Time Domain Decomposition Methods. Springer, 2015, pp Yalla & Engquist Parallel in Time Algorithms April 17, / 20
5 The Parareal Algorithm Consider the dynamical system: u (t) = f(u(t)), t (t 0, t f ) u(t 0 ) = u 0 Divide the time domain (t 0, t f ) into N equal subdomains Ω n = (T n, T n+1 ). u n(t) = f(u n (t)) t (T n, T n+1 ) u(t n ) = U n This is simply a shooting method in time, and it equivalent to solving: U 0 u 0 U 1 φ T0 (U 0 ). = 0 (3) U N 1 φ TN 2 (U N 2 )) where φ Tn (U n ) is solution of (1) with initial condition U n after time T n. Yalla & Engquist Parallel in Time Algorithms April 17, / 20 (1) (2)
6 The Parareal Algorithm Consider the dynamical system: u (t) = f(u(t)), t (t 0, t f ) u(t 0 ) = u 0 Divide the time domain (t 0, t f ) into N equal subdomains Ω n = (T n, T n+1 ). Parareal Algorithm (1) U k+1 0 = u 0 U k+1 n = G(U k+1 n 1 ) + [F(U k n 1) G(U k n 1) ] Un k is the solution u(t n) on the k th iteration of the algorithm. G denotes a coarse solver F denotes a fine solver Yalla & Engquist Parallel in Time Algorithms April 17, / 20
7 U k 0 = u 0 U k n = G(U k n 1 ) + [ F(Un 1 k 1 ] k 1 ) G(Un 1 ) Video Source: Yalla & Engquist Parallel in Time Algorithms April 17, / 20
8 U k 0 = u 0 U k n = G(U k n 1 ) + [ F(Un 1 k 1 ] k 1 ) G(Un 1 ) Efficiency = Iterations for Convergence / Total Number of Steps Video Source: Yalla & Engquist Parallel in Time Algorithms April 17, / 20
9 Drawbacks of the Parareal Algorithm For Hamiltonian systems, very high accuracy is required for the coarse solver! 2 A good coarse approximation is needed. Typical coarse integrators fail to predict phase correctly, leading to blow up of error during correction step with fine solver. 2 Martin J Gander and Ernst Hairer. Analysis for parareal algorithms applied to Hamiltonian differential equations. In: Journal of Computational and Applied Mathematics 259 (2014), pp Yalla & Engquist Parallel in Time Algorithms April 17, / 20
10 Phase Plane Map Consider a set of M initial conditions (training points) {u i 0 }M i=1 and a set of M corresponding (target points) {v i } M i=1 defined by, v i = F ( u i ) 0 i = 1,..., M The set {u i 0 vi } M i=1 acts as a type of look-up table for future initials conditions 3. Definition Define a phase plane map G map ( { u i 0 vi} M i=1, U n) to be a coarse solver that uses the information contained in { u i 0 vi} M i=1 and a point U n at time T n, to determine the solution U n+1 at time T n+1. G map can be defined through, e.g., Interpolation or a Neural Network. 3 Jürg Nievergelt. Parallel methods for integrating ordinary differential equations. In: Communications of the ACM 7.12 (1964), pp Yalla & Engquist Parallel in Time Algorithms April 17, / 20
11 Phase Plane Map v i = F ( u i ) 0 U n+1 G map ( { u i 0 vi} M i=1, U n) Yalla & Engquist Parallel in Time Algorithms April 17, / 20
12 Phase Plane Map The computation of v i = F(u i 0 ) is embarrassingly parallel. Yalla & Engquist Parallel in Time Algorithms April 17, / 20
13 Phase Plane Map For autonomous ODEs, G map can be applied over each subdomain Ω n. Else, generate {u i j vi } M i=1 for each Ω j, j = 0,..., N 1. Yalla & Engquist Parallel in Time Algorithms April 17, / 20
14 Phase Plane Map 1 How to select {u i 0 }M i=1? 2 How to select M? 3 Which method is best for G map? Yalla & Engquist Parallel in Time Algorithms April 17, / 20
15 Numerical Examples Traditional parareal algorithm = G = RK4 with step size T/N Modified parareal algorithm = G = phase plane map with step size T/N For both, the fine solver is an adaptive RK45 solver. T = Total Simulation Time. N = Number of time-subdomains Ω n. Yalla & Engquist Parallel in Time Algorithms April 17, / 20
16 Harmonic Oscillator Consider a simple harmonic oscillator: with ε = 0.01 T = 70 (1000 revolutions) x(0) = [ 2/9, 2/3] T x 1 = 1 ε x 2 x 2 = 1 ε x 1 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
17 Harmonic Oscillator N = 100, M = 2 Iterations: 1 N = 10, 000 Iterations: 42 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
18 Lennard-Jones System r = 2 m F (r) F (r) = V (r) [ (σ ) 12 ( σ ) ] 6 V (r) = 4ε r r ε = well depth σ = zero distance r min = 1.12 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
19 Lennard-Jones System Iterations: 10 Iterations: 53 T = 50 N = 400 r 0 = [1.2, 0] T M=50 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
20 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
21 ( ) r = 1.12 M Iterations ( ) r = 1.2 M Iterations ( ) r = 1.4 M Iterations Yalla & Engquist Parallel in Time Algorithms April 17, / 20
22 High Dimensional Harmonic Oscillator ẋ = 1 ε Ax where x = [x 1, v 1, x 2, v 2, x 3, v 3, x 4, v 4 ] T, ε = 0.01, and A R 8 8 is defined as, 0 a a b A = 0 0 b c c d d 0 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
23 High Dimensional Harmonic Oscillator ẋ = 1 ε Ax where x = [x 1, v 1, x 2, v 2, x 3, v 3, x 4, v 4 ] T, ε = 0.01, and A R 8 8 is defined as, 0 a a b A = 0 0 b c c d d 0 *Neural network defined with logistic activation function, size , with a limited memory BFGS solver for optimization. Yalla & Engquist Parallel in Time Algorithms April 17, / 20
24 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
25 Localized Multiscale Problem 1 Type A: Problems that contain local defects or singularities. Microscale model needed only near defects. 2 Type B: Microscale model needed everywhere as a supplement to macroscale model, e.g., resolve constitutive model. 4. Type A Example: An N-body Problem: m i q i = n j i m i m j (q j q i ) q j q i 3 m i denotes the mass of the i th body. q i denotes the position of the i th body. Let n = 2 and assume m 2 >> m 1. 4 Weinan E. Principles of multiscale modeling. Cambridge University Press, Yalla & Engquist Parallel in Time Algorithms April 17, / 20
26 Apply a traditional parareal algorithm with RK45 as fine scale solver with adaptive step size. RK4 as coarse solver with step size T/N. T = 0.5 N = 5 Iterations = 2 Yalla & Engquist Parallel in Time Algorithms April 17, / 20
27 Conclusion & Future Work Parallel in time algorithms needed for massively parallel simulations of time dependent dynamical systems. Results for parareal + phase plane map show promise. Demonstrate a proof of concept. Future Work Potential Improvements include use of modern sparse grid and adaptive methods for interpolation or optimizing neural network approach. More realistic examples in the realm of molecular dynamics. Scaling/speedup results on modern supercomputing platforms. Yalla & Engquist Parallel in Time Algorithms April 17, / 20
28 References Weinan E. Principles of multiscale modeling. Cambridge University Press, Martin J Gander. 50 years of time parallel time integration. In: Multiple Shooting and Time Domain Decomposition Methods. Springer, 2015, pp Martin J Gander and Ernst Hairer. Analysis for parareal algorithms applied to Hamiltonian differential equations. In: Journal of Computational and Applied Mathematics 259 (2014), pp Jacques-Louis Lions, Yvon Maday, et al. A parareal method in time discretization of pde s. In: Comptes Rendus de l Académie des Sciences-Series I-Mathematics (2001), pp Jürg Nievergelt. Parallel methods for integrating ordinary differential equations. In: Communications of the ACM 7.12 (1964), pp *Code will be available at Yalla & Engquist Parallel in Time Algorithms April 17, / 20
Multiscale dynamical systems and parareal algorithms
Colloquium du CERMICS Multiscale dynamical systems and parareal algorithms Bjorn Engquist (University of Texas at Austin) 7juin2018 Multiscale dynamical systems and parareal algorithms Bjorn Engquist,
More informationAsynchronous Parareal in time discretization for partial differential equations
Asynchronous Parareal in time discretization for partial differential equations Frédéric Magoulès, Guillaume Gbikpi-Benissan April 7, 2016 CentraleSupélec IRT SystemX Outline of the presentation 01 Introduction
More informationParareal Methods. Scott Field. December 2, 2009
December 2, 2009 Outline Matlab example : Advantages, disadvantages, and survey of usage in literature Overview Overview Notation Algorithm Model equation example Proposed by Lions, Maday, Turinici in
More informationParallel-in-time integrators for Hamiltonian systems
Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationHeterogeneous Multiscale Methods: A Review
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 2, No. 3, pp. 367-45 Commun. Comput. Phys. June 27 REVIEW ARTICLE Heterogeneous Multiscale Methods: A Review Weinan E 1,, Bjorn Engquist 2, Xiantao Li 3, Weiqing
More informationSpacetime Computing. A Dynamical Systems Approach to Turbulence. Bruce M. Boghosian. Department of Mathematics, Tufts University
Spacetime Computing A Dynamical Systems Approach to Turbulence Bruce M. Boghosian, Tufts University Presented at 3rd NA-HPC Roadmap Workshop, Royal Society, London January 26, 2009 Acknowledgements: Hui
More informationAn Adjoint Approach for Stabilizing the Parareal Method
An Adjoint Approach for Stabilizing the Parareal Method Feng Chen a, Jan S Hesthaven b, Yvon Maday c, Allan Nielsen b a Department of Mathematics, Baruch College, New York, NY 10010, USA b Mathematics
More informationAn asymptotic-parallel-in-time method for weather and climate prediction
An asymptotic-parallel-in-time method for weather and climate prediction Beth Wingate University of Exeter Terry Haut Los Alamos National Laboratory Jared Whitehead Brigham Young University Outline o A
More informationRATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS
RATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS Nabil R. Nassif Noha Makhoul Karam Yeran Soukiassian Mathematics Department, American University of Beirut IRISA, Campus Beaulieu,
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 informationc 2007 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 29, No. 2, pp. 556 578 c 27 Society for Industrial and Applied Mathematics ANALYSIS OF THE PARAREAL TIME-PARALLEL TIME-INTEGRATION METHOD MARTIN J. GANDER AND STEFAN VANDEWALLE
More informationPARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr
More informationHETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING
Proceedings of ALGORITMY 2009 pp. 219 225 HETEROGENEOUS MULTISCALE METHOD IN EDDY CURRENTS MODELING JÁN BUŠA, JR. AND VALDEMAR MELICHER Abstract. The induction of eddy currents in a conductive piece is
More informationMultiple time scale numerical methods for the inverted pendulum problem
Multiple time scale numerical methods for the inverted pendulum problem Richard Sharp 1, Yen-Hsi Tsai, and Bjorn Engquist 1 Program in Applied and Computational Mathematics, Princeton University, NJ 08544
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 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 informationDeep Learning for Partial Differential Equations (PDEs)
Deep Learning for Partial Differential Equations (PDEs) Kailai Xu kailaix@stanford.edu Bella Shi bshi@stanford.edu Shuyi Yin syin3@stanford.edu Abstract Partial differential equations (PDEs) have been
More informationPoisson Equation in 2D
A Parallel Strategy Department of Mathematics and Statistics McMaster University March 31, 2010 Outline Introduction 1 Introduction Motivation Discretization Iterative Methods 2 Additive Schwarz Method
More informationAcceleration of Time Integration
Acceleration of Time Integration edited version, with extra images removed Rick Archibald, John Drake, Kate Evans, Doug Kothe, Trey White, Pat Worley Research sponsored by the Laboratory Directed Research
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
More 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 informationOverlapping Schwarz Waveform Relaxation for Parabolic Problems
Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03038-7 Overlapping Schwarz Waveform Relaxation for Parabolic Problems Martin J. Gander 1. Introduction We analyze a new domain decomposition algorithm
More informationA GENERAL STRATEGY FOR DESIGNING SEAMLESS MULTISCALE METHODS
A GENERAL STRATEGY FOR DESIGNING SEAMLESS MULTISCALE METHODS WEINAN E, WEIQING REN, AND ERIC VANDEN-EIJNDEN Abstract. We present a new general framework for designing multiscale methods. Compared with
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 informationParallel Numerical Methods for Ordinary Differential Equations: a Survey
Parallel Numerical Methods for Ordinary Differential Equations: a Survey Svyatoslav I. Solodushkin 1,2 and Irina F. Iumanova 1 1 Ural Federal University, Yekaterinburg, Russia 2 Krasovskii Institute of
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
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
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 informationCoupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations
Coupling of Multi fidelity Models Applications to PNP cdft and local nonlocal Poisson equations P. Bochev, J. Cheung, M. D Elia, A. Frishknecht, K. Kim, M. Parks, M. Perego CM4 summer school, Stanford,
More informationParallel Simulation of Subsurface Fluid Flow
Parallel Simulation of Subsurface Fluid Flow Scientific Achievement A new mortar domain decomposition method was devised to compute accurate velocities of underground fluids efficiently using massively
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
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 informationDraft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov
TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL Nikolai Shegunov, Ivan Hristov March 9, Seminar in IMI,BAS,Soa 1/25 Mathematical model We consider the Hamiltonian :
More informationMultilevel 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 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 informationIntroduction to Multiscale Modeling
Introduction to Multiscale Modeling Bjorn Engquist University of Texas at Austin Summer School on Multiscale Modeling and Analysis, University of Texas at Austin, July 21 - August 8, 2008 Outline I Overview
More informationMéthodes de raffinement espace-temps
Méthodes de raffinement espace-temps Caroline Japhet Université Paris 13 & University of Maryland Laurence Halpern Université Paris 13 Jérémie Szeftel Ecole Normale Supérieure, Paris Pascal Omnes CEA,
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationSolving 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 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 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 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 informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationA Time-Parallel Implicit Method for Accelerating the Solution of Nonlinear Structural Dynamics Problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2008; 0:1 25 [Version: 2000/01/19 v2.0] A Time-Parallel Implicit Method for Accelerating the Solution of Nonlinear
More informationA Survey of Computational High Frequency Wave Propagation I
A Survey of Computational High Frequency Wave Propagation I Bjorn Engquist University of Texas at Austin CSCAMM Workshop on High Frequency Wave Propagation, University of Maryland, September 19-22, 2005
More informationOptimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits
Optimized Schwarz Waveform Relaxation: Roots, Blossoms and Fruits Laurence Halpern 1 LAGA, Université Paris XIII, 99 Avenue J-B Clément, 93430 Villetaneuse, France, halpern@math.univ-paris13.fr 1 Introduction:
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 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 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 informationIntroduction to Wavelet Based Numerical Homogenization
Introduction to Wavelet Based Numerical Homogenization Olof Runborg Numerical Analysis, School of Computer Science and Communication, KTH RTG Summer School on Multiscale Modeling and Analysis University
More informationAbsorbing Boundary Conditions for Nonlinear Wave Equations
Absorbing Boundary Conditions for Nonlinear Wave Equations R. B. Gibson December, 2010 1 1 Linear Wave Equations 1.1 Boundary Condition 1 To begin, we wish to solve where with the boundary condition u
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
More informationMultigrid Interpretations of the Parareal Algorithm Leading to an Overlapping Variant and MGRIT
Noname manuscript No (will be inserted by the editor) Multigrid Interpretations of the Parareal Algorithm Leading to an Overlapping Variant and MGRIT Martin J Gander Felix Kwok Hui Zhang Received: date
More informationTime Parallel Time Integration
Time Inegraion marin.gander@unige.ch Universiy of Geneva SSSTC Symposium, Geneva, Augus 17 Solving Evoluion Problems in Parallel? The Lorenz equaions: d d dy d dz d = σ +σy = z +r y = y bz 4 4 15 1 5 5
More informationA Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation
A Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation Andreas Frommer, Karsten Kahl, Stefan Krieg, Björn Leder and Matthias Rottmann Bergische Universität Wuppertal April 11, 2013
More informationPARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS
PARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS GUILLAUME BAL Abstract. This paper analyzes some properties of the parareal algorithm, which can be used to parallelize the time
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 informationCIMPA Summer School on Current Research in Finite Element Methods
CIMPA Summer School on Current Research in Finite Element Methods Local zooming techniques for the finite element method Alexei Lozinski Laboratoire de mathématiques de Besançon Université de Franche-Comté,
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
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 informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
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 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 informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part II Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More informationNonparametric density estimation for elliptic problems with random perturbations
Nonparametric density estimation for elliptic problems with random perturbations, DqF Workshop, Stockholm, Sweden, 28--2 p. /2 Nonparametric density estimation for elliptic problems with random perturbations
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 informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
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 informationNUMERICAL JUSTIFICATION OF ASYMPTOTIC EDDY CURRENTS MODEL FOR HETEROGENEOUS MATERIALS
Proceedings of ALGORITMY 2009 pp. 212 218 NUMERICAL JUSTIFICATION OF ASYMPTOTIC EDDY CURRENTS MODEL FOR HETEROGENEOUS MATERIALS V. MELICHER AND J.BUŠA JR. Abstract. We study electromagnetic properties
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 informationICES REPORT Parareal methods for highly oscillatory ordinary differential equations
ICES REPORT 15-6 March 215 Parareal methods for highly oscillatory ordinary differential equations by G. Ariel, S. Kim, and R. Tsai The Institute for Computational Engineering and Sciences The University
More informationRESEARCH PROPOSAL This research proposal is devoted to the development and analysis of mathematical and numerical methods, as well as computer
RESEARCH PROPOSAL This research proposal is devoted to the development and analysis of mathematical and numerical methods, as well as computer implementation in quantum relativistic physics. With the advent
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
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 informationCFD with OpenSource software
CFD with OpenSource software The Harmonic Balance Method in foam-extend Gregor Cvijetić Supervisor: Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture University of Zagreb gregor.cvijetic@fsb.hr
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationc 2004 Society for Industrial and Applied Mathematics
MULTISCALE MODEL. SIMUL. Vol. 3, No., pp. 65 88 c 4 Society for Industrial and Applied Mathematics ON WAVELET-BASED NUMERICAL HOMOGENIZATION ALINA CHERTOCK AND DORON LEVY Abstract. Recently, a wavelet-based
More informationParallel-in-Time Optimization with the General-Purpose XBraid Package
Parallel-in-Time Optimization with the General-Purpose XBraid Package Stefanie Günther, Jacob Schroder, Robert Falgout, Nicolas Gauger 7th Workshop on Parallel-in-Time methods Wednesday, May 3 rd, 2018
More informationSpace-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries
Space-Time Nonconforming Optimized Schwarz Waveform Relaxation for Heterogeneous Problems and General Geometries Laurence Halpern, Caroline Japhet, and Jérémie Szeftel 3 LAGA, Université Paris XIII, Villetaneuse
More informationCache Oblivious Stencil Computations
Cache Oblivious Stencil Computations S. HUNOLD J. L. TRÄFF F. VERSACI Lectures on High Performance Computing 13 April 2015 F. Versaci (TU Wien) Cache Oblivious Stencil Computations 13 April 2015 1 / 19
More informationOpleiding Wiskunde & Informatica
Opleiding Wiskunde & Informatica Numerics and continuation for Reaction-Diffusion equations Renzo Baasdam s1524054 Supervisors: Martina Chirilus-Bruckner (MI) & Michael Emmerich (LIACS) BACHELOR THESIS
More informationRandom walks and anisotropic interpolation on graphs. Filip Malmberg
Random walks and anisotropic interpolation on graphs. Filip Malmberg Interpolation of missing data Assume that we have a graph where we have defined some (real) values for a subset of the nodes, and that
More informationSpace-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure
Space-time domain decomposition for a nonlinear parabolic equation with discontinuous capillary pressure Elyes Ahmed, Caroline Japhet, Michel Kern INRIA Paris Maison de la Simulation ENPC University Paris
More informationJournal of Computational Physics
ournal of Computational Physics 9 () 397 3987 Contents lists available at ScienceDirect ournal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A multiscale coupling method for the
More informationMassively parallel electronic structure calculations with Python software. Jussi Enkovaara Software Engineering CSC the finnish IT center for science
Massively parallel electronic structure calculations with Python software Jussi Enkovaara Software Engineering CSC the finnish IT center for science GPAW Software package for electronic structure calculations
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 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 informationEfficient solution of stationary Euler flows with critical points and shocks
Efficient solution of stationary Euler flows with critical points and shocks Hans De Sterck Department of Applied Mathematics University of Waterloo 1. Introduction consider stationary solutions of hyperbolic
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 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 informationParallelizing large scale time domain electromagnetic inverse problem
Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,
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 informationA PARALLEL SPACE-TIME ALGORITHM
SIAM J. SCI. COMPUT. Vol. 34, No. 5, pp. C233 C248 c 2012 Society for Industrial and Applied Mathematics A PARALLEL SPACE-TIME ALGORITHM ANDREW J. CHRISTLIEB, RONALD D. HAYNES, AND BENJAMIN W. ONG Abstract.
More informationarxiv: v2 [math.na] 21 May 2018
SHORT-MT: Optimal Solution of Linear Ordinary Differential Equations by Conjugate Gradient Method arxiv:1805.01085v2 [math.na] 21 May 2018 Wenqiang Yang 1, Wenyuan Wu 1, and Robert M. Corless 2 1 Chongqing
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More information