code: def rhochange(self):
          
  
          if self.opt['HighMemSolve']:
              self.c = sl.solvedbi_sm_c(self.gDf, np.conj(self.gDf), self.rho,
                                        self.cri.axisM)
docs: Updated cached c array when rho changes.