def mpi_max(value): global_max = np.zeros(1, dtype='float64') local_max = np.max(value).astype('float64') MPI.COMM_WORLD.Reduce(local_max, global_max, op=MPI.MAX) return global_max[0]