[Numpy-discussion] Matrix square root

Andrew Nesbit alnesbit at students.cs.mu.oz.au
Sun Sep 7 03:50:07 EDT 2003

Bernard Frankpitt <frankpit at erols.com> writes:

[snip... finding matrix square root through eigenvalue decomposition.]

> The only problem with this approach is that you have to find all the
> eigenvalues and eigenvectors to get the spectral decomposition.  For
> large, poorly conditioned matrices this can be a source of significant
> errors, and a lot of work. An alternative, and often faster and more
> accurate approach is to use a series expansion for the square
> root. Evaluating the series expansion only requires multiplication.

This technique sounds like it's definitely worth checking out.

Thankyou for the tip.


More information about the NumPy-Discussion mailing list