code: def relax_AX(self):
          
  
          # We need to keep the non-relaxed version of AX since it is
          # required for computation of primal residual r
          self.AXnr = self.cnst_A(self.X, self.Xf)
          if self.rlx == 1.0:
              # If RelaxParam option is 1.0 there is no relaxation
              self.AX = self.AXnr
          else:
              # Avoid calling cnst_c() more than once in case it is expensive
              # (e.g. due to allocation of a large block of memory)
              if not hasattr(self, '_cnst_c'):
                  self._cnst_c = self.cnst_c()
              # Compute relaxed version of AX
              alpha = self.rlx
              self.AX = alpha*self.AXnr - (1-alpha)*(self.cnst_B(self.Y) -
                                                     self._cnst_c)
docs: Implement relaxation if option ``RelaxParam`` != 1.0.