global mesh.hmax()?

Question

I noticed that mesh.hmax() is a local operation in that every process returns the maximum diameter in its part of the mesh. At least the documentation should be clarified here.

Anyways, is there a way to form the global maximum from within Python?

Answer 1

hmax = MPI.max(mesh.hmax())

Comments

In DOLFIN 1.4, this will become

hmax = MPI.max(mesh.mpi_comm(), mesh.hmax())