def U(self): ''' The the upper triangular matrix. ''' if self.order=='LU': if self.is_scalar: dl=ones(self.n) else: dl=[identity(self.p)]*self.n elif self.order=='UL' or self.order=='UDL' or self.order=='LDU': dl=self.hl U=build_sr(dl=dl,ul=self.ul) return U
def U(self): ''' The the upper triangular matrix. ''' if self.order == 'LU': if self.is_scalar: dl = ones(self.n) else: dl = [identity(self.p)] * self.n elif self.order == 'UL' or self.order == 'UDL' or self.order == 'LDU': dl = self.hl U = build_sr(dl=dl, ul=self.ul) return U
def L(self): ''' The the lower triangular matrix. ''' if self.order=='UL': if self.is_scalar: dl=ones(self.n) else: dl=[identity(self.p)]*n elif self.order=='LU' or self.order=='UDL' or self.order=='LDU': dl=self.hl L=build_sr(dl=dl,ll=self.ll) return L
def L(self): ''' The the lower triangular matrix. ''' if self.order == 'UL': if self.is_scalar: dl = ones(self.n) else: dl = [identity(self.p)] * n elif self.order == 'LU' or self.order == 'UDL' or self.order == 'LDU': dl = self.hl L = build_sr(dl=dl, ll=self.ll) return L