| Contents |
| Propagation of waves in 2 dimensions, or in the presence of damping, c.f. viscoelastodynamics, follows the weak Huygen's principle. This makes long time computations expensive when using time-domain boundary integral equations as the complete history needs to be stored . We will show how the smoothness of this tail can be exploited to perform fast computations. In particular, we will show that the late time propagation is goverened by a parabolic operator and that consequently oblivious quadrature can be applied. Oblivious quadrature requires the storage to increase logarithmically rather than linearly in the number of time-steps. The method will be illustrated by numerical examples. |
|