Abstract:The computation of electron repulsion integrals (ERIs) over Gaussian-type orbitals (GTOs) is a challenging problem in quantum-mechanics-based atomistic simulations. In practical simulations, several trillions of ERIs may have to be computed for every time step. In this work, we investigate FPGAs as accelerators for the ERI computation. We use template parameters, here within the Intel oneAPI tool flow, to create customized designs for 256 different ERI quartet classes, based on their orbitals. To maximize data reuse, all intermediates are buffered in FPGA on-chip memory with customized layout. The pre-calculation of intermediates also helps to overcome data dependencies caused by multi-dimensional recurrence relations. The involved loop structures are partially or even fully unrolled for high throughput of FPGA kernels. Furthermore, a lossy compression algorithm utilizing arbitrary bitwidth integers is integrated in the FPGA kernels. To our best knowledge, this is the first work on ERI computation on FPGAs that supports more than just the single most basic quartet class. Also, the integration of ERI computation and compression it a novelty that is not even covered by CPU or GPU libraries so far. Our evaluation shows that using 16-bit integer for the ERI compression, the fastest FPGA kernels exceed the performance of 10 GERIS ($10 \times 10^9$ ERIs per second) on one Intel Stratix 10 GX 2800 FPGA, with maximum absolute errors around $10^{-7}$ - $10^{-5}$ Hartree. The measured throughput can be accurately explained by a performance model. The FPGA kernels deployed on 2 FPGAs outperform similar computations using the widely used libint reference on a two-socket server with 40 Xeon Gold 6148 CPU cores of the same process technology by factors up to 6.0x and on a new two-socket server with 128 EPYC 7713 CPU cores by up to 1.9x.

A computation of D(9) using FPGA Supercomputing

A computation of the ninth Dedekind Number

Solving systems of equations on antichains for the computation of the ninth Dedekind Number

Computation of a 30750-Bit Binary Field Discrete Logarithm

Fast computation of generalized dedekind sums

Computing and Compressing Electron Repulsion Integrals on FPGAs

Numerical Verification of the Ternary Goldbach Conjecture up to 8.875e30

Highly optimized quantum circuits synthesized via data-flow engines

Efficient FPGA-based multipliers for F_{3^97} and F_{3^{6*97}}

High-Efficiency Compressor Trees for Latest AMD FPGAs

Error Characteristics of Fiber Distributed Data Interface (FDDI)

Efficient Ab-Initio Molecular Dynamic Simulations by Offloading Fast Fourier Transformations to FPGAs

Towards Lattice Quantum Chromodynamics on FPGA devices

Efficient Fpga Emulation Of Quantum Fourier Transform

High Performance Scalable FPGA Accelerator for Deep Neural Networks

An algorithm and computation to verify Legendre's Conjecture up to $3.33\cdot10^{13}$

Optimizing FPGA-Based DNN Accelerator with Shared Exponential Floating-Point Format

Scalable and Parallel Optimization of the Number Theoretic Transform Based on FPGA

GH CORDIC-Based Architecture for Computing $N$ Th Root of Single-Precision Floating-Point Number

Specific Hardware Implementation for Cofactorization in Gnfs