Weiran Yuan1, Pu Chen1,2, Kaishin Liu1,3
CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.2, pp. 115-134, 2007, DOI:10.3970/cmes.2007.017.115
Abstract In this paper we propose a direct solution method for the quasi-unsymmetric sparse matrix (QUSM) arising in the Meshless Local Petrov-Galerkin method (MLPG). QUSM, which is conventionally treated as a general unsymmetric matrix, is unsymmetric in its numerical values, but nearly symmetric in its nonzero distribution of upper and lower triangular portions. MLPG employs trial and test functions in different functional spaces in the local domain weak form of governing equations. Consequently the stiffness matrix of the resultant linear system is a QUSM. The new solver for QUSM conducts a two-level unrolling technique for LDU factorization More >