06f40aff5fafd3db599ae361d8d7aa33f3142038,pymanopt/manifolds/grassmann.py,Grassmann,retr,#Grassmann#,90

Before Change

            // Calculate "thin" qr decomposition of X + G
            q, r = np.linalg.qr(X + G)
            // Unflip any flipped signs
            XNew = np.dot(q, np.diag(np.sign(np.sign(np.diag(r))+.5)))
            return XNew
        else:
            XNew = X + G
            for i in xrange(self._k):

After Change

            // columns. Compare this with the Stiefel manifold.

            // Compute the polar factorization of Y = X+G
            u, s, v = svd(X + G, full_matrices=False)
            XNew = u.dot(v.T);
        else:
            XNew = np.zeros(np.shape(X))

Instances