Research Article

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

Table 2

GPU run-time comparisons for FEM forward solver computations of 80 source distributions in cylindrical meshes of different node densities. Real-valued problems were solved with a conjugate gradient solver, complex problems with a biconjugate gradient stabilised solver. Values in parentheses are CPU solution times.

Node number Runtime [s]
Real single Complex single Real double Complex double

8987 11.07 13.43 11.1 14.4
(8.49) (26.15) (3.74) (10.46)
82517 23.24 60.03 26.85 75.42
(193.94) (678.65) (96.21) (271.9)
245917 82.56 433.89 117.41 523.76
(1789.7) (12996.8) (837.39) (2351.25)
327617 127.19 1060.42 189.06 909.45
(3258.46) (−) (1509.22) (4699.71)