 			Real Wd = W3(d);
 			Real3 velMasB = FETCH(sortedVelMas, j);
 			//			deltaVDenom += mR4(
-			//					velMasB.w / rhoPresMuB.x * mR3(velMasB) * Wd,
-			//					velMasB.w / rhoPresMuB.x * Wd);
+			//					paramsD.markerMass / rhoPresMuB.x * velMasB * Wd,
+			//					paramsD.markerMass / rhoPresMuB.x * Wd);
 			//			deltaVDenom += mR4(
-			//					mR3(velMasB) * Wd,
+			//					velMasB * Wd,
 			//					Wd);
 
 			if (rhoPresMuB.w < -.1) { // only fluid pressure is used to update BCE pressure see Eq 27 of Adami, 2012 paper

				// colagrossi (2003), the other phase (i.e. rigid) should not be considered)
				Real multRho = 2.0f / (rhoPresMuA.x + rhoPresMuB.x);
				Real3 velMasB = FETCH(sortedVelMas, j);
				deltaV += paramsD.markerMass * velMasB - velMasA * W3(d) * multRho;
			}
		}
	}

	// write new velocity back to original unsorted location
	//  uint originalIndex = gridMarkerIndex[index];

	// Real newDensity = densityShare + paramsD.markerMass * W3(0); //?$ include the particle in its summation as well
	// if (rhoPreMuA.w < -.1) { rhoPreMuA.x = newDensity; }
	// rhoPreMuA.y = Eos(rhoPreMuA.x, rhoPreMuA.w);
	//   oldRhoPreMu[originalIndex] = rhoPreMuA;
