×

An application of partition method for solving 3D Stokes equation. (English) Zbl 1443.65123

Summary: In our previous work we have studied the performance of a parallel algorithm, based on a direction splitting approach, for solving of time dependent Stokes equation. We used a rectangular uniform mesh, combined with a central difference scheme for the second derivatives. Hence, the proposed algorithm required only solution of tridiagonal linear systems. In our work, we are targeting massively parallel computers, as well as clusters of multi-core nodes. The somehow slower (experimentally-established) performance of the proposed approach was observed when using all cores on a single node of a cluster. To remedy this problem, we tried to use LAPACK subroutines from the multi-threaded layer library, but the parallel performance of the code (while improved) was still not satisfactory on a single (multi-core) node.
Our current work considers hybrid parallelization based on the MPI and OpenMP standards. It is motivated by the need to maximize the parallel efficiency of our implementation of the proposed algorithm. Essential improvements of the parallel algorithm are achieved by introducing two levels of parallelism: (i) between-node parallelism based on the MPI and (ii) inside-node parallelism based on the OpenMP. The implementation was tested on Linux clusters with Intel processors and on the IBM supercomputer.

MSC:

65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65Y05 Parallel numerical computation
76D07 Stokes and related (Oseen, etc.) flows

Software:

MPI; LAPACK
Full Text: DOI

References:

[2] Hayryan, E. A.; Busa, J.; Donets, E. E.; Pokorny, I.; Strel’tsova, O. I., Numerical studies of perturbed static solutions decay in the coupled system of Yang-Mills-dilaton equations with use of MPI technology, Tech. Rep. (2004), Laboratory of Information Technologies, Joint Institute for Nuclear Research: Laboratory of Information Technologies, Joint Institute for Nuclear Research Dubna (Russian Federation)
[3] Wang, H., A parallel method for tridiagonal equations, ACM Trans. Math. Software (TOMS), 7, 2, 170-183 (1981) · Zbl 0473.65010
[5] Hockney, R. W., A fast direct solution of Poissons equation using Fourier analysis, J. ACM, 12, 1, 95-113 (1965) · Zbl 0139.10902
[6] Stone, H. S., An efficient parallel algorithm for the solution of a tridiagonal linear system of equations, J. ACM (JACM), 20, 1, 27-38 (1973) · Zbl 0269.65018
[7] Walshaw, C. H., Diagonal dominance in the parallel partition method for tridiagonal systems, SIAM J. Matrix Anal. Appl., 16, 4, 1086-1099 (1995) · Zbl 0841.65017
[8] Yalamov, P.; Pavlov, V., On the stability of a partitioning algorithm for tridiagonal systems, SIAM J. Matrix Anal. Appl., 20, 1, 159-181 (1998) · Zbl 0924.65019
[9] Amodio, P.; Brugnano, L., Parallel factorizations and parallel solvers for tridiagonal linear systems, Linear Algebra Appl., 172, 347-364 (1992) · Zbl 0756.65030
[10] Amodio, P.; Brugnano, L.; Politi, T., Parallel factorizations for tridiagonal matrices, SIAM J. Numer. Anal., 30, 3, 813-823 (1993) · Zbl 0809.65015
[11] Sun, X.-H.; Sun, H. Z.; Ni, L. M., Parallel algorithms for solution of tridiagonal systems on multicomputers, (Proceedings of the 3rd International Conference on Supercomputing (1989), ACM), 303-312
[12] Bondeli, S., Divide and conquer: A parallel algorithm for the solution of a tridiagonal linear system of equations, Parallel Comput., 17, 4, 419-434 (1991) · Zbl 0739.65016
[13] Hofhaus, J.; Van de Velde, E. F., Alternating-direction line-relaxation methods on multicomputers, SIAM J. Sci. Comput., 17, 2, 454-478 (1996) · Zbl 0851.65065
[14] Guermond, J.-L.; Minev, P., A new class of fractional step techniques for the incompressible Navier-Stokes equations using direction splitting, C. R. Math., 348, 9-10, 581-585 (2010) · Zbl 1352.76090
[15] Guermond, J.-L.; Minev, P., A new class of massively parallel direction splitting for the incompressible Navier-Stokes equations, Comput. Methods Appl. Mech. Engrg., 200, 23, 2083-2093 (2011) · Zbl 1228.76095
[16] Ganzha, M.; Georgiev, K.; Lirkov, I.; Margenov, S.; Paprzycki, M., Highly parallel alternating directions algorithm for time dependent problems, (Christov, C.; Todorov, M., Applications of Mathematics in Technical and Natural Sciences, AMiTaNS 2011. Applications of Mathematics in Technical and Natural Sciences, AMiTaNS 2011, AIP Conference Proceedings, vol. 1404 (2011)), 210-217
[17] Lirkov, I.; Paprzycki, M.; Ganzha, M., Performance analysis of parallel alternating directions algorithm for time dependent problems, (Wyrzykowski, R.; Dongarra, J.; Karczewski, K.; Waśniewski, J., 9th international Conference on Parallel Processing and Applied Mathematics, PPAM 2011, Part I. 9th international Conference on Parallel Processing and Applied Mathematics, PPAM 2011, Part I, Lecture Notes in Computer Science, vol. 7203 (2012), Springer), 173-182
[18] Lirkov, I.; Paprzycki, M.; Ganzha, M.; Gepner, P., Performance evaluation of MPI/OpenMP algorithm for 3D time dependent problems, (Ganzha, M.; Maciaszek, L.; Paprzycki, M., Preprints of Position Papers of the Federated Conference on Computer Science and Information Systems. Preprints of Position Papers of the Federated Conference on Computer Science and Information Systems, Annals of Computer Science and Information Systems, vol. 2 (2013)), 27-32
[19] Snir, M.; Otto, S.; Huss-Lederman, S.; Walker, D.; Dongarra, J., (MPI: The Complete Reference. MPI: The Complete Reference, Scientific and Engineering Computation Series (1997), The MIT Press: The MIT Press Cambridge, Massachusetts)
[20] Walker, D.; Dongarra, J., MPI: a standard Message Passing Interface, Supercomputer, 63, 56-68 (1996)
[21] Chandra, R.; Menon, R.; Dagum, L.; Kohr, D.; Maydan, D.; McDonald, J., Parallel Programming in OpenMP (2000), Morgan Kaufmann
[22] Chapman, B.; Jost, G.; Van Der Pas, R., Using OpenMP: Portable Shared Memory Parallel Programming, Vol. 10 (2008), MIT Press
[23] Anderson, E.; Bai, Z.; Bischof, C.; Blackford, S.; Demmel, J.; Dongarra, J.; Croz, J. D.; Greenbaum, A.; Hammarling, S.; McKenney, A.; Sorensen, D., LAPACK Users’ Guide (1999), SIAM: SIAM Philadelphia · Zbl 0934.65030
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.