There is a significant reduction of processing time and speedup of performance in computer graphics with the emergence of Graphic Processing Units (GPUs). GPUs have been developed to surpass Central Processing Unit (CPU) in terms of performance and processing speed. This evolution has opened up a new area in computing and researches where highly parallel GPU has been used for non-graphical algorithms. Physical or phenomenal simulations and modelling can be accelerated through General Purpose Graphic Processing Units (GPGPU) and Compute Unified Device Architecture (CUDA) implementations. These phenomena can be represented with mathematical models in the form of Ordinary Differential Equations (ODEs) which encompasses the gist of change rate between independent and dependent variables. ODEs are numerically integrated over time in order to simulate these behaviours. The classical Runge-Kutta (RK) scheme is the common method used to numerically solve ODEs. The Runge Kutta Fehlberg (RKF) scheme has been specially developed to provide an estimate of the principal local truncation error at each step, known as embedding estimate technique. This paper delves into the implementation of RKF scheme for GPU devices and compares its result with Dorman Prince method. A pseudo code is developed to show the implementation in detail. Hence, practitioners will be able to understand the data allocation in GPU, formation of RKF kernels and the flow of data to/from GPU-CPU upon RKF kernel evaluation. The pseudo code is then written in C Language and two ODE models are executed to show the achievable speedup as compared to CPU implementation. The accuracy and efficiency of the proposed implementation method is discussed in the final section of this paper.
Skip Nav Destination
Article navigation
10 July 2014
PROCEEDINGS OF THE 21ST NATIONAL SYMPOSIUM ON MATHEMATICAL SCIENCES (SKSM21): Germination of Mathematical Sciences Education and Research towards Global Sustainability
6–8 November 2013
Penang, Malaysia
Research Article|
July 10 2014
GPU acceleration of Runge Kutta-Fehlberg and its comparison with Dormand-Prince method
Wo Mei Seen;
Wo Mei Seen
School of Informatics & Applied Mathematics, University Malaysia Terengganu, 20130 Kuala Terengganu, Terengganu,
Malaysia
Search for other works by this author on:
R. U. Gobithaasan;
R. U. Gobithaasan
School of Informatics & Applied Mathematics, University Malaysia Terengganu, 20130 Kuala Terengganu, Terengganu,
Malaysia
Search for other works by this author on:
Kenjiro T. Miura
Kenjiro T. Miura
Graduate School of Science & Technology, Shizuoka University, Shizuoka Prefecture,
Japan
Search for other works by this author on:
AIP Conf. Proc. 1605, 16–21 (2014)
Citation
Wo Mei Seen, R. U. Gobithaasan, Kenjiro T. Miura; GPU acceleration of Runge Kutta-Fehlberg and its comparison with Dormand-Prince method. AIP Conf. Proc. 10 July 2014; 1605 (1): 16–21. https://doi.org/10.1063/1.4887558
Download citation file:
Citing articles via
The effect of a balanced diet on improving the quality of life in malignant neoplasms
Yu. N. Melikova, A. S. Kuryndina, et al.
Animal intrusion detection system using Mask RCNN
C. Vijayakumaran, Dakshata, et al.
Recognition of cat ras of face and body using convolutional neural networks
Akhmad Wahyu Aji, Esmeralda Contessa Djamal, et al.
Related Content
A comparative study of Taylor method, fourth order Runge-Kutta method and Runge-Kutta Fehlberg method to solve ordinary differential equations
AIP Conf. Proc. (March 2024)
Asymptotic solution for heat convection-radiation equation
AIP Conference Proceedings (July 2014)
Variational iteration method for solving sea-air oscillator of the ENSO model
AIP Conference Proceedings (October 2015)
Marangoni mixed convection boundary layer flow with suction and injection
AIP Conference Proceedings (April 2013)
Numerical solutions of the stagnation-point flow and heat transfer towards an exponentially stretching/ shrinking sheet with constant heat flux
AIP Conference Proceedings (February 2015)