ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers
|
|
|
- Kevin Norris
- 5 years ago
- Views:
Transcription
1 ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers Victor Yu and the ELSI team Department of Mechanical Engineering & Materials Science Duke University
2 Kohn-Sham Density-Functional Theory density DGDFT NWCHEM Hc = εsc L.M. Ghiringhelli et al., Towards a Common Format for Computational Materials Science Data, Psi-k Scientific Highlight of the Month 131, July wave function 1/13
3 Kohn-Sham Density-Functional Theory density DGDFT NWCHEM density matrix Hc = εsc L.M. Ghiringhelli et al., Towards a Common Format for Computational Materials Science Data, Psi-k Scientific Highlight of the Month 131, July wave function 1/13
4 Kohn-Sham Density-Functional Theory density DGDFT NWCHEM density matrix Hc = εsc L.M. Ghiringhelli et al., Towards a Common Format for Computational Materials Science Data, Psi-k Scientific Highlight of the Month 131, July /13
5 Kohn-Sham Density-Functional Theory density DGDFT NWCHEM density matrix Hc = εsc L.M. Ghiringhelli et al., Towards a Common Format for Computational Materials Science Data, Psi-k Scientific Highlight of the Month 131, July wave function 1/13
6 The Cubic Wall in Kohn-Sham Density-Functional Theory DFT-PBE with FHI-aims Graphene monolayer (2D) From 4,050 to 7,200 atoms 14 numeric atom-centered orbitals per carbon atom The O(N 3 ) Kohn-Sham eigenproblem must be addressed by all electronic structure codes based on KS-DFT O(N) O(N 3 ) Edison (Intel Ivy Bridge) 80 nodes x 24 CPU cores per node 2/13
7 Kohn-Sham Electronic Structure Solvers ELPA (eigensolution) Dense, direct eigensolver O(N 3 ) libomm (density matrix) Orbital minimization O(N 3 ) reduced prefactor PEXSI (density matrix) Pole expansion, selected inversion O(N 2 ) SIPs (eigensolution) Sparse, iterative eigensolver O(N 2 ) CheSS (density matrix) Chebyshev expansion O(N) And many others, e.g. iterative solvers and ScaLAPACK A unified software interface and solver comparison platform 3/13
8 ELSI Timeline ELSI Kick-Off Meeting, Duke University, Durham, 2017 Aug 2015
9 ELSI: Connection between DFT Codes and KS Solvers Create Hamiltonian and overlap matrices Different matrix storage strategies Different application programming interfaces Different programming languages Compute eigenvalues, eigenvectors, density matrix 4/13
10 ELSI: Connection between DFT Codes and KS Solvers Create Hamiltonian and overlap matrices Convert data format and distribution automatically Analyze the problem: Basis set size, sparsity, desired output, available computational resource, (ongoing) Suggest a solver with reasonable settings (future) Compute eigenvalues, eigenvectors, density matrix 4/13
11 ELSI Timeline Electronic Structure Library Coding Workshop, ZCAM, Zaragoza, 2016 Aug 2015 Jun 2016
12 Philosophy of the ELSI Interface Easy to implement Designed to be flexible for rapid integration into a variety of KS-DFT codes Versatile Includes density matrix and eigensystem formalisms on equal footing Flexible All technical settings in solver libraries are adjustable for experienced users API of ELSI interface elsi_init! BEGIN SCF LOOP elsi_set elsi_{dm ev}! END SCF LOOP elsi_finalize 5/13
13 ELSI Timeline First Release of ELSI Aug 2015 Jun 2016 May 2017
14 First Release of ELSI: May Software Package - ELSI interface - Current ELPA, libomm, PEXSI Programming Language - Fortran C/C++ interface Functionalities - Eigenvalues and eigenvectors - Density matrix - Chemical potential - Occupation numbers Portability - GNU, Intel, PGI, IBM, Cray - Cori, Edison, Mira, Theta, Titan Flexibility - 2 matrix I/O formats - 2 parallel modes - Concurrent instances - Tested in DGDFT, FHI-aims, NWChem (Global Arrays), SIESTA Scalability - (Hundreds of) thousands of MPI tasks 6/13
15 Supported Matrix Formats a) 2D block-cyclic distributed dense storage (BLACS standard) b) 1D block distributed compressed column storage (PEXSI standard) c) 1D block distributed compressed row storage supported implicitly due to symmetry Auto-conversion between (a) and (b) when necessary Example: BLACS input è PEXSI format è PEXSI solver è BLACS output 7/13
16 Performance of the Solvers (ELSI May 2017) DFT-PBE with FHI-aims Graphene monolayer (2D) From 1,800 to 7,200 atoms V. Yu et al., ELSI: A Unified Software Interface for Kohn-Sham Electronic Structure Solvers, arxiv: [physics.comp-ph]. Edison (Intel Ivy Bridge) 80 nodes x 24 CPU cores per node 14 numeric atom-centered orbitals per carbon atom PEXSI starts outperforming ELPA and libomm at 3,000 and 7,000 atoms, respectively 8/13
17 Performance of the Solvers (Today) DFT-PBE with FHI-aims Graphene monolayer (2D) From 800 to 7,200 atoms Edison (Intel Ivy Bridge) 80 nodes x 24 CPU cores per node 14 numeric atom-centered orbitals per carbon atom PEXSI starts outperforming ELPA and libomm at 3, and 7,000 3,200 atoms, respectively! Up-to-date 9/13
18 Enhanced Feature: Energy-Weighted Density Matrix ELPA libomm PEXSI DM * C () C () W(2I / S / )W * Im 4 w 6 () H z 6 + μ S EDM * C =() C =() W[ 2I / S / H / ]W * Im 4 w 6 =() : 6;< : 6;< H z 6 + μ S Pulay force (localized orbitals) Available with all three solvers Only computed when calling elsi_get_edm Geometry optimization 10/13
19 ELSI Timeline Electronic Structure Library Coding Workshop, ICTP, Trieste, 2017 Aug 2015 Jun 2016 May 2017 Jul 2017
20 Enhanced Feature: K-Point and Spin Calculations elsi_init(handle,...) New MPI communicator for communications among all k-points and spins! Set MPI communicator elsi_set_mpi(handle,my_mpi_comm) elsi_set_mpi_all(handle,mpi_comm_all)! Set spin and/or k-point elsi_set_spin(handle,n_spin,my_spin) elsi_set_kpt(handle,n_kpt,my_kpt,my_weight)! Compute density matrices elsi_dm_real(handle,my_h,my_s,my_dm,energy) elsi_finalize(handle) 11/13
21 Enhanced Feature: K-Point and Spin Calculations elsi_init(handle,...)! Set MPI communicator elsi_set_mpi(handle,my_mpi_comm) elsi_set_mpi_all(handle,mpi_comm_all) New Number of spins and my spin index! Set spin and/or k-point elsi_set_spin(handle,n_spin,my_spin) elsi_set_kpt(handle,n_kpt,my_kpt,my_weight)! Compute density matrices elsi_dm_real(handle,my_h,my_s,my_dm,energy) elsi_finalize(handle) 11/13
22 Enhanced Feature: K-Point and Spin Calculations elsi_init(handle,...)! Set MPI communicator elsi_set_mpi(handle,my_mpi_comm) elsi_set_mpi_all(handle,mpi_comm_all) New Number of k-points, my k-point index, and my k-point weight! Set spin and/or k-point elsi_set_spin(handle,n_spin,my_spin) elsi_set_kpt(handle,n_kpt,my_kpt,my_weight)! Compute density matrices elsi_dm_real(handle,my_h,my_s,my_dm,energy) elsi_finalize(handle) 11/13
23 Enhanced Feature: K-Point and Spin Calculations Minimal additions to the existing ELSI interface Capability to handle k-point and/or spin in parallel No overhead if not interested in k-point and/or spin elsi_init(handle,...)! Set MPI communicator elsi_set_mpi(handle,my_mpi_comm) elsi_set_mpi_all(handle,mpi_comm_all)! Set spin and/or k-point elsi_set_spin(handle,n_spin,my_spin) elsi_set_kpt(handle,n_kpt,my_kpt,my_weight) Spin K-point! Compute density matrices elsi_dm_real(handle,my_h,my_s,my_dm,energy) elsi_finalize(handle) 11/13
24 New Feature: CheSS and SIPs Solvers CheSS (density matrix) Fermi operator expansion with Chebyshev polynomials O(N) Both solvers connected to ELSI - SOLVER = 4 (CheSS) - SOLVER = 5 (SIPs) Flexibility of the ELSI infrastructure: Given compatible matrix formats, a new solver can be quickly added to ELSI in hours SIPs (eigensolution) Sparse, iterative eigensolver based on spectrum slicing O(N 2 ) Both are giving correct results Todo - Performance optimization - Benchmark O(N) solver Spectrum slicing 12/13
25 ELSI Timeline Aug 2015 Jun 2016 May 2017 Jul 2017 Future
26 Future Plans & Acknowledgments Performance-optimized, up-to-date solver libraries Iterative solvers, e.g. Davidson and CG solvers from Quantum Espresso ELSI decision layer for automatic, optimal choice of solvers More matrix formats used in KS- DFT codes to be supported Integration of ELSI into KS-DFT codes for production-level usage ELSI workshop in 2018 ELSI is an NSF SI2-SSI supported software infrastructure project under Grant Number Any opinions, findings, and conclusions or recommendations expressed here are those of the author(s) and do not necessarily reflect the views of NSF. 13/13
Recent Developments in the ELSI Infrastructure for Large-Scale Electronic Structure Theory
elsi-interchange.org MolSSI Workshop and ELSI Conference 2018, August 15, 2018, Richmond, VA Recent Developments in the ELSI Infrastructure for Large-Scale Electronic Structure Theory Victor Yu 1, William
DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations
DGDFT: A Massively Parallel Method for Large Scale Density Functional Theory Calculations The recently developed discontinuous Galerkin density functional theory (DGDFT)[21] aims at reducing the number
A knowledge-based approach to high-performance computing in ab initio simulations.
Mitglied der Helmholtz-Gemeinschaft A knowledge-based approach to high-performance computing in ab initio simulations. AICES Advisory Board Meeting. July 14th 2014 Edoardo Di Napoli Academic background
Comparing the Efficiency of Iterative Eigenvalue Solvers: the Quantum ESPRESSO experience
Comparing the Efficiency of Iterative Eigenvalue Solvers: the Quantum ESPRESSO experience Stefano de Gironcoli Scuola Internazionale Superiore di Studi Avanzati Trieste-Italy 0 Diagonalization of the Kohn-Sham
Verbundprojekt ELPA-AEO. Eigenwert-Löser für Petaflop-Anwendungen Algorithmische Erweiterungen und Optimierungen
Verbundprojekt ELPA-AEO http://elpa-aeo.mpcdf.mpg.de Eigenwert-Löser für Petaflop-Anwendungen Algorithmische Erweiterungen und Optimierungen BMBF Projekt 01IH15001 Feb 2016 - Jan 2019 7. HPC-Statustagung,
Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
Parallel Eigensolver Performance on the HPCx System
Parallel Eigensolver Performance on the HPCx System Andrew Sunderland, Elena Breitmoser Terascaling Applications Group CCLRC Daresbury Laboratory EPCC, University of Edinburgh Outline 1. Brief Introduction
Experiences with Self-Consistent Tight Binding and ELSI
Experiences with Self-Consistent Tight Binding and ELSI Ben Hourahine benjamin.hourahine@strath.ac.uk Motivation via computational cost Method Eform. (ev) Time Orthogonal TB 3.2 1 Self-Consistent Orthogonal
Parallelization of the Molecular Orbital Program MOS-F
Parallelization of the Molecular Orbital Program MOS-F Akira Asato, Satoshi Onodera, Yoshie Inada, Elena Akhmatskaya, Ross Nobes, Azuma Matsuura, Atsuya Takahashi November 2003 Fujitsu Laboratories of
CP2K. New Frontiers. ab initio Molecular Dynamics
CP2K New Frontiers in ab initio Molecular Dynamics Jürg Hutter, Joost VandeVondele, Valery Weber Physical-Chemistry Institute, University of Zurich Ab Initio Molecular Dynamics Molecular Dynamics Sampling
FEAST eigenvalue algorithm and solver: review and perspectives
FEAST eigenvalue algorithm and solver: review and perspectives Eric Polizzi Department of Electrical and Computer Engineering University of Masachusetts, Amherst, USA Sparse Days, CERFACS, June 25, 2012
Recent developments of Discontinuous Galerkin Density functional theory
Lin Lin Discontinuous Galerkin DFT 1 Recent developments of Discontinuous Galerkin Density functional theory Lin Lin Department of Mathematics, UC Berkeley; Computational Research Division, LBNL Numerical
Large Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
ESLW_Drivers July 2017
ESLW_Drivers 10-21 July 2017 Volker Blum - ELSI Viktor Yu - ELSI William Huhn - ELSI David Lopez - Siesta Yann Pouillon - Abinit Micael Oliveira Octopus & Abinit Fabiano Corsetti Siesta & Onetep Paolo
INITIAL INTEGRATION AND EVALUATION
INITIAL INTEGRATION AND EVALUATION OF SLATE PARALLEL BLAS IN LATTE Marc Cawkwell, Danny Perez, Arthur Voter Asim YarKhan, Gerald Ragghianti, Jack Dongarra, Introduction The aim of the joint milestone STMS10-52
References. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
ELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS
FROM RESEARCH TO INDUSTRY 32 ème forum ORAP 10 octobre 2013 Maison de la Simulation, Saclay, France ELECTRONIC STRUCTURE CALCULATIONS FOR THE SOLID STATE PHYSICS APPLICATION ON HPC, BLOCKING POINTS, Marc
CRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel?
CRYSTAL in parallel: replicated and distributed (MPP) data Roberto Orlando Dipartimento di Chimica Università di Torino Via Pietro Giuria 5, 10125 Torino (Italy) roberto.orlando@unito.it 1 Why parallel?
Making electronic structure methods scale: Large systems and (massively) parallel computing
AB Making electronic structure methods scale: Large systems and (massively) parallel computing Ville Havu Department of Applied Physics Helsinki University of Technology - TKK Ville.Havu@tkk.fi 1 Outline
The ELPA Library Scalable Parallel Eigenvalue Solutions for Electronic Structure Theory and Computational Science
TOPICAL REVIEW The ELPA Library Scalable Parallel Eigenvalue Solutions for Electronic Structure Theory and Computational Science Andreas Marek 1, Volker Blum 2,3, Rainer Johanni 1,2 ( ), Ville Havu 4,
The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and
Home Search Collections Journals About Contact us My IOPscience The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and computational science This content has been
Before we start: Important setup of your Computer
Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann
Parallel Eigensolver Performance on High Performance Computers
Parallel Eigensolver Performance on High Performance Computers Andrew Sunderland Advanced Research Computing Group STFC Daresbury Laboratory CUG 2008 Helsinki 1 Summary (Briefly) Introduce parallel diagonalization
PFEAST: A High Performance Sparse Eigenvalue Solver Using Distributed-Memory Linear Solvers
PFEAST: A High Performance Sparse Eigenvalue Solver Using Distributed-Memory Linear Solvers James Kestyn, Vasileios Kalantzis, Eric Polizzi, Yousef Saad Electrical and Computer Engineering Department,
DFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
Efficient implementation of the overlap operator on multi-gpus
Efficient implementation of the overlap operator on multi-gpus Andrei Alexandru Mike Lujan, Craig Pelissier, Ben Gamari, Frank Lee SAAHPC 2011 - University of Tennessee Outline Motivation Overlap operator
Atomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
Cyclops Tensor Framework
Cyclops Tensor Framework Edgar Solomonik Department of EECS, Computer Science Division, UC Berkeley March 17, 2014 1 / 29 Edgar Solomonik Cyclops Tensor Framework 1/ 29 Definition of a tensor A rank r
COMPUTATIONAL TOOL. Fig. 4.1 Opening screen of w2web
CHAPTER -4 COMPUTATIONAL TOOL Ph.D. Thesis: J. Maibam CHAPTER: 4 4.1 The WIEN2k code In this work, all the calculations presented are performed using the WIEN2k software package (Blaha et al., 2001). The
Massively 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
A Reconfigurable Quantum Computer
A Reconfigurable Quantum Computer David Moehring CEO, IonQ, Inc. College Park, MD Quantum Computing for Business 4-6 December 2017, Mountain View, CA IonQ Highlights Full Stack Quantum Computing Company
ILAS 2017 July 25, 2017
Polynomial and rational filtering for eigenvalue problems and the EVSL project Yousef Saad Department of Computer Science and Engineering University of Minnesota ILAS 217 July 25, 217 Large eigenvalue
CP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
Basic introduction of NWChem software
Basic introduction of NWChem software Background! NWChem is part of the Molecular Science Software Suite! Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
4th year Project demo presentation
4th year Project demo presentation Colm Ó héigeartaigh CASE4-99387212 coheig-case4@computing.dcu.ie 4th year Project demo presentation p. 1/23 Table of Contents An Introduction to Quantum Computing The
DFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
ab initio Electronic Structure Calculations
ab initio Electronic Structure Calculations New scalability frontiers using the BG/L Supercomputer C. Bekas, A. Curioni and W. Andreoni IBM, Zurich Research Laboratory Rueschlikon 8803, Switzerland ab
Improvements for Implicit Linear Equation Solvers
Improvements for Implicit Linear Equation Solvers Roger Grimes, Bob Lucas, Clement Weisbecker Livermore Software Technology Corporation Abstract Solving large sparse linear systems of equations is often
Parallel Eigensolver Performance on High Performance Computers 1
Parallel Eigensolver Performance on High Performance Computers 1 Andrew Sunderland STFC Daresbury Laboratory, Warrington, UK Abstract Eigenvalue and eigenvector computations arise in a wide range of scientific
Opportunities for ELPA to Accelerate the Solution of the Bethe-Salpeter Eigenvalue Problem
Opportunities for ELPA to Accelerate the Solution of the Bethe-Salpeter Eigenvalue Problem Peter Benner, Andreas Marek, Carolin Penke August 16, 2018 ELSI Workshop 2018 Partners: The Problem The Bethe-Salpeter
Efficient algorithms for symmetric tensor contractions
Efficient algorithms for symmetric tensor contractions Edgar Solomonik 1 Department of EECS, UC Berkeley Oct 22, 2013 1 / 42 Edgar Solomonik Symmetric tensor contractions 1/ 42 Motivation The goal is to
Large-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
Algorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part II Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
Matrix Eigensystem Tutorial For Parallel Computation
Matrix Eigensystem Tutorial For Parallel Computation High Performance Computing Center (HPC) http://www.hpc.unm.edu 5/21/2003 1 Topic Outline Slide Main purpose of this tutorial 5 The assumptions made
Mitglied der Helmholtz-Gemeinschaft. Linear algebra tasks in Materials Science: optimization and portability
Mitglied der Helmholtz-Gemeinschaft Linear algebra tasks in Materials Science: optimization and portability ADAC Workshop, July 17-19 2017 Edoardo Di Napoli Outline Jülich Supercomputing Center Chebyshev
DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
The Augmented Spherical Wave Method
Introduction Institut für Physik, Universität Augsburg Electronic Structure in a Nutshell Outline 1 Fundamentals Generations 2 Outline 1 Fundamentals Generations 2 Outline Fundamentals Generations 1 Fundamentals
OPENATOM for GW calculations
OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW
Improving the performance of applied science numerical simulations: an application to Density Functional Theory
Improving the performance of applied science numerical simulations: an application to Density Functional Theory Edoardo Di Napoli Jülich Supercomputing Center - Institute for Advanced Simulation Forschungszentrum
TR 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
Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale
Fortran Expo: 15 Jun 2012 Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale Lianheng Tong Overview Overview of Conquest project Brief Introduction
Presentation of XLIFE++
Presentation of XLIFE++ Eigenvalues Solver & OpenMP Manh-Ha NGUYEN Unité de Mathématiques Appliquées, ENSTA - Paristech 25 Juin 2014 Ha. NGUYEN Presentation of XLIFE++ 25 Juin 2014 1/19 EigenSolver 1 EigenSolver
Metal Conquest. Lianheng Tong. 2 nd March 2011
Metal Conquest Lianheng Tong nd March 011 Abstract This report describes the work done in the one year Distributed Computational Science and Engineering (dcse) project aimed to develop an ab initio Density
Performance optimization of WEST and Qbox on Intel Knights Landing
Performance optimization of WEST and Qbox on Intel Knights Landing Huihuo Zheng 1, Christopher Knight 1, Giulia Galli 1,2, Marco Govoni 1,2, and Francois Gygi 3 1 Argonne National Laboratory 2 University
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
Sakurai-Sugiura algorithm based eigenvalue solver for Siesta. Georg Huhs
Sakurai-Sugiura algorithm based eigenvalue solver for Siesta Georg Huhs Motivation Timing analysis for one SCF-loop iteration: left: CNT/Graphene, right: DNA Siesta Specifics High fraction of EVs needed
Introduction to Benchmark Test for Multi-scale Computational Materials Software
Introduction to Benchmark Test for Multi-scale Computational Materials Software Shun Xu*, Jian Zhang, Zhong Jin xushun@sccas.cn Computer Network Information Center Chinese Academy of Sciences (IPCC member)
Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
Eigenvalue Problems CHAPTER 1 : PRELIMINARIES
Eigenvalue Problems CHAPTER 1 : PRELIMINARIES Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Preliminaries Eigenvalue problems 2012 1 / 14
An Optimized and Scalable Eigensolver for Sequences of Eigenvalue Problems
An Optimized and Scalable Eigensolver for Sequences of Eigenvalue Problems Mario Berljafa Daniel Wortmann Edoardo Di Napoli arxiv:1404.4161v2 [cs.ms] 6 Jul 2014 July 8, 2014 Abstract In many scientific
Divide and conquer algorithms for large eigenvalue problems Yousef Saad Department of Computer Science and Engineering University of Minnesota
Divide and conquer algorithms for large eigenvalue problems Yousef Saad Department of Computer Science and Engineering University of Minnesota PMAA 14 Lugano, July 4, 2014 Collaborators: Joint work with:
Scalable Hybrid Programming and Performance for SuperLU Sparse Direct Solver
Scalable Hybrid Programming and Performance for SuperLU Sparse Direct Solver Sherry Li Lawrence Berkeley National Laboratory Piyush Sao Rich Vuduc Georgia Institute of Technology CUG 14, May 4-8, 14, Lugano,
Leveraging sparsity and symmetry in parallel tensor contractions
Leveraging sparsity and symmetry in parallel tensor contractions Edgar Solomonik Department of Computer Science University of Illinois at Urbana-Champaign CCQ Tensor Network Workshop Flatiron Institute,
The Quantum ESPRESSO Software Distribution
The Quantum ESPRESSO Software Distribution The DEMOCRITOS center of Italian INFM is dedicated to atomistic simulations of materials, with a strong emphasis on the development of high-quality scientific
Time-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
Contents. 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
Density Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
A Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters
A Quantum Chemistry Domain-Specific Language for Heterogeneous Clusters ANTONINO TUMEO, ORESTE VILLA Collaborators: Karol Kowalski, Sriram Krishnamoorthy, Wenjing Ma, Simone Secchi May 15, 2012 1 Outline!
Basic introduction of NWChem software
Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
Preliminary Results of GRAPES Helmholtz solver using GCR and PETSc tools
Preliminary Results of GRAPES Helmholtz solver using GCR and PETSc tools Xiangjun Wu (1),Lilun Zhang (2),Junqiang Song (2) and Dehui Chen (1) (1) Center for Numerical Weather Prediction, CMA (2) School
2.1 Introduction: The many-body problem
Chapter 2 Smeagol: Density Functional Theory and NEGF s 2.1 Introduction: The many-body problem In solid state physics one is interested in systems comprising many atoms, and consequently many electrons.
Strassen s Algorithm for Tensor Contraction
Strassen s Algorithm for Tensor Contraction Jianyu Huang, Devin A. Matthews, Robert A. van de Geijn The University of Texas at Austin September 14-15, 2017 Tensor Computation Workshop Flatiron Institute,
Accelerating Linear Algebra on Heterogeneous Architectures of Multicore and GPUs using MAGMA and DPLASMA and StarPU Schedulers
UT College of Engineering Tutorial Accelerating Linear Algebra on Heterogeneous Architectures of Multicore and GPUs using MAGMA and DPLASMA and StarPU Schedulers Stan Tomov 1, George Bosilca 1, and Cédric
Solving the Inverse Toeplitz Eigenproblem Using ScaLAPACK and MPI *
Solving the Inverse Toeplitz Eigenproblem Using ScaLAPACK and MPI * J.M. Badía and A.M. Vidal Dpto. Informática., Univ Jaume I. 07, Castellón, Spain. badia@inf.uji.es Dpto. Sistemas Informáticos y Computación.
Porting a sphere optimization program from LAPACK to ScaLAPACK
Porting a sphere optimization program from LAPACK to ScaLAPACK Mathematical Sciences Institute, Australian National University. For presentation at Computational Techniques and Applications Conference
VASP: running on HPC resources. University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria
VASP: running on HPC resources University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria The Many-Body Schrödinger equation 0 @ 1 2 X i i + X i Ĥ (r 1,...,r
1 Overview. 2 Adapting to computing system evolution. 11 th European LS-DYNA Conference 2017, Salzburg, Austria
1 Overview Improving LSTC s Multifrontal Linear Solver Roger Grimes 3, Robert Lucas 3, Nick Meng 2, Francois-Henry Rouet 3, Clement Weisbecker 3, and Ting-Ting Zhu 1 1 Cray Incorporated 2 Intel Corporation
Domain specific libraries. Material science codes on innovative HPC architectures Anton Kozhevnikov, CSCS December 5, 2016
Domain specific libraries Material science codes on innovative HPC architectures Anton Kozhevnikov, CSCS December 5, 2016 Part 1: Introduction Kohn-Shame equations 1 2 Eigen-value problem + v eff (r) j(r)
Arnoldi Methods in SLEPc
Scalable Library for Eigenvalue Problem Computations SLEPc Technical Report STR-4 Available at http://slepc.upv.es Arnoldi Methods in SLEPc V. Hernández J. E. Román A. Tomás V. Vidal Last update: October,
Preconditioned Eigenvalue Solvers for electronic structure calculations. Andrew V. Knyazev. Householder Symposium XVI May 26, 2005
1 Preconditioned Eigenvalue Solvers for electronic structure calculations Andrew V. Knyazev Department of Mathematics and Center for Computational Mathematics University of Colorado at Denver Householder
Opportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems
Seminar CENTRE FOR PREDICTIVE MODELLING, WARWICK Opportunities from Accurate and Efficient Density Functional Theory Calculations for Large Systems Luigi Genovese L_Sim CEA Grenoble October 30, 2017 http://bigdft.org
Scalability Programme at ECMWF
Scalability Programme at ECMWF Picture: Stan Tomov, ICL, University of Tennessee, Knoxville Peter Bauer, Mike Hawkins, George Mozdzynski, Tiago Quintino, Deborah Salmond, Stephan Siemen, Yannick Trémolet
All-electron density functional theory on Intel MIC: Elk
All-electron density functional theory on Intel MIC: Elk W. Scott Thornton, R.J. Harrison Abstract We present the results of the porting of the full potential linear augmented plane-wave solver, Elk [1],
MODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
Weather Research and Forecasting (WRF) Performance Benchmark and Profiling. July 2012
Weather Research and Forecasting (WRF) Performance Benchmark and Profiling July 2012 Note The following research was performed under the HPC Advisory Council activities Participating vendors: Intel, Dell,
J S Parker (QUB), Martin Plummer (STFC), H W van der Hart (QUB) Version 1.0, September 29, 2015
Report on ecse project Performance enhancement in R-matrix with time-dependence (RMT) codes in preparation for application to circular polarised light fields J S Parker (QUB), Martin Plummer (STFC), H
Susumu YAMADA 1,3 Toshiyuki IMAMURA 2,3, Masahiko MACHIDA 1,3
Dynamical Variation of Eigenvalue Problems in Density-Matrix Renormalization-Group Code PP12, Feb. 15, 2012 1 Center for Computational Science and e-systems, Japan Atomic Energy Agency 2 The University
Density-Matrix-Based Algorithms for Solving Eingenvalue Problems
University of Massachusetts - Amherst From the SelectedWorks of Eric Polizzi 2009 Density-Matrix-Based Algorithms for Solving Eingenvalue Problems Eric Polizzi, University of Massachusetts - Amherst Available
PASC17 - Lugano June 27, 2017
Polynomial and rational filtering for eigenvalue problems and the EVSL project Yousef Saad Department of Computer Science and Engineering University of Minnesota PASC17 - Lugano June 27, 217 Large eigenvalue
H ψ = E ψ. Introduction to Exact Diagonalization. Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden
H ψ = E ψ Introduction to Exact Diagonalization Andreas Läuchli, New states of quantum matter MPI für Physik komplexer Systeme - Dresden http://www.pks.mpg.de/~aml laeuchli@comp-phys.org Simulations of
CP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
Introduction 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
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on Matthew R. Norman Scientific Computing Group National Center for Computational Sciences Oak Ridge
AMS526: 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
Symmetric Pivoting in ScaLAPACK Craig Lucas University of Manchester Cray User Group 8 May 2006, Lugano
Symmetric Pivoting in ScaLAPACK Craig Lucas University of Manchester Cray User Group 8 May 2006, Lugano Introduction Introduction We wanted to parallelize a serial algorithm for the pivoted Cholesky factorization
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
Accelerating and Scaling Lanczos Diagonalization with GPGPU
Accelerating and Scaling Lanczos Diagonalization with GPGPU Bill Brouwer, Filippo Spiga, Pierre-Yves Taunay, Sreejith GJ Nvidia GTC 2013 Outline Introduction Applications QE FQHE Theory Diagonalization
The Plane-wave Pseudopotential Method
The Plane-wave Pseudopotential Method k(r) = X G c k,g e i(g+k) r Chris J Pickard Electrons in a Solid Nearly Free Electrons Nearly Free Electrons Nearly Free Electrons Electronic Structures Methods Empirical
Performance Analysis of Lattice QCD Application with APGAS Programming Model
Performance Analysis of Lattice QCD Application with APGAS Programming Model Koichi Shirahata 1, Jun Doi 2, Mikio Takeuchi 2 1: Tokyo Institute of Technology 2: IBM Research - Tokyo Programming Models
MAGMA MIC 1.0: Linear Algebra Library for Intel Xeon Phi Coprocessors
MAGMA MIC 1.0: Linear Algebra Library for Intel Xeon Phi Coprocessors J. Dongarra, M. Gates, A. Haidar, Y. Jia, K. Kabir, P. Luszczek, and S. Tomov University of Tennessee, Knoxville 05 / 03 / 2013 MAGMA: