VectorHelper.h

Go to the documentation of this file.

/****************************************************************************/
// A simple vector of SUMOReals
/****************************************************************************/
// 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.
//
/****************************************************************************/
#ifndef VectorHelper_h
#define VectorHelper_h


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

#include <vector>
#include <algorithm>
#include <iostream>


// ===========================================================================
// definitions
// ===========================================================================
typedef std::vector<bool> BoolVector;
typedef std::vector<SUMOReal> DoubleVector;
typedef std::vector<int> IntVector;


// ===========================================================================
// class definitions
// ===========================================================================
template<class T>
class VectorHelper {
public:
    static T sum(const std::vector<T> &v) {
        T sum = 0;
        for (typename std::vector<T>::const_iterator i=v.begin(); i!=v.end(); i++) {
            sum += *i;
        }
        return sum;
    }

    static void normalise(std::vector<T> &v, T msum=1.0) {
        if (msum==0) {
            // is an error; do nothing
            return;
        }
        T rsum = sum(v);
        if (rsum==0) {
            set(v, (T) 1.0/msum/(T) v.size());
            return;
        }
        div(v, rsum/msum);
    }

    static void div(std::vector<T> &v, T by) {
        for (typename std::vector<T>::iterator i=v.begin(); i!=v.end(); i++) {
            *i /= by;
        }
    }

    static void removeDouble(std::vector<T> &v) {
        typename std::vector<T>::iterator i=v.begin();
        while (i!=v.end()) {
            for (typename std::vector<T>::iterator j=i+1; j!=v.end();) {
                if (*i==*j) {
                    j = v.erase(j);
                } else {
                    j++;
                }
            }
            i++;
        }
    }


    static void set(std::vector<T> &v, T to) {
        for (typename std::vector<T>::iterator i=v.begin(); i!=v.end(); i++) {
            *i = to;
        }
    }

    static T maxValue(const std::vector<T> &v) {
        SUMOReal m = *(v.begin());
        for (typename std::vector<T>::const_iterator j=v.begin()+1; j!=v.end(); j++) {
            if ((*j)>m) {
                m = *j;
            }
        }
        return m;
    }

    static T minValue(const std::vector<T> &v) {
        SUMOReal m = *(v.begin());
        for (typename std::vector<T>::const_iterator j=v.begin()+1; j!=v.end(); j++) {
            if ((*j)<m) {
                m = *j;
            }
        }
        return m;
    }

    static void remove_smaller_than(std::vector<T> &v, T swell) {
        for (typename std::vector<T>::iterator j=v.begin(); j!=v.end();) {
            if ((*j)<swell) {
                j = v.erase(j);
            } else {
                j++;
            }
        }
    }

    static void remove_larger_than(std::vector<T> &v, T swell) {
        for (typename std::vector<T>::iterator j=v.begin(); j!=v.end();) {
            if ((*j)>swell) {
                j = v.erase(j);
            } else {
                j++;
            }
        }
    }

    static void add2All(std::vector<T> &v, T what) {
        for (typename std::vector<T>::iterator j=v.begin(); j!=v.end(); j++) {
            (*j) += what;
        }
    }

    static bool subSetExists(const std::vector<T> &v1, const std::vector<T> &v2) {
        for (typename std::vector<T>::const_iterator i=v1.begin(); i!=v1.end(); i++) {
            int val1 = (*i);
            if (find(v2.begin(), v2.end(), val1)!=v2.end()) {
                return true;
            }
        }
        return false;
    }



};

template<class T>
std::ostream &operator<<(std::ostream &os, const std::vector<T> &v) {
    for (typename std::vector<T>::const_iterator i=v.begin(); i!=v.end(); i++) {
        if (i!=v.begin()) {
            os << ", ";
        }
        os << (*i);
    }
    return os;
}



#endif

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