Respect solvability options

Author: YingboMa

src/structural_transformation/utils.jl

@@ -231,7 +231,7 @@ function find_solvables!(state::TearingState; kwargs...)
     return nothing
 end
 
-function linear_subsys_adjmat!(state::TransformationState)
+function linear_subsys_adjmat!(state::TransformationState; kwargs...)
     graph = state.structure.graph
     if state.structure.solvable_graph === nothing
         state.structure.solvable_graph = BipartiteGraph(nsrcs(graph), ndsts(graph))
@@ -243,7 +243,7 @@ function linear_subsys_adjmat!(state::TransformationState)
     coeffs = Int[]
     to_rm = Int[]
     for i in eachindex(eqs)
-        all_int_vars, rhs = find_eq_solvables!(state, i, to_rm, coeffs)
+        all_int_vars, rhs = find_eq_solvables!(state, i, to_rm, coeffs; kwargs...)
 
         # Check if all states in the equation is both linear and homogeneous,
         # i.e. it is in the form of

src/systems/abstractsystem.jl

@@ -1035,7 +1035,7 @@ function structural_simplify(sys::AbstractSystem, io = nothing; simplify = false
     has_io = io !== nothing
     has_io && markio!(state, io...)
     state, input_idxs = inputs_to_parameters!(state, io)
-    sys, ag = alias_elimination!(state)
+    sys, ag = alias_elimination!(state; kwargs...)
     #ag = AliasGraph(length(ag))
     # TODO: avoid construct `TearingState` again.
     #state = TearingState(sys)

src/systems/alias_elimination.jl

@@ -4,8 +4,8 @@ using Graphs.Experimental.Traversals
 
 const KEEP = typemin(Int)
 
-function alias_eliminate_graph!(state::TransformationState)
-    mm = linear_subsys_adjmat!(state)
+function alias_eliminate_graph!(state::TransformationState; kwargs...)
+    mm = linear_subsys_adjmat!(state; kwargs...)
     if size(mm, 1) == 0
         ag = AliasGraph(ndsts(state.structure.graph))
         return ag, mm, ag, mm, BitSet() # No linear subsystems
@@ -48,11 +48,12 @@ function extreme_var(var_to_diff, v, level = nothing, ::Val{descend} = Val(true)
 end
 
 alias_elimination(sys) = alias_elimination!(TearingState(sys))[1]
-function alias_elimination!(state::TearingState)
+function alias_elimination!(state::TearingState; kwargs...)
     sys = state.sys
     complete!(state.structure)
     graph_orig = copy(state.structure.graph)
-    ag, mm, complete_ag, complete_mm, updated_diff_vars = alias_eliminate_graph!(state)
+    ag, mm, complete_ag, complete_mm, updated_diff_vars = alias_eliminate_graph!(state;
+                                                                                 kwargs...)
     isempty(ag) && return sys, ag
 
     fullvars = state.fullvars