                         if (d > SuppRadii)
                             continue;
                         Real4 rhoPresMuB = sortedRhoPreMu[j];
-                        if (rhoPresMuA.w > -.1 && rhoPresMuB.w > -.1) {  // no rigid-rigid force
+                        if (rhoPresMuA.w > -0.5 && rhoPresMuB.w > -0.5) {  // no rigid-rigid force
                             continue;
                         }
                         Real invd = 1.0 / d;

                                       Real4 rhoPresMuA,
                                       Real4 rhoPresMuB,
                                       Real multViscosity) {
    if (rhoPresMuA.w > -1 && rhoPresMuB.w > -1)
        return mR4(0.0);

    Real3 gradW = GradWh(dist3, (posRadA.w + posRadB.w) * 0.5);

        if (d > SuppRadii)
            continue;
        Real4 rhoPresMuB = sortedRhoPreMu[j];
        if (rhoPresMuA.w > -0.5 && rhoPresMuB.w > -0.5) {  
            continue; // No BCE-BCE interaction
        }
        Real invd = 1.0 / d;
