Maximum of a function from DG space

As the title says, I need to find the maximum of a field interpolated to functions from discontinuous galerkin space. I see some old threads where the problem was addressed for continuous functions. Is there a similar solution worked out for discontinuous functions?

Maybe this function from scifem would be useful?

Thank you for the suggestion. The function does not seem to be available in 0.9 which I must use due to dependence on periodic mesh (is that ported to the latest release?) But hope I will be able to adapt it.

What specific periodic mesh function are you using?
Are you using: Parallel periodic mesh · GitHub ?
If so, I just updated it to 0.11

Yes, I am using parallel version of the periodic mesh. Thanks a lot for referring me to the update and of course for making it happen. I will try it right away.