Anti-Gauss cubature rules with applications to Fredholm integral equations on the square
Patricia Diaz de Alba, Luisa Fermo, Giuseppe Rodriguez
TL;DR
The paper addresses accurate numerical integration on the square $\mathcal{S}$ for functions with endpoint algebraic singularities by factoring the integrand with a Jacobi-type weight and introduces anti-Gauss cubature rules as a direct multivariate construction. It develops a weighted Nyström method built from Gauss and anti-Gauss cubatures, proves that the two Nyström interpolants bracket the true solution and shows that their average improves accuracy, with convergence tied to the best polynomial approximation in weighted spaces. The authors formulate efficient linear solvers for the resulting dense systems, including GMRES-FM for general tensors and GMRES-SK for separable kernels, and establish stability with discretization-size-independent condition numbers. Numerical experiments confirm the opposite-sign errors of Gauss and anti-Gauss rules and demonstrate substantial accuracy gains and computational savings from the averaged Nyström interpolant, including favorable performance for large, low-regularity problems and separable kernels. These results provide a robust framework for solving weighted Fredholm integral equations on the square with boundary singularities, with clear avenues for extension to additional averaged cubature schemes.
Abstract
The purpose of this paper is to develop the anti-Gauss cubature rule for approximating integrals defined on the square whose integrand function may have algebraic singularities at the boundaries. An application of such a rule to the numerical solution of Fredholm integral equations of the second-kind is also explored. The stability, convergence, and conditioning of the proposed Nyström-type method are studied. The numerical solution of the resulting dense linear system is also investigated and several numerical tests are presented.
