These are already computed when doing anisotropic scattering, there has been interest in having them output