High-precision quadrature via local Fourier extension: analytic integration, uniform sampling, and correction for piecewise smooth integrands

Abstract: We propose a high-precision numerical quadrature framework based on local Fourier extension (LFE) approximations. The method constructs, on each subinterval, a truncated-SVD stabilized local Fourier continuation of the integrand on an extended periodic domain, and then evaluates the integral \emph{analytically} from the resulting Fourier coefficients. Under uniform sampling, the discrete LFE matrix and its TSVD factors are precomputed once and reused across all windows, yielding an efficient offline/online implementation that remains compatible with classical composite rules. We provide an error bound that reduces the quadrature error to the LFE approximation error and derive algebraic convergence rates for Sobolev-regular integrands. Numerical experiments demonstrate that, on smooth functions, the proposed quadrature reaches near machine precision with substantially fewer nodes than the composite Simpson rule. The advantage persists for oscillatory and variable-frequency integrands and becomes more pronounced for nonuniform phase structures. For continuous piecewise smooth integrands, we develop a correction strategy driven by coefficient-energy outliers to identify singularity-containing windows, followed by a localized procedure that brackets the singular point within one grid cell and corrects only the affected window contribution. The corrected quadrature restores near-spectral accuracy in the reported tests, including cases where the singularity is not aligned with the window endpoints.