Re: [ublas] How to do sqrt on matrix

From: Georg Baum (Georg.Baum_at_[hidden])
Date: 2007-09-10 10:27:52

Am Montag 10 September 2007 16:21 schrieb Paul C. Leopardi:

> The second is a matrix square root, ie. M1 == sqrt(M2) <=> M1 * M1 == M2.
> This corresponds to the Matlab sqrtm function.
> See
> http://www.mathworks.com/access/helpdesk/help/techdoc/ref/sqrtm.html

This is actually the more interesting one :-) If you need to implement this, 
have a look at Nick Higham's papers on matrix funtions: 
http://www.maths.manchester.ac.uk/~higham/pap-mf.html

The matlab sqrtm is (or was, I don't know about recent versions) based on 
his code.

Georg