Jacobians work!

SciML · Nov 21, 2017 · 1634906 · 1634906
1 parent 5cd03d3
commit 1634906
Show file tree

Hide file tree

Showing 3 changed files with 136 additions and 119 deletions.
diff --git a/src/common.jl b/src/common.jl
@@ -35,10 +35,6 @@ function solve{uType,duType,tType,isinplace,LinearSolver}(
                 warn("Explicit t-gradient given to this stiff solver is ignored.")
                 warned = true
             end
-            if has_jac(prob.f)
-                warn("Explicit Jacobian given to this stiff solver is ignored.")
-                warned = true
-            end
         end
         warned && warn_compat()
     end
@@ -146,7 +142,8 @@ function solve{uType,duType,tType,isinplace,LinearSolver}(
     res = DASKR.res_c(f!)
     rt = Int32[0]
     if has_jac(f!)
-      jac = jac_c((t, y, yp, pd, cj) -> f(Val{:jac},t,y,dy,cj,pd))
+      jac = common_jac_c(f!)
+      info[5] = 1 # Enables Jacobian
     else
       jac = Int32[0]
     end