A symmetric Hessian propagation is proposed for explicit or implicit dynamics. The approach is presented for discrete- and continuous-time problems. A three-sweep propagation scheme can reduce the memory requirements. An open-source implementation is presented as part of the ACADO Toolkit.