I find a small bug in momentTemplates3D.h

At the line 488 \palabos-master\src\latticeBoltzmann\momentTemplates3D.h, I think the line compute_PiNeq(f, rhoBar, j, PiNeq, rhoBar) should be compute_PiNeq(f, rhoBar, j, PiNeq, invRho).
what about you?

1 Like

In this file, there are three same functions “compute_rhoBar_thetaBar_j_PiNeq”. These functions call compute_PiNeq in a wrong way: the last arg “incompr” is bool, which should be T.

Oh. that may be the case indeed. Could you please open an issue on the gitlab repo please?

ok, i can. At first ,I register an account on gitlab.