See these for the ideas
1d and 2d cartesian: https://doi.org/10.1090/mcom/3998
On triangles: https://arxiv.org/abs/2409.09620
They use smoothness indicator based on jumps in solution and its derivatives. Computing higher order derivatives is expensive on unstructured grids.
Implement this idea but using only jump in solution and first derivative for defining smoothness indicator.
Do this in system_lagrange_mpi version of the code.
See these for the ideas
1d and 2d cartesian: https://doi.org/10.1090/mcom/3998
On triangles: https://arxiv.org/abs/2409.09620
They use smoothness indicator based on jumps in solution and its derivatives. Computing higher order derivatives is expensive on unstructured grids.
Implement this idea but using only jump in solution and first derivative for defining smoothness indicator.
Do this in
system_lagrange_mpiversion of the code.