Computational-Complexity on Numerical Integrals

XH WANG,DF HAN

1991-01-01

Abstract:This paper studies the computational complexity on numerical integrals initiated by S. Smale. For a general quadrature functional with k - 1 algebraic precision, the average error of a complexification quadrature formula in the Sobolev space H(k) is obtained. The approach applied in this paper may be also available to other linear approximate problems.

What problem does this paper attempt to address?

The Iyengar Inequality, Revisited

X. H. Wang,C. C. Xie,S. J. Yang

DOI: https://doi.org/10.1007/s10543-007-0153-9

2007-01-01

Abstract:Detailed analysis shows that the famous Iyengar inequality actually says that the Trapezoidal formula is a central algorithm for approximating integrals over an appropriate interval for the class of functions whose derivatives are bounded by a positive number K in L ∞ -sense. The inherent nonlinearity from central algorithms reflects the importance of the Iyengar inequality and thus makes familiar linear methods malfunction when one tries to generalize it. It is shown that the generalization depends on a nonlinear system of equations satisfied by a set of free nodes of a perfect spline. Explicit constructions are obtained in the spirit of the Iyengar inequality for the class of functions whose r th ( r ≤4) derivatives are bounded by a positive number K in L ∞ -sense because a closed solution to the nonlinear system is only available for r ≤4. Connections with computational mathematics, especially with best interpolation and best quadrature, are discussed. Numerical experiments are also included.
Error asymptotic expansion for the numerical approximation of logarithmic-kernel integral equations on closed curves

Hu Li,Jin Huang

DOI: https://doi.org/10.1007/s40314-023-02361-3

2023-06-25

Computational and Applied Mathematics

Abstract:We present a modified quadrature method to solve the logarithmic-kernel integral equations, because of its importance in mathematical physics. At first, the quadrature method is constructed by the trapezoidal rule, which is used to discrete the integral equation. In addition, when the integral curve can be represented by a periodic function, the error asymptotic expansion with single parameter will be proved, which shows that the quadrature method has high accuracy order , and can improve the accuracy of numerical solutions to order accuracy by extrapolation algorithm once. Three numerical examples show our theoretical analysis.

mathematics, applied
Numerical computation of real or complex elliptic integrals

B. C. Carlson

DOI: https://doi.org/10.1007/bf02198293

IF: 2.37

1995-03-01

Numerical Algorithms

Abstract:Algorithms for numerical computation of symmetric elliptic integrals of all three kinds are improved in several ways and extended to complex values of the variables (with some restrictions in the case of the integral of the third kind). Numerical check values, consistency checks, and relations to Legendre's integrals and Bulirsch's integrals are included.

mathematics, applied
Some Simpson- and Ostrowski-Type Integral Inequalities for Generalized Convex Functions in Multiplicative Calculus with Their Computational Analysis

Xinlin Zhan,Abdul Mateen,Muhammad Toseef,Muhammad Aamir Aamir Ali

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

IF: 2.4

2024-06-02

Mathematics

Abstract:Integral inequalities are very useful in finding the error bounds for numerical integration formulas. In this paper, we prove some multiplicative integral inequalities for first-time differentiable s-convex functions. These new inequalities help in finding the error bounds for different numerical integration formulas in multiplicative calculus. The use of s-convex function extends the results for convex functions and covers a large class of functions, which is the main motivation for using s-convexity. To prove the inequalities, we derive two different integral identities for multiplicative differentiable functions in the setting of multiplicative calculus. Then, with the help of these integral identities, we prove some integral inequalities of the Simpson and Ostrowski types for multiplicative generalized convex functions. Moreover, we provide some numerical examples and computational analysis of these newly established inequalities, to show the validity of the results for multiplicative s-convex functions. We also give some applications to quadrature formula and special means of real numbers within the framework of multiplicative calculus.

mathematics
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
Singularity Swapping Method for Nearly Singular Integrals Based on Trapezoidal Rule

Gang Bao,Wenmao Hua,Jun Lai,Jinrui Zhang

DOI: https://doi.org/10.1137/23m1571666

IF: 3.039

2024-04-09

SIAM Journal on Numerical Analysis

Abstract:SIAM Journal on Numerical Analysis, Volume 62, Issue 2, Page 974-997, April 2024. Accurate evaluation of nearly singular integrals plays an important role in many boundary integral equation based numerical methods. In this paper, we propose a variant of singularity swapping method to accurately evaluate the layer potentials for arbitrarily close targets. Our method is based on the global trapezoidal rule and trigonometric interpolation, resulting in an explicit quadrature formula. The method achieves spectral accuracy for nearly singular integrals on closed analytic curves. In order to extract the singularity from the complexified distance function, an efficient root finding method is proposed based on contour integration. Through the change of variables, we also extend the quadrature method to integrals on the piecewise analytic curves. Numerical examples for Laplace and Helmholtz equations show that high-order accuracy can be achieved for arbitrarily close field evaluation.

mathematics, applied
Improving the Accuracy of Exponentially Converging Quadratures

Khokhlachev, V. S.

DOI: https://doi.org/10.1134/s0965542524010020

2024-03-21

Computational Mathematics and Mathematical Physics

Abstract:Evaluation of one-dimensional integrals arises in many problems in physics and technology. This is most often done using simple quadratures of midpoints, trapezoids and Simpson on a uniform grid. For integrals of periodic functions over the full period, the convergence of these quadratures drastically accelerates and depends on the number of grid steps according to an exponential law. In this paper, new asymptotically accurate estimates of the error of such quadratures are obtained. They take into account the location and multiplicity of the poles of the integrand in the complex plane. A generalization of these estimates is constructed for the case when there is no a priori information about the poles of the integrand. An error extrapolation procedure is described that drastically accelerates the convergence of quadratures.

mathematics, applied,physics, mathematical
Numerical Algorithms for Approximation of Fractional Integrals and Derivatives Based on Quintic Spline Interpolation

Mariusz Ciesielski

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

2024-02-18

Symmetry

Abstract:Numerical algorithms for calculating the left- and right-sided Riemann–Liouville fractional integrals and the left- and right-sided fractional derivatives in the Caputo sense using spline interpolation techniques are derived. The spline of the fifth degree (the so-called quintic spline) is mainly taken into account, but the linear and cubic splines are also considered to compare the quality of the developed method and numerical calculations. The estimation of errors for the derived approximation algorithms is presented. Examples of the numerical evaluation of the fractional integrals and derivatives are executed using 128-bit floating-point numbers and arithmetic routines. For each derived algorithm, the experimental orders of convergence are calculated. Also, an illustrative computational example showing the action of the considered fractional operators on the symmetric function in the interval is presented.

multidisciplinary sciences
On Newton–Cotes Formula-Type Inequalities for Multiplicative Generalized Convex Functions via Riemann–Liouville Fractional Integrals with Applications to Quadrature Formulas and Computational Analysis

Abdul Mateen,Serap Özcan,Zhiyue Zhang,Bandar Bin-Mohsin

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

IF: 3.577

2024-09-19

Fractal and Fractional

Abstract:In this article, we develop multiplicative fractional versions of Simpson's and Newton's formula-type inequalities for differentiable generalized convex functions with the help of established identities. The main motivation for using generalized convex functions lies in their ability to extend results beyond traditional convex functions, encompassing a broader class of functions, and providing optimal approximations for both lower and upper bounds. These inequalities are very useful in finding the error bounds for the numerical integration formulas in multiplicative calculus. Applying these results to the Quadrature formulas demonstrates their practical utility in numerical integration. Furthermore, numerical analysis provides empirical evidence of the effectiveness of the derived findings. It is also demonstrated that the newly proven inequalities extend certain existing results in the literature.

mathematics, interdisciplinary applications
Numerical calculation method for function integration on submanifolds of $\mathbb{R}^n$ or compact Riemannian manifolds

Fusheng Deng,Gang Huang,Yingyi Wu

2024-09-21

Abstract:In this paper, we present a method for digitally representing the "volume element" and calculating the integral of a function on compact hypersurfaces with or without boundary, and low-dimensional submanifolds in $\mathbb{R}^n$. We also extend such calculation to hypersurfaces in compact Riemannnian manifolds.

Numerical Analysis,Differential Geometry
Numerical analysis of a nonsmooth quasilinear elliptic control problem: II. Finite element discretization and error estimates

Christian Clason,Vu Huu Nhu,Arnd Rösch

2024-02-23

Abstract:In this paper, we carry out the numerical analysis of a nonsmooth quasilinear elliptic optimal control problem, where the coefficient in the divergence term of the corresponding state equation is not differentiable with respect to the state variable. Despite the lack of differentiability of the nonlinearity in the quasilinear elliptic equation, the corresponding control-to-state operator is of class $C^1$ but not of class $C^2$. Analogously, the discrete control-to-state operators associated with the approximated control problems are proven to be of class $C^1$ only. By using an explicit second-order sufficient optimality condition, we prove a priori error estimates for a variational approximation, a piecewise constant approximation, and a continuous piecewise linear approximation of the continuous optimal control problem. The numerical tests confirm these error estimates.

Numerical Analysis,Optimization and Control
On quadrature for singular integral operators with complex symmetric quadratic forms

Jeremy Hoskins,Manas Rachh,Bowei Wu

2023-12-12

Abstract:This paper describes a trapezoidal quadrature method for the discretization of weakly singular, singular and hypersingular boundary integral operators with complex symmetric quadratic forms. Such integral operators naturally arise when complex coordinate methods or complexified contour methods are used for the solution of time-harmonic acoustic and electromagnetic interface problems in three dimensions. The quadrature is an extension of a locally corrected punctured trapezoidal rule in parameter space wherein the correction weights are determined by fitting moments of error in the punctured trapezoidal rule, which is known analytically in terms of the Epstein zeta function. In this work, we analyze the analytic continuation of the Epstein zeta function and the generalized Wigner limits to complex quadratic forms; this analysis is essential to apply the fitting procedure for computing the correction weights. We illustrate the high-order convergence of this approach through several numerical examples.

Numerical Analysis
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
Computing singular and near-singular integrals over curved boundary elements: The strongly singular case

Hadrien Montanelli,Francis Collino,Houssem Haddar

2024-06-21

Abstract:We present algorithms for computing strongly singular and near-singular surface integrals over curved triangular patches, based on singularity subtraction, the continuation approach, and transplanted Gauss quadrature. We demonstrate the accuracy and robustness of our method for quadratic basis functions and quadratic triangles by integrating it into a boundary element code and solving several scattering problems in 3D. We also give numerical evidence that the utilization of curved boundary elements enhances computational efficiency compared to conventional planar elements.

Numerical Analysis
High-order numerical integration on self-affine sets

Patrick Joly,Maryna Kachanovska,Zoïs Moitier

2024-10-01

Abstract:We construct an interpolatory high-order cubature rule to compute integrals of smooth functions over self-affine sets with respect to an invariant measure. The main difficulty is the computation of the cubature weights, which we characterize algebraically, by exploiting a self-similarity property of the integral. We propose an $h$-version and a $p$-version of the cubature, present an error analysis and conduct numerical experiments.

Numerical Analysis
On the numerical solution of nonlinear integral equations on non-rectangular domains utilizing thin plate spline collocation method

Pouria Assari,Mehdi Dehghan,Pouria Assari,Mehdi Dehghan

DOI: https://doi.org/10.1007/s12044-019-0511-y

2019-08-09

Abstract:This article investigates an approximate scheme to solve nonlinear Fredholm integral equations of the second kind on non-rectangular domains. The integral equations considered in the current paper are considered together with either smooth or weakly singular kernels. The offered method utilizes thin plate splines as a basis in the discrete collocation method. We can regard thin plate splines as a type of the free shape parameter radial basis functions. These basis functions establish an accurate and stable technique to estimate an unknown function by using a set of scattered points on the solution domains. Since the thin plate splines have limited smoothness, the integrals appeared in the scheme cannot be estimated by classical integration rules. Therefore, we introduce a special precise quadrature formula on non-rectangular domains to compute these integrals. The proposed scheme does not require any mesh generations, so it is meshless and does not depend on the domain form. Error analysis is also provided for the method. The performance and convergence of the new approach are tested on four two-dimensional integral equations given on the wing, mushroom, pentagon and fish-like domains.
Two Numerical Schemes Via an Auxiliary Function for Nonlocal Problems With Integrable Kernels

Yutong Li,Xiaobo Yin

DOI: https://doi.org/10.1002/num.23153

2024-10-11

Numerical Methods for Partial Differential Equations

Abstract:In this study, we explore volume constrained nonlocal problem with integrable kernels. A central challenge addressed is the discontinuity of true solutions across the boundary, posing difficulties for designing efficient numerical methods. To address this issue, we introduce an auxiliary function which serves as the solution of an auxiliary problem. This approach effectively transforms the original problem to a new nonlocal problem with a true solution which is continuous across the boundary. We propose two numerical schemes to solve the new nonlocal problem: one is exclusively applicable to one‐dimensional case, while the other is suitable for general dimensional cases. Experimental results show that both schemes exhibit optimal convergence order under mild condition.

mathematics, applied
On quadrature of highly oscillatory Bessel function via asymptotic analysis of simplex integrals

Yongxiong Zhou,Ruyun Chen

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

IF: 2.872

2024-08-31

Journal of Computational and Applied Mathematics

Abstract:In this article, two methods for evaluating highly oscillatory Bessel integrals are explored. Firstly, a polynomial is analyzed as an effective approximation of the simplex integral of a highly oscillatory Bessel function based on Laplace transform, and its error rapidly decreases as the frequency increases. Furthermore, the inner product of f and highly oscillatory Bessel function can be approximated by two other forms of inner product by which one depends on a polynomial and the higher derivatives of f , another depends on Bessel function and the interpolation polynomial of f . In addition, three issues related to highly oscillatory Bessel integrals have also been discussed: inequalities for the convergence rate of Filon-type methods, evaluation of Cauchy principal values, and simplified evaluation on infinite intervals. Through some preliminary numerical experiments, our theoretical analysis has been preliminarily confirmed, and the proposed numerical method is accurate and effective.

mathematics, applied
Numerical approximations for a hyperbolic integrodifferential equation with a non-positive variable-sign kernel and nonlinear-nonlocal damping

Wenlin Qiu,Xiangcheng Zheng,Kassem Mustapha

2024-12-10

Abstract:This work considers the Galerkin approximation and analysis for a hyperbolic integrodifferential equation, where the non-positive variable-sign kernel and nonlinear-nonlocal damping with both the weak and viscous damping effects are involved. We derive the long-time stability of the solution and its finite-time uniqueness. For the semi-discrete-in-space Galerkin scheme, we derive the long-time stability of the semi-discrete numerical solution and its finite-time error estimate by technical splitting of intricate terms. Then we further apply the centering difference method and the interpolating quadrature to construct a fully discrete Galerkin scheme and prove the long-time stability of the numerical solution and its finite-time error estimate by designing a new semi-norm. Numerical experiments are performed to verify the theoretical findings.

Numerical Analysis

Computational-Complexity on Numerical Integrals

The Iyengar Inequality, Revisited

Error asymptotic expansion for the numerical approximation of logarithmic-kernel integral equations on closed curves

Numerical computation of real or complex elliptic integrals

Some Simpson- and Ostrowski-Type Integral Inequalities for Generalized Convex Functions in Multiplicative Calculus with Their Computational Analysis

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

Singularity Swapping Method for Nearly Singular Integrals Based on Trapezoidal Rule

Improving the Accuracy of Exponentially Converging Quadratures

Numerical Algorithms for Approximation of Fractional Integrals and Derivatives Based on Quintic Spline Interpolation

On Newton–Cotes Formula-Type Inequalities for Multiplicative Generalized Convex Functions via Riemann–Liouville Fractional Integrals with Applications to Quadrature Formulas and Computational Analysis

Numerical calculation method for function integration on submanifolds of $\mathbb{R}^n$ or compact Riemannian manifolds

Numerical analysis of a nonsmooth quasilinear elliptic control problem: II. Finite element discretization and error estimates

On quadrature for singular integral operators with complex symmetric quadratic forms

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

Computing singular and near-singular integrals over curved boundary elements: The strongly singular case

High-order numerical integration on self-affine sets

On the numerical solution of nonlinear integral equations on non-rectangular domains utilizing thin plate spline collocation method

Two Numerical Schemes Via an Auxiliary Function for Nonlocal Problems With Integrable Kernels

On quadrature of highly oscillatory Bessel function via asymptotic analysis of simplex integrals

Numerical approximations for a hyperbolic integrodifferential equation with a non-positive variable-sign kernel and nonlinear-nonlocal damping