Hello! I currently have a working Helmholtz solution operator (A) on my nice custom mesh. I’m now looking for the most efficient way to create a CoefficientFunction on the form
Norm(A(f(x,y)))
where f for each x,y in my mesh is itself a different CoefficientFunction (e.g. f could be a delta-Dirac or thin Gaussian centered in (x,y)). I’m then interested in visualising the norm of each solution.
I found an old thread that seems to want the same thing - Integral over product space - Kunena - but it is a few years old, and I was wondering if there is a modern, canonical way to do this, ideally without having to learn C++?
Best,
Christian