Using graphics processing units to accelerate perturbation Monte Carlo simulation in a turbid medium
2012 (English)In: Journal of Biomedical Optics, ISSN 1083-3668, Vol. 17, no 4, 040502- p.Article in journal (Refereed) Published
We report a fast perturbation Monte Carlo (PMC) algorithm accelerated by graphics processing units (CPU). The two-step PMC simulation [Opt. Lett. 36, 2095 (2011)] is performed by storing the seeds instead of the photon's trajectory, and thus the requirement in computer random-access memory (RAM) becomes minimal. The two-step PMC is extremely suitable for implementation onto CPU. In a standard simulation of spatially-resolved photon migration in the turbid media, the acceleration ratio between using GPU and using conventional CPU is about 1000. Furthermore, since in the two-step PMC algorithm one records the effective seeds, which is associated to the photon that reaches a region of interest in this letter, and then re-run the MC simulation based on the recorded effective seeds, radiative transfer equation (RTE) can be solved by two-step PMC not only with an arbitrary change in the absorption coefficient, but also with large change in the scattering coefficient.
Place, publisher, year, edition, pages
2012. Vol. 17, no 4, 040502- p.
biomedical optics, Monte Carlo, scattering
Radiology, Nuclear Medicine and Medical Imaging Biochemistry and Molecular Biology
IdentifiersURN: urn:nbn:se:kth:diva-98337DOI: 10.1117/1.JBO.17.4.040502ISI: 000304382900002OAI: oai:DiVA.org:kth-98337DiVA: diva2:536942
QC 201206252012-06-252012-06-252013-11-19Bibliographically approved