solve_for->solve, some docs

Author: shashi

src/solve.jl

@@ -56,8 +56,16 @@ function A_b(eqs, vars)
     A, b
 end
 
-function solve_for(eqs, vars)
+"""
+    solve(eqs::Vector, vars::Vector)
 
+Solve the vector of equations `eqs` for a set of variables `vars`.
+
+Assumes `length(eqs) == length(vars)`
+
+Currently only works if all equations are linear.
+"""
+function solve(eqs, vars)
     A, b = A_b(eqs, vars)
     A = SymbolicUtils.simplify.(to_symbolic.(A), polynorm=true)
     b = SymbolicUtils.simplify.(to_symbolic.(b), polynorm=true)

src/systems/reduction.jl

@@ -54,7 +54,7 @@ function alias_elimination(sys::ODESystem)
     eliminate = setdiff(convert.(Variable, all_vars), newstates)
 
     vars = map(x->x(sys.iv()), eliminate)
-    outputs = solve_for(eqs[alg_idxs], vars)
+    outputs = solve(eqs[alg_idxs], vars)
 
     diffeqs = eqs[setdiff(1:length(eqs), alg_idxs)]