On the implementation of large-scale integral operators with modern HPC solutions - Application to 3D Marchenko imaging by least-squares inversion.