code: def setdict(self, D=None):
          
  
          if D is not None:
              self.D = np.asarray(D, dtype=self.dtype)
          self.Df = sl.rfftn(self.D, self.cri.Nv, self.cri.axisN)
          # Compute D^H S
          self.DSf = np.conj(self.Df) * self.Sf
          if self.cri.Cd > 1:
              self.DSf = np.sum(self.DSf, axis=self.cri.axisC, keepdims=True)
          if self.opt['HighMemSolve'] and self.cri.Cd == 1:
              self.c = sl.solvedbd_sm_c(
                  self.Df, np.conj(self.Df), self.mu*self.GHGf + self.rho,
                  self.cri.axisM)
          else:
              self.c = None
docs: Set dictionary array.