We present a fast method for numerically solving the inhomogeneous Helmholtz equation. Our iterative method is based on the Born series, which we modified to achieve convergence for scattering media of arbitrary size and scattering strength. Compared to pseudospectral time-domain simulations, our modified Born approach is two orders of magnitude faster and nine orders of magnitude more accurate in benchmark tests in 1, 2, and 3-dimensional systems.

The ability to implement reconfigurable linear optical circuits is a fundamental building block for the implementation of scalable quantum technologies. Here, we implement such circuits in a multimode fiber by harnessing its complex mixing using wavefront shaping techniques. We program linear transformations involving spatial and polarization modes of the fiber and experimentally demonstrate their accuracy and robustness using two-photon quantum states. In particular, we illustrate the reconfigurability of our platform by emulating a tunable coherent absorption experiment, where output probabilities of single-and two-photon survivals can be controlled. By demonstrating complex, reprogrammable, reliable, linear transformations, with the potential to scale, our results highlight the potential of complex media driven by wavefront shaping for quantum information processing.