:returns: the matrix exponential
D, V = scipy.linalg.eigh(Ci)
D = numpy.diag(numpy.exp(D))
Out = numpy.dot(numpy.dot(V, D), V.T)
return Out
After Change
:returns: the matrix exponential
return _matrix_operator(Ci, numpy.exp)
def invsqrtm(Ci):
Return the inverse matrix square root of a covariance matrix defined by :
