Before revision:  A numerical implementation of these algorithms in C/C++ is benchmarked showing computation times on the order of 10-20 \mus for the computation of the Coriolis matrix and 40-120 \mus for the computation of the Christoffel symbols for systems with 20 degrees of freedom.
After revision:  A numerical implementation of these algorithms in C/C++ is benchmarked showing computation times on the order of 10-20 \mus for the computation of the Coriolis matrix and 40-120 \mus for the Christoffel symbols on systems with 20 degrees of freedom.