VECTORIZED STRONGLY IMPLICIT SOLVING PROCEDURE FOR A SEVEN DIAGONAL COEFFICIENT MATRIX英文文献.pdfVIP

INT. J . NUM. METH. HEAT FLUID FLOW, VOL. 4 , 159-17 2 (1994) VECTORIZED STRONGLY IMPLICIT SOLVING PROCEDURE FOR A SEVEN-DIAGONAL COEFFICIENT MATRIX H.-J . LEISTER AND M . PERIĆ Lehrstuhl für Strömungsmechanik, Universitāt Erlangen-Nürnberg, Cauerstr. 4, D-8520 Erlangen, Germany ABSTRACT 1 The paper presents an extension of Stones strongly implicit procedure for solving linear equation systems resulting from the discretization of partial differential equations to three-dimensional problems. The solver is applicable to seven-diagonal coefficient matrices, as are obtained when central-difference approximations are used for discretization. The algorithm is implemented in a way which allows vector processing on modern supercomputers, in spite of its recursive structure. Other solvers, using incomplete lower-upper decomposition (ILU), can be vectorized in the same way. Test calculations show solver performance of about 150 Mflops on CRAY-YMP and over 200 Mflops on FUJITSU-VP200 computers. A listing of the FORTRAN code is provided. KEY WORDS Implicit matrix solver Vectorization ILU decomposition INTRODUCTION Solution of large linear equation systems uses most of computing time in computational fluid dynamics (CFD) codes . Due to the non-linearity and coupling of partial differential equations , outer iterations must be used to update the coefficients and source terms in the linear equation systems. These a