HelpersHBEFA.h

Go to the documentation of this file.

/****************************************************************************/
// Helper methods for HBEFA-based emission computation
/****************************************************************************/
// SUMO, Simulation of Urban MObility; see http://sumo.sourceforge.net/
// Copyright (C) 2001-2012 DLR (http://www.dlr.de/) and contributors
/****************************************************************************/
//
//   This file is part of SUMO.
//   SUMO is free software: you can redistribute it and/or modify
//   it under the terms of the GNU General Public License as published by
//   the Free Software Foundation, either version 3 of the License, or
//   (at your option) any later version.
//
/****************************************************************************/
#ifndef HelpersHBEFA_h
#define HelpersHBEFA_h


// ===========================================================================
// included modules
// ===========================================================================
#ifdef _MSC_VER
#include <windows_config.h>
#else
#include <config.h>
#endif

#include <vector>
#include <cassert>
#include "StdDefs.h"
#include "SUMOVehicleClass.h"
#include <limits>
#include <cmath>


// ===========================================================================
// definitions
// ===========================================================================
#ifndef PI
#define PI 3.1415926535897932384626433832795
#endif


// ===========================================================================
// class definitions
// ===========================================================================
class HelpersHBEFA {
public:
    static SUMOReal computeCO(SUMOEmissionClass c, double v, double a) ;


    static SUMOReal computeCO2(SUMOEmissionClass c, double v, double a) ;


    static SUMOReal computeHC(SUMOEmissionClass c, double v, double a) ;


    static SUMOReal computeNOx(SUMOEmissionClass c, double v, double a) ;


    static SUMOReal computePMx(SUMOEmissionClass c, double v, double a) ;


    static SUMOReal computeFuel(SUMOEmissionClass c, double v, double a) ;


private:
    static inline SUMOReal compute(SUMOEmissionClass c, const int offset, double v, const double a) {
        switch (c) {
            case SVE_ZERO_EMISSIONS:
                return 0.;
            case SVE_UNKNOWN:
                c = SVE_P_LDV_7_7;
                break;
            default:
                break;
        }
        v *= 3.6;
        if (c > 42) {
            const double* f = myFunctionParameter[c - 42] + offset;
            return (SUMOReal) MAX2((f[0] + f[3] * v + f[4] * v * v + f[5] * v * v * v) / 3.6, 0.);
        }
        if (a < 0.) {
            return 0.;
        }
        const double* f = myFunctionParameter[c] + offset;
        const double alpha = asin(a / 9.81) * 180. / PI;
        return (SUMOReal) MAX2((f[0] + f[1] * alpha * v + f[2] * alpha * alpha * v + f[3] * v + f[4] * v * v + f[5] * v * v * v) / 3.6, 0.);
    }


private:
    static double myFunctionParameter[42][36];

};


#endif

/****************************************************************************/