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?