Abstract:Large linear systems play an important role in high-energy theory, appearing in amplitude bootstraps and during integral reduction. This paper introduces FiniteFieldSolve , a general-purpose toolkit for exactly solving large linear systems over the rationals. The solver interfaces directly with Mathematica, is straightforward to install, and seamlessly replaces Mathematica's native solvers. In testing, FiniteFieldSolve is approximately two orders of magnitude faster than Mathematica and uses an order of magnitude less memory. The package also compares favorably against other public solvers in FiniteFieldSolve 's intended use cases. As the name of the package suggests, solutions are obtained via well-known finite field methods. These methods suffer from introducing an inordinate number of modulo (or integer division) operations with respect to different primes. By automatically recompiling itself for each prime, FiniteFieldSolve converts the division operations into much faster combinations of instructions, dramatically improving performance. The technique of compiling the prime can be applied to any finite field solver, where the time savings will be solver dependent. The operation of the package is illustrated through a detailed example of an amplitude bootstrap. Program summary Program Title: FiniteFieldSolve CPC Library link to program files: https://doi.org/10.17632/ntxvp58mjg.1 Developer's repository link: https://github.com/jfmangan/FiniteFieldSolve Licensing provisions: GPLv3 Programming language: Mathematica, C++ Nature of problem: Exactly solving large linear systems over the rationals occurs in various settings in high-energy theory, for example when performing integral reduction or bootstrapping an amplitude. Solution method: The linear system is solved by repeatedly row reducing over different finite fields (see Ref [1] and references therein). Finite fields avoid the intermediary expression swell inherent to arbitrary precision rationals and bypass roundoff errors from floating point numbers. A downside to using modular arithmetic is that it introduces a tremendous number of integer divisions, but this can be mitigated by compiling the divisions down to simpler instructions. The solver is designed to handle arbitrarily dense systems such as those that appear in certain amplitudes bootstraps. References [1] M. Kauers, Fast solvers for dense linear systems, Nucl. Phys. B Proc. Suppl. 183, 245–250 (2008), 10.1016/j.nuclphysbps.2008.09.111

Requirements Scaling Properties in Large Scale Computing

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

Parallel Fast Algorithm for Large-Scale Electromagnetic Scattering

Parallel Multilevel Fast Multipole Method for Solving Large-Scale Problems

Docker-Enabled Scalable Parallel MLFMA System for RCS Evaluation

Towards a Scalable Hierarchical High-order CFD Solver

An efficient parallel multilevel fast multipole algorithm for large-scale scattering problems

Scalable Hierarchical Parallel Algorithm for the Solution of Super Large-Scale Sparse Linear Equations

A Two-Level Parallel Decomposition Approach for Transient Stability Constrained Optimal Power Flow

Scaling of ab-initio nuclear physics calculations on multicore computer architectures

10-Million Atoms Simulation of First-Principle Package LS3DF

Review Of Large Scale Computing In Electromagnetics With Fast Integral Equation Solvers

Imfssc: An In-Memory Distributed File System Framework For Super Computing

FiniteFieldSolve: Exactly solving large linear systems in high-energy theory

Parallel Multi-Level Fast Multipole Algorithm with Combined Field Integral Equation on Distributed Memory System

A parallel solver framework for fully implicit monolithic fluid-structure interaction

FMM-LU: A Fast Direct Solver for Multiscale Boundary Integral Equations in Three Dimensions

Experience and Analysis of Scalable High-Fidelity Computational Fluid Dynamics on Modular Supercomputing Architectures

Scalable Linear Solvers for Sparse Linear Systems from Large-Scale Numerical Simulations.

LPFSC: A Light Weight Parallel Framework for Super Computing

Performance analysis of electronic structure codes on HPC systems: A case study of SIESTA