![](https://secure.gravatar.com/avatar/d41fa6e1fe29e6c5c5821b5a3f31f190.jpg?s=120&d=mm&r=g)
Stuart Murdock wrote:
Hi
I am trying to calculate the matrix square root (lets say B) of a Matrix (A) which is an array.
A = array([[1,2],[3,4]])
such that
A = matrixmultiply(B,B)
are there any functions in SciPy to do that ?
Thanks
Stuart
AFAIK there is no direct support for sqrtm but you may use funm.
a = array(([1,2],[3,4])) b = linalg.funm(a,sqrt) b array([[ 0.55368857+0.46439416j, 0.80696073-0.21242648j], [ 1.21044109-0.31863972j, 1.76412966+0.14575444j]]) dot(b,b)-a array([[ 2.22044605e-16 -3.19839641e-17j, 0.00000000e+00 -2.83518868e-17j], [ 4.44089210e-16 -1.06821019e-16j, 0.00000000e+00 -1.69677640e-17j]])
HTH, Nils -- Dr.-Ing. Nils Wagner Institut A für Mechanik Universität Stuttgart Pfaffenwaldring 9 D-70550 Stuttgart Tel.: (+49) 0711 685 6262 Fax.: (+49) 0711 685 6282 E-mail: nwagner@mecha.uni-stuttgart.de URL : http://www.mecha.uni-stuttgart.de