Distribution_Points.cpp

Go to the documentation of this file.

/****************************************************************************/
// The description of a distribution by a curve
/****************************************************************************/
// SUMO, Simulation of Urban MObility; see http://sumo.sourceforge.net/
// Copyright 2001-2010 DLR (http://www.dlr.de/) and contributors
/****************************************************************************/
//
//   This program 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 2 of the License, or
//   (at your option) any later version.
//
/****************************************************************************/


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

#include <cassert>
#include "Distribution.h"
#include <utils/geom/Position2DVector.h>
#include "Distribution_Points.h"
#include <utils/common/StdDefs.h>

#ifdef CHECK_MEMORY_LEAKS
#include <foreign/nvwa/debug_new.h>
#endif // CHECK_MEMORY_LEAKS


// ===========================================================================
// method definitions
// ===========================================================================
Distribution_Points::Distribution_Points(const std::string &id,
        const Position2DVector &points,
        bool interpolating) throw()
        : Distribution(id), myPoints(points), myProbabilitiesAreComputed(false),
        myInterpolateDist(interpolating) {}


Distribution_Points::~Distribution_Points() throw() {}


SUMOReal
Distribution_Points::getMax() const {
    assert(myPoints.size()>0);
    const Position2D &p = myPoints[-1];
    return p.x();
}


size_t
Distribution_Points::getAreaNo() const {
    return myPoints.size()-1;
}


SUMOReal
Distribution_Points::getAreaBegin(size_t index) const {
    return myPoints[(int) index].x();
}


SUMOReal
Distribution_Points::getAreaEnd(size_t index) const {
    return myPoints[(int) index+1].x();
}


SUMOReal
Distribution_Points::getAreaPerc(size_t index) const {
    if (!myProbabilitiesAreComputed) {
        SUMOReal sum = 0;
        size_t i;
        if (myInterpolateDist) {
            for (i=0; i<myPoints.size()-1; i++) {
                SUMOReal width = getAreaEnd(i) - getAreaBegin(i);
                SUMOReal minval = MIN2(myPoints[(int) i].y(), myPoints[(int) i].y());
                SUMOReal maxval = MAX2(myPoints[(int) i].y(), myPoints[(int) i].y());
                SUMOReal amount = minval * width + (maxval-minval) * width / (SUMOReal) 2.;
                myProbabilities.push_back(amount);
                sum += amount;
            }
        } else {
            for (i=0; i<myPoints.size()-1; i++) {
                myProbabilities.push_back(myPoints[(int) i].y());
                sum += myPoints[(int) i].y();
            }
        }
        // normalize
        if (myInterpolateDist) {
            for (i=0; i<myPoints.size()-1; i++) {
                myProbabilities[i] = myProbabilities[i] / sum;
            }
        } else {
            for (i=0; i<myPoints.size()-1; i++) {
                myProbabilities[i] = myProbabilities[i] / sum;
            }
        }
        myProbabilitiesAreComputed = true;
    }
    return myProbabilities[index];
}



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