Before revision:  We suggest simple effective iterative methods based on FFT to implement the schemes whose convergence rate, under the stability condition, is fast and independent on both the meshes and variable sound speed .
After revision:  We suggest simple effective iterative methods based on FFT to implement the schemes .