Optimal Control of the Schrödinger Equation on Many-Core Architectures
|
|
- Jordan Quinn
- 5 years ago
- Views:
Transcription
1 Optimal Control of the Schrödinger Equation on Many-Core Architectures Manfred Liebmann Institute for Mathematics and Scientific Computing University of Graz July 15, 2014
2 Motivation We investigate efficient parallel algorithms for the optimal control problem of the timedepent Schrödinger equation on modern many-core architectures. We consider controls for multi-component Schrödinger systems with external laser fields: { min(1 φ, ψ(t ) 2 ) + γ 1 u 2 + γ 2 u 2 subject to i t ψ = ( 1 2m + V + u(t)b)ψ, ψ(0) = ψ 0 Where φ L 2 (R d, C N ) is the desired target state to be achieved at time T. The external laser field amplitude u : [0, T ] R is the control. The regularizing terms, which penalize overall field energy or overall field fluctuation, are measured in suitable norms and weighted by the regularization parameters γ 1, γ 2 > 0. The time-depent Schrödinger equation describes the effective quantum motion of nuclei with average mass 0 < m 1 and general initial conditions ψ 0 L 2 (R d, C N ). The dipole operator B initiates the transfer between the potential energy surfaces of the matrix-valued potential V : R d C N N. (1) Optimal Control of the Schrödinger Equation on Many-Core Architectures 1
3 Principal Investigators The international Research Training Group (IGDK 1754) Optimization and Numerical Analysis for Partial Differential Equations with Nonsmooth Structures is a joint project of Technische Universität München, Universität der Bundeswehr München, Karl-Franzens- Universität Graz, and Technische Universität Graz. P 12 Optimal control for Schrödinger dynamics on multiple potential energy surfaces Graduate Students Felix Henneke David Sattlegger Principal Investigators Gero Friesecke Caroline Lasser Karl Kunisch Optimal Control of the Schrödinger Equation on Many-Core Architectures 2
4 Introduction Numerical integrators for time-depent Hamiltonians in the Schrödinger equation i dψ dt = H(t)ψ, ψ(t 0) = ψ 0 (2) can be constructed as a combination of a generalized Suzuki-Trotter (GSTM) approximation scheme for the exponential operator and the Magnus expansion to capture the timedepence of the Hamiltonian. The general initial value problem ψ = Aψ for a time-indepent linear operator A with initial state ψ(0) = ψ 0 has the solution: ψ(t) = exp(ta)ψ 0 (3) Optimal Control of the Schrödinger Equation on Many-Core Architectures 3
5 Construction Theorem Theorem 1. Let B be a Banach algebra and A B an operator. Let Q 1 (t) B be an analytic operator function and an approximation of first order to the exponential operator exp(ta), exp(ta) Q 1 (t) = o(t) (t 0) (4) and let r > 1 be a fixed integer number. Then an analytic approximation of order m > 1 to the exponential operator is given by Q m (t) = r j=1 Q m 1 (p m,j t). (5) Where the coefficients p m,j for 1 j r are restricted by the two constraints r p m,j = 1, j=1 r p m m,j = 0. (6) j=1 Optimal Control of the Schrödinger Equation on Many-Core Architectures 4
6 Explicit Approximation Scheme We use the construction theorem to build a fractal approximation scheme for the operator function Q 1 (t) = 1 + ta and with r = 2. The coefficients p m,1, p m,2 are determined by the two restrictions The solution of this nonlinear system is given by p m,1 + p m,2 = 1, p m m,1 + p m m,2 = 0. (7) p m,1 = e iπ/m, p m,2 = 1. (8) 1 + eiπ/m This can also be written as p m = i 2 tan(π/2m), p m,1 = p m, p m,2 = p m (9) Optimal Control of the Schrödinger Equation on Many-Core Architectures 5
7 The complete approximation scheme is given as an operator product Q 2 (t) = Q 1 (p 2 t)q 1 ( p 2 t) (10) Q 3 (t) = Q 1 (p 3 p 2 t)q 1 (p 3 p 2 t)q 1 ( p 3 p 2 t)q 1 ( p 3 p 2 t) (11) Q m (t) = Q 1 (p m p 3 p 2 t)q 1 (p m p 3 p 2 t) Q 1 ( p m p 3 p 2 t) (12) and the complete formula for the approximation scheme reads: Q m (t) = z m,n = 2 m 1 n=1 m k=2 (1 + tz m,n A) (13) e iπb k,n/k 1 + e iπ/k (14) b k,n = (n 1)/2 k 2 mod 2 (15) Optimal Control of the Schrödinger Equation on Many-Core Architectures 6
8 Numerical Solution of the Optimal Control Problem The discretized optimal control problem reads min jh (u h ) = 1 u h U h 2 ψh N(u h ), O h ψn(u h h ) H h + α 2 uh 2 U h, (16) with the solution of the state equation calculated by the GSTM ψ h n(u h ) = n p=1 Q p (u h )ψ h 0, n = 1,..., N. (17) Numerical solution methods for the minimazation Newton Trust Region Barzilai-Borwein Optimal Control of the Schrödinger Equation on Many-Core Architectures 7
9 Data: u h 0 Result: optimal u h u h u h 0; ρ 1; r 1; j h functional(u h ); v h gradient(u h ); while v h TOL do δu h subproblem(u h, v h, r); j h new functional(u h + δu h ); w h hessian(u h, δu h ); ρ (j h j h new)/ δu h, v h 0.5w h U h; if ρ < η then r 0.25r; if ρ 0 then u h u h + δu h ; j h j h new; v h gradient(u h ); if ρ η + then r 2r; Algorithm 1: Newton trust region method Optimal Control of the Schrödinger Equation on Many-Core Architectures 8
10 Data: u h, v h, r Result: δu h p 0; res v h ; d res; while res > TOL 2 do w h hessian(u h, d); s d, w h U h; if s 0 then find τ such that p + τd = r; p p + τd; break; α res /s; if p + αd > r then find τ 0 such that p + τd = r; p p + τd; break; p p + αd; β res αw h 2 / res 2 ; res res αw h ; d res + βd; δu h p Algorithm 2: Newton trust region subproblem (Steihaug CG) Optimal Control of the Schrödinger Equation on Many-Core Architectures 9
11 Functional, Gradient and Hessian Evaluation Data: u h Result: j h ψ h ψ h 0 ; for n = 1,..., N do ψ h Q n (u h ) ψ h ; j h 1 2 ψh, O h ψ h H h + α 2 uh, u h U h Algorithm 3: Functional evaluation Optimal Control of the Schrödinger Equation on Many-Core Architectures 10
12 Data: u h Result: v h ψ h ψ h 0 ; X 0; for n = 1,..., N do ψ h Q n (u h )ψ h ; ϕ h O h ψ h ; for n = N..., 1 do assemble rhs gradient(u h, n, X) solve for Z: MZ = X; v h K k=1 Z kh k + αu h ; Algorithm 4: Gradient evaluation Optimal Control of the Schrödinger Equation on Many-Core Architectures 11
13 Data: u h, n, X Result: X ψ h Q n (u h )ψ h ; χ ( h ) ϕ h ; ( ϕ h Q n (u h ) 0 χ h ) ( ϕ h Q n (u h )(δu h ) Q n (u h ) χ h χ h ϕ h ; δx χ h, ψ h H h; X n X n + δx; X n 1 X n 1 + δx; Algorithm 5: assemble rhs gradient for DTO piecewise linear χ h ) ; Optimal Control of the Schrödinger Equation on Many-Core Architectures 12
14 Data: u h, δu Result: ( ) w h ( ) ψ h ψ h ψ h 0 ; 0 Y 0; for ( n = 1 )..., N do ψ h ( ϕ h ϕ h ψ h ) ( Qn (u h ) 0 Q n(u h )(δu h ) Q n (u h ) ( O h ψ h ) ; ) ( ψ h ψ h O h ψ h for n = N..., 1 do assemble rhs hessian(u h, δu h, n, Y ) solve for Z: MZ = Y ; w h K k=1 Z kh k + αδu h ; Algorithm 6: Hessian evaluation ) ; Optimal Control of the Schrödinger Equation on Many-Core Architectures 13
15 Data: u h, δu h, n, Y ( Result: ) Y ( ) ( ψ h Qn (u ψ h h ) 0 ψ h Q n (u h )(δu h ) Q n (u h ) ψ h χ h 2 ϕ h ; χ ( h 3 ) ϕ h ; ( ) ( ϕ h Q χ h n (u h ) 0 ϕ h 2 Q n (u h )(h n ) Q n (u h ) χ h 2 χ h 2 χ h 2 ϕ h ; ϕ h χ h 3; χ h 1 ϕ h ; ϕ h ϕ h χ h 1 χ h 3 Q n (u h ) Q n (u h )(δu h ) Q n (u h ) 0 0 Q n (u h )(h n ) 0 Q n (u h ) 0 Q n (u h )(δu h, h n ) Q n (u h )(h n ) Q n (u h )(δu h ) Q n (u h ) χ h 1 χ h 1 ϕ h ; χ h 3 χ h 3 ϕ h χ h 2; δy χ h 1, ψ h H h + χ h 2, ψ h H h + χ h 3, ψ h H h; Y n Y n + δy ; Y n 1 Y n 1 + δy ; Algorithm 7: assemble rhs hessian for DTO piecewise linear ) ; ) ; Optimal Control of the Schrödinger Equation on Many-Core Architectures 14 ϕ h ϕ h χ h 1 χ h 3 ;
16 Time Optimal Control 1D Model Problem Figure 1: Two potential energy surfaces V 11, V 22, V 12 = V 21 = 0 on the left. Components of the magnetic dipole field B 11 = B 22, B 12 = B 21 on the right. Spatial grid size is 256. Optimal Control of the Schrödinger Equation on Many-Core Architectures 15
17 Newton Trust Region Algorithm Figure 2: Initial control function for the laser field u on the left. Control function calculated by the Newton trust region algorithm after 30 iterations on the right. The functional was reduced from to Temporal grid size is Compute time seconds on Intel Core2Duo 2.8GHz. Optimal Control of the Schrödinger Equation on Many-Core Architectures 16
18 Manfred Liebmann July 15, 2014 Figure 3: Visualization of the complex wavefunctions ψ1, ψ2 and the corresponding probability distribution in a space time diagram. Optimal Control of the Schro dinger Equation on Many-Core Architectures 17
19 Parallelization of the Optimal Control Problem One-dimensional problems are too small for effective parallelization although the compute time is already quite high. Two-dimensional test problem on a spacial grid and 2048 point temporal grid with a similar simulation setup to the 1D problem. Cores GNU g++ Intel icpc compact KMP AFFINITY 4KB page mm malloc More parallel code Table 1: Parallel scaling with OpenMP parallelization on Intel Xeon 2.00GHz for 10 iterations of the Barzilai-Borwein method. Optimal Control of the Schrödinger Equation on Many-Core Architectures 18
20 Many-Core Parallelization for the Optimal Control Problem 2x Intel Xeon E Intel Xeon Phi 5110P Nvidia Tesla K20 16 Cores 60 Cores / 240 Threads 2496 Cuda Cores Table 2: Benchmark for Newton trust region method for 10 iterations on a space-time grid. Timings are given in seconds. Optimal Control of the Schrödinger Equation on Many-Core Architectures 19
21 Conclusions Optimal Control of the Schrödinger Equation Optimal control problems are well suited for many-core architectures Algorithms are essentially bound by the performance of the PDE solver. Flexible C++ framework for different state equations and solution methods Future work: Complete framework for MPI/OpenMP/SSE/AVX/GPU/PHI Optimal Control of the Schrödinger Equation on Many-Core Architectures 20
22 Thank you! Optimal Control of the Schrödinger Equation on Many-Core Architectures 21
A Massively Parallel Eigenvalue Solver for Small Matrices on Multicore and Manycore Architectures
A Massively Parallel Eigenvalue Solver for Small Matrices on Multicore and Manycore Architectures Manfred Liebmann Technische Universität München Chair of Optimal Control Center for Mathematical Sciences,
More informationParallel Algorithms for Semiclassical Quantum Dynamics
Parallel Algorithms for Semiclassical Quantum Dynamics Manfred Liebmann David Sattlegger May, 015 Institute for Mathematics and Scientific Computing, University of Graz, Heinrichstraße 36, 8010 Graz, Austria
More informationSimulation based optimization
SimBOpt p.1/52 Simulation based optimization Feb 2005 Eldad Haber haber@mathcs.emory.edu Emory University SimBOpt p.2/52 Outline Introduction A few words about discretization The unconstrained framework
More informationAlgorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method
Algorithm for Sparse Approximate Inverse Preconditioners in the Conjugate Gradient Method Ilya B. Labutin A.A. Trofimuk Institute of Petroleum Geology and Geophysics SB RAS, 3, acad. Koptyug Ave., Novosibirsk
More informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationOPER 627: Nonlinear Optimization Lecture 9: Trust-region methods
OPER 627: Nonlinear Optimization Lecture 9: Trust-region methods Department of Statistical Sciences and Operations Research Virginia Commonwealth University Sept 25, 2013 (Lecture 9) Nonlinear Optimization
More informationA robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme
A robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme International Conference On Preconditioning Techniques For Scientific And Industrial Applications 17-19
More informationTraining Schro dinger s Cat
Training Schro dinger s Cat Quadratically Converging algorithms for Optimal Control of Quantum Systems David L. Goodwin & Ilya Kuprov d.goodwin@soton.ac.uk comp-chem@southampton, Wednesday 15th July 2015
More informationDesign of optimal RF pulses for NMR as a discrete-valued control problem
Design of optimal RF pulses for NMR as a discrete-valued control problem Christian Clason Faculty of Mathematics, Universität Duisburg-Essen joint work with Carla Tameling (Göttingen) and Benedikt Wirth
More informationAPPLICATION OF CUDA TECHNOLOGY FOR CALCULATION OF GROUND STATES OF FEW-BODY NUCLEI BY FEYNMAN'S CONTINUAL INTEGRALS METHOD
APPLICATION OF CUDA TECHNOLOGY FOR CALCULATION OF GROUND STATES OF FEW-BODY NUCLEI BY FEYNMAN'S CONTINUAL INTEGRALS METHOD M.A. Naumenko, V.V. Samarin Joint Institute for Nuclear Research, Dubna, Russia
More informationNumerical Modeling of Methane Hydrate Evolution
Numerical Modeling of Methane Hydrate Evolution Nathan L. Gibson Joint work with F. P. Medina, M. Peszynska, R. E. Showalter Department of Mathematics SIAM Annual Meeting 2013 Friday, July 12 This work
More informationIntroduction to the GRAPE Algorithm
June 8, 2010 Reference: J. Mag. Res. 172, 296 (2005) Journal of Magnetic Resonance 172 (2005) 296 305 www.elsevier.com/locate/jmr Optimal control of coupled spin dynamics: design of NMR pulse sequences
More informationSparse Control of Quantum Systems
SpezialForschungsBereich F 3 Karl Franzens Universit at Graz Technische Universit at Graz Medizinische Universit at Graz Sparse Control of Quantum Systems G. Friesecke F. Henneke K. Kunisch SFB-Report
More informationIntroduction to numerical computations on the GPU
Introduction to numerical computations on the GPU Lucian Covaci http://lucian.covaci.org/cuda.pdf Tuesday 1 November 11 1 2 Outline: NVIDIA Tesla and Geforce video cards: architecture CUDA - C: programming
More informationLecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms)
Lecture 1: Numerical Issues from Inverse Problems (Parameter Estimation, Regularization Theory, and Parallel Algorithms) Youzuo Lin 1 Joint work with: Rosemary A. Renaut 2 Brendt Wohlberg 1 Hongbin Guo
More informationSemismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations
Semismooth Newton Methods for an Optimal Boundary Control Problem of Wave Equations Axel Kröner 1 Karl Kunisch 2 and Boris Vexler 3 1 Lehrstuhl für Mathematische Optimierung Technische Universität München
More informationHigher-Order Methods
Higher-Order Methods Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. PCMI, July 2016 Stephen Wright (UW-Madison) Higher-Order Methods PCMI, July 2016 1 / 25 Smooth
More informationatoms and light. Chapter Goal: To understand the structure and properties of atoms.
Quantum mechanics provides us with an understanding of atomic structure and atomic properties. Lasers are one of the most important applications of the quantummechanical properties of atoms and light.
More informationequations in the semiclassical regime
Adaptive time splitting for nonlinear Schrödinger equations in the semiclassical regime Local error representation and a posteriori error estimator Workshop, Numerical Analysis of Evolution Equations Vill,
More informationQuantum mechanics is a physical science dealing with the behavior of matter and energy on the scale of atoms and subatomic particles / waves.
Quantum mechanics is a physical science dealing with the behavior of matter and energy on the scale of atoms and subatomic particles / waves. It also forms the basis for the contemporary understanding
More informationA Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations
A Newton-Galerkin-ADI Method for Large-Scale Algebraic Riccati Equations Peter Benner Max-Planck-Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory
More informationMore NMR Relaxation. Longitudinal Relaxation. Transverse Relaxation
More NMR Relaxation Longitudinal Relaxation Transverse Relaxation Copyright Peter F. Flynn 2017 Experimental Determination of T1 Gated Inversion Recovery Experiment The gated inversion recovery pulse sequence
More information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
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 informationBound and Scattering Solutions for a Delta Potential
Physics 342 Lecture 11 Bound and Scattering Solutions for a Delta Potential Lecture 11 Physics 342 Quantum Mechanics I Wednesday, February 20th, 2008 We understand that free particle solutions are meant
More informationLarge-scale Electronic Structure Simulations with MVAPICH2 on Intel Knights Landing Manycore Processors
Large-scale Electronic Structure Simulations with MVAPICH2 on Intel Knights Landing Manycore Processors Hoon Ryu, Ph.D. (E: elec1020@kisti.re.kr) Principal Researcher / Korea Institute of Science and Technology
More informationExact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation
Exact Solutions to the Focusing Discrete Nonlinear Schrödinger Equation Francesco Demontis (based on a joint work with C. van der Mee) Università degli Studi di Cagliari Dipartimento di Matematica e Informatica
More informationPhysics 342 Lecture 30. Solids. Lecture 30. Physics 342 Quantum Mechanics I
Physics 342 Lecture 30 Solids Lecture 30 Physics 342 Quantum Mechanics I Friday, April 18th, 2008 We can consider simple models of solids these highlight some special techniques. 30.1 An Electron in a
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationIntroduction. New Nonsmooth Trust Region Method for Unconstraint Locally Lipschitz Optimization Problems
New Nonsmooth Trust Region Method for Unconstraint Locally Lipschitz Optimization Problems Z. Akbari 1, R. Yousefpour 2, M. R. Peyghami 3 1 Department of Mathematics, K.N. Toosi University of Technology,
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationFinite Elements for Magnetohydrodynamics and its Optimal Control
Finite Elements for Magnetohydrodynamics and its Karl Kunisch Marco Discacciati (RICAM) FEM Symposium Chemnitz September 25 27, 2006 Overview 1 2 3 What is Magnetohydrodynamics? Magnetohydrodynamics (MHD)
More informationControl of bilinear systems: multiple systems and perturbations
Control of bilinear systems: multiple systems and perturbations Gabriel Turinici CEREMADE, Université Paris Dauphine ESF OPTPDE Workshop InterDyn2013, Modeling and Control of Large Interacting Dynamical
More informationThe Particle in a Box
Page 324 Lecture 17: Relation of Particle in a Box Eigenstates to Position and Momentum Eigenstates General Considerations on Bound States and Quantization Continuity Equation for Probability Date Given:
More informationScaled gradient projection methods in image deblurring and denoising
Scaled gradient projection methods in image deblurring and denoising Mario Bertero 1 Patrizia Boccacci 1 Silvia Bonettini 2 Riccardo Zanella 3 Luca Zanni 3 1 Dipartmento di Matematica, Università di Genova
More informationThe Quantum Theory of Atoms and Molecules
The Quantum Theory of Atoms and Molecules The postulates of quantum mechanics Dr Grant Ritchie The postulates.. 1. Associated with any particle moving in a conservative field of force is a wave function,
More informationInterpolation with Radial Basis Functions on GPGPUs using CUDA
Interpolation with Radial Basis Functions on GPGPUs using CUDA Gundolf Haase in coop. with: Dirk Martin [VRV Vienna] and Günter Offner [AVL Graz] Institute for Mathematics and Scientific Computing University
More informationMotion Estimation (I) Ce Liu Microsoft Research New England
Motion Estimation (I) Ce Liu celiu@microsoft.com Microsoft Research New England We live in a moving world Perceiving, understanding and predicting motion is an important part of our daily lives Motion
More informationLecture 16: Relaxation methods
Lecture 16: Relaxation methods Clever technique which begins with a first guess of the trajectory across the entire interval Break the interval into M small steps: x 1 =0, x 2,..x M =L Form a grid of points,
More informationTheoretical and numerical aspects of quantum control problems
Theoretical and numerical aspects of quantum control problems Dipartimento e Facoltà di Ingegneria, Università degli Studi del Sannio, Italy Institute for Mathematics and Scientific Computing Karl-Franzens-University,
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationPage 684. Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02
Page 684 Lecture 40: Coordinate Transformations: Time Transformations Date Revised: 2009/02/02 Date Given: 2009/02/02 Time Transformations Section 12.5 Symmetries: Time Transformations Page 685 Time Translation
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationHPMPC - A new software package with efficient solvers for Model Predictive Control
- A new software package with efficient solvers for Model Predictive Control Technical University of Denmark CITIES Second General Consortium Meeting, DTU, Lyngby Campus, 26-27 May 2015 Introduction Model
More informationPerm State University Research-Education Center Parallel and Distributed Computing
Perm State University Research-Education Center Parallel and Distributed Computing A 25-minute Talk (S4493) at the GPU Technology Conference (GTC) 2014 MARCH 24-27, 2014 SAN JOSE, CA GPU-accelerated modeling
More informationB O S Z. - Boussinesq Ocean & Surf Zone model - International Research Institute of Disaster Science (IRIDeS), Tohoku University, JAPAN
B O S Z - Boussinesq Ocean & Surf Zone model - Volker Roeber 1 Troy W. Heitmann 2, Kwok Fai Cheung 2, Gabriel C. David 3, Jeremy D. Bricker 1 1 International Research Institute of Disaster Science (IRIDeS),
More informationA classification of gapped Hamiltonians in d = 1
A classification of gapped Hamiltonians in d = 1 Sven Bachmann Mathematisches Institut Ludwig-Maximilians-Universität München Joint work with Yoshiko Ogata NSF-CBMS school on quantum spin systems Sven
More informationReactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)
Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li
More informationLogo. A Massively-Parallel Multicore Acceleration of a Point Contact Solid Mechanics Simulation DRAFT
Paper 1 Logo Civil-Comp Press, 2017 Proceedings of the Fifth International Conference on Parallel, Distributed, Grid and Cloud Computing for Engineering, P. Iványi, B.H.V Topping and G. Várady (Editors)
More informationINTRODUCTION TO FINITE ELEMENT METHODS ON ELLIPTIC EQUATIONS LONG CHEN
INTROUCTION TO FINITE ELEMENT METHOS ON ELLIPTIC EQUATIONS LONG CHEN CONTENTS 1. Poisson Equation 1 2. Outline of Topics 3 2.1. Finite ifference Method 3 2.2. Finite Element Method 3 2.3. Finite Volume
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationNumerical Methods for Large-Scale Nonlinear Systems
Numerical Methods for Large-Scale Nonlinear Systems Handouts by Ronald H.W. Hoppe following the monograph P. Deuflhard Newton Methods for Nonlinear Problems Springer, Berlin-Heidelberg-New York, 2004 Num.
More informationMODELING MATTER AT NANOSCALES. 4. Introduction to quantum treatments Outline of the principles and the method of quantum mechanics
MODELING MATTER AT NANOSCALES 4. Introduction to quantum treatments 4.01. Outline of the principles and the method of quantum mechanics 1 Why quantum mechanics? Physics and sizes in universe Knowledge
More informationRegularization of An Optimal Control Problem With BV-Functions
Regularization of An Optimal Control Problem With BV-Functions Dominik Hafemeyer 1, Florian Kruse 2 1 Technische Universität Müchen, Chair for Optimal Control, 2 Universität Graz, Institute of Mathematics
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationPerturbation Theory and Numerical Modeling of Quantum Logic Operations with a Large Number of Qubits
Contemporary Mathematics Perturbation Theory and Numerical Modeling of Quantum Logic Operations with a Large Number of Qubits G. P. Berman, G. D. Doolen, D. I. Kamenev, G. V. López, and V. I. Tsifrinovich
More informationCombined systems in PT-symmetric quantum mechanics
Combined systems in PT-symmetric quantum mechanics Brunel University London 15th International Workshop on May 18-23, 2015, University of Palermo, Italy - 1 - Combined systems in PT-symmetric quantum
More informationComputational Numerical Integration for Spherical Quadratures. Verified by the Boltzmann Equation
Computational Numerical Integration for Spherical Quadratures Verified by the Boltzmann Equation Huston Rogers. 1 Glenn Brook, Mentor 2 Greg Peterson, Mentor 2 1 The University of Alabama 2 Joint Institute
More informationTowards a highly-parallel PDE-Solver using Adaptive Sparse Grids on Compute Clusters
Towards a highly-parallel PDE-Solver using Adaptive Sparse Grids on Compute Clusters HIM - Workshop on Sparse Grids and Applications Alexander Heinecke Chair of Scientific Computing May 18 th 2011 HIM
More informationAspects of Two- and Three-Flavor Chiral Phase Transitions
Aspects of Two- and Three-Flavor Chiral Phase Transitions Mario Karl-Franzens-Universität Graz Institut für Physik Fachbereich Theoretische Physik Kyoto, September 6, 211 Table of Contents 1 Motivation
More informationTrust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization
Trust-Region SQP Methods with Inexact Linear System Solves for Large-Scale Optimization Denis Ridzal Department of Computational and Applied Mathematics Rice University, Houston, Texas dridzal@caam.rice.edu
More informationAffine covariant Semi-smooth Newton in function space
Affine covariant Semi-smooth Newton in function space Anton Schiela March 14, 2018 These are lecture notes of my talks given for the Winter School Modern Methods in Nonsmooth Optimization that was held
More informationA minimum effort optimal control problem for the wave equation
www.oeaw.ac.at A minimum effort optimal control problem for the wave equation A. Kröner, K. Kunisch RICAM-Report 2013-02 www.ricam.oeaw.ac.at Noname manuscript No. (will be inserted by the editor) A minimum
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationNumerical simulation of the Gross-Pitaevskii equation by pseudo-spectral and finite element methods comparison of GPS code and FreeFem++
. p.1/13 10 ec. 2014, LJLL, Paris FreeFem++ workshop : BECASIM session Numerical simulation of the Gross-Pitaevskii equation by pseudo-spectral and finite element methods comparison of GPS code and FreeFem++
More informationOn Lagrange multipliers of trust region subproblems
On Lagrange multipliers of trust region subproblems Ladislav Lukšan, Ctirad Matonoha, Jan Vlček Institute of Computer Science AS CR, Prague Applied Linear Algebra April 28-30, 2008 Novi Sad, Serbia Outline
More informationIntroduction into Implementation of Optimization problems with PDEs: Sheet 3
Technische Universität München Center for Mathematical Sciences, M17 Lucas Bonifacius, Korbinian Singhammer www-m17.ma.tum.de/lehrstuhl/lehresose16labcourseoptpde Summer 216 Introduction into Implementation
More informationMathematical Tripos Part IB Michaelmas Term Example Sheet 1. Values of some physical constants are given on the supplementary sheet
Mathematical Tripos Part IB Michaelmas Term 2015 Quantum Mechanics Dr. J.M. Evans Example Sheet 1 Values of some physical constants are given on the supplementary sheet 1. Whenasampleofpotassiumisilluminatedwithlightofwavelength3
More informationPerformance of the fusion code GYRO on three four generations of Crays. Mark Fahey University of Tennessee, Knoxville
Performance of the fusion code GYRO on three four generations of Crays Mark Fahey mfahey@utk.edu University of Tennessee, Knoxville Contents Introduction GYRO Overview Benchmark Problem Test Platforms
More informationNonlinear Model Reduction for Rubber Components in Vehicle Engineering
Nonlinear Model Reduction for Rubber Components in Vehicle Engineering Dr. Sabrina Herkt, Dr. Klaus Dreßler Fraunhofer Institut für Techno- und Wirtschaftsmathematik Kaiserslautern Prof. Rene Pinnau Universität
More informationChapter 2 The Group U(1) and its Representations
Chapter 2 The Group U(1) and its Representations The simplest example of a Lie group is the group of rotations of the plane, with elements parametrized by a single number, the angle of rotation θ. It is
More informationMEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS
MEASURE VALUED DIRECTIONAL SPARSITY FOR PARABOLIC OPTIMAL CONTROL PROBLEMS KARL KUNISCH, KONSTANTIN PIEPER, AND BORIS VEXLER Abstract. A directional sparsity framework allowing for measure valued controls
More informationComputation of Local ISS Lyapunov Function Via Linear Programming
Computation of Local ISS Lyapunov Function Via Linear Programming Huijuan Li Joint work with Robert Baier, Lars Grüne, Sigurdur F. Hafstein and Fabian Wirth Institut für Mathematik, Universität Würzburg
More informationNumerical Integration of the Wavefunction in FEM
Numerical Integration of the Wavefunction in FEM Ian Korey Eaves ike26@drexel.edu December 14, 2013 Abstract Although numerous techniques for calculating stationary states of the schrödinger equation are
More informationSPACE MAPPING FOR OPTIMAL CONTROL OF PARTIAL DIFFERENTIAL EQUATIONS
SPACE MAPPING FOR OPTIMAL CONTROL OF PARTIAL DIFFERENTIAL EQUATIONS MICHAEL HINTERMÜLLER AND LUíS N. VICENTE Abstract. Solving optimal control problems for nonlinear partial differential equations represents
More informationInverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 2007 Technische Universiteit Eindh ove n University of Technology
Inverse problems Total Variation Regularization Mark van Kraaij Casa seminar 23 May 27 Introduction Fredholm first kind integral equation of convolution type in one space dimension: g(x) = 1 k(x x )f(x
More informationChebychev Propagator for Inhomogeneous Schrödinger Equations
Chebychev Propagator for Inhomogeneous Schrödinger Equations Michael Goerz May 16, 2011 Solving the Schrödinger Equation Schrödinger Equation Ψ(t) = Ĥ Ψ(t) ; e.g. Ĥ = t ( V1 (R) ) µɛ(t) µɛ(t) V 2 (R) Solving
More informationRobust Principal Component Pursuit via Alternating Minimization Scheme on Matrix Manifolds
Robust Principal Component Pursuit via Alternating Minimization Scheme on Matrix Manifolds Tao Wu Institute for Mathematics and Scientific Computing Karl-Franzens-University of Graz joint work with Prof.
More informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
More informationTight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction
16/05/2017, CTU in Prague Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej
More informationDense Arithmetic over Finite Fields with CUMODP
Dense Arithmetic over Finite Fields with CUMODP Sardar Anisul Haque 1 Xin Li 2 Farnam Mansouri 1 Marc Moreno Maza 1 Wei Pan 3 Ning Xie 1 1 University of Western Ontario, Canada 2 Universidad Carlos III,
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationInfinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems
Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems D. Georges, Control Systems Dept - Gipsa-lab, Grenoble INP Workshop on Irrigation Channels and Related Problems,
More informationarxiv: v1 [physics.comp-ph] 10 Aug 2017
A QUANTUM KINETIC MONTE CARLO METHOD FOR QUANTUM MANY-BODY SPIN DYNAMICS arxiv:79.386v [physics.comp-ph] Aug 27 ZHENNING CAI AND JIANFENG LU Abstract. We propose a general framework of quantum kinetic
More informationParallel Transposition of Sparse Data Structures
Parallel Transposition of Sparse Data Structures Hao Wang, Weifeng Liu, Kaixi Hou, Wu-chun Feng Department of Computer Science, Virginia Tech Niels Bohr Institute, University of Copenhagen Scientific Computing
More informationMultiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU
Multiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU Khramtsov D.P., Nekrasov D.A., Pokusaev B.G. Department of Thermodynamics, Thermal Engineering and Energy Saving Technologies,
More informationAlgorithms for PDE-Constrained Optimization
GAMM-Mitteilungen, 31 January 2014 Algorithms for PDE-Constrained Optimization Roland Herzog 1 and Karl Kunisch 2 1 Chemnitz University of Technology, Faculty of Mathematics, Reichenhainer Straße 41, D
More informationTR A Comparison of the Performance of SaP::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems
TR-0-07 A Comparison of the Performance of ::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems Ang Li, Omkar Deshmukh, Radu Serban, Dan Negrut May, 0 Abstract ::GPU is a
More information(A) Opening Problem Newton s Law of Cooling
Lesson 55 Numerical Solutions to Differential Equations Euler's Method IBHL - 1 (A) Opening Problem Newton s Law of Cooling! Newton s Law of Cooling states that the temperature of a body changes at a rate
More informationExplicit approximate controllability of the Schrödinger equation with a polarizability term.
Explicit approximate controllability of the Schrödinger equation with a polarizability term. Morgan MORANCEY CMLA, ENS Cachan Sept. 2011 Control of dispersive equations, Benasque. Morgan MORANCEY (CMLA,
More informationA Randomized Algorithm for the Approximation of Matrices
A Randomized Algorithm for the Approximation of Matrices Per-Gunnar Martinsson, Vladimir Rokhlin, and Mark Tygert Technical Report YALEU/DCS/TR-36 June 29, 2006 Abstract Given an m n matrix A and a positive
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
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 informationFrom the N-body problem to the cubic NLS equation
From the N-body problem to the cubic NLS equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Los Alamos CNLS, January 26th, 2005 Formal derivation by N.N. Bogolyubov
More informationOn the computation of the reciprocal of floating point expansions using an adapted Newton-Raphson iteration
On the computation of the reciprocal of floating point expansions using an adapted Newton-Raphson iteration Mioara Joldes, Valentina Popescu, Jean-Michel Muller ASAP 2014 1 / 10 Motivation Numerical problems
More informationExp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 +
Exp. 4. Quantum Chemical calculation: The potential energy curves and the orbitals of H2 + 1. Objectives Quantum chemical solvers are used to obtain the energy and the orbitals of the simplest molecules
More information2 Nonlinear least squares algorithms
1 Introduction Notes for 2017-05-01 We briefly discussed nonlinear least squares problems in a previous lecture, when we described the historical path leading to trust region methods starting from the
More information