L^∞ Norm

Hi Everyone,

is there a good way to compute the (element-wise) L^∞-norm of a GridFunction? I need this as part of a stabilising bilinear form.

Best wishes,
Henry

Hi Henry,
I think that point evaluation and some finger crossing is the easiest way. For efficient point evaluation have a look at this thread:
need help to implement the bound limiter (again...)
If you need more accuracy you could combine this with an optimizer from scipy

Best
Christopher

Hi Christopher,
Thank you very much! I will go and have a look at that !
Best wishes,
Henry