Direct Numerical Computation of the Stationary Radiative Transport Equation on Modern Parallel Architectures

Abstract

We study direct numerical computations of the stationary radiative transport equation for light propagation in tissue with a three dimensional MR image and bio-optical parameters. We employ an upwind finite difference scheme which has mathematical reliability. The computational times in modern three parallel architectures are compared. Although its discretization is a large scale problem with several billions unknowns, our result shows practicality of direct numerical computation from viewpoints of computational time and modern computer resources.