The eciency of the frontal method for the solution of ยฎnite-element problems depends on the order in which the elements are assembled. This paper looks at using variants of Sloan's algorithm to reorder the elements. Both direct and indirect reordering algorithms are considered and are used in combin
A new row ordering strategy for frontal solvers
โ Scribed by Jennifer A. Scott
- Publisher
- John Wiley and Sons
- Year
- 1999
- Tongue
- English
- Weight
- 99 KB
- Volume
- 6
- Category
- Article
- ISSN
- 1070-5325
No coin nor oath required. For personal study only.
โฆ Synopsis
The frontal method is a variant of Gaussian elimination that has been widely used since the mid 1970s. In the innermost loop of the computation the method exploits dense linear algebra kernels, which are straightforward to vectorize and parallelize. This makes the method attractive for modern computer architectures. However, unless the matrix can be ordered so that the front is never very large, frontal methods can require many more floating-point operations for factorization than other approaches. We are interested in matrices that have a highly asymmetric structure. We use the idea of a row graph of an unsymmetric matrix combined with a variant of Sloan's profile reduction algorithm to reorder the rows. We also look at applying the spectral method to the row graph. Numerical experiments performed on a range of practical problems illustrate that our proposed MSRO and hybrid MSRO row ordering algorithms yield substantial reductions in the front sizes and, when used with a frontal solver, significantly enhance its performance both in terms of the factorization time and storage requirements.
๐ SIMILAR VOLUMES
The e ciency of the row-by-row frontal method for the solution of unsymmetric sparse linear systems of equations is dependent on the row ordering used. Numerical experience has shown us that it can be advantageous to reverse a given row ordering. We present two results on invariances under the rever
In recent years there have been a number of important developments in frontal algorithms for solving the large sparse linear systems of equations that arise from รฟnite-element problems. We report on the design of a new fully portable and e cient frontal solver for large-scale real and complex unsymm