Higher-order derivatives of the QR and of the real symmetric eigenvalue decomposition in forward and reverse mode algorithmic differentiation

We address the task of higher-order derivative evaluation of computer programs that contain QR decompositions and real symmetric eigenvalue decompositions. The approach is a combination of univariate Taylor polynomial arithmetic and matrix calculus in the (combined) forward/reverse mode of Algorithmic Differentiation (AD). Explicit algorithms are derived and presented in an accessible form. The approach is illustrated via examples.