code: def ccmod_ystep():
      
  
      mAXU = np.mean(mp_D_X + mp_D_U, axis=0)
      mp_D_Y[:] = mp_dprox(mAXU)
docs: Do the Y step of the ccmod stage. There are no parameters
      or return values because all inputs and outputs are from and to
      global variables.