Source author record

Jae-Seok Huh

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA

Catalog footprint

Research overview

2works

1topics

1close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2012arXiv

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

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^{λr}/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 $λ$ . 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.

preprint2012arXiv

Singular quadrature rules and fast convolutions for Fourier spectral methods

We present a generic scheme to construct corrected trapezoidal rules with spectral accuracy for integral operators with weakly singular kernels in arbitrary dimensions. We assume that the kernel factorization of the form, $K=αϕ+\widetilde{K}$ with smooth $α$ and $\widetilde{K}$ , is available so that the operations on the smooth factors can be performed accurately on the basis of standard Fourier spectral methods. To achieve high precision results, our approach utilizes the exact evaluation of the Fourier coefficients of the radial singularity $ϕ$ , which can be obtained in arbitrary dimensions by the singularity isolation/truncation described in this article. We provide a complete set of formulas for singularities of the type: $\log(r)$ and $r^{-ν}$ . Convergence analysis shows that the constructed quadrature rules exhibit almost identical rate of convergence to the trapezoidal rule applied for non-singular integrands. Especially, for smooth data, the corrected trapezoidal rules converge super-algebraically.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.