get them all

Author: ChrisRackauckas

src/discretization/generate_ic_defaults.jl

@@ -11,7 +11,7 @@ function PDEBase.generate_ic_defaults(tconds, s::DiscreteSpace,
             indexmap = Dict([args[i] => i for i in 1:length(args)])
             D = ic.order == 0 ? identity : (Differential(t)^ic.order)
             defaultvars = D.(s.discvars[depvar(ic.u, s)])
-            broadcastable_rhs = [solve_for(ic.eq, D(ic.u))]
+            broadcastable_rhs = [symbolic_linear_solve(ic.eq, D(ic.u))]
             out = substitute.(
                 broadcastable_rhs, valmaps(s, depvar(ic.u, s), ic.depvars, indexmap))
             vec(defaultvars .=>

src/system_parsing/pde_system_transformation.jl

@@ -271,7 +271,7 @@ function generate_bc_rules(bcs, v)
         deriv = bc.order == 0 ? identity : (Differential(bc.x)^bc.order)
         bcrule_lhs = deriv(operation(bc.u)(v.args[operation(bc.u)]...))
         bcterm = deriv(bc.u)
-        rhs = solve_for(bc.eq, bcterm)
+        rhs = symbolic_linear_solve(bc.eq, bcterm)
         bcrule_lhs => rhs
     end
 end