Line data    Source code

       1              : /****************************************************************************/
       2              : // Eclipse SUMO, Simulation of Urban MObility; see https://eclipse.dev/sumo
       3              : // Copyright (C) 2001-2024 German Aerospace Center (DLR) and others.
       4              : // This program and the accompanying materials are made available under the
       5              : // terms of the Eclipse Public License 2.0 which is available at
       6              : // https://www.eclipse.org/legal/epl-2.0/
       7              : // This Source Code may also be made available under the following Secondary
       8              : // Licenses when the conditions for such availability set forth in the Eclipse
       9              : // Public License 2.0 are satisfied: GNU General Public License, version 2
      10              : // or later which is available at
      11              : // https://www.gnu.org/licenses/old-licenses/gpl-2.0-standalone.html
      12              : // SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
      13              : /****************************************************************************/
      14              : /// @file    MSCFModel_KraussX.cpp
      15              : /// @author  Tobias Mayer
      16              : /// @author  Daniel Krajzewicz
      17              : /// @author  Jakob Erdmann
      18              : /// @author  Michael Behrisch
      19              : /// @author  Laura Bieker
      20              : /// @date    Mon, 04 Aug 2009
      21              : ///
      22              : // Krauss car-following model, changing accel and speed by slope
      23              : /****************************************************************************/
      24              : #include <config.h>
      25              : 
      26              : #include <microsim/lcmodels/MSAbstractLaneChangeModel.h>
      27              : #include <microsim/MSVehicle.h>
      28              : #include <microsim/MSNet.h>
      29              : #include "MSCFModel_KraussX.h"
      30              : 
      31              : 
      32              : #define OVERBRAKING_THRESHOLD -3
      33              : 
      34              : // ===========================================================================
      35              : // method definitions
      36              : // ===========================================================================
      37           12 : MSCFModel_KraussX::MSCFModel_KraussX(const MSVehicleType* vtype):
      38              :     MSCFModel_Krauss(vtype),
      39           12 :     myTmp1(vtype->getParameter().getCFParam(SUMO_ATTR_TMP1, 0.0)),
      40           12 :     myTmp2(vtype->getParameter().getCFParam(SUMO_ATTR_TMP2, 0.0)) {
      41           12 : }
      42              : 
      43              : 
      44           24 : MSCFModel_KraussX::~MSCFModel_KraussX() {}
      45              : 
      46              : 
      47              : MSCFModel*
      48            0 : MSCFModel_KraussX::duplicate(const MSVehicleType* vtype) const {
      49            0 :     return new MSCFModel_KraussX(vtype);
      50              : }
      51              : 
      52              : 
      53              : double
      54      3708488 : MSCFModel_KraussX::patchSpeedBeforeLC(const MSVehicle* veh, double vMin, double vMax) const {
      55      3708488 :     return dawdleX(veh->getSpeed(), vMin, vMax, veh->getRNG());
      56              : }
      57              : 
      58              : 
      59              : double
      60      3708488 : MSCFModel_KraussX::dawdleX(double vOld, double vMin, double vMax, SumoRNG* rng) const {
      61              :     double speed = vMax;
      62      3708488 :     if (!MSGlobals::gSemiImplicitEulerUpdate) {
      63              :         // in case of the ballistic update, negative speeds indicate
      64              :         // a desired stop before the completion of the next timestep.
      65              :         // We do not allow dawdling to overwrite this indication
      66            0 :         if (speed < 0) {
      67              :             return speed;
      68              :         }
      69              :     }
      70              :     // extra slow to start
      71      3708488 :     if (vOld < myAccel) {
      72       162560 :         speed -= ACCEL2SPEED(myTmp1 * myAccel);
      73              :     }
      74      3708488 :     const double random = RandHelper::rand(rng);
      75      3708488 :     speed -= ACCEL2SPEED(myDawdle * myAccel * random);
      76              :     // overbraking
      77      3708488 :     if (vOld > vMax) {
      78       296528 :         speed -= ACCEL2SPEED(myTmp2 * myAccel * random);
      79              :         //std::cout << " vMin=" << vMin << " vMax=" << vMax << "speed=" << speed << " d1=" << ACCEL2SPEED(myDawdle * myAccel * random) << " d2=" << ACCEL2SPEED(myTmp2 * myAccel * random) << " unexpectedDecel=" << (speed < vMin) << "\n";
      80       296528 :         if (MSGlobals::gSemiImplicitEulerUpdate) {
      81              :             speed = MAX2(0.0, speed);
      82              :         }
      83              :     }
      84              :     speed = MAX2(vMin, speed);
      85              :     return speed;
      86              : }
      87              : 
      88              : 
      89              : /****************************************************************************/