code: def U(self):
          
          if self._U is None:
              sinv = N.diag(1/self.singular_values)
              self._U = dot(self.arr,self.V.T,sinv)
          return self._U
docs: Property to support lazy evaluation of residuals