Parallel Explicit Group Domain Decomposition Methods for the Telegraph Equation

In a previous work, we presented the numerical solution of the two dimensional second order telegraph partial differential equation discretized by the centred and rotated five-point finite difference discretizations, namely the explicit group (EG) and explicit decoupled group (EDG) iterative methods, respectively. In this paper, we utilize a domain decomposition algorithm on these group schemes to divide the tasks involved in solving the same equation. The objective of this study is to describe the development of the parallel group iterative schemes under OpenMP programming environment as a way to reduce the computational costs of the solution processes using multicore technologies. A detailed performance analysis of the parallel implementations of points and group iterative schemes will be reported and discussed.

Authors:

• Kew Lee Ming
• Norhashidah Hj. Mohd. Ali

Keywords:

• Telegraph equation
• explicit group iterative scheme
• domain decomposition algorithm
• parallelization.

References:

[1] Kew, L.M. and Ali, N.H.M. "Explicit Group Iterative Methods for the
Solution of Telegraph Equations", Lecture Notes In Engineering and
Computer Science, World Congress On Engineering 2010, The 2010
International Conference of Applied and Engineering Mathematics, 30th
June - 2nd July 2010, Imperial College, London, 2010, pp.1770- 1775.
[2] Saad, Y. "Iterative Methods for Sparse Linear Systems". 2nd Edition.
pp. 382-421.
[3] Abdullah, A.R. "The Four Point Explicit Decoupled Group EDG
Method: A Fast Poisson Solver". International Journal of Computer
Mathematics, 38, 1991, pp.61-70.
[4] Yousif, W.S. and Evans, D.J. "Explicit De-Coupled Group Iterative
Methods and Their Parallel Implementations". Parallel Algorithms and
Applications, 7, 1995, pp.53-71.
[5] Evans, D.J. "Group Explicit Methods for the Numerical Solution of
Partial Differential Equations". Loughborough University of
Technology, UK, Gordon and Breach Science Publisher, The
Netherlands, 1997, pp.106-162.
[6] Dryja, M. and Widlund, O.B. "Some Domain Decomposition
Algorithms for Elliptic Problems", in: L. Hayes, D. Kincaid (Eds.),
Iterative Methods for Large Linear Systems, Academic Press, San
Diego, CA. 1989.
[7] Cai, X.-C. and Widlund, O.B. "Domain Decomposition Algorithms for
Indefinite Elliptic Problems", SIAM J. Sci. Statist. Comput., 13, 1992,
pp. 243-258.
[8] Dawson, C.N., Du, Q., and Dupont, T.F. "A Finite Difference Domain
Decomposition Algorithm for Numerical Solution of the Heat
Equation". Mathematics of Computation, 57(195), 1991, pp. 63-71.
[9] Jun, Y. and Mai, T.Z. (2006). "IPIC Domain Decomposition Algorithm
for Parabolic Problems". Applied Mathematics and Computation, 177,
2006, pp. 352-364.
[10] Ali, A. "Introduction to Hybrid MPI/OpenMP Programming". Presented
at ICTP Advanced School in High Performance and Grid Computing,
11-22 April 2011, Abdus Salam ICTP - Trieste.
[11] Dehghan, M. and Shokri, A., "A Meshless Method for Numerical
Solution of a Linear Hyperbolic Equation with Variable Coefficients in
Two Space Dimensions." Numerical Methods for Partial Differential
Equations, 2008, pp.494-506.