Solution Optimization of RBF Interpolation for Implicit Modeling of Orebody

Abstract

To interpolate large drillhole data sets efficiently in implicit modeling of orebody, we optimize the solution of the RBF equation based on the two-level domain decomposition method. The solution performance is improved by balancing the two goals of the convergence rate and the iteration efficiency. The solution method converts the large domain into small subdomains which can be solved directly. The optimum division of subdomains can improve both the convergence of iteration and the efficiency of a single iteration, whereas the subdivision of a trivial subdomain leads the iteration failed to converge. For the efficiency of a single iteration, the kernel independent fast multipole method is used to calculate the matrix-vector product efficiently. Moreover, the initial solution of iteration and center reduction strategy are optimized for the dynamic updating of implicit surface. Combined with the preconditioned Krylov subspace method, the optimized solution method was implemented. The experimental results of several drillhole data sets show that the optimized solution method converges rapidly within a small number of iterations.

Keywords

• Radial basis functions
• implicit modeling
• domain decomposition method
• implicit surface
• surface reconstruction