Before revision:  For the arbitrary high order DG method introduced in H. Liu and Z. Wang, J. Comput. Phys., 328, 2017 , solution positivity is only numerically enforced.
After revision:  For the arbitrary high order DG method introduced in  , solution positivity is only numerically enforced.