 	)
 {
 	real_ l = ANCF_Beam_LengthD[i];
-	printf("l, %f\n", l);
 	int2 nodesPortion = ANCF_ReferenceArrayNodesOnBeamsD[i];
 	int numNodes = nodesPortion.y - nodesPortion.x;
 	int numElements = numNodes - 1;

	cudaThreadSynchronize();
	CUT_CHECK_ERROR("Kernel execution failed: CalcElasticForcesD");
#else
	CalcElasticForcesD __KERNEL__(nBlock_FlexBodies, nThreads_FlexBodies)(
			R3CAST(flex_FSI_NodesForcesH1), R3CAST(flex_FSI_NodesForcesH2),
			R3CAST(ANCF_NodesH), R3CAST(ANCF_SlopesH), R3CAST(ANCF_NodesVelH), R3CAST(ANCF_SlopesVelH),
			I2CAST(ANCF_ReferenceArrayNodesOnBeamsH), R1CAST(ANCF_Beam_LengthH));
#endif

	int totalNumberOfFlexNodes = flex_FSI_NodesForcesD1.size();
	thrust::device_vector<real3> ANCF_NodesAccD(totalNumberOfFlexNodes);
	thrust::device_vector<real3> ANCF_SlopesAccD(totalNumberOfFlexNodes);

	cudaThreadSynchronize();
	CUT_CHECK_ERROR("Kernel execution failed: SolveAndIntegrateInTime");
#else
	SolveForAccD(
			R3CAST(ANCF_NodesAccH), R3CAST(ANCF_SlopesAccH),
			R3CAST(flex_FSI_NodesForcesH1), R3CAST(flex_FSI_NodesForcesH2),
			I2CAST(ANCF_ReferenceArrayNodesOnBeamsH), R1CAST(ANCF_Beam_LengthH), BCAST(ANCF_IsCantileverH));
#endif

	printf("a2\n");


#if flexGPU
	IntegrateInTimeD __KERNEL__(nBlock_FlexBodies, nThreads_FlexBodies)(
