 		Zero12(e_ee);
 		// Elastic Force, 1/2: tension force, GQ 5th order. Maybe 4th order is enough as well.
 		for (int k = 0; k < 5; k ++) {
-			real_ gqPoint = (lE - 0) / 2 * GQD.GQ5_p[k] + (lE + 0) / 2;
+			real_ gqPoint = (lE - 0) / 2 * GQ.GQ5_p[k] + (lE + 0) / 2;
 			eps_eps_e(e_ee, gqPoint, lE, e);
-			SumElem12(f_e, e_ee, flexParamsD.E * flexParamsD.A * (lE - 0) / 2 * GQD.GQ5_w[k]);
+			SumElem12(f_e, e_ee, flexParams.E * flexParams.A * (lE - 0) / 2 * GQ.GQ5_w[k]);
 		}
 		// Elastic Force, 2/2: bending force, GQ 3rd order.
 		Elem12 k_ke;
 		Zero12(k_ke);
 		for (int k = 0; k < 3; k ++) {
-			real_ gqPoint = (lE - 0) / 2 * GQD.GQ3_p[k] + (lE + 0) / 2;
+			real_ gqPoint = (lE - 0) / 2 * GQ.GQ3_p[k] + (lE + 0) / 2;
 			kappa_kappa_e(k_ke, gqPoint, lE, e);
-			SumElem12(f_e, k_ke, flexParamsD.E * flexParamsD.I  * (lE - 0) / 2 * GQD.GQ3_w[k]);
+			SumElem12(f_e, k_ke, flexParams.E * flexParams.I  * (lE - 0) / 2 * GQ.GQ3_w[k]);
 		}
 		// Gravitational Foce
 		Elem12 f_g;
 		Zero12(f_g);
-		gravitational_force(f_g, lE, flexParamsD.rho, flexParamsD.A, flexParamsD.gravity);
+		gravitational_force(f_g, lE, flexParams.rho, flexParams.A, flexParams.gravity);
 		SumElem12(f_e, f_g, -1);
 		// Add element forces to associated nodes
 		Add_f_ToForces(flex_FSI_NodesForcesD1, flex_FSI_NodesForcesD2, -1, f_e, nodeIdx);
 	}
 }
 //------------------------------------------------------------------------------
-__global__ void CalcElasticForcesD(
+__GLOBAL__ void CalcElasticForcesD(
 		real3 * flex_FSI_NodesForcesD1,
 		real3 * flex_FSI_NodesForcesD2,
 		const real3 * ANCF_NodesD,

		const int2 * ANCF_ReferenceArrayNodesOnBeamsD,
		const real_ * ANCF_Beam_LengthD,
		const bool * ANCF_IsCantileverD,
		uint i
	)
{

		CalcElasticForcesKernel(
				flex_FSI_NodesForcesD1, flex_FSI_NodesForcesD2,
				ANCF_NodesD, ANCF_SlopesD, ANCF_NodesVelD, ANCF_SlopesVelD,
				ANCF_ReferenceArrayNodesOnBeamsD, ANCF_Beam_LengthD, flexParamsD, GQD, i);
#if !flexGPU
	}
#endif
