“…However, as discussed in Section 4.1, much larger step sizes may be used with the implicit algorithm without compromising on accuracy, as is expected from an implicit integration method in a path-dependent problem. This proves useful in performance critical scenarios, as considered e.g., by Wölfle et al [4], in which the additional effort of a Newton iteration might thereby be balanced against the need to compute the thermoelastic input solution on a much finer discretization, which certainly comes at a cost, especially if double convolution integrals are evaluated numerically over the whole domain for each grid point as is commonly done [9,11,[13][14][15][16][17]23].…”