Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
10346879 | Computers & Mathematics with Applications | 2005 | 7 Pages |
Abstract
where the domain D is a cube containing the point of singularity (0,0,0) and v is a Câ function defined on â3. The procedure employed here is a generalization to 3-D of the method of central corrections for logarithmic singularities [1] in one dimension, and [2] in two dimensions. As in one and two dimensions, the correction coefficients for high-order trapezoidal rules for J(v) are independent of the number of sampling points used to discretize the cube D. When v is compactly supported in D, the approximation is the trapezoidal rule plus a local weighted sum of the values of v around the point of singularity. These quadrature rules provide an efficient, stable and accurate way of approximating J(v). We demonstrate the performance of these quadratures of orders up to 17 for highly oscillatory functions v. These type of integrals appear in scattering calculations in 3-D.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science (General)
Authors
J.C. Aguilar, Yu Chen,