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:
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.