Fix tests

Author: YingboMa

src/structural_transformation/StructuralTransformations.jl

@@ -22,7 +22,7 @@ using ModelingToolkit: ODESystem, AbstractSystem, var_from_nested_derivative, Di
                        get_postprocess_fbody, vars!,
                        IncrementalCycleTracker, add_edge_checked!, topological_sort,
                        invalidate_cache!, Substitutions, get_or_construct_tearing_state,
-                       AliasGraph
+                       AliasGraph, filter_kwargs
 
 using ModelingToolkit.BipartiteGraphs
 import .BipartiteGraphs: invview

src/structural_transformation/codegen.jl

@@ -384,7 +384,6 @@ function build_observed_function(state, ts, var_eq_matching, var_sccs,
 
     fullvars = state.fullvars
     s = state.structure
-    graph = s.graph
     solver_states = fullvars[is_solver_state_idxs]
     algvars = fullvars[.!is_solver_state_idxs]
 
@@ -525,6 +524,7 @@ function ODAEProblem{iip}(sys,
     has_difference = any(isdifferenceeq, eqs)
     cbs = process_events(sys; callback, has_difference, kwargs...)
 
+    kwargs = filter_kwargs(kwargs)
     if cbs === nothing
         ODEProblem{iip}(fun, u0, tspan, p; kwargs...)
     else

src/structural_transformation/symbolics_tearing.jl

@@ -591,12 +591,12 @@ function tearing_reassemble(state::TearingState, var_eq_matching; simplify = fal
     end
     deleteat!(oldobs, sort!(removed_obs))
     @set! sys.substitutions = Substitutions(subeqs, deps)
-    @set! state.sys = sys
-    @set! sys.tearing_state = state
 
     der2expr = Dict(eq.lhs => eq.rhs for eq in equations(sys) if isdiffeq(eq))
     obs = substitute.([oldobs; subeqs], (der2expr,))
     @set! sys.observed = obs
+    @set! state.sys = sys
+    @set! sys.tearing_state = state
     return invalidate_cache!(sys)
 end