Solving Unsymmetric Sparse Systems of Linear Equations with PARDISO
Informazioni aggiuntive
Autori
Schenk O.,
Gaertner K.
Tipo
Articolo pubblicato in rivista scientifica
Anno
2002
Lingua
Inglese
Sommario
Supernode partitioning for unsymmetric matrices together with complete block diagonal supernode pivoting and asynchronous computation can achieve high gigaflop rates for parallel sparse LU factorization on shared memory parallel computers. The progress in weighted graph matching algorithms helps to extend these concepts further and unsymmetric prepermutation of rows is used to place large matrix entries on the diagonal. Complete block diagonal supernode pivoting allows dynamical interchanges of columns and rows during the factorization process. The level-3 BLAS efficiency is retained and an advanced two-level left–right looking scheduling scheme results in good speedup on SMP machines. These algorithms have been integrated into the recent unsymmetric version of the PARDISO solver. Experiments demonstrate that a wide set of unsymmetric linear systems can be solved and high performance is consistently achieved for large sparse unsymmetric matrices from real world applications
Parole chiave
Computational sciences, Direct solver, Numerical linear algebra, Unsymmetric linear systems
Periodico
Future Generation Computer Systems
Pagine (o numero dell’articolo)
475-487
