A stable non-deterministic parallel algortihm for general unsymmetric sparse Lu factorization:
Gespeichert in:
Bibliographische Detailangaben
Beteiligte Personen: Davis, Timothy A. (VerfasserIn), Yew, Pen-Chung (VerfasserIn)
Format: Buch
Sprache:Englisch
Veröffentlicht: Urbana, Ill. 1990
Schriftenreihe:Center for Supercomputing Research and Development <Urbana, Ill.>: CSRD report 908
Schlagwörter:
Abstract:Abstract: "We present a parallel algorithm for the direct LU factorization of general unsymmetric sparse matrices. The algorithm, D2, is based on a new nondeterministic parallel pivot search that finds a compatible pivot set S of size m, followed by a parallel rank-m update. These two steps alternate until switching to dense matrix code or until the matrix is factored. The algorithm is based on a shared-memory MIMD model and takes advantage of both concurrency and (gather-scatter) vectorization. The detection of parallelism due to sparsity is based on Markowitz's strategy, an unsymmetric ordering method
As a result, D2 finds more potential parallelism for matrices with highly asymmetric nonzero patterns than algorithms that construct an elimination tree using a symmetric ordering method (minimum degree or nested dissection, for example) applied to the symmetric pattern of A + A[superscript T] or A[superscript T]A. The pivot search exploits more parallelism than previous algorithms that are based on unsymmetric ordering methods. Possible extensions to the D2 algorithm are discussed, including the use of dense matrix kernels and a software combining tree to enhance parallelsim in the pivot search.
Umfang:32 S.