utilize Wfact overloads

Author: ChrisRackauckas

src/ode/ode_simple_nonlinear_prob.jl

@@ -52,7 +52,7 @@ eqs = [D(x) ~ μ*((1-x^2)*y - x),
        D(y) ~ y]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-#ModelingToolkit.generate_factorized_W(de)
+ModelingToolkit.calculate_factorized_W(de)
 van = ODEFunction(de, [x,y], [μ])
 
 """
@@ -95,7 +95,7 @@ eqs = [D(y₁) ~ -k₁*y₁+k₃*y₂*y₃,
        D(y₃) ~ k₂*y₂^2]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-#ModelingToolkit.generate_factorized_W(de)
+ModelingToolkit.calculate_factorized_W(de)
 rober = ODEFunction(de, [y₁,y₂,y₃], [k₁,k₂,k₃])
 
 """
@@ -163,7 +163,7 @@ eqs = [D(y₁) ~ I₁*y₂*y₃,
        D(y₃) ~ I₃*y₁*y₂]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-#ModelingToolkit.generate_factorized_W(de)
+ModelingToolkit.calculate_factorized_W(de)
 rigid = ODEFunction(de, [y₁,y₂,y₃], [I₁,I₂,I₃])
 
 """
@@ -292,7 +292,7 @@ eqs = [D(y1) ~ -p1*y1 + p2*y2 + p3*y3 + p4,
        D(y8) ~ -p10*y6*y8 + p12*y7]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-#ModelingToolkit.generate_factorized_W(de)
+ModelingToolkit.calculate_factorized_W(de)
 hires = ODEFunction(de, [y1,y2,y3,y4,y5,y6,y7,y8],
                         [p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12])
 
@@ -337,7 +337,7 @@ eqs = [D(y1) ~ p1*(y2+y1*(1-p2*y1-y2)),
        D(y3) ~ p3*(y1-y3)]
 de = ODESystem(eqs)
 jac = calculate_jacobian(de)
-#ModelingToolkit.generate_factorized_W(de)
+ModelingToolkit.calculate_factorized_W(de)
 orego = ODEFunction(de, [y1,y2,y3], [p1,p2,p3])
 
 """