Suppose you have a curve gamma on a riemannian manifold (M,g). Can you check if gamma is a geodesic?
I tried
nabla =g.connection()
nabla(gamma.nabla.tangent_vector_fields())
But it seems to be impossible to compute covariant derivatives of vectors field along curves. Is there any way to do this?
