 		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;

		real3 * ANCF_NodesAccD,
		real3 * ANCF_SlopesAccD,

		const real3 * flex_FSI_NodesForcesD1,
		const real3 * flex_FSI_NodesForcesD2,

		const int2 * ANCF_ReferenceArrayNodesOnBeamsD,
		const real_ * ANCF_Beam_LengthD,

		real3 * ANCF_NodesAccD,
		real3 * ANCF_SlopesAccD,

		const real3 * flex_NodesForcesD1,
		const real3 * flex_NodesForcesD2,

		const int2 * ANCF_ReferenceArrayNodesOnBeamsD,
		const real_ * ANCF_Beam_LengthD,
