 		thrust::device_vector<real3> vel_XSPH_D(numAllMarkers);
 				//rigids
 		thrust::device_vector<real3> posRigidD2 = posRigidD;
-		thrust::device_vector<real3> posRadRigidCumulativeD2 = posRigidCumulativeD;
+		thrust::device_vector<real3> posRigidCumulativeD2 = posRigidCumulativeD;
 		thrust::device_vector<real4> velMassRigidD2 = velMassRigidD;
 		thrust::device_vector<real3> omegaLRF_D2 = omegaLRF_D;
 				thrust::device_vector<real3> AD1_2 = AD1;

	}
}
//--------------------------------------------------------------------------------------------------------------------------------
__global__ void SumSurfaceInteractionForces(real3 * totalForcesRigid3, real4 * totalSurfaceInteractionRigid4, real4 * velMassRigidD) {
	uint rigidSphereA = blockIdx.x * blockDim.x + threadIdx.x;
	if (rigidSphereA >= numRigidBodiesD) {
		return;

//	//TODO: add paramsH.gravity to Flex objects
//	thrust::device_vector<real3> gravityForces3(numRigidBodies);
//	thrust::fill(gravityForces3.begin(), gravityForces3.end(), paramsH.gravity);
//	thrust::transform(totalAccRigid3.begin(), totalAccRigid3.end(), gravityForces3.begin(), totalAccRigid3.begin(), thrust::plus<real3>());
//	gravityForces3.clear();
//	//

