Ok, this no trivial,
because min, max are not integral, and le max of P2 function is on the node.
So the compute this first build the sub-domain with trunc., when build a finite element function on this mesh interpolate and take the min, max of dof array.