Gauss quadrature for integrals and sums

A. D. Alhaidari

DOI: https://doi.org/10.51483/IJPAMR.3.1.2023.1-10

2023-05-29

Abstract:Gauss quadrature integral approximation is extended to include integrals with a measure consisting of continuous as well as discrete components. That is, we give an approximation for the integral of a function plus its sum over a discrete weighted set.

Numerical Analysis

What problem does this paper attempt to address?

On Lebesgue Integral Quadrature

Vladislav Gennadievich Malyshkin

DOI: https://doi.org/10.2139/ssrn.3229363

2020-02-25

Abstract:A new type of quadrature is developed. The Gaussian quadrature, for a given measure, finds optimal values of a function's argument (nodes) and the corresponding weights. In contrast, the Lebesgue quadrature developed in this paper, finds optimal values of function (value-nodes) and the corresponding weights. The Gaussian quadrature groups sums by function argument; it can be viewed as a $n$-point discrete measure, producing the Riemann integral. The Lebesgue quadrature groups sums by function value; it can be viewed as a $n$-point discrete distribution, producing the Lebesgue integral. Mathematically, the problem is reduced to a generalized eigenvalue problem: Lebesgue quadrature value-nodes are the eigenvalues and the corresponding weights are the square of the averaged eigenvectors. A numerical estimation of an integral as the Lebesgue integral is especially advantageous when analyzing irregular and stochastic processes. The approach separates the outcome (value-nodes) and the probability of the outcome (weight). For this reason, it is especially well-suited for the study of non-Gaussian processes. The software implementing the theory is available from the authors.

Numerical Analysis,Machine Learning
A Note on the Gauss-Jacobi Quadrature Formulae for Singular Integral Equations of the Second Kind

Lifeng Ma,Alexander M. Korsunsky

DOI: https://doi.org/10.1023/b:frac.0000031158.62052.3c

IF: 2.635

2004-01-01

International Journal of Fracture

Abstract:A fast and efficient numerical method based on the Gauss-Jacobi quadrature is described that is suitable for solving Fredholm singular integral equations of the second kind that are frequently encountered in fracture and contact mechanics. Here we concentrate on the case when the unknown function is singular at both ends of the interval. Quadrature formulae involve fixed nodal points and provide exact results for polynomials of degree 2n − 1, where n is the number of nodes. Finally, an application of the method to a plane problem involving complete contact is presented.
Optimal numerical integration and approximation of functions on $\mathbb{R}^d$ equipped with Gaussian measure

Dinh Dũng,Van Kien Nguyen

2023-06-17

Abstract:We investigate the numerical approximation of integrals over $\mathbb{R}^d$ equipped with the standard Gaussian measure $\gamma$ for integrands belonging to the Gaussian-weighted Sobolev spaces $W^\alpha_p(\mathbb{R}^d, \gamma)$ of mixed smoothness $\alpha \in \mathbb{N}$ for $1 < p < \infty$. We prove the asymptotic order of the convergence of optimal quadratures based on $n$ integration nodes and propose a novel method for constructing asymptotically optimal quadratures. As for related problems, we establish by a similar technique the asymptotic order of the linear, Kolmogorov and sampling $n$-widths in the Gaussian-weighted space $L_q(\mathbb{R}^d, \gamma)$ of the unit ball of $W^\alpha_p(\mathbb{R}^d, \gamma)$ for $1 \leq q < p < \infty$ and $q=p=2$.

Numerical Analysis
A note on the accuracy of adaptive Gauss–Hermite quadrature

Shaobo Jin,Björn Andersson

DOI: https://doi.org/10.1093/biomet/asz080

IF: 3.0279

2020-02-10

Biometrika

Abstract:Summary Numerical quadrature methods are needed for many models in order to approximate integrals in the likelihood function. In this note, we correct the error rate given by Liu & Pierce (1994) for integrals approximated with adaptive Gauss–Hermite quadrature and show that the approximation is less accurate than previously thought. We discuss the relationship between the error rates of adaptive Gauss–Hermite quadrature and Laplace approximation, and provide a theoretical explanation of simulation results obtained in previous studies regarding the accuracy of adaptive Gauss–Hermite quadrature.

statistics & probability,mathematical & computational biology,biology
Monotonicity, bounds and extrapolation of Block-Gauss and Gauss-Radau quadrature for computing $B^T ϕ(A) B$

Jörn Zimmerling,Vladimir Druskin,Valeria Simoncini

2024-08-11

Abstract:In this paper, we explore quadratures for the evaluation of $B^T \phi(A) B$ where $A$ is a symmetric nonnegative-definite matrix in $\mathbb{R}^{n \times n}$, $B$ is a tall matrix in $\mathbb{R}^{n \times p}$, and $\phi(\cdot)$ represents a matrix function that is regular enough in the neighborhood of $A$'s spectrum, e.g., a Stieltjes or exponential function. These formulations, for example, commonly arise in the computation of multiple-input multiple-output (MIMO) transfer functions for diffusion PDEs. We propose an approximation scheme for $B^T \phi(A) B$ leveraging the block Lanczos algorithm and its equivalent representation through Stieltjes matrix continued fractions. We extend the notion of Gauss-Radau quadrature to the block case, facilitating the derivation of easily computable error bounds. For problems stemming from the discretization of self-adjoint operators with a continuous spectrum, we obtain sharp estimates grounded in potential theory for Padé approximations and justify extrapolation algorithms at no added computational cost. The obtained results are illustrated on large-scale examples of 2D diffusion and 3D Maxwell's equations as well as a graph from the SNAP repository. We also present promising experimental results on convergence acceleration via random enrichment of the initial block $B$.

Numerical Analysis
Evaluation of the Gauss integral

Dmitri Martila,Stefan Groote

DOI: https://doi.org/10.3390/stats5020032

2022-02-22

Abstract:The normal or Gaussian distribution plays a prominent role in almost all fields of science. However, it is well known that the Gauss (or Euler--Poisson) integral over a finite boundary, as it is necessary for instance for the error function or the cumulative distribution of the normal distribution, cannot be expressed by analytic functions. This is proven by the Risch algorithm. Still, there are proposals for approximate solutions. In this paper, we give a new solution in terms of normal distributions by applying a geometric procedure iteratively to the problem.

Data Analysis; Statistics and Probability,Probability,Numerical Analysis
Gaussian quadrature rules for composite highly oscillatory integrals

Menghan Wu,Haiyong Wang

2023-06-08

Abstract:Highly oscillatory integrals of composite type arise in electronic engineering and their calculations is a challenging problem. In this paper, we propose two Gaussian quadrature rules for computing such integrals. The first one is constructed based on the classical theory of orthogonal polynomials and its nodes and weights can be computed efficiently by using tools of numerical linear algebra. We show that the rate of convergence of this rule depends solely on the regularity of the non-oscillatory part of the integrand. The second one is constructed with respect to a sign-changing function and the classical theory of Gaussian quadrature can not be used anymore. We explore theoretical properties of this Gaussian quadrature, including the trajectories of the quadrature nodes and the convergence rate of these nodes to the endpoints of the integration interval, and prove its asymptotic error estimate under suitable hypotheses. Numerical experiments are presented to demonstrate the performance of the proposed methods.

Numerical Analysis
Gauss $q$-ed from Heine cubed

P.L. Robinson

DOI: https://doi.org/10.48550/arXiv.2201.07049

2022-01-18

Abstract:We consider $q$-analytic derivations of the $q$-Gauss summation formula for a $\, _2\phi _1$ that respect the symmetry in its upper parameters.

Classical Analysis and ODEs
Equivariant Gauss sum of finite quadratic forms

Shouhei Ma

DOI: https://doi.org/10.48550/arXiv.1703.07504

2017-03-22

Abstract:The classical quadratic Gauss sum can be thought of as an exponential sum attached to a quadratic form on a cyclic group. We introduce an equivariant version of Gauss sum for arbitrary finite quadratic forms, which is an exponential sum twisted by the action of the orthogonal group. We prove that simple arithmetic formulae hold for some basic classes of quadratic forms. In application, such invariant appears in the dimension formula for certain vector-valued modular forms.

Number Theory
Gaussian integrals depending by a quantum parameter in finite dimension

Simone Camosso

DOI: https://doi.org/10.57016/TM-ZOVB7625

2021-07-15

Abstract:A common theme in mathematics is the evaluation of Gauss integrals. This, coupled with the fact that they are used in different branches of science, makes the topic always actual and interesting. In these notes we shall analyze a particular class of Gaussian integrals that depends by the quantum parameter $\hbar$. Starting from classical results, we will present an overview on methods, examples and analogies regarding the practice of solving quantum Gaussian integrals.

General Mathematics
On Weighted Generalized Gauss Quadratures for Müntz Systems

Huaijin Wang,Chuanju Xu

DOI: https://doi.org/10.48550/arXiv.2310.13534

2023-10-20

Abstract:A novel recurrence formula for moments with respect to Müntz-Legendre polynomials is proposed and applied to construct a numerical method for solving generalized Gauss quadratures with power function weight for Müntz systems. These quadrature rules exhibit several properties similar to the classical Gaussian quadratures for polynomial systems, including positive weights, rapid convergence, and others. They are applicable to a wide range of functions, including smooth functions and functions with endpoint singularities, commonly found in integral equations with singular kernels, complex analysis, potential theory, and other areas.

Numerical Analysis
Spherical configurations and quadrature methods for integral equations of the second kind

Congpei An,Hao-Ning Wu

2024-08-27

Abstract:In this paper, we propose and analyze a product integration method for the second-kind integral equation with weakly singular and continuous kernels on the unit sphere $\mathbb{S}^2$. We employ quadrature rules that satisfy the Marcinkiewicz--Zygmund property to construct hyperinterpolation for approximating the product of the continuous kernel and the solution, in terms of spherical harmonics. By leveraging this property, we significantly expand the family of candidate quadrature rules and establish a connection between the geometrical information of the quadrature points and the error analysis of the method. We then utilize product integral rules to evaluate the singular integral with the integrand being the product of the singular kernel and each spherical harmonic. We derive a practical $L^{\infty}$ error bound, which consists of two terms: one controlled by the best approximation of the product of the continuous kernel and the solution, and the other characterized by the Marcinkiewicz--Zygmund property and the best approximation polynomial of this product. Numerical examples validate our numerical analysis.

Numerical Analysis
Gauss Relations in Feynman Integrals

Tai-Fu Feng,Yang Zhou,Hai-Bin Zhang

2024-07-15

Abstract:Embedding Feynman integrals in Grassmannians, we express Feynman integrals as linear combinations of generalized hypergeometric functions. Here we present general methods to obtain Gauss relations among those generalized hypergeometric functions. The hypergeometric expressions of Feynman integral are analytically continued from some connected component to another by the Gauss inverse relations, then continued to the whole domain of definition by the Gauss-Kummer relations. The Laurant series of the Feynman integral around the time-space dimension $D=4$ is obtained by the Gauss adjacent relations where the coefficients of powers of $D-4$ are given as some finite linear combinations of hypergeometric functions with integer parameters. As an example, we illustrate how to use the method to obtain the analytic expression of the Feynman integral of one-loop self energy in its whole domain of definition.

High Energy Physics - Theory
Error bounds for Gauss–Lobatto quadrature of analytic functions on an ellipse

Hiroshi Sugiura,Takemitsu Hasegawa

DOI: https://doi.org/10.1016/j.cam.2024.116326

IF: 2.872

2024-10-12

Journal of Computational and Applied Mathematics

Abstract:For the ( n+2 )-point Gauss–Jacobi–Lobatto quadrature to integrals with the Jacobi weight function (1−t)α(1+t)β ( α>−1 , β>−1 ) over the interval [−1,1] , we estimate the location where the kernel of the error functional for functions analytic on an ellipse and its interior in the complex plane attains its maximum modulus. As in our previous work on the Gauss–Jacobi rule, when α≠β , the location is the intersection point of the ellipse with the real axis in the complex plane. When α=β (the Gegenbauer weight), it is the intersection points with the real axis for −1<α≤0 or for 0<α≤52 and 1≤n≤n+(α) , and with the imaginary axis for 0<α≤52 and n>n+(α) or for α>52 . Here, n+(α) ( >0 ) is a monotonously decreasing function for α>0 with n+(52)=1 and limα→+0n+(α)=∞ . Some numerical results are given to confirm the locations.

mathematics, applied
A Golub-Welsch version for simultaneous Gaussian quadrature

Walter Van Assche

DOI: https://doi.org/10.1007/s11075-024-01767-2

2024-02-05

Abstract:The zeros of type II multiple orthogonal polynomials can be used for quadrature formulas that approximate $r$ integrals of the same function $f$ with respect to $r$ measures $\mu_1,\ldots,\mu_r$ in the spirit of Gaussian quadrature. This was first suggested by Borges in 1994, even though he does not mention multiple orthogonality. We give a method to compute the quadrature nodes and the quadrature weights which extends the Golub-Welsch approach using the eigenvalues and left and right eigenvectors of a banded Hessenberg matrix. This method was already described by Coussement and Van Assche in 2005 but it seems to have gone unnoticed. We describe the result in detail for $r=2$ and give some examples.

Numerical Analysis,Classical Analysis and ODEs
Evaluations of a Weighted Average of Gauss Sums

Wen-Kai Shao,Yuan He

DOI: https://doi.org/10.1155/2021/6206329

2021-01-01

Journal of Mathematics

Abstract:In this paper, we perform a further investigation for a weighted average of Gauss sums. By making use of some properties of the cotangent function and the Bernoulli polynomials, we explicitly evaluate the weighted average of Gauss sums in terms of the special values of Dirichlet L -functions at positive integers.
Generalized quadrature for finite temperature Green’s function methods

Jie Gu,Jia Chen,Yang Wang,X.-G. Zhang

DOI: https://doi.org/10.1016/j.cpc.2020.107178

IF: 4.717

2020-08-01

Computer Physics Communications

Abstract:<p>In electronic structure and quantum transport calculations, many physical quantities are integrations over the electron energy, weighted by the Fermi-Dirac distribution function. Green's function based approaches commonly circumvent the numerically difficult real energy integration by extending the integrand analytically into the complex energy plane, and using a Gaussian quadrature integration over a complex energy contour for zero temperature. For finite temperatures, a much slower convergent sum over the Matsubara frequencies is necessary. We present a generalized quadrature method that uses orthogonal polynomials on the manifold of Matsubara frequencies to enable rapid convergence. Both Gaussian quadrature integration and Matsubara frequency summation methods are shown to be limiting cases of the generalized method. Tests on an all-electron <em>ab initio</em> code and total energy calculation of interacting Anderson impurity model show convergence with a small number of energy mesh points.</p>

physics, mathematical,computer science, interdisciplinary applications
Quadrature Methods for Highly Oscillatory Singular Integrals

Jing Gao,Marissa Condon,Arieh Iserles,Benjamin Gilvey,Jon Trevelyan

DOI: https://doi.org/10.4208/jcm.1911-m2019-0044

2020-01-01

Journal of Computational Mathematics

Abstract:We address the evaluation of highly oscillatory integrals, with power-law and logarithmic singularities. Such problems arise in numerical methods in engineering. Notably, the evaluation of oscillatory integrals dominates the run-time for wave-enriched boundary integral formulations for wave scattering, and many of these exhibit singularities. We show that the asymptotic behaviour of the integral depends on the integrand and its derivatives at the singular point of the integrand, the stationary points and the endpoints of the integral. A truncated asymptotic expansion achieves an error that decays faster for increasing frequency. Based on the asymptotic analysis, a Filon-type method is constructed to approximate the integral. Unlike an asymptotic expansion, the Filon method achieves high accuracy for both small and large frequency. Complex-valued quadrature involves interpolation at the zeros of polynomials orthogonal to a complex weight function. Numerical results indicate that the complex-valued Gaussian quadrature achieves the highest accuracy when the three methods are compared. However, while it achieves higher accuracy for the same number of function evaluations, it requires significant additional cost of computation of orthogonal polynomials and their zeros.
High-Order Numerical Integration on Domains Bounded by Intersecting Level Sets

Lauritz Beck,Florian Kummer

2023-08-21

Abstract:We present a high-order method that provides numerical integration on volumes, surfaces, and lines defined implicitly by two smooth intersecting level sets. To approximate the integrals, the method maps quadrature rules defined on hypercubes to the curved domains of the integrals. This enables the numerical integration of a wide range of integrands since integration on hypercubes is a well known problem. The mappings are constructed by treating the isocontours of the level sets as graphs of height functions. Numerical experiments with smooth integrands indicate a high-order of convergence for transformed Gauss quadrature rules on domains defined by polynomial, rational, and trigonometric level sets. We show that the approach we have used can be combined readily with adaptive quadrature methods. Moreover, we apply the approach to numerically integrate on difficult geometries without requiring a low-order fallback method.

Numerical Analysis

Gauss quadrature for integrals and sums

On Lebesgue Integral Quadrature

A Note on the Gauss-Jacobi Quadrature Formulae for Singular Integral Equations of the Second Kind

Optimal numerical integration and approximation of functions on $\mathbb{R}^d$ equipped with Gaussian measure

A note on the accuracy of adaptive Gauss–Hermite quadrature

Monotonicity, bounds and extrapolation of Block-Gauss and Gauss-Radau quadrature for computing $B^T ϕ(A) B$

Evaluation of the Gauss integral

Gaussian quadrature rules for composite highly oscillatory integrals

Gauss $q$-ed from Heine cubed

Equivariant Gauss sum of finite quadratic forms

Gaussian integrals depending by a quantum parameter in finite dimension

On Weighted Generalized Gauss Quadratures for Müntz Systems

Spherical configurations and quadrature methods for integral equations of the second kind

Gauss Relations in Feynman Integrals

Error bounds for Gauss–Lobatto quadrature of analytic functions on an ellipse

A Golub-Welsch version for simultaneous Gaussian quadrature

Evaluations of a Weighted Average of Gauss Sums

Generalized quadrature for finite temperature Green’s function methods

Quadrature Methods for Highly Oscillatory Singular Integrals

High-Order Numerical Integration on Domains Bounded by Intersecting Level Sets