Here is the comparison for Zeldovich test using different reconstruction methods after the shock is developed, all cases use HLLC and VL, the Dual Energy parameters are:

Here using PPMP:

PPMC Reconstruction:

From the CHOLLA paper:

Our implementation of this method closely follows that outlined in Stone et al. (2008).

But from Stone et al. (2008) equation 46:

The signs are different. Changing the signs in ppmc_cuda.cu to match Stone (2008), I get the next result:

Seems similar to the evolution using PPMP, looks like there is an inconsistency between reconstruction methods.