Good question and I would also like to see any hints to improvements (although that’s probably a really hard problem). We experienced something similar in our hyperbolic PDE solver framework Trixi.jl, cf. this Discourse thread.
1 Like