Very large scale finite differences in modeling of seismic waves