Raise if we try to construct a reduction over an axis >= ndim
Tests for nelements over 3d arrays
Implement general nelements reduction on numpy arrays
We already raise in construction