On the implementation of a modified Sag–Szekeres quadrature method

We describe a modified Sag–Szekeres multidimensional quadrature algorithm and discuss its implementation as a general-purpose library procedure on serial and parallel architectures. Examples illustrate its effectiveness for both smooth and singular integrands.