Efficiently Finding Maximum Function Value in a Mesh with FEniCSx

Thank you for sharing the updated link! It’s really helpful for understanding the finite element definitions in more detail.

I’ve been thinking about this and trying to start the implementation, but I find it really difficult to even begin the process—I feel stuck at the very first step. From my understanding, the first step involves combining the coefficients with the basis functions on the reference element to form the target function. Is that correct? I also tried searching the community for similar issues and came across some potentially helpful posts, like Accessing the basis function of a Finite Element, but I’m not sure if they are relevant or necessary for this optimization problem.

Could I ask for your further guidance on this? Thanks again for your patience and support—I really appreciate it.