Affiliation:
1. Univ. of Florida, Gainesville
2. Rutherford Appleton Laboratory, Chilton, Didcot, UK
Abstract
We discuss the organization of frontal matrices in multifrontal methods for the solution of large sparse sets of unsymmetric linear equations. In the multifrontal method, work on a frontal matrix can be suspended, the frontal matrix can be stored for later reuse, and a new frontal matrix can be generated. There are thus several frontal matrices stored during the factorization, and one or more of these are assembled (summed) when creating a new frontal matrix. Although this means that arbitrary sparsity patterns can be handled efficiently, extra work is required to sum the frontal matrices together and can be costly because indirect addressing is requred. The (uni)frontal method avoids this extra work by factorizing the matrix with a single frontal matrix. Rows and columns are added to the frontal matrix, and pivot rows and columns are removed. Data movement is simpler, but higher fill-in can result if the matrix cannot be permuted into a variable-band form with small profile. We consider a combined unifrontal/multifrontal algorithm to enable general fill-in reduction orderings to be applied without the data movement of previous multifrontal approaches. We discuss this technique in the context of a code designed for the solution of sparse systems with unsymmetric pattern.
Publisher
Association for Computing Machinery (ACM)
Subject
Applied Mathematics,Software
Reference31 articles.
1. Vectorization of a multiprocessor multifrontal code;AMESTOY P. R.;Int. J. Supercomput. Appl. High Perform. Eng.,1989
2. An Approximate Minimum Degree Ordering Algorithm
3. Solving sparse linear systems with sparse backward error;ARIOLI M.;SIAM J. Matrix Anal. Appl.,1989
4. A linear time implementation of the reverse Cuthill- Mckee algorithm;CHAN W. M.;BIT,1980
Cited by
207 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献