“…This algorithm has also been used for a wide range of other calculations, such as solving linear systems [20,25], computing the action of the matrix exponential on a vector, see, e.g. [23,27], and even solving differential equations using exponentially-fitted methods, e.g., [16,17,12]. In exact arithmetic, the algorithm can be recast as a simple three-term recurrence, namely,…”