 							ANCF_IsCantilever,
 							referenceArray,
 
-							paramsH,
 							flexParams,
 							numObjects,
 							float(tStep)/stepEnd,


					//			/* post_process for Segre-Silberberg */ goes before ApplyBoundary
					//			if(tStep >= 0) {
					//				real2 channelCenter = .5 * R2(paramsH.cMax.y + paramsH.cMin.y, paramsH.cMax.z + paramsH.cMin.z);
					//				FindPassesFromTheEnd(posRigidD, distributionD, numRigidBodies, channelCenter, channelRadius, numberOfSections);
					//			}


	paramsH_B.bodyForce4 = R4(0);
	paramsH_B.gravity = R3(0);

	SimParams currentParamsH = paramsH;

	//for (int tStep = 0; tStep < 0; tStep ++) {
	for (int tStep = 0; tStep < stepEnd + 1; tStep++) {
		//edit  since yu deleted cyliderRotOmegaJD
