OK, I know this is a damn tiny corner cases... If we already have inf in this kind of function, the argmax may return the index to the NaN instead of the inf. Getting around it is likely tedious (would require getting the `mask.all(axis=axis)` by looking at the value of the mask `ind` points to instead), but should we note it in the docs?