Fix calculate_jacobian

Author: HarrisonGrodin

src/differentials.jl

@@ -79,6 +79,6 @@ macro derivatives(x...)
     esc(_differential_macro(x))
 end
 
-function calculate_jacobian(eqs, dvs, iv)
-    Expression[Differential(dv(iv()))(eq) for eq ∈ eqs, dv ∈ dvs]
+function calculate_jacobian(eqs, dvs)
+    Expression[Differential(dv)(eq) for eq ∈ eqs, dv ∈ dvs]
 end

src/systems/diffeqs/diffeqsystem.jl

@@ -62,7 +62,10 @@ function calculate_jacobian(sys::ODESystem)
     isempty(sys.jac[]) || return sys.jac[]  # use cached Jacobian, if possible
     rhs = [eq.rhs for eq ∈ sys.eqs]
 
-    jac = expand_derivatives.(calculate_jacobian(rhs, sys.dvs, sys.iv))
+    iv = sys.iv()
+    dvs = [dv(iv) for dv ∈ sys.dvs]
+
+    jac = expand_derivatives.(calculate_jacobian(rhs, dvs))
     sys.jac[] = jac  # cache Jacobian
     return jac
 end