Is there an easy method to compute the maximum of u[] in parallel code?
u[]
fespace Vh(Th, P2); Vh u;
fespace Vh(Th, P2);
Vh u;
Yes using mpiMAX, cf. Parallelization.
mpiMAX