Circumventing the ill-conditioning problem with multiquadric radial basis functions: Applications to elliptic partial differential equations

Madych and Nelson [1] proved multiquadric (MQ) mesh-independent radial basi s functions (RBFs) enjoy exponential convergence. The primary disadvantage of the MQ scheme is that it is global, hence, the coefficient matrices obta ined from this discretization scheme are full. Full matrices tend to become progressively more ill-conditioned as the rank increases. In this paper, w e explore several techniques, each of which improves the conditioning of th e coefficient matrix and the solution accuracy. The methods that were inves tigated are (1) replacement of global solvers by block partitioning, LU decomposition s chemes, (2) matrix preconditioners, (3) variable Mg shape parameters based upon the local radius of curvature o f the function being solved, (4) a truncated MQ basis function having a finite, rather than a full band- width, (5) multizone methods for large simulation problems, and (6) knot adaptivity that minimizes the total number of knots required in a simulation problem. The hybrid combination of these methods contribute to very accurate solutio ns. Even though FEM gives rise to sparse coefficient matrices; these matrices i n practice can become very ill-conditioned. We recommend using what has bee n learned from the FEM practitioners and combining their methods with what has been learned in RBF simulations to form a flexible, hybrid approach to solve complex multidimensional problems. (C) 2000 Elsevier Science Ltd. All rights reserved.