[SciPy-User] raising a matrix to float power

Alexey Brazhe brazhe at gmail.com
Sat Jul 10 19:57:02 EDT 2010


>Sure, M**0.5 is cho_factor(M). For other non-integers I am not sure what
>matrix exponentiation could possibly mean.

>Are you sure you don't mean array exponentiation?

Indeed, I needed to raise a matrix (not array) to power 1/2 (in fact, -1/2).
More precisely, I need to compute W(W^TW)^(-1/2).
cho_factor fails with "matrix not positive definite", and I don't know how
to avoid that

A.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20100711/2b59f339/attachment.html>


More information about the SciPy-User mailing list