L^∞ Norm

Hi Everyone,

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

Best wishes,

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


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