code: def relax_AX(self):
          
  
          self.AXnr = self.cnst_A(self.X, self.Xf)
          if self.rlx == 1.0:
              self.AX = self.AXnr
          else:
              if not hasattr(self, 'c0'):
                  self.c0 = self.cnst_c0()
              if not hasattr(self, 'c1'):
                  self.c1 = self.cnst_c1()
              alpha = self.rlx
              self.AX = alpha*self.AXnr + (1-alpha)*self.block_cat(
                  self.var_y0() + self.c0, self.var_y1() + self.c1)
docs: Implement relaxation if option ``RelaxParam`` != 1.0.