This paper focuses on designing and implementing parallel Moving Least Squares (MLS) interpolation algorithms by exploiting the Graphics Processing Unit (GPU) for the usage in Meshfree methods. The MLS method is an approach for scattered points' approximation / interpolation, which is commonly employed as the shape functions in various Meshfree methods. To improve the computational efficiency in building stiffness matrices in Meshfree methods, we are specifically interested in parallelizing the MLS interpolation on the GPU. The accelerated Meshfree methods can be employed to numerically analyze large deformations of soil and rock masses such as the landslides. In this paper, we first introduce our previously proposed method for finding the k nearest neighboring points located within the local region of the interest point in MLS. We then develop one sequential and three parallel implementations for each of three variations of the MLS interpolation, including the serial implementation, the parallel implementation on the multi-core CPU, the parallel implementation on a single GPU, and the parallel implementation on multi-GPUs.To evaluate the computational performance of our GPU implementations, five groups of benchmark tests are conducted in both two-dimensions and three-dimensions. We observe that our GPU implementations can achieve satisfied speedups over the sequential CPU implementation for varied sizes of testing data. To benefit the community, all source code and testing data related to the presented parallel MLS interpolation are publicly available.
The building of large-scale Digital Elevation Models (DEMs) using various interpolation algorithms is one of the key issues in geographic information science. Different choices of interpolation algorithms may trigger significant differences in interpolation accuracy and computational efficiency, and a proper interpolation algorithm needs to be carefully used based on the specific characteristics of the scene of interpolation. In this paper, we comparatively investigate the performance of parallel Radial Basis Function (RBF)-based, Moving Least Square (MLS)-based, and Shepard’s interpolation algorithms for building DEMs by evaluating the influence of terrain type, raw data density, and distribution patterns on the interpolation accuracy and computational efficiency. The drawn conclusions may help select a suitable interpolation algorithm in a specific scene to build large-scale DEMs.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.