[SciPy-user] Matrix square root

Pearu Peterson pearu at scipy.org
Mon Oct 4 10:59:33 CDT 2004



On Mon, 4 Oct 2004, Nils Wagner wrote:

> Pearu Peterson wrote:
>> On Mon, 4 Oct 2004, Nils Wagner wrote:
>> 
>>> AFAIK there is no direct support for sqrtm but you may use funm.
>> 
>> I just discovered that there might be a bug in funm. I have added
>> a test to scipy.linalg.test_matfuncs that demonstrates the bug.
>> So, be aware to check the results of funm(..,sqrt).
>> 
> That might give a reason to implement sqrtm seperately.
> What is your opinion ?

The fact that funm is buggy does not give a good reason to implement 
sqrtm.

The following tasks are independent: funm must be fixed and sqrtm must be 
implemented.

Pearu



More information about the SciPy-user mailing list