GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography

We introduce a GPU-accelerated finite element forward solver for the computation of light transport in scattering media. The forward model is the computationally most expensive component of iterative methods for image reconstruction in diffuse optical tomography, and performance optimisation of the...

Full description

Saved in:
Bibliographic Details
Main Author: Martin Schweiger
Format: Article
Language:English
Published: Wiley 2011-01-01
Series:International Journal of Biomedical Imaging
Online Access:http://dx.doi.org/10.1155/2011/403892
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1832553373633609728
author Martin Schweiger
author_facet Martin Schweiger
author_sort Martin Schweiger
collection DOAJ
description We introduce a GPU-accelerated finite element forward solver for the computation of light transport in scattering media. The forward model is the computationally most expensive component of iterative methods for image reconstruction in diffuse optical tomography, and performance optimisation of the forward solver is therefore crucial for improving the efficiency of the solution of the inverse problem. The GPU forward solver uses a CUDA implementation that evaluates on the graphics hardware the sparse linear system arising in the finite element formulation of the diffusion equation. We present solutions for both time-domain and frequency-domain problems. A comparison with a CPU-based implementation shows significant performance gains of the graphics accelerated solution, with improvements of approximately a factor of 10 for double-precision computations, and factors beyond 20 for single-precision computations. The gains are also shown to be dependent on the mesh complexity, where the largest gains are achieved for high mesh resolutions.
format Article
id doaj-art-ffed561103864003ac4ea892b6143f48
institution Kabale University
issn 1687-4188
1687-4196
language English
publishDate 2011-01-01
publisher Wiley
record_format Article
series International Journal of Biomedical Imaging
spelling doaj-art-ffed561103864003ac4ea892b6143f482025-02-03T05:54:21ZengWileyInternational Journal of Biomedical Imaging1687-41881687-41962011-01-01201110.1155/2011/403892403892GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical TomographyMartin Schweiger0Department of Computer Science, University College London, Gower Street, London WC1E 6BT, UKWe introduce a GPU-accelerated finite element forward solver for the computation of light transport in scattering media. The forward model is the computationally most expensive component of iterative methods for image reconstruction in diffuse optical tomography, and performance optimisation of the forward solver is therefore crucial for improving the efficiency of the solution of the inverse problem. The GPU forward solver uses a CUDA implementation that evaluates on the graphics hardware the sparse linear system arising in the finite element formulation of the diffusion equation. We present solutions for both time-domain and frequency-domain problems. A comparison with a CPU-based implementation shows significant performance gains of the graphics accelerated solution, with improvements of approximately a factor of 10 for double-precision computations, and factors beyond 20 for single-precision computations. The gains are also shown to be dependent on the mesh complexity, where the largest gains are achieved for high mesh resolutions.http://dx.doi.org/10.1155/2011/403892
spellingShingle Martin Schweiger
GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
International Journal of Biomedical Imaging
title GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
title_full GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
title_fullStr GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
title_full_unstemmed GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
title_short GPU-Accelerated Finite Element Method for Modelling Light Transport in Diffuse Optical Tomography
title_sort gpu accelerated finite element method for modelling light transport in diffuse optical tomography
url http://dx.doi.org/10.1155/2011/403892
work_keys_str_mv AT martinschweiger gpuacceleratedfiniteelementmethodformodellinglighttransportindiffuseopticaltomography