“…Recently the efficiency of TDIE techniques has increased dramatically with the advent of acceleration schemes such as the plane wave time domain algorithm [11] and the fast Fourier transform accelerated schemes [12]. Furthermore, the marching-on-in-time methods for solving TDIE, that have been shown prone to late time instabilities, have also been successfully treated for various scenarios through the use of spatial and/or temporal filtering techniques [13][14][15][16], and the implementation of accurate spatial interpolation rules with implicit time-stepping schemes [17][18][19][20]. In addition, a backward finite difference approximation, resulting in an implicit technique, to minimise these late-time oscillations has also been proposed [21,22].…”