Multiscale computation with interpolating wavelets

From MaRDI portal



Abstract: Multiresolution analyses based upon interpolets, interpolating scaling functions introduced by Deslauriers and Dubuc, are particularly well-suited to physical applications because they allow exact recovery of the multiresolution representation of a function from its sample values on a finite set of points in space. We present a detailed study of the application of wavelet concepts to physical problems expressed in such bases. The manuscript describes algorithms for the associated transforms which, for properly constructed grids of variable resolution, compute correctly without having to introduce extra grid points. We demonstrate that for the application of local homogeneous operators in such bases, the non-standard multiply of Beylkin, Coifman and Rokhlin also proceeds exactly for inhomogeneous grids of appropriate form. To obtain less stringent conditions on the grids, we generalize the non-standard multiply so that communication may proceed between non-adjacent levels. The manuscript concludes with timing comparisons against naive algorithms and an illustration of the scale-independence of the convergence rate of the conjugate gradient solution of Poisson's equation using a simple preconditioning, suggesting that this approach leads to an O(n) solution of this equation.


The authors study the application of interpolatory scaling functions/wavelets to physical problems expressed in such bases. They describe how nonuniform bases can be conceptualized in the framework of interpolet expansions and use their results to develop algorithms for the associated transforms which for properly constructed grids compute correctly without having to introduce extra grid points. In particular, they focus on the nonstandard form of operators introduced by \textit{G. Beylkin}, \textit{R. Coifman} and \textit{V. Rokhlin} [Commun. Pure Appl. Math. 44, No. 2, 141-183 (1991; Zbl 0722.65022)]. As detailed example (including presentation of implementation features and numerical computations) they consider the solution of Poisson's equation for the full three-dimensional electrostatic potential arising from the nuclei in the nitrogen molecule.











This page was built for publication: Multiscale computation with interpolating wavelets

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