Abstract
Scattered data interpolation using conditionally positive definite radial basis functions (RBFs) requires the solution of a symmetric saddle-point system. Based on an approximation of the system matrix as a hierarchical matrix, we solve the system iteratively using the GMRes algorithm and a domain decomposition preconditioner. The novelty of our work lies in the proposed solution of the subdomain problems using the nullspace method with an orthogonal basis represented as a sequence of Householder reflectors. The resulting positive definite subdomain systems are solved either directly or using an inner GMRes iteration with H-Cholesky preconditioning. Numerical tests demonstrate the effectiveness of this solution process for up to N=160000 centers in two and three dimensions.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
