I was interested so I tried this and in my case it just crashed Python. I think that the crash might be related to an issue that posted earlier related to numpy 1.7.1. In fact any of the linalg functions in scipy 0.12.0 seem to cause Python to crash on my computer.

On Saturday, April 13, 2013 8:47:25 AM UTC-7, Vivek Kulkarni wrote:

Hi,
I am implementing spectral clustering for my course work and am using the sqrtm function to find the square root of a matrix . But its far too slow, my matrix is of size (1258,1258). Any suggestions on how I can speed things up or some other function that scipy supports which can be used.

I am implementing the algorithm as described here:
http://books.nips.cc/papers/files/nips14/AA35.pdf

Finding the square root of D^(-1) is way too slow for D^-1 of size (1258,1258).

Thanks.
Vivek.