Previous Work on Multiple Scattering
; calculate scattered light at each voxel
; gather scattered light along each ray
direction bins of constant intensity;
Phase function is ignored in the 2nd path
limited to layered surfaces
LU-decomposition of matrix equation;