(adjoint) Jacobi fields for non mutating manifolds

[kellertuer]

For manifolds that do not do mutating operations like Circle (adjoint) jacobi fields currently fail.

Example

using Manifolds, ManifoldsDiff
M = Circle()
p = π/2
q = 0.
X = 1.
Y = ManifoldDiff.adjoint_Jacobi_field(M, p, q, 0.5, X, ManifoldDiff.βdifferential_shortest_geodesic_startpoint)