Hi,<br>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. <br><br>I am implementing the algorithm as described here:<br>http://books.nips.cc/papers/files/nips14/AA35.pdf<br><br>Finding the square root of D^(-1) is way too slow for D^-1 of size (1258,1258).<br><br><br>Thanks.<br>Vivek.<br>