             Real3 rij = Distance(posRadA, posRadB);
             // Real d = length(rij);
             Real dd = rij.x*rij.x + rij.y*rij.y + rij.z*rij.z;
-            if (dd > SqRadii || sortedRhoPreMu[j].w <= -2)
+            if (dd > SqRadii || sortedRhoPreMu[j].w < -1.5)
                 continue;
             Real3 grad_i_wij = GradWh(rij, hA);
             Real3 grw_vj = grad_i_wij * paramsD.volume0;

                        Real3 posRadB = mR3(sortedPosRad[j]);
                        Real3 rij = Distance(posRadA, posRadB);
                        Real d = length(rij);
                        if (d > SuppRadii || sortedRhoPreMu[j].w <= -2)
                            continue;
                        Real3 eij = rij / d;


                            continue;
                        Real4 rhoPresMuB = sortedRhoPreMu[j];

                        if (rhoPresMuB.w > -0.5 || rhoPresMuB.w < -1.5)
                            continue;
                        Real3 velMasB = sortedVelMas[j];
                        Real rho_bar = 0.5 * (rhoPreMuA.x + rhoPresMuB.x);
