Andrew Winters

As a followup to the discussion today I threw together a quick loop that computes the shortest element side length of a mesh given only the `mesh` and using its...

> For a real iregular unstructured mesh, I wonder if the length of the sides is a good measure for the size "h" of the element. This is somewhat related...

> If possible, it would of course be nice to have something that's somewhat "dimension independent". The Verdict library (and HOHQMesh) also has hex element statistics that we can experiment...

I figured this would be expensive and would likely only be used for debugging purposes. I like the idea of a callback but see that dissecting the stage information from...

Thinking about it more it might also be good to print a solution `h5` file once such an abort is encountered. Then the user can also visualize what is going...

Here is an augmented version of the suggested workaround given above by @ranocha; however, it is changed to operate with `P4estMesh{2}` ```julia function unstable_check(checks, u_ode, semi, t) u = Trixi.wrap_array(u_ode,...

Hi James, These divergence errors are quite large (especially the L∞). The main issue is that there is a bug in the computation of div(B) on the unstructured mesh type...

> That's much better! I've only dealt with Cartesian finite difference grids with divergencelessness baked into the numerical scheme (i.e. like constrained transport) or with codes that use the vector...

The handling of the nonconservative terms in the volume / boundary will be the subject of a future PR. It will probably involve: 1. Writing an alternative `split_form_kernel!` that uses...

I need to look back at some experimenting I did on this subject for the shallow water equations. From what I remember, using the standard `derivative_matrix` would allow us to...