  The Cayley-Hamilton theorem is used to implement an iterative process for the
efficient numerical computation of matrix power series and their differentials.
In addition to straight-forward applications in lattice gauge theory
simulations e.g. to reduce the computational cost of smearing, the method can
also be used to simplify the evaluation of SU(N) one-link integrals or the
computation of SU(N) matrix logarithms.
