baggepinnen
commented
1 month ago This is a bug / poor error message in structural simplify
BoundsError: attempt to access 302-element Vector{Vector{Int64}} at index [303]
[2] 𝑑neighbors
@ C:\Users\Bernt\.julia\packages\ModelingToolkit\7em8Z\src\bipartite_graph.jl:369 [inlined]
[3] 𝑑neighbors
@ C:\Users\Bernt\.julia\packages\ModelingToolkit\7em8Z\src\bipartite_graph.jl:368 [inlined]
[4] check_consistency(state::TearingState{ODESystem}, orig_inputs::Set{Any})looks similar to https://github.com/SciML/ModelingToolkit.jl/issues/2727
This issue can probably be closed in favor of other similar ones, this one does not have a reproducer
B-LIE
ChrisRackauckas
pmc4
When I try to linearize my model, I get an "BoundsError: attempt to access 302-element Vector{Vector{Int64}} at index [303]" error🐞
Here is my Julia/MTK statement with partial error message:
Expected behavior
I expected Julia to respond with a collection of linearized matrices A, B, C, D, etc.
Minimal Reproducible Example 👇
I have not been able to create a minimal reproducible example (MRE). I tried to make a couple of similar examples, but in both of those cases, my code works.
Still, to illustrate what I do, I include a similar but slightly simpler example. In this simplified example, I consider a Telegraph Line module:
which I then cascade into a Telegraph Line system:
I then instantiate this model + run the simulation successfully:
This works fine.
Next, I want to linearize the model. To do so, without using
ModelingToolkitStandardLibraryand Analysis Points, I create a copy of the system model where I comment out the definition of input(s), i.e., I create a version with "no inputs" (_noi):which I then instantiate and do the linearization:
This works (as expected).
My problem is that when I do the exactly same thing for my other, slightly more complex model, (i) doing simulation with the balanced model works, but (ii) when creating a model with inputs commented out, linearization leads to the error message with BoundsError.
Error & Stacktrace ⚠️
Environment (please complete the following information):
using Pkg; Pkg.status()using Pkg; Pkg.status(; mode = PKGMODE_MANIFEST)versioninfo()Additional context
The problem where I get an error message is somewhat too complex to include here. Also, it is part of a research project, and I'd prefer to not include it in public. However, if some MTK developers need to get access to my code, I can provide it via e-mail.