code: def Ge(self):
          
  
          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
docs: Result of US from the SVD decomposition G = USVᵀ.