All, On a 2D array, numpy.core.umath.maximum.reduce returns a 1D array (the default axis is 0). An axis=None is not recognized as a valid argument for numpy.core.umath.maximum, unfortunately... Is this the wanted behavior ? Thanks a lot in advance P.