A robust solver for elliptic PDEs in 3D complex geometries