[Nipy-devel] Affine orientation ordering - please review and test

Jonathan Taylor jonathan.taylor@stanford....
Tue Jan 19 14:15:57 CST 2010


>
> -    RS = RZS / zooms
> +    RS = np.dot(RZS, np.diag(zooms))
>
>
Yes, it is easy to explain -- it's a bug and should be reversed. Sorry about
that.
It is meant to be np.diag(1/zooms).

I used the np.dot
just because I wasn't sure that in the non-square case broadcasting would
work (it does, I've now checked).

Interesting that none of the tests picked up this bug, I suppose one would
have to exaggerate one row severely.


-- 
Jonathan Taylor
Dept. of Statistics
Sequoia Hall, 137
390 Serra Mall
Stanford, CA 94305
Tel:   650.723.9230
Fax:   650.725.8977
Web: http://www-stat.stanford.edu/~jtaylo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.scipy.org/pipermail/nipy-devel/attachments/20100119/a2a7feda/attachment.html 


More information about the Nipy-devel mailing list