geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
217 stars 232 forks source link

Enable limiter by default for bilinear least squares #5939

Open gassmoeller opened 1 week ago

gassmoeller commented 1 week ago

This is one of the default changes I would like to happen for #5568. The least squares interpolators are prone to over- and undershoots and for most applications (i.e. if particles track composition) this is unwanted. The quadratic least squares interpolator already enables the limiter by default. It is still possible to disable the limiter in case the user is sure the tracked property is smooth and accuracy is more important than boundedness of the interpolated property.