SINGINT: Automatic numerical integration of singular integrands Original Research Article

We explore the combination of deterministic and Monte Carlo methods to facilitate efficient automatic numerical computation of multidimensional integrals with singular integrands. Two adaptive algorithms are presented that employ recursion and are runtime and memory optimized, respectively. SINGINT, a C implementation of the algorithms, is introduced and its utilization in the calculation of particle scattering amplitudes is exemplified.