Use interface defined in SymbolicIndexingInterface

- `symbolic_indexing_interface.jl` moved to its own package
- Update tests accordingly

Author: AayushSabharwal

Project.toml

@@ -16,6 +16,7 @@ LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 RecipesBase = "3cdcf5f2-1ef4-517c-9805-6587b60abb01"
 StaticArraysCore = "1e83bf80-4336-4d27-bf5d-d5a4f845583c"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
+SymbolicIndexingInterface = "2efcf032-c050-4f8e-a9bb-153293bab1f5"
 Tables = "bd369af6-aec1-5ad0-b16a-f7cc5008161c"
 ZygoteRules = "700de1a5-db45-46bc-99cf-38207098b444"
 
@@ -30,6 +31,7 @@ GPUArraysCore = "0.1"
 IteratorInterfaceExtensions = "1"
 RecipesBase = "0.7, 0.8, 1.0"
 StaticArraysCore = "1.1"
+SymbolicIndexingInterface = "0.1"
 Tables = "1"
 ZygoteRules = "0.2"
 julia = "1.6"

src/RecursiveArrayTools.jl

@@ -7,6 +7,7 @@ module RecursiveArrayTools
 using DocStringExtensions
 using RecipesBase, StaticArraysCore, Statistics,
       ArrayInterfaceCore, LinearAlgebra
+using SymbolicIndexingInterface
 
 import ChainRulesCore
 import ChainRulesCore: NoTangent
@@ -37,9 +38,6 @@ import GPUArraysCore
 Base.convert(T::Type{<:GPUArraysCore.AbstractGPUArray}, VA::AbstractVectorOfArray) = T(VA)
 ChainRulesCore.rrule(T::Type{<:GPUArraysCore.AbstractGPUArray}, xs::AbstractVectorOfArray) = T(xs), ȳ -> (NoTangent(),ȳ)
 
-export independent_variables, is_indep_sym, states, state_sym_to_index, is_state_sym,
-       parameters, param_sym_to_index, is_param_sym, SymbolCache
-
 export VectorOfArray, DiffEqArray, AbstractVectorOfArray, AbstractDiffEqArray,
        AllObserved, vecarr_to_arr, vecarr_to_vectors, tuples
 

src/symbolic_indexing_interface.jl

@@ -1,32 +1,5 @@
 using RecursiveArrayTools, Test
 
-sc = SymbolCache(nothing, nothing, nothing)
-@test isempty(independent_variables(sc))
-@test !is_indep_sym(sc, :a)
-@test isempty(states(sc))
-@test isnothing(state_sym_to_index(sc, :a))
-@test !is_state_sym(sc, :a)
-@test isempty(parameters(sc))
-@test isnothing(param_sym_to_index(sc, :a))
-@test !is_param_sym(sc, :a)
-
-sc = SymbolCache([:a, :b], [:t], [:c, :d])
-@test independent_variables(sc) == [:t]
-@test is_indep_sym(sc, :t)
-@test !is_indep_sym(sc, :a)
-@test states(sc) == [:a, :b]
-@test state_sym_to_index(sc, :a) == 1
-@test state_sym_to_index(sc, :b) == 2
-@test isnothing(state_sym_to_index(sc, :t))
-@test all(is_state_sym.((sc,), [:a, :b]))
-@test !is_state_sym(sc, :c)
-@test parameters(sc) == [:c, :d]
-@test param_sym_to_index(sc, :c) == 1
-@test param_sym_to_index(sc, :d) == 2
-@test isnothing(param_sym_to_index(sc, :a))
-@test all(is_param_sym.((sc,), [:c, :d]))
-@test !is_param_sym(sc, :b)
-
 t = 0.0:0.1:1.0
 f(x) = 2x
 f2(x) = 3x