Parallel Direct Integration Variable Step Block Method for Solving Large System of Higher Order Ordinary Differential Equations

The aim of this paper is to investigate the performance of the developed two point block method designed for two processors for solving directly non stiff large systems of higher order ordinary differential equations (ODEs). The method calculates the numerical solution at two points simultaneously and produces two new equally spaced solution values within a block and it is possible to assign the computational tasks at each time step to a single processor. The algorithm of the method was developed in C language and the parallel computation was done on a parallel shared memory environment. Numerical results are given to compare the efficiency of the developed method to the sequential timing. For large problems, the parallel implementation produced 1.95 speed-up and 98% efficiency for the two processors.

Authors:

• Zanariah Abdul Majid
• Mohamed Suleiman

Keywords:

• Numerical methods
• parallel method
• block method
• higher order ODEs.

References:

[1] P.C. Chakravarti, P.B. Worland. A class of self starting methods for the
numerical solution of yÔÇ▓ÔÇ▓ = f (x, y), BIT 11, pp 368-383, 1971.
[2] M.B.Suleiman, Solving higher order odes directly by the direct
integration method, Applied Mathematics and Computation 33, pp 197-
219., 1989.
[3] S O.Fatunla, Block methods for second order odes. Intern. J. Computer
Math 40, pp 55-63, .1990.
[4] Z. Omar, M. Suleiman, Parallel two-point explicit block method for
solving high-order ordinary differential equations. Int. J. of Simulation
and Process Modelling. Vol. 2, No.3/4 pp. 227 - 231, 2006.
[5] Z. Omar and M. Suleiman, Parallel r-point implicit block method for
solving higher order ordinary differential equations directly, Journal of
ICT, 3(1), pp 53-66, 2004.
[6] N.H. Cong, K. Strehmel, R. Weiner, and H. Podhaisky, Runge-Kutta-
Nystrom-type parallel block predictor-corrector methods. Advances in
Computational Mathematics 10, pp 115-133, 1999.
[7] Z. A. Majid and M. Suleiman, Two point block direct integration
implicit variable steps method for solving higher order systems of
ordinary differential equations. International Conference of Applied and
Engineering Mathematics. WCE (London). Proceeding of the World
Congress on Engineering 2007, WCE 2007, Volume II, pp 812-815,
2007.
[8] Z. Omar, Developing parallel block methods for solving higher order
odes directly, Ph.D. Thesis, University Putra Malaysia, Malaysia, 1999.
[9] Hairer.E., Norsett. S.P. and Wanner. G., Solving Ordinary Differential
Equations I: Nonstiff Problems. Berlin: Springer-Verlag. pp 26, 1993.
[10] Cong, N.H., Podhaisky, H. and Weiner, R., Performance of explicit
pseudo two-step RKN methods on a shared memory computer, 2001.
(http://www.mathematik.uni-halle.de/reports/rep-num.html)