Merge pull request #1349 from SciML/myb/fix

Make sure `x in AliasGraphKeySet` is in bounds

Author: YingboMa

.github/workflows/Downstream.yml

@@ -24,6 +24,7 @@ jobs:
           - {user: SciML, repo: NeuralPDE.jl, group: NNPDE}
           - {user: SciML, repo: DataDrivenDiffEq.jl, group: Standard}
           - {user: SciML, repo: StructuralIdentifiability.jl, group: All}
+          - {user: SciML, repo: ModelingToolkitStandardLibrary.jl}
     steps:
       - uses: actions/checkout@v2
       - uses: julia-actions/setup-julia@v1

Project.toml

@@ -1,7 +1,7 @@
 name = "ModelingToolkit"
 uuid = "961ee093-0014-501f-94e3-6117800e7a78"
 authors = ["Chris Rackauckas <[email protected]>"]
-version = "7.1.1"
+version = "7.1.2"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"

src/systems/alias_elimination.jl

@@ -181,7 +181,10 @@ struct AliasGraphKeySet <: AbstractSet{Int}
 end
 Base.keys(ag::AliasGraph) = AliasGraphKeySet(ag)
 Base.iterate(agk::AliasGraphKeySet, state...) = Base.iterate(agk.ag.eliminated, state...)
-Base.in(i::Int, agk::AliasGraphKeySet) = agk.ag.aliasto[i] !== nothing
+function Base.in(i::Int, agk::AliasGraphKeySet)
+    aliasto = agk.ag.aliasto
+    1 <= i <= length(aliasto) && aliasto[i] !== nothing
+end
 
 count_nonzeros(a::AbstractArray) = count(!iszero, a)