## JMI2011C-10 Image segmentation using CUDA implementations of the Runge-Kutta-Merson and GMRES methods (pp.73-79)

Author(s)： Tomáš Oberhuber, Atsushi Suzuki, Jan Vacata and Vítězslav Žabka

J. Math-for-Ind. **3C** (2011) 73-79.

- File： JMI2011C-10.pdf (334KB)

Abstract

Modern GPUs are well suited for performing image processing tasks. We utilize their high computational performance and memory bandwidth for image segmentation purposes. We segment cardiac MRI data by means of numerical solution of an anisotropic partial differential equation of the Allen-Cahn type. We implement two different algorithms for solving the equation on the CUDA architecture. One of them is based on the Runge-Kutta-Merson method for the approximation of solutions of ordinary differential equations, the other uses the GMRES method for the numerical solution of systems of linear equations. In our experiments, the CUDA implementations of both algorithms are about 3--9 times faster than corresponding 12-threaded OpenMP implementations.

Keyword(s). CUDA, image segmentation, Allen-Cahn equation, Runge-Kutta-Merson, GMRES