fix the other jacobian.

JuliaManifolds · Sep 23, 2021 · 378bbba · 378bbba
1 parent a7d1edb
commit 378bbba
Showing 1 changed file with 16 additions and 8 deletions.
diff --git a/src/manifolds/ConnectionManifold.jl b/src/manifolds/ConnectionManifold.jl
@@ -106,15 +106,23 @@ function christoffel_symbols_second_jacobian(
     M::AbstractManifold,
     p,
     B::AbstractBasis;
-    backend::AbstractDiffBackend=diff_backend(),
+    backend::AbstractDiffBackend=default_differential_backend(),
 )
-    n = size(p, 1)
+    d = manifold_dimension(M)
     ∂Γ = reshape(
-        _jacobian(q -> christoffel_symbols_second(M, q, B; backend=backend), p, backend),
-        n,
-        n,
-        n,
-        n,
+        _jacobian(
+            c -> christoffel_symbols_second(
+                M,
+                exp(M, p, get_vector(M, p, c, B)),
+                B;
+                backend=backend
+            ),
+            p,
+            backend),
+        d,
+        d,
+        d,
+        d,
     )
     return ∂Γ
 end
@@ -236,7 +244,7 @@ function riemann_tensor(
     M::AbstractManifold,
     p,
     B::AbstractBasis;
-    backend::AbstractDiffBackend=diff_backend(),
+    backend::AbstractDiffBackend=default_differential_backend(),
 )
     n = size(p, 1)
     Γ = christoffel_symbols_second(M, p, B; backend=backend)