An accurate and quadrature-free evaluation of multipole expansion of functions represented by multiwavelets

arXiv1211.6076MaRDI QIDQ6237490FDOQ6237490

Authors: Jae-Seok Huh 

Publication date: 26 November 2012

Abstract: We present formulas for accurate numerical conversion between functions represented by multiwavelets and their multipole/local expansions with respect to the kernel of the form, $e^{lambdar}/r$. The conversion is essential for the application of fast multipole methods for functions represented by multiwavelets. The corresponding separated kernels exhibit near-singular behaviors at large $lambda$. Moreover, a multiwavelet basis function oscillates more wildly as its degree increases. These characteristics in combination render any brute-force approach based on numerical quadratures impractical. Our approach utilizes the series expansions of the modified spherical Bessel functions and the Cartesian expansions of solid harmonics so that the multipole-multiwavelet conversion matrix can be evaluated like a special function. The result is a quadrature-free, fast, reliable, and machine precision accurate scheme to compute the conversion matrix with predictable sparsity patterns.

This page was built for publication: An accurate and quadrature-free evaluation of multipole expansion of functions represented by multiwavelets

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6237490)