add tests for Jacobians of LabelledArrays

Author: ChrisRackauckas

test/labelledarrays.jl

@@ -13,14 +13,20 @@ eqs = [D(x) ~ σ*(y-x),
        D(z) ~ x*y - β*z]
 
 de = ODESystem(eqs)
-f = ODEFunction(de, [x,y,z], [σ,ρ,β])
+ff = ODEFunction(de, [x,y,z], [σ,ρ,β], jac=true)
 
 a = @SVector [1.0,2.0,3.0]
 b = SLVector(x=1.0,y=2.0,z=3.0)
 c = [1.0,2.0,3.0]
 p = SLVector(σ=10.0,ρ=26.0,β=8/3)
-@test f(a,p,0.0) isa SVector
-@test typeof(f(b,p,0.0)) <: SLArray
-@test f(c,p,0.0) isa Vector
-@test f(a,p,0.0) == f(b,p,0.0)
-@test f(a,p,0.0) == f(c,p,0.0)
+@test ff(a,p,0.0) isa SVector
+@test typeof(ff(b,p,0.0)) <: SLArray
+@test ff(c,p,0.0) isa Vector
+@test ff(a,p,0.0) == ff(b,p,0.0)
+@test ff(a,p,0.0) == ff(c,p,0.0)
+
+@test ff.jac(a,p,0.0) isa SMatrix
+@test typeof(ff.jac(b,p,0.0)) <: SMatrix
+@test ff.jac(c,p,0.0) isa Matrix
+@test ff.jac(a,p,0.0) == ff.jac(b,p,0.0)
+@test ff.jac(a,p,0.0) == ff.jac(c,p,0.0)