a1be7bf2f0926707e25b02570c17ba373878f768,pymanopt/manifolds/psd.py,PositiveDefinite,exp,#PositiveDefinite#Any#Any#,110

Before Change

        // TODO: Check which method is faster depending on n, k.
        if self._k == 1:
            // Use manopt method
            return x.dot(sp.linalg.expm(sp.linalg.solve(x, u, sym_pos=True)))
        else:
            c = la.cholesky(x)
            c_inv = la.inv(c)

After Change

        if self._k > 1:
            e = np.zeros(np.shape(x))
            for i in range(self._k):
                e[i] = sp.linalg.expm(x_inv_u[i])
        else:
            e = sp.linalg.expm(x_inv_u)
        return multiprod(x, e)
        // This alternative implementation is sometimes faster though less

Instances