def Ge(self):
        """
        Result of US from the SVD decomposition G = USVᵀ.
        """

        from scipy.linalg import svd
        from numpy_sugar.linalg import ddot

        U, S, _ = svd(self._G, full_matrices=False, check_finite=False)
        if U.shape[1] < self._G.shape[1]:
            return ddot(U, S)
        return self._G