code: def ystep(self):
          
  
          super(ConvCnstrMODMaskDcpl_Consensus, self).ystep()
          AXU1 = self.AX1 + self.U1
          self.Y1 = self.rho*(AXU1 - self.S) / (self.W**2 + self.rho)
docs: The parent class ystep method is overridden to allow also
          performing the ystep for the additional variables introduced in
          the modification to the baseline algorithm.