TY - GEN
AB - CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular and biological systems. It is especially aimed at massively-parallel and linear-scaling electronic structure methods and state-of-the-art ab-initio molecular dynamics simulations. Excellent performance for electronic structure calculations is achieved using novel algorithms implemented for modern high-performance computing systems. This review revisits the main capabilities of CP2k to perform efficient and accurate electronic structure simulations. The emphasis is put on density functional theory and multiple post-Hartree-Fock methods using the Gaussian and plane wave approach and its augmented all-electron extension.
AU - Kühne, Thomas
AU - Iannuzzi, Marcella
AU - Ben, Mauro Del
AU - Rybkin, Vladimir V.
AU - Seewald, Patrick
AU - Stein, Frederick
AU - Laino, Teodoro
AU - Khaliullin, Rustam Z.
AU - Schütt, Ole
AU - Schiffmann, Florian
AU - Golze, Dorothea
AU - Wilhelm, Jan
AU - Chulkov, Sergey
AU - Mohammad Hossein Bani-Hashemian, Mohammad Hossein Bani-Hashemian
AU - Weber, Valéry
AU - Borstnik, Urban
AU - Taillefumier, Mathieu
AU - Jakobovits, Alice Shoshana
AU - Lazzaro, Alfio
AU - Pabst, Hans
AU - Müller, Tiziano
AU - Schade, Robert
AU - Guidon, Manuel
AU - Andermatt, Samuel
AU - Holmberg, Nico
AU - Schenter, Gregory K.
AU - Hehn, Anna
AU - Bussy, Augustin
AU - Belleflamme, Fabian
AU - Tabacchi, Gloria
AU - Glöß, Andreas
AU - Lass, Michael
AU - Bethune, Iain
AU - Mundy, Christopher J.
AU - Plessl, Christian
AU - Watkins, Matt
AU - VandeVondele, Joost
AU - Krack, Matthias
AU - Hutter, Jürg
ID - 16277
TI - CP2K: An Electronic Structure and Molecular Dynamics Software Package -- Quickstep: Efficient and Accurate Electronic Structure Calculations
ER -
TY - GEN
AB - In scientific computing, the acceleration of atomistic computer simulations
by means of custom hardware is finding ever growing application. A major
limitation, however, is that the high efficiency in terms of performance and
low power consumption entails the massive usage of low-precision computing
units. Here, based on the approximate computing paradigm, we present an
algorithmic method to rigorously compensate for numerical inaccuracies due to
low-accuracy arithmetic operations, yet still obtaining exact expectation
values using a properly modified Langevin-type equation.
AU - Rengaraj, Varadarajan
AU - Lass, Michael
AU - Plessl, Christian
AU - Kühne, Thomas
ID - 12878
T2 - arXiv:1907.08497
TI - Accurate Sampling with Noisy Forces from Approximate Computing
ER -
TY - CONF
AB - Stratix 10 FPGA cards have a good potential for the acceleration of HPC workloads since the Stratix 10 product line introduces devices with a large number of DSP and memory blocks. The high level synthesis of OpenCL codes can play a fundamental role for FPGAs in HPC, because it allows to implement different designs with lower development effort compared to hand optimized HDL. However, Stratix 10 cards are still hard to fully exploit using the Intel FPGA SDK for OpenCL. The implementation of designs with thousands of concurrent arithmetic operations often suffers from place and route problems that limit the maximum frequency or entirely prevent a successful synthesis. In order to overcome these issues for the implementation of the matrix multiplication, we formulate Cannon's matrix multiplication algorithm with regard to its efficient synthesis within the FPGA logic. We obtain a two-level block algorithm, where the lower level sub-matrices are multiplied using our Cannon's algorithm implementation. Following this design approach with multiple compute units, we are able to get maximum frequencies close to and above 300 MHz with high utilization of DSP and memory blocks. This allows for performance results above 1 TeraFLOPS.
AU - Gorlani, Paolo
AU - Kenter, Tobias
AU - Plessl, Christian
ID - 15478
T2 - Proceedings of the International Conference on Field-Programmable Technology (FPT)
TI - OpenCL Implementation of Cannon's Matrix Multiplication Algorithm on Intel Stratix 10 FPGAs
ER -
TY - CONF
AB - This paper describes a data structure and a heuristic to plan and map arbitrary resources in complex combinations while applying time dependent constraints. The approach is used in the planning based workload manager OpenCCS at the Paderborn Center for Parallel Computing (PC\(^2\)) to operate heterogeneous clusters with up to 10000 cores. We also show performance results derived from four years of operation.
AU - Keller, Axel
ED - Klusáček, D.
ED - Cirne, W.
ED - Desai, N.
ID - 22
KW - Scheduling Planning Mapping Workload management
SN - 978-3-319-77398-8
T2 - Proc. Workshop on Job Scheduling Strategies for Parallel Processing (JSSPP)
TI - A Data Structure for Planning Based Workload Management of Heterogeneous HPC Systems
VL - 10773
ER -
TY - CONF
AB - The exploration of FPGAs as accelerators for scientific simulations has so far mostly been focused on small kernels of methods working on regular data structures, for example in the form of stencil computations for finite difference methods. In computational sciences, often more advanced methods are employed that promise better stability, convergence, locality and scaling. Unstructured meshes are shown to be more effective and more accurate, compared to regular grids, in representing computation domains of various shapes. Using unstructured meshes, the discontinuous Galerkin method preserves the ability to perform explicit local update operations for simulations in the time domain. In this work, we investigate FPGAs as target platform for an implementation of the nodal discontinuous Galerkin method to find time-domain solutions of Maxwell's equations in an unstructured mesh. When maximizing data reuse and fitting constant coefficients into suitably partitioned on-chip memory, high computational intensity allows us to implement and feed wide data paths with hundreds of floating point operators. By decoupling off-chip memory accesses from the computations, high memory bandwidth can be sustained, even for the irregular access pattern required by parts of the application. Using the Intel/Altera OpenCL SDK for FPGAs, we present different implementation variants for different polynomial orders of the method. In different phases of the algorithm, either computational or bandwidth limits of the Arria 10 platform are almost reached, thus outperforming a highly multithreaded CPU implementation by around 2x.
AU - Kenter, Tobias
AU - Mahale, Gopinath
AU - Alhaddad, Samer
AU - Grynko, Yevgen
AU - Schmitt, Christian
AU - Afzal, Ayesha
AU - Hannig, Frank
AU - Förstner, Jens
AU - Plessl, Christian
ID - 1588
KW - tet_topic_hpc
T2 - Proc. Int. Symp. on Field-Programmable Custom Computing Machines (FCCM)
TI - OpenCL-based FPGA Design to Accelerate the Nodal Discontinuous Galerkin Method for Unstructured Meshes
ER -
TY - CONF
AU - Riebler, Heinrich
AU - Vaz, Gavin Francis
AU - Kenter, Tobias
AU - Plessl, Christian
ID - 1204
KW - htrop
SN - 9781450349826
T2 - Proc. ACM SIGPLAN Symposium on Principles and Practice of Parallel Programming (PPoPP)
TI - Automated Code Acceleration Targeting Heterogeneous OpenCL Devices
ER -
TY - CONF
AB - We present the submatrix method, a highly parallelizable method for the approximate calculation of inverse p-th roots of large sparse symmetric matrices which are required in different scientific applications. Following the idea of Approximate Computing, we allow imprecision in the final result in order to utilize the sparsity of the input matrix and to allow massively parallel execution. For an n x n matrix, the proposed algorithm allows to distribute the calculations over n nodes with only little communication overhead. The result matrix exhibits the same sparsity pattern as the input matrix, allowing for efficient reuse of allocated data structures.
We evaluate the algorithm with respect to the error that it introduces into calculated results, as well as its performance and scalability. We demonstrate that the error is relatively limited for well-conditioned matrices and that results are still valuable for error-resilient applications like preconditioning even for ill-conditioned matrices. We discuss the execution time and scaling of the algorithm on a theoretical level and present a distributed implementation of the algorithm using MPI and OpenMP. We demonstrate the scalability of this implementation by running it on a high-performance compute cluster comprised of 1024 CPU cores, showing a speedup of 665x compared to single-threaded execution.
AU - Lass, Michael
AU - Mohr, Stephan
AU - Wiebeler, Hendrik
AU - Kühne, Thomas
AU - Plessl, Christian
ID - 1590
KW - approximate computing
KW - linear algebra
KW - matrix inversion
KW - matrix p-th roots
KW - numeric algorithm
KW - parallel computing
SN - 978-1-4503-5891-0/18/07
T2 - Proc. Platform for Advanced Scientific Computing (PASC) Conference
TI - A Massively Parallel Algorithm for the Approximate Calculation of Inverse p-th Roots of Large Sparse Matrices
ER -
TY - CONF
AU - Kenter, Tobias
AU - Förstner, Jens
AU - Plessl, Christian
ID - 1592
KW - tet_topic_hpc
T2 - Proc. Int. Conf. on Field Programmable Logic and Applications (FPL)
TI - Flexible FPGA design for FDTD using OpenCL
ER -
TY - CONF
AU - Kenter, Tobias
AU - Plessl, Christian
ID - 24
T2 - Proc. Workshop on Heterogeneous High-performance Reconfigurable Computing (H2RC)
TI - Microdisk Cavity FDTD Simulation on FPGA using OpenCL
ER -
TY - CONF
AU - Riebler, Heinrich
AU - Vaz, Gavin Francis
AU - Plessl, Christian
AU - Trainiti, Ettore M. G.
AU - Durelli, Gianluca C.
AU - Bolchini, Cristiana
ID - 31
T2 - Proc. HiPEAC Workshop on Reonfigurable Computing (WRC)
TI - Using Just-in-Time Code Generation for Transparent Resource Management in Heterogeneous Systems
ER -