reefactor: update solve_for to symbolic_linear_solve

Author: SebastianM-C

src/ModelingToolkit.jl

@@ -59,7 +59,7 @@ using RuntimeGeneratedFunctions: drop_expr
 
 using Symbolics: degree
 using Symbolics: _parse_vars, value, @derivatives, get_variables,
-                 exprs_occur_in, solve_for, build_expr, unwrap, wrap,
+                 exprs_occur_in, symbolic_linear_solve, build_expr, unwrap, wrap,
                  VariableSource, getname, variable, Connection, connect,
                  NAMESPACE_SEPARATOR, set_scalar_metadata, setdefaultval,
                  initial_state, transition, activeState, entry, hasnode,

src/structural_transformation/StructuralTransformations.jl

@@ -12,7 +12,7 @@ using SymbolicUtils: maketerm, iscall
 using ModelingToolkit
 using ModelingToolkit: ODESystem, AbstractSystem, var_from_nested_derivative, Differential,
                        unknowns, equations, vars, Symbolic, diff2term, value,
-                       operation, arguments, Sym, Term, simplify, solve_for,
+                       operation, arguments, Sym, Term, simplify, symbolic_linear_solve,
                        isdiffeq, isdifferential, isirreducible,
                        empty_substitutions, get_substitutions,
                        get_tearing_state, get_iv, independent_variables,

src/structural_transformation/symbolics_tearing.jl

@@ -141,7 +141,7 @@ function tearing_assignments(sys::AbstractSystem)
 end
 
 function solve_equation(eq, var, simplify)
-    rhs = value(solve_for(eq, var; simplify = simplify, check = false))
+    rhs = value(symbolic_linear_solve(eq, var; simplify = simplify, check = false))
     occursin(var, rhs) && throw(EquationSolveErrors(eq, var, rhs))
     var ~ rhs
 end
@@ -457,7 +457,7 @@ function tearing_reassemble(state::TearingState, var_eq_matching,
                 dx = D(simplify_shifts(lower_varname_withshift(
                     fullvars[lv], idep, order - 1)))
                 eq = dx ~ simplify_shifts(Symbolics.fixpoint_sub(
-                    Symbolics.solve_for(neweqs[ieq],
+                    Symbolics.symbolic_linear_solve(neweqs[ieq],
                         fullvars[iv]),
                     total_sub; operator = ModelingToolkit.Shift))
                 for e in 𝑑neighbors(graph, iv)