On evaluation of the solution of Gaussian impulse propagation

We consider the solution of the 2D wave equation with the initial data prescribing a Gaussian pulse. This solution has integral representations. Previously proposed algorithm of its evaluation is optimized by the use of more precise error bounds of numerical quadratures.