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)
  
          self.GDf = self.Gf * (self.Wtv * self.Df)[..., np.newaxis]
  
          # 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)
docs: Set dictionary array.