numpy 1.5 does not support the  `out=` keyword argument to ufuncs. I have [this PR](https://github.com/scipy/scipy/pull/3163) stalled because of that. @WarrenWeckesser commented there that if you leave the `out=` out and make it a positional argument it works, but I haven't been able to test it.