A Semi-Lagrangian Scheme with Radial Basis Approximation for Surface Reconstruction

We propose a Semi-Lagrangian scheme coupled with Radial Basis Function interpolation for approximating a curvature-related level set model, which has been proposed by Zhao et al. in \cite{ZOMK} to reconstruct unknown surfaces from sparse, possibly noisy data sets. The main advantages of the proposed scheme are the possibility to solve the level set method on unstructured grids, as well as to concentrate the reconstruction points in the neighbourhood of the data set, with a consequent reduction of the computational effort. Moreover, the scheme is explicit. Numerical tests show the accuracy and robustness of our approach to reconstruct curves and surfaces from relatively sparse data sets.Comment: 14 pages, 26 figure

On the role of particles and radial basis functions in a finite element level set method for bubble dynamics

The aim of this presentation is to highlight the role that Particle-based simulations and Radial Basis Functions (RBFs) have played in the development of a computationally efficient, level-set, Finite Element method for the simulation of Newtonian and non-Newtonian interface flows. First, we introduce the mathemat- ical formulation and the interface-capturing technique used in the simulation of multiphase flows, underscoring the influence of marker particles on the enhanced definition of the interface. Then, we explore the effect of adding polymer parti- cles to the domain to perform Brownian Dynamics Simulations of polymer flows. Finally, we leverage RBFs to reconstruct, in an almost free-independent way the polymer stress tensor retrieved from the polymer particles. Numerical simulations of pure advection flows and bubble dynamics simulations of complex flows on two-dimensional configurations emphasize the improvements offered by this hybrid, Finite Element/RBF/Particle-based method

Rapid evaluation of radial basis functions

Over the past decade, the radial basis function method has been shown to produce high quality solutions to the multivariate scattered data interpolation problem. However, this method has been associated with very high computational cost, as compared to alternative methods such as finite element or multivariate spline interpolation. For example. the direct evaluation at M locations of a radial basis function interpolant with N centres requires O(M N) floating-point operations. In this paper we introduce a fast evaluation method based on the Fast Gauss Transform and suitable quadrature rules. This method has been applied to the Hardy multiquadric, the inverse multiquadric and the thin-plate spline to reduce the computational complexity of the interpolant evaluation to O(M + N) floating point operations. By using certain localisation properties of conditionally negative definite functions this method has several performance advantages against traditional hierarchical rapid summation methods which we discuss in detail

Perturbation of the Eigenvectors of the Graph Laplacian: Application to Image Denoising

The original contributions of this paper are twofold: a new understanding of the influence of noise on the eigenvectors of the graph Laplacian of a set of image patches, and an algorithm to estimate a denoised set of patches from a noisy image. The algorithm relies on the following two observations: (1) the low-index eigenvectors of the diffusion, or graph Laplacian, operators are very robust to random perturbations of the weights and random changes in the connections of the patch-graph; and (2) patches extracted from smooth regions of the image are organized along smooth low-dimensional structures in the patch-set, and therefore can be reconstructed with few eigenvectors. Experiments demonstrate that our denoising algorithm outperforms the denoising gold-standards