 
 
 	//------------------------ delete stuff
-//	flex_FSI_NodesForces1.clear();
-//	flex_FSI_NodesForces2.clear();
+	flex_FSI_NodesForces1.clear();
+	flex_FSI_NodesForces2.clear();
 }
 //##############################################################################################################################################
 // the main function, which updates the particles and implements BC

	thrust::fill(flex_FSI_NodesForces1.begin(), flex_FSI_NodesForces1.end(),R3(0));
	thrust::fill(flex_FSI_NodesForces2.begin(), flex_FSI_NodesForces2.end(),R3(0));

//	thrust::device_vector<real3> flexNodesForcesAllMarkers1(totalNumberOfFlexMultNodes);
//	thrust::device_vector<real3> flexNodesForcesAllMarkers2(totalNumberOfFlexMultNodes);
//
//	uint nBlocks_numFlex_SphMarkers;
//	uint nThreads_SphMarkers;
//	computeGridSize(numFlex_SphMarkers, 256, nBlocks_numFlex_SphMarkers, nThreads_SphMarkers);


						//ff1 for Radu
					//	for (int i = 0; i < ANCF_Beam_LengthD.size(); i++) {

//	//

	//################################################### update rigid body things
	computeGridSize(numFlex_SphMarkers, 256, nBlocks_numFlex_SphMarkers, nThreads_SphMarkers);
	UpdateFlexMarkersPosition<<<nBlocks_numFlex_SphMarkers, nThreads_SphMarkers>>>(
			R3CAST(posRadD), R4CAST(velMasD),
			I1CAST(flexIdentifierD),
			R3CAST(flexSPH_MeshPos_LRF_D),
			R3CAST(flexSPH_MeshSlope_Initial_D),
			R1CAST(flexParametricDistD),
			R1CAST(ANCF_Beam_LengthD),
			I2CAST(ANCF_ReferenceArrayNodesOnBeamsD),
			R3CAST(ANCF_NodesD2),
			R3CAST(ANCF_SlopesD2),
			R3CAST(ANCF_NodesVelD2),
			R3CAST(ANCF_SlopesVelD2)
			);

	cudaThreadSynchronize();

