Wavelet based Adaptive RBF Method for Nearly Singular Poisson-Type Problems on Irregular Domains

We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs over irregularly shaped domains. For a problem defined over Ω ∈ ℜd, the boundary of an irregularly shaped domain, Γ, is defined as a boundary curve that is a product of a Heaviside function along the normal direction and a piecewise continuous tangential curve. The link between the original wavelet based adaptive method presented in Libre, Emdadi, Kansa, Shekarchi, and Rahimian (2008, 2009) or LEKSR method and the generalized one is given through the use of simple Heaviside masking procedure. In addition level dependent thresholding were introduced to improve the efficiency and convergence rate of the solution. We will show how the generalized wavelet based adaptive method can be applied for detecting nearly singularities in Poisson type PDEs over irregular domains. The numerical examples have illustrated that the proposed method is powerful to analyze the Poisson type PDEs with rapid changes in gradients and nearly singularities

Adaptive Scattered Data Fitting with Tensor Product Spline-Wavelets

The core of the work we present here is an algorithm that constructs a least squares approximation to a given set of unorganized points. The approximation is expressed as a linear combination of particular B-spline wavelets. It implies a multiresolution setting which constructs a hierarchy of approximations to the data with increasing level of detail, proceeding from coarsest to finest scales. It allows for an efficient selection of the degrees of freedom of the problem and avoids the introduction of an artificial uniform grid. In fact, an analysis of the data can be done at each of the scales of the hierarchy, which can be used to select adaptively a set of wavelets that can represent economically the characteristics of the cloud of points in the next level of detail. The data adaption of our method is twofold, as it takes into account both horizontal distribution and vertical irregularities of data. This strategy can lead to a striking reduction of the problem complexity. Furthermore, among the possible ways to achieve a multiscale formulation, the wavelet approach shows additional advantages, based on good conditioning properties and level-wise orthogonality. We exploit these features to enhance the efficiency of iterative solution methods for the system of normal equations of the problem. The combination of multiresolution adaptivity with the numerical properties of the wavelet basis gives rise to an algorithm well suited to cope with problems requiring fast solution methods. We illustrate this by means of numerical experiments that compare the performance of the method on various data sets working with different multi-resolution bases. Afterwards, we use the equivalence relation between wavelets and Besov spaces to formulate the problem of data fitting with regularization. We find that the multiscale formulation allows for a flexible and efficient treatment of some aspects of this problem. Moreover, we study the problem known as robust fitting, in which the data is assumed to be corrupted by wrong measurements or outliers. We compare classical methods based on re-weighting of residuals to our setting in which the wavelet representation of the data computed by our algorithm is used to locate the outliers. As a final application that couples two of the main applications of wavelets (data analysis and operator equations), we propose the use of this least squares data fitting method to evaluate the non-linear term in the wavelet-Galerkin formulation of non-linear PDE problems. At the end of this thesis we discuss efficient implementation issues, with a special interest in the interplay between solution methods and data structures

Multiresolution Analysis on Irregular Surface Meshes

International audienceWavelet-based methods have proven their efficiency for the visu-alization at different levels of detail, progressive transmission, and compression of large data sets. The required core of all wavelet-based methods is a hierarchy of meshes that satisfies subdivision-connectivity: this hierarchy has to be the result of a subdivision process starting from a base mesh. Examples include quadtree uniform 2D meshes, octree uniform 3D meshes, or 4-to-1 split triangular meshes. In particular, the necessity of subdivision-connectivity prevents the application of wavelet-based methods on irregular triangular meshes. In this paper a " wavelet-like " decomposition is introduced , that works on piecewise constant data sets over irregular triangular surface meshes. The decomposition/reconstruction algorithms are based on an extension of wavelet-theory allowing hierarchical meshes without subdivision-connectivity property. Among others, this approach has the following features: it allows exact reconstruction of the data set, even for non-regular triangulations, it extents previous results on Haar-wavelets over 4-to-1 split triangulations