consolidate coasting/pseudoFriction

[namdre]

sumo/src/microsim/MSVehicle.cpp

Lines 4138 to 4143 in 3937085

	MSVehicle::setBrakingSignals(double vNext) {
	// To avoid casual blinking brake lights at high speeds due to dawdling of the
	// leading vehicle, we don't show brake lights when the deceleration could be caused
	// by frictional forces and air resistance (i.e. proportional to v^2, coefficient could be adapted further)
	double pseudoFriction = (0.05 + 0.005 * getSpeed()) * getSpeed();
	bool brakelightsOn = vNext < getSpeed() - ACCEL2SPEED(pseudoFriction);

vs

sumo/src/utils/emissions/PollutantsInterface.cpp

Lines 196 to 210 in 3937085

	PollutantsInterface::Helper::getCoastingDecel(const SUMOEmissionClass c, const double v, const double a, const double slope, const EnergyParams* param) const {
	// the interpolation for small v is basically the same as in PHEMlightdllV5::CEP::GetDecelCoast
	if (v < PHEMlightdllV5::Constants::SPEED_DCEL_MIN) {
	return v / PHEMlightdllV5::Constants::SPEED_DCEL_MIN * getCoastingDecel(c, PHEMlightdllV5::Constants::SPEED_DCEL_MIN, a, slope, param);
	}
	if (param == nullptr) {
	param = EnergyParams::getDefault();
	}
	// the magic numbers below come from a linear interpolation with http://ts-sim-service-ba/svn/simo/trunk/projects/sumo/data/emissions/linear.py
	const double mass = param->getDouble(SUMO_ATTR_MASS);
	const double area = param->getDouble(SUMO_ATTR_WIDTH) * param->getDouble(SUMO_ATTR_HEIGHT) * M_PI / 4.;
	const double incl = area / mass * -9.05337017 + -0.00017774;
	const double grad = PHEMlightdllV5::Constants::GRAVITY_CONST * slope / 100.;
	return MIN2(0., incl * v + 0.00001066 * mass + -0.38347107 - 20.0 * incl - grad);
	}

[namdre]

also, there should be a way to have the cfModel override the coasting decelerating (trains compute this based on their resistance curve in MSCFModel_Rail::maxNextSpeed and MSCFModel_Rail::minNextSpeed)