fix: use is_diff_equation instead of isdiffeq when finding algeeqs

Author: vyudu

src/systems/callbacks.jl

@@ -295,7 +295,6 @@ function make_affect(affect::Vector{Equation}; iv = nothing, algeeqs::Vector{Equ
     for u in dvs
         aff_map[u] = u
     end
-    @show explicit
 
     AffectSystem(affectsys, collect(dvs), params, discretes, aff_map, explicit)
 end
@@ -842,7 +841,6 @@ Compile an affect defined by a set of equations. Systems with algebraic equation
 function compile_equational_affect(aff::Union{AffectSystem, Vector{Equation}}, sys; reset_jumps = false, kwargs...)
     if aff isa AbstractVector
         aff = make_affect(aff; iv = get_iv(sys))
-        @show is_explicit(aff)
     end
     affsys = system(aff)
     ps_to_update = discretes(aff)

src/systems/diffeqs/odesystem.jl

@@ -305,7 +305,7 @@ function ODESystem(deqs::AbstractVector{<:Equation}, iv, dvs, ps;
         throw(ArgumentError("System names must be unique."))
     end
 
-    algeeqs = filter(eq -> eq.lhs isa Union{Symbolic, Number} && !isdiffeq(eq), deqs)
+    algeeqs = filter(eq -> eq.lhs isa Union{Symbolic, Number} && !is_diff_equation(eq), deqs)
     cont_callbacks = SymbolicContinuousCallbacks(continuous_events; algeeqs, iv)
     disc_callbacks = SymbolicDiscreteCallbacks(discrete_events; algeeqs, iv)
 

src/systems/diffeqs/sdesystem.jl

@@ -264,7 +264,7 @@ function SDESystem(deqs::AbstractVector{<:Equation}, neqs::AbstractArray, iv, dv
     Wfact = RefValue(EMPTY_JAC)
     Wfact_t = RefValue(EMPTY_JAC)
 
-    algeeqs = filter(eq -> eq.lhs isa Union{Symbolic, Number} && !isdiffeq(eq), deqs)
+    algeeqs = filter(eq -> eq.lhs isa Union{Symbolic, Number} && !is_diff_equation(eq), deqs)
     cont_callbacks = SymbolicContinuousCallbacks(continuous_events; algeeqs, iv)
     disc_callbacks = SymbolicDiscreteCallbacks(discrete_events; algeeqs, iv)
     if is_dde === nothing

src/systems/model_parsing.jl

@@ -1115,7 +1115,6 @@ function parse_equations!(exprs, eqs, dict, body)
 end
 
 function parse_continuous_events!(c_evts, dict, body)
-    @show body
     dict[:continuous_events] = []
     Base.remove_linenums!(body)
     for arg in body.args
@@ -1125,7 +1124,6 @@ function parse_continuous_events!(c_evts, dict, body)
 end
 
 function parse_discrete_events!(d_evts, dict, body)
-    @show body
     dict[:discrete_events] = []
     Base.remove_linenums!(body)
     for arg in body.args

test/initializationsystem.jl

@@ -1280,9 +1280,9 @@ end
     @parameters β γ S0
     @variables S(t)=S0 I(t) R(t)
     rate₁ = β * S * I
-    affect₁ = [S ~ S - 1, I ~ I + 1]
+    affect₁ = [S ~ Pre(S) - 1, I ~ Pre(I) + 1]
     rate₂ = γ * I
-    affect₂ = [I ~ I - 1, R ~ R + 1]
+    affect₂ = [I ~ Pre(I) - 1, R ~ Pre(R) + 1]
     j₁ = ConstantRateJump(rate₁, affect₁)
     j₂ = ConstantRateJump(rate₂, affect₂)
     j₃ = MassActionJump(2 * β + γ, [R => 1], [S => 1, R => -1])