jwf, IDL from posersurface approximates inverse square fall-off for each bounce, given that the
surfaces are large (not point sources), hence the math may be complex, as mentioned by bill. 
some of these things are not fully documented yet.

not complex as in complex number (real and imaginary component).