Accurate numerical methods for micromagnetics simulations with general geometries
Journal of Computational Physics
A wideband fast multipole method for the Helmholtz equation in three dimensions
Journal of Computational Physics
High performance BLAS formulation of the multipole-to-local operator in the fast multipole method
Journal of Computational Physics
Journal of Computational Physics
Computers & Mathematics with Applications
Efficient discretization of Laplace boundary integral equations on polygonal domains
Journal of Computational Physics
A Nonlinear Optimization Procedure for Generalized Gaussian Quadratures
SIAM Journal on Scientific Computing
A new Fast Multipole formulation for the elastodynamic half-space Green's tensor
Journal of Computational Physics
Advances in Computational Mathematics
Hi-index | 0.03 |
We present a nonlinear optimization procedure for the design of generalized Gaussian quadratures for a fairly broad class of functions. While some of the components of the algorithm have been published previously, we introduce an improved procedure for the determination of an acceptable initial point for the continuation scheme that stabilizes the Newton-type process used to find the quadratures. The resulting procedure never failed when applied to Chebyshev systems (for which the existence and uniqueness of generalized Gaussian quadratures are well known); it also worked for many non-Chebyshev systems, for which the generalized Gaussian quadratures are not guaranteed to exist. The performance of the algorithm is illustrated with several numerical examples; some of the presented quadratures integrate efficiently large classes of singular functions.