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

Authors: Zanariah Abdul Majid, Mohamed Suleiman

Abstract:

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.

Keywords: Numerical methods, parallel method, block method, higher order ODEs.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1082501

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)