In a problem of getting the maximum of a multi-variable function, I would use mpi::reduce and provide it with a function object that compares the evaluations of this function.

Each process evaluates a part of the space, ie a number of n-tuples of the n-dimension space.

How do I get the n-tuple that yields the max of f(), given that MPI_MAXLOC is not available in boost mpi?

Thanks