@@ -4,26 +4,25 @@ mutable struct DiffEqSystem <: AbstractSystem
44 dvs:: Vector{Variable}
55 ps:: Vector{Variable}
66 jac:: Matrix{Expression}
7- function DiffEqSystem (eqs, ivs, dvs, ps, jac )
7+ function DiffEqSystem (eqs, ivs, dvs, ps)
88 all (! isintermediate, eqs) ||
99 throw (ArgumentError (" no intermediate equations permitted in DiffEqSystem" ))
1010
11+ jac = Matrix {Expression} (undef, 0 , 0 )
1112 new (eqs, ivs, dvs, ps, jac)
1213 end
1314end
1415
15- DiffEqSystem (eqs, ivs, dvs, ps) = DiffEqSystem (eqs, ivs, dvs, ps, Matrix {Expression} (undef,0 ,0 ))
16-
1716function DiffEqSystem (eqs)
1817 dvs, = extract_elements (eqs, [_is_dependent])
1918 ivs = unique (vcat ((dv. dependents for dv ∈ dvs). .. ))
2019 ps, = extract_elements (eqs, [_is_parameter (ivs)])
21- DiffEqSystem (eqs, ivs, dvs, ps, Matrix {Expression} (undef, 0 , 0 ) )
20+ DiffEqSystem (eqs, ivs, dvs, ps)
2221end
2322
2423function DiffEqSystem (eqs, ivs)
2524 dvs, ps = extract_elements (eqs, [_is_dependent, _is_parameter (ivs)])
26- DiffEqSystem (eqs, ivs, dvs, ps, Matrix {Expression} (undef, 0 , 0 ) )
25+ DiffEqSystem (eqs, ivs, dvs, ps)
2726end
2827
2928isintermediate (eq:: Equation ) = ! (isa (eq. lhs, Operation) && isa (eq. lhs. op, Differential))
0 commit comments