@@ -10,7 +10,7 @@ function calculate_tgrad(sys::AbstractODESystem;
10
10
xs = states (sys)
11
11
rule = Dict (map ((x, xt) -> xt=> x, detime_dvs .(xs), xs))
12
12
rhs = substitute .(rhs, Ref (rule))
13
- tgrad = [expand_derivatives (ModelingToolkit . Differential (iv)(r), simplify) for r in rhs]
13
+ tgrad = [expand_derivatives (Differential (iv)(r), simplify) for r in rhs]
14
14
reverse_rule = Dict (map ((x, xt) -> x=> xt, detime_dvs .(xs), xs))
15
15
tgrad = Num .(substitute .(tgrad, Ref (reverse_rule)))
16
16
get_tgrad (sys)[] = tgrad
@@ -102,8 +102,8 @@ function generate_function(
102
102
p = map (x-> time_varying_as_func (value (x), sys), ps)
103
103
t = get_iv (sys)
104
104
105
- if has_preface (sys)
106
- pre = ex -> Let (preface (sys) , ex)
105
+ if has_preface (sys) && (pre = preface (sys); pre != = nothing )
106
+ pre = ex -> Let (pre , ex)
107
107
else
108
108
pre = ex -> ex
109
109
end
@@ -415,13 +415,13 @@ function ODEFunctionExpr{iip}(sys::AbstractODESystem, dvs = states(sys),
415
415
=#
416
416
417
417
fsym = gensym (:f )
418
- _f = :($ fsym = ModelingToolkit . ODEFunctionClosure ($ f_oop, $ f_iip))
418
+ _f = :($ fsym = $ ODEFunctionClosure ($ f_oop, $ f_iip))
419
419
tgradsym = gensym (:tgrad )
420
420
if tgrad
421
421
tgrad_oop, tgrad_iip = generate_tgrad (sys, dvs, ps;
422
422
simplify= simplify,
423
423
expression= Val{true }, kwargs... )
424
- _tgrad = :($ tgradsym = ModelingToolkit . ODEFunctionClosure ($ tgrad_oop, $ tgrad_iip))
424
+ _tgrad = :($ tgradsym = $ ODEFunctionClosure ($ tgrad_oop, $ tgrad_iip))
425
425
else
426
426
_tgrad = :($ tgradsym = nothing )
427
427
end
@@ -431,7 +431,7 @@ function ODEFunctionExpr{iip}(sys::AbstractODESystem, dvs = states(sys),
431
431
jac_oop,jac_iip = generate_jacobian (sys, dvs, ps;
432
432
sparse= sparse, simplify= simplify,
433
433
expression= Val{true }, kwargs... )
434
- _jac = :($ jacsym = ModelingToolkit . ODEFunctionClosure ($ jac_oop, $ jac_iip))
434
+ _jac = :($ jacsym = $ ODEFunctionClosure ($ jac_oop, $ jac_iip))
435
435
else
436
436
_jac = :($ jacsym = nothing )
437
437
end
0 commit comments