Played with the back transform from rho, j, phi to populations, but did not manage to fix the problem. The back conversion is definitely the culprit, in both the LB and the LB_GPU. Also removed some warnings and made some minor style changes.