> Because using |add.reduce(w, axis=axis)| is generally faster than 
|w.sum(axis)|

This is clearly a bug in Numpy. The code in Numpy does call 
`add.reduce`, but this call goes through a callback back to Python code 
(_methods.py)!