reimplement the scattering calculation with the eigen matrix