Beam propagation method based on the iterated Crank-Nicolson scheme for solving large-scale wave propagation problems