more symbolics removal

Author: ChrisRackauckas

Project.toml

@@ -22,7 +22,6 @@ SparseArrays = "2f01184e-e22b-5df5-ae63-d93ebab69eaf"
 SparseDiffTools = "47a9eef4-7e08-11e9-0b38-333d64bd3804"
 StaticArrays = "90137ffa-7385-5640-81b9-e52037218182"
 SuiteSparse = "4607b0f0-06f3-5cda-b6b1-a6196a1729e9"
-SymbolicUtils = "d1185830-fcd6-423d-90d6-eec64667417b"
 
 [compat]
 BandedMatrices = "0.15.11, 0.16"
@@ -40,7 +39,6 @@ RuntimeGeneratedFunctions = "0.4, 0.5"
 SciMLBase = "1.11"
 SparseDiffTools = "1.17"
 StaticArrays = "0.10, 0.11, 0.12, 1.0"
-SymbolicUtils = "0.11, 0.12, 0.13, 0.14, 0.15, 0.16, 0.17, 0.18"
 julia = "1.6"
 
 [extras]

src/utils.jl

@@ -1,5 +1,3 @@
-using SymbolicUtils
-
 """
 A function that creates a tuple of CartesianIndices of unit length and `N` dimensions, one pointing along each dimension.
 """
@@ -22,61 +20,6 @@ function cartesian_to_linear(I::CartesianIndex, s)  #Not sure if there is a buil
     return out
 end
 
-# Counts the Differential operators for given variable x. This is used to determine
-# the order of a PDE.
-function count_differentials(term, x::Symbolics.Symbolic)
-    S = Symbolics
-    SU = SymbolicUtils
-    if !S.istree(term)
-        return 0
-    else
-        op = SU.operation(term)
-        count_children = sum(map(arg -> count_differentials(arg, x), SU.arguments(term)))
-        if op isa Differential && op.x === x
-            return 1 + count_children
-        end
-        return count_children
-    end
-end
-
-# return list of differential orders in the equation
-function differential_order(eq, x::Symbolics.Symbolic)
-    S = Symbolics
-    SU = SymbolicUtils
-    orders = Set{Int}()
-    if S.istree(eq)
-        op = SU.operation(eq)
-        if op isa Differential
-            push!(orders, count_differentials(eq, x))
-        else
-            for o in map(ch -> differential_order(ch, x), SU.arguments(eq))
-                union!(orders, o)
-            end
-        end
-    end
-    return filter(!iszero, orders)
-end
-
-# find all the dependent variables given by depvar_ops in an expression
-function get_depvars(eq,depvar_ops)
-    S = Symbolics
-    SU = SymbolicUtils
-    depvars = Set()
-    if eq isa Num
-       eq = eq.val
-    end
-    if S.istree(eq)
-        if eq isa Term && any(u->isequal(operation(eq),u),depvar_ops)
-              push!(depvars, eq)
-        else
-            for o in map(x->get_depvars(x,depvar_ops), SU.arguments(eq))
-                union!(depvars, o)
-            end
-        end
-    end
-    depvars
-end
-
 add_dims(A::AbstractArray, n::Int; dims::Int = 1) = cat(ndims(A) + n, A, dims = dims)
 
 ""