Merge pull request #2993 from SebastianM-C/depwarn

Fix depwarns from `solve_for`

Author: ChrisRackauckas

Project.toml

@@ -119,7 +119,7 @@ SpecialFunctions = "0.7, 0.8, 0.9, 0.10, 1.0, 2"
 StaticArrays = "0.10, 0.11, 0.12, 1.0"
 SymbolicIndexingInterface = "0.3.28"
 SymbolicUtils = "3.2"
-Symbolics = "6"
+Symbolics = "6.3"
 URIs = "1"
 UnPack = "0.1, 1.0"
 Unitful = "1.1"

docs/src/examples/perturbation.md

@@ -26,7 +26,7 @@ function solve_coef(eqs, ps)
 
     for i in 1:length(ps)
         eq = substitute(eqs[i], vals)
-        vals[ps[i]] = Symbolics.solve_for(eq ~ 0, ps[i])
+        vals[ps[i]] = Symbolics.symbolic_linear_solve(eq ~ 0, ps[i])
     end
     vals
 end

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)