Numerical implementation of diagrammatic scattering theory for few-photon transport through Bose-Hubbard lattices