×

A parallel frontal solver for finite element applications. (English) Zbl 1008.76043

Summary: In finite element simulations, the overall computing time is dominated by the time needed to solve large sparse linear systems of equations. We report on the design and development of a parallel frontal code that can significantly reduce the wallclock time needed for the solution of these systems. The algorithm used is based on dividing the finite element domain into subdomains and applying the frontal method to each subdomain in parallel. The so-called multiple front approach is shown to reduce the amount of work and memory required compared with the frontal method and, when run on a small number of processes, achieves good speedups. The code, \(\text{HSL}_- \text{MP42}\), has been developed for the Harwell Subroutine Library (http://www.numerical.rl.ac.uk/hsl). It is written in Fortran 90 and, by using MPI for message passing, achieves portability across a wide range of modern computer architectures.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
76S05 Flows in porous media; filtration; seepage
65Y05 Parallel numerical computation
86A05 Hydrology, hydrography, oceanography

Software:

HSL; HSL_MP42; MA32; MA42
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] MA42?a new frontal code for solving sparse unsymmetric systems. Technical Report RAL-93-064, Rutherford Appleton Laboratory, 1993.
[2] Duff, ACM Transactions of the Mathematical Software 22 pp 30– (1996) · Zbl 0884.65018
[3] MA32?a package for solving sparse unsymmetric systems using the frontal method. Report AERE R10079, Her Majesty’s Stationery Office, London, 1981.
[4] HSL. A Collection of Fortran Codes for Large Scale Scientific Computation, 2000. Full details from http://www.numerical.rl.ac.uk/hsl.
[5] Dongarra, ACM Transactions of the Mathematical Software 16 pp 1– (1990) · Zbl 0900.65115
[6] The use of multiple fronts in Gaussian elimination. Technical Report RAL-94-040, Rutherford Appleton Laboratory, 1994. · Zbl 0819.65023
[7] The use of multiple fronts in Gaussian elimination. In: Proceedings of the 5th SIAM Conference Applied Linear Algebra, (ed). SIAM, Philadelphia, PA, 1994; 567-571. · Zbl 0819.65023
[8] Irons, International Journal for Numerical Methods in Engineering 2 pp 5– (1970) · Zbl 0252.73050
[9] Duff, SIAM Journal of Scientific and Statistical Computing 5 pp 270– (1984) · Zbl 0541.65017
[10] Duff, International Journal for Numerical Methods in Engineering 28 pp 2555– (1989) · Zbl 0725.65045
[11] Benner, International Journal of Supercomputer Applications 1 pp 26– (1987)
[12] Zang, Computers and Structures 38 pp 202– (1991)
[13] Direct solution of two-dimensional finite element equations on distributed memory parallel computers. In: High Performance Computing, (eds). Elsevier: Amsterdam, 1991.
[14] Mallya, Computers in Chemical Engineering 21 pp s439– (1997)
[15] The design of a parallel frontal solver. Technical Report RAL-TR-99-075, Rutherford Appleton Laboratory, 1999.
[16] Scott, Communications in Numerical Methods in Engineering 15 pp 309– (1999) · Zbl 0941.65042
[17] Exploiting zeros in frontal solvers. Technical Report RAL-TR-98-041, Rutherford Appleton Laboratory, 1997.
[18] The Chaco user’s guide: Version 2.0. Technical Report SAND94-2692, Sandia National Laboratories, Albuquerque, NM, 1995.
[19] Cliffe, International Journal for Numerical Methods in Engineering 42 pp 127– (1998) · Zbl 0909.65019
[20] Mallya, A.I.Ch.E. Journal 43 pp 1032– (1997)
[21] Amestoy, International Journal of Supercomputer Applications 3 pp 41– (1989)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.