-
-
Notifications
You must be signed in to change notification settings - Fork 57
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Broadcast over VectorOfArray
doesn't respect dimensions of non-vector parent arrays
#373
Comments
I think that's the proper way for it to be handled. |
Closed by #374 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Is your feature request related to a problem? Please describe.
VectorOfArray
was extended to multi-dimensional parent arrays in #357. However, these arrays don't work with all solvers in OrdinaryDiffEq.jl (for example,SSPRK43
).The error happens during the call to
integrator = init(ode, SSPRK43())
, specifically this line: https://github.com/SciML/OrdinaryDiffEq.jl/blob/1f2e058c14d2a5c997bf7d5f055cebfb406cbc1b/src/initdt.jl#L51. An MWE isThe manually truncated stacktrace is
Describe the solution you’d like
The error above is due to the fact that for
u::VectorOfArray
,zero.(u)
returns aVectorOfArray
whose parent array is aVector
. ForVectorOfArray
with multi-dimensional parents, the broadcast should ideally respect the dimension of the parent array.Additional context
I'm looking at this now, but I'm not the most experienced with broadcast. If anyone has any suggestions on where to start and/or how to implement, I'd welcome them.
The text was updated successfully, but these errors were encountered: