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