 __device__ inline real4 DifVelocityRho(
 		const real3 & dist3,
 		const real_ & d,
-		const real_ & rSPH,
 		const real4 & velMasA,
 		const real3 & vel_XSPH_A,
 		const real4 & velMasB,

				real3 posRadB = FETCH(sortedPosRad, j);
				real3 dist3Alpha = posRadA - posRadB;
				real3 dist3 = Distance(posRadA, posRadB);
				real_ rSPH = paramsD.HSML;
				real_ d = length(dist3);
				if (d > RESOLUTION_LENGTH_MULT * paramsD.HSML) continue;
				real4 velMasB = FETCH(sortedVelMas, j);
				real4 rhoPresMuB = FETCH(sortedRhoPreMu, j);


__device__ inline real3 DifVelocity_SSI_Lubrication(
				const real3 & dist3,
				const real_ & d,
				const real4 & velMasA,
				const real4 & velMasB) {
//printf("** DifVelocity_SSI_DEM\n");
