                         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;

                        Real3 posRadB = mR3(sortedPosRad[j]);
                        Real3 rij = Distance(posRadA, posRadB);
                        Real d = length(rij);
                        if (d > SuppRadii || sortedRhoPreMu[j].w <= -2)
                            continue;
                        Real h_j = sortedPosRad[j].w;
                        Real h_ij = 0.5 * (h_j + h_i);

                        Real d = length(dist3);
                        if (d > SuppRadii)
                            continue;
                        if (sortedRhoPreMu_old[j].w > -1.5 && sortedRhoPreMu_old[j].w < -0.5) {
                            Real h_j = sortedPosRad[j].w;
                            Real m_j = paramsD.markerMass;
                            Real W3 = W3h(d, 0.5 * (h_j + h_i));
