 	rhoPresMu.y = Eos(rho2, rhoPresMu.w);
 	rhoPresMu.x = rho2;
 	if (!(isfinite(rhoPresMu.x) && isfinite(rhoPresMu.y) && isfinite(rhoPresMu.z) && isfinite(rhoPresMu.w))) {
-		printf("Error! particle rho pressure is NAN: thrown from SDKCollisionSystem.cu, UpdateFluidD !\n");
+		printf("Error! particle rho pressure is NAN: thrown from SDKCollisionSystem.cu, UpdateFluidDKernel !\n");
+		*isErrorD = true;
 		return;
 	}
 	rhoPresMuD[index] = rhoPresMu;  // rhoPresMuD updated

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

}

/**
 * @brief Wrapper function for collide
 * @details
 * 		See SDKCollisionSystem.cuh for informaton on collide
 */
