Adapt calculation of macroscopic quantities

We discussed the calculation of macroscopic quantities in our VirtualFluids Stammtisch 25.10.23.

As a result, we decided to use the conditioned version implemented in the "old" K17Collision Kernel:

real drho = ((((f_PPP + f_MMM) + (f_MPM + f_PMP)) + ((f_MPP + f_PMM) + (f_MMP + f_PPM))) + (((f_0MP + f_0PM) + (f_0MM + f_0PP)) + ((f_M0P + f_P0M) + (f_M0M + f_P0P)) + ((f_MP0 + f_PM0) + (f_MM0 + f_PP0))) + ((f_M00 + f_P00) + (f_0M0 + f_0P0) + (f_00M + f_00P))) + f_000;

We also discussed if we can exchange the last two summations which results in this:

real drho = ((((f_PPP + f_MMM) + (f_MPM + f_PMP)) + ((f_MPP + f_PMM) + (f_MMP + f_PPM))) + (((f_0MP + f_0PM) + (f_0MM + f_0PP)) + ((f_M0P + f_P0M) + (f_M0M + f_P0P)) + ((f_MP0 + f_PM0) + (f_MM0 + f_PP0))) + f_000) + ((f_M00 + f_P00) + (f_0M0 + f_0P0) + (f_00M + f_00P));

here is the resulted deviation: conditioning_rho.txt

(DrivenCavity_new dataset contains the second version above)