 
         CopySortedToOriginal_NonInvasive_R3(fsiGeneralData->vel_XSPH_D, vel_XSPH_Sorted_D, 
             markersProximityD->gridMarkerIndexD);
-        // CopySortedToOriginal_NonInvasive_R4(sphMarkersD->posRadD, sortedSphMarkersD->posRadD, 
-        //     markersProximityD->gridMarkerIndexD);
     }
 
     if (density_initialization % paramsH->densityReinit == 0)

                                       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);

                        Real3 posRadB = mR3(sortedPosRad[j]);
                        Real3 rij = Distance(posRadA, posRadB);
                        Real d = length(rij);
                        if (d > SuppRadii || sortedRhoPreMu[j].w < -1.5)
                            continue;
                        Real3 grad_i_wij = GradWh(rij, h_i);
                        Real3 grw_vj = grad_i_wij * paramsD.volume0;
