Abstract:Fast and accurate numerical integration always represented a bottleneck in high-performance computational physics, especially in large and multiscale industrial simulations involving Finite (FEM) and Boundary Element Methods (BEM). The computational demand escalates significantly in problems modelled by irregular or endpoint singular behaviours which can be approximated with generalised polynomials of real degree. This is due to both the practical limitations of finite-arithmetic computations and the inefficient samples distribution of traditional Gaussian quadrature rules. We developed a non-iterative mathematical software implementing an innovative numerical quadrature which largely enhances the precision of Gauss-Legendre formulae (G-L) for integrands modelled as generalised polynomial with the optimal amount of nodes and weights capable of guaranteeing the required numerical precision. This methodology avoids to resort to more computationally expensive techniques such as adaptive or composite quadrature rules. From a theoretical point of view, the numerical method underlying this work was preliminary presented in [1] by constructing the monomial transformation itself and providing all the necessary conditions to ensure the numerical stability and exactness of the quadrature up to machine precision. The novel contribution of this work concerns the optimal implementation of said method, the extension of its applicability at run-time with different type of inputs, the provision of additional insights on its functionalities and its straightforward implementation, in particular FEM applications or other mathematical software either as an external tool or embedded suite. The open-source, cross-platform C++ library Monomial Transformation Quadrature Rule (MTQR) has been designed to be highly portable, fast and easy to integrate in larger codebases. Numerical examples in multiple physical applications showcase the improved efficiency and accuracy when compared to traditional schemes. Program summary Program title: MTQR CPC Library link to program files: https://doi.org/10.17632/276f78wzsw.1 Developer's repository link: https://github.com/MTQR/MTQR Licensing provisions: GNU General Public License 3 Programming language: C++ (C++17 standard) Supplementary material: User manual (for installation and execution) Nature of problem: Accuracy and time of execution of the current implementations of high-precision numerical integration routines for singular and irregular integrands modelled by generalised polynomials are restricted by: limitations of the floating-point (f.p.) finite-arithmetic of the machine; inability of classical Gaussian quadrature rules to efficiently capture irregular behaviours or end-point singularities using an optimal number of samples; relying on significantly expensive techniques as adaptive or composite quadrature rules that severely increase the number of steps necessary to converge to the desired accuracy threshold. However by precisely manipulating the G-L samples using an ad-hoc monomial transformation we achieve a one-shot, non-iterative, machine-precise quadrature rule with straightforward scalability in higher dimensions. The advantages brought by a non-adaptive technique are greatly emphasized whenever the problem on hand is characterised by a numerous set of integrand functions that behave similarly to sets of generalised polynomials. Solution method: The underlying algorithm is an application of a monomial transformation of (real) order to the nmin∈N quadrature samples (nodes and weights) of the G-L rule. The order of the transformation depends on the infimum (λmin) and supremum (λmax) of the Müntz sequence of real degrees of the integrands modelled by generalised polynomials. With an a-priori full analysis of the set of integrands to be integrated, the design and the action of such transformation ensures that the bounding monomials are integrated with machine-epsilon precision in double floating-point (f.p.) format [2] without resorting to less efficient schemes as adaptive or composite quadrature rules, singularity subtraction and cancellation methods with limited and uncontrolled precision. The effects of the action of the monomial transformation onto the original G-L samples is clearly visible by the clustering of the nodes around the endpoint singularity of the integrand function (see Figure 1.2 of the user manual shipped with the source code and located in the MTQR/doc sub-directory as indic -Abstract Truncated-

High-performance evaluation of high angular momentum 4-center Gaussian integrals on modern accelerated processors

3-center and 4-center 2-particle Gaussian AO integrals on modern accelerated processors

High-performance GPU-accelerated evaluation of electron repulsion integrals

Parallelized Implementation of the Finite Particle Method for Explicit Dynamics in GPU

Libcint: an Efficient General Integral Library for Gaussian Basis Functions

Efficient evaluation of two-center Gaussian integrals in periodic systems

Numerical Performance and Throughput Benchmark for Electronic Structure Calculations in Pc-Linux Systems with New Architectures, Updated Compilers, and Libraries

Optimizing two-electron repulsion integral calculations with McMurchie–Davidson method on graphic processing unit

Extending GPU-Accelerated Gaussian Integrals in the TeraChem Software Package to f Type Orbitals: Implementation and Applications

A High-Performance Implementation of Atomistic Spin Dynamics Simulations on x86 CPUs

Multi-GPU RI-HF Energies and Analytic Gradients $-$ Towards High Throughput Ab Initio Molecular Dynamics

Performance engineering for the Lattice Boltzmann method on GPGPUs: Architectural requirements and performance results

Quadrature of functions with endpoint singular and generalised polynomial behaviour in computational physics

Accelerating Relativistic Exact-Two-Component Density Functional Theory Calculations with Graphical Processing Units

Efficient molecular dynamics simulations with many-body potentials on graphics processing units

A High Performance Block Eigensolver for Nuclear Configuration Interaction Calculations.

GPGPU Acceleration of All-Electron Electronic Structure Theory Using Localized Numeric Atom-Centered Basis Functions

Kernel Fusion in Atomistic Spin Dynamics Simulations on Nvidia GPUs using Tensor Core

Computationally Efficient Molecular Integrals of Solid Harmonic Gaussian Orbitals Using Quantum Entanglement of Angular Momentum

GPU Acceleration of Numerical Atomic Orbitals-Based Density Functional Theory Algorithms within the ABACUS package

Double Exponential Transformation For Computing Three-Centre Nuclear Attraction Integrals