A parallel and vectorial implementation of basic linear algebra subroutines in iterative solving of large sparse linear systems of equations

摘要

Electromagnetic field analysis by finite element methods, which involve the solution of large sparse systems of linear equations, is discussed. Though no discernible structure for the distribution of nonzero elements can be found (e.g. multidiagonal structures), subsets of independent equations can be determined. Equations that are in the same subset are then solved in parallel. A good choice for the storage scheme of sparse matrices is also very important to speed up the resolution by vectorization. The modifications made to data structures are presented, and the possibility of using other schemes is discussed.