 
 //------------------------------------------------------------------------------
 __device__ __host__ inline void CalcElasticForcesKernel(
-		real3 * flex_FSI_NodesForcesD1,
-		real3 * flex_FSI_NodesForcesD2,
+		real3 * flex_NodesForcesD1,
+		real3 * flex_NodesForcesD2,
 		const real3 * ANCF_NodesD,
 		const real3 * ANCF_SlopesD,
 		const real3 * ANCF_NodesVelD,

	thrust::copy(ANCF_SlopesVelH2.begin(), ANCF_SlopesVelH2.end(), ANCF_SlopesVelD2.begin());
	//---------------------------------------------

	ANCF_NodesH2.clear();
	ANCF_SlopesH2.clear();
	ANCF_NodesVelH2.clear();

		const ANCF_Params & flexParams,
		real_ dT)
{
	thrust::device_vector<real3> flex_NodesForcesD1 = flex_FSI_NodesForcesD1;
	thrust::device_vector<real3> flex_NodesForcesD2 = flex_FSI_NodesForcesD2;

	thrust::host_vector<real3> ANCF_NodesH2 = ANCF_NodesD2;
	thrust::host_vector<real3> ANCF_SlopesH2 = ANCF_SlopesD2;
	thrust::host_vector<real3> ANCF_NodesVelH2 = ANCF_NodesVelD2;
