Computation of three-dimensional, pulsed, nonlinear acoustic wavefields from medical phased array transducers in very large domains