A Numerical Method for Solving Ill-Conditioned Equation Systems Arising from Radial Basis Functions  

作　　者：Edward J. Kansa Edward J. Kansa(Convergent Solutions, Livermore, CA, USA)

机构地区：[1]Convergent Solutions, Livermore, CA, USA

出　　处：《American Journal of Computational Mathematics》2023年第2期356-370,共15页美国计算数学期刊（英文）

摘　　要：Continuously differentiable radial basis functions (C^∞-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C^∞-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C^∞-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.Continuously differentiable radial basis functions (C^∞-RBFs), while being theoretically exponentially convergent are considered impractical computationally because the coefficient matrices are full and can become very ill- conditioned. Similarly, the Hilbert and Vandermonde have full matrices and become ill-conditioned. The difference between a coefficient matrix generated by C^∞-RBFs for partial differential or integral equations and Hilbert and Vandermonde systems is that C^∞-RBFs are very sensitive to small changes in the adjustable parameters. These parameters affect the condition number and solution accuracy. The error terrain has many local and global maxima and minima. To find stable and accurate numerical solutions for full linear equation systems, this study proposes a hybrid combination of block Gaussian elimination (BGE) combined with arbitrary precision arithmetic (APA) to minimize the accumulation of rounding errors. In the future, this algorithm can execute faster using preconditioners and implemented on massively parallel computers.

关 键 词：Continuously Differentiable Radial Basis Functions Global Maxima and Minima Solutions of Ill-Conditioned Linear Equations Block Gaussian Elimination Arbitrary Precision Arithmetic 

分 类 号：O17[理学—数学]