code: def relax_AX(self):
          
  
          self.AXnr = self.X
          if self.rlx == 1.0:
              self.AX = self.X
          else:
              alpha = self.rlx
              self.AX = alpha*self.X + (1 - alpha)*self.Y[..., np.newaxis]
docs: Implement relaxation if option ``RelaxParam`` != 1.0.