1. I need to be efficient
Implementation corrections:
•
S
1sc
does not vary with θ within ΔΩ (Bθ,
θ
bθ constant along transport line)
bθ evaluated recursively as a moving
average of previous and current step
1
θ
2
•
S
1sc
generated per unit r is constant
within the same voxel (scerma does not
θ
3
vary considerably within voxels)
in high scerma gradient regions scerma is
estimated piecewise from a log linear
-
interpolation over r




