code: def setdict(self, D):
          
  
          self.D = np.asarray(D, dtype=self.dtype)
          self.DTS = self.D.T.dot(self.S)
          # Factorise dictionary for efficient solves
          self.lu, self.piv = sl.cho_factor(self.D, self.rho)
          self.lu = np.asarray(self.lu, dtype=self.dtype)
docs: Set dictionary array.