Computational aspects of the choice of operator and sampling interval for numerical differentiation in large-scale simulation of wave phenomena