TY - EJOU
AU - Liu, Chein-Shan
AU - Yeih, Weichung
AU - Atluri, Satya N.
TI - On Solving the Ill-Conditioned System Ax=b: General-Purpose Conditioners Obtained From the Boundary-Collocation Solution of the Laplace Equation, Using Trefftz Expansions With Multiple Length Scales
T2 - Computer Modeling in Engineering \& Sciences
PY - 2009
VL - 44
IS - 3
SN - 1526-1506
AB - Here we develop a general purpose pre/post conditionerT, to solve an ill-posed system of linear equations,Ax=b. The conditionerTis obtained in the course of the solution of the Laplace equation, through a boundary-collocation Trefftz method, leading to:Ty=x, whereyis the vector of coefficients in the Trefftz expansion, andxis the boundary data at the discrete points on a unit circle. We show that the quality of the conditionerTis greatly enhanced by using multiple characteristic lengths (Multiple Length Scales) in the Trefftz expansion. We further show thatTcan be multiplicatively decomposed into a dilationTDand a rotationTR. For an odd-orderedA, we develop four conditioners based on the solution of the Laplace equation for Dirichlet boundary conditions, while for an even-orderedAwe develop four conditioners employing the Neumann boundary conditions. All these conditioners are well-behaved and easily invertible. Several examples involving ill-conditionedA, such as the Hilbert matrices, those arising from the Method of Fundamental Solutions, those arising from very-high order polynomial interpolations, and those resulting from the solution of the first-kind Fredholm integral equations, are presented. The results demonstrate that the presently proposed conditioners result in very high computational efficiency and accuracy, whenAx=bis highly ill-conditioned, andbis noisy.
KW - Ill-posed linear equations
KW - Multi-Scale Trefftz Method (MSTM)
KW - Multi-Scale Trefftz-Collocation Laplacian Conditioner (MSTCLC)
KW - Transformation matrix
KW - Dilation matrix
KW - Rotation matrix
DO - 10.3970/cmes.2009.044.281