 }
 
 /**
- * @brief Wrapper function for collideD
+ * @brief Wrapper function for collide
  * @details
  * 		See SDKCollisionSystem.cuh for informaton on collide
  */

			if (paramsD.enableAggressiveTweak) {
				vel_XSPH = mR3(0);
			} else {
				printf("Error! particle vel_XSPH is NAN: thrown from SDKCollisionSystem.cu, UpdateFluidD !\n");
				return;
			}
		}

	// write new velocity back to original unsorted location
	// *** let's tweak a little bit :)
	if (!(isfinite(derivVelRho.x) && isfinite(derivVelRho.y)
			&& isfinite(derivVelRho.z) )) {
		printf("Error! particle derivVel is NAN: thrown from SDKCollisionSystem.cu, collideD !\n");
		*isErrorD = true;
	}
	if (!(isfinite(derivVelRho.w))) {
		printf("Error! particle derivRho is NAN: thrown from SDKCollisionSystem.cu, collideD !\n");
		*isErrorD = true;
	}
	sortedDerivVelRho_fsi_D[index] = derivVelRho;
}
