4346ff930104b0d9b6735090a1ce4c485d40b981,pymanopt/manifolds/grassmann.py,Grassmann,retr,#Grassmann#,105

Before Change

            u, s, vt = svd(X + G, full_matrices=False)
            XNew = u.dot(vt)
        else:
            XNew = np.zeros(np.shape(X))
            for i in xrange(self._k):
                // q, r = np.linalg.qr(Y[i])
                // XNew[i] = np.dot(q, np.diag(np.sign(np.sign(np.diag(r))+.5)))
                u, s, vt = svd(X[i] + G[i], full_matrices=False)
                XNew[i] = u.dot(vt)
        return XNew

    def norm(self, X, G):
        // Norm on the tangent space is simply the Euclidean norm.

After Change

        // Compute the polar factorization of Y = X+G
        u, s, vt = svd(X + G, full_matrices=False)
        return multiprod(u, vt)

    def norm(self, X, G):
        // Norm on the tangent space is simply the Euclidean norm.
        return np.linalg.norm(G)

Instances