Abstract
The results so far appear to indicate that promising alternatives to standard backward and forward eliminations in a parallel programming environment do exist. At present a lack of available hardware prohibits the author from proving their worth further here, but at least the timing models themselves reflect how comparatively inefficient orthodox parallel triangular solvers are. Even though not every banded matrix of this type can be expected to reduce to just one sparse factor, solutions based on several sparse factors would still appear to be competetive with the existing elimination process. The 2 colour, or red/black reordering strategy also merits further investigation. This method is particularly appealing because formulating the actual reordering is fairly trivial once the nx, ny, nz parameters of the problem are known. The main disadvantage comes with the effect both reorderings may have on the convergence rate of the problem [4]. Initial investigations have found a three fold increase in the number of iterations when preconditioned TFQMR is applied to a the 2 colour reordered problem. Whether such a degredation in performance is acceptable within the context of the overall solution, or is even advisable regarding the quality of the solution, is yet to be seen. With this in mind, we may find ourselves limited to strategies (5) and (8), which only reorder when solving the triangular factorized systems in the preconditioning stage.
This work was funded by the U.K.'s Department of Trade and Industry as part of the EUREKA Project EU 638, PARSIM.
Preview
Unable to display preview. Download preview PDF.
References
F.L. Alvarado and R. Schreiber. Optimal parallel solution of sparse triangular systems. SIAM Journal of Scientific Computation, 14: 446–459, 1993.
C. Addison. An update into sparse iterative solvers for oil reservoir simulation. PARSIM Deliverable D3.3-1, August 1994.
C. Addison, T. Christensen, J. Larsen, T. Oliver and A. Sunderland. Supporting an oil reservoir simulator in a distributed memory environment. Proceedings of HPCN'94, Munich, Germany, April 1994.
I.S. Duff, G.A. Meurant. The effect of ordering on preconditioned conjugate gradients. BIT, 29: 635–657, 1989.
R. Schreiber, W.P. Tang. Vectorizing the conjugate gradient method. Unpublished Information, Dept. of Computer Science, Stanford University, 1982.
M.T. Jones, P.E. Plassmann.Scalable iterative solution of sparse linear systems. Parallel Computing, 20: 753–773, 1994.
O. Axelsson. Iterative Solution Methods. Cambridge University Press, 1994.
R. W. Freund, A Transpose-Free Quasi-Minimal Residual Algorithm for non-Hermitian Linear Systems, SIAM J. Sci. Comput, 14: No. 2, 470–482, 1993.
Author information
Authors and Affiliations
Editor information
Rights and permissions
Copyright information
© 1995 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Sunderland, A. (1995). Parallel solution strategies for triangular systems arising from oil reservoir simulations. In: Hertzberger, B., Serazzi, G. (eds) High-Performance Computing and Networking. HPCN-Europe 1995. Lecture Notes in Computer Science, vol 919. Springer, Berlin, Heidelberg. https://doi.org/10.1007/BFb0046623
Download citation
DOI: https://doi.org/10.1007/BFb0046623
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-59393-5
Online ISBN: 978-3-540-49242-9
eBook Packages: Springer Book Archive