Jesse Chan

A note - I've run into memory errors related to this limitation when the size of inner(0) > 256.

no worries. I don't need it at the moment, but would it be useful to just add a warning flag during the OKL build?

Note to self: the issue disappears when specifying an initial `dt` and `adaptive=false`. If an initial `dt` isn't specified, the error is in https://github.com/SciML/OrdinaryDiffEq.jl/blob/b15f3c53601c274a69216ffea242a46561bbe121/src/initdt.jl#L23 ```julia @.. broadcast=false sk = abstol+internalnorm(u0,t)*reltol...

Added it, but same error. Should that have done something, or were you just making a note to yourself?

The tests really take _way_ too long to run IMO. It's so far been 15 minutes and my tests are still running. It would be nice to have a minimal...

One possibility would be to create a testset intended for "local" testing, which could exclude some of the CI tests.

`ODESolution` is just OrdinaryDiffEq.jl's solution type, so it doesn't know anything about the spatial coordinates `x`. That information is usually stored in `semi.cache` (I believe specifically in `semi.cache.elements.node_coordinates`). Trixi also...

I should note that the isosurface looks significantly fancier even if it's just colored according to a random quantity, such as the `z`-coordinate. Right now, the isosurfaces are all just...

Hm...this is probably due to some default "convert on setindex" behavior introduced by Tim Holy (see https://github.com/JuliaArrays/StructArrays.jl/pull/227 and https://github.com/JuliaArrays/StructArrays.jl/issues/216 for an earlier discussion). My guess is that `PlotData2D` creates a...

We should be able to allow any number of variables (just call `solution_variables` on a single element of `u` to get the number of outputs and initialize `uplot` accordingly). We'd...