Jesse Chan

@DanielDoehring @sloede parabolic mortars are on me and @apey236's list, but we were going to focus on MPI parallelization first. If you are interested in adding mortars for parabolic terms,...

Hi @DanielDoehring - yes, I implemented this originally in `DGMulti`, with @sloede and @ranocha's help for `TreeMesh`. Let me take a look at that PR and see if I can...

Can you explain briefly what makes it hard to use SplitODEProblem with AMR?

Bumping to add that we should do the same for `DGMultiMesh` constructors. For example, `cells_per_dimension` is a keyword argument here ```julia DGMultiMesh(dg::DGMulti; cells_per_dimension, coordinates_min=(-1.0, -1.0), coordinates_max=(1.0, 1.0), is_on_boundary=nothing, periodicity=ntuple(_ ->...

In case it's still helpful, we now have (somewhat experimental) support for parabolic terms. See https://trixi-framework.github.io/Trixi.jl/stable/tutorials/adding_new_parabolic_terms/ and the tutorial before it.

Does `Threads.@threads :static` help at all?

As a for loop annotation - e.g., `Threads.@threads :static for ...`. It seemed to do pretty well compared to Polyester in some of my experiments.

I probably won't have time to finish this before then, so I'm happy to wait until the next release.

Yes - we should be able to wrap a array of arrays `u` with `VectorOfArray(u)` to get around OrdinaryDiffEq.jl dropping support for arrays of arrays. We can unwrap `u` via...

There's one more issue that needs to be resolved upstream related to broadcasting over `VectorOfArray` with multidimensional parent arrays https://github.com/SciML/RecursiveArrayTools.jl/issues/373