double Simulation::getTotalRate() 
{

    double buffer = 0;
    for(Event e : _events)
    {
        buffer += e.getRate() * e.getNumberOfTransitions();
    }
    return buffer;
}

#include <memory>
#include <iostream>
#include <vector>

using namespace std;

double fRand(double fMin, double fMax)
{
    double f = (double)rand() / RAND_MAX;
    return fMin + f * (fMax - fMin);
}

class test
{
public:
    test(double rate, int items) {
        _rate = rate;
        _n_items = items;
    };
    double getTotalRate() {
        return _rate * _n_items;
    }

private:
    int _n_items;
    double _rate;
};

double calcNumbers(vector<test> a)
{
    double buffer = 0;
    for(auto t : a)
        buffer += t.getTotalRate();
    return buffer;
}

int main()
{
    vector<test> vec;
    double running_val = 0;

    for(int i = 0; i < 1000; ++i)
    {
        double r = fRand(0, 100);
        int    n = rand() % 100; 

        test t(r, n);

        vec.push_back(t);

        for(int j = 0; j < n; ++j)
            running_val += r;
    }

    cout << running_val - calcNumbers(vec) << endl;
}

vec3 getGerstnerPos(float wL, float amp, vec2 dir)
{
   vec3 gPos = vec3(.1,.1,.1);
   float w = 2.*3.1415926/wL;
   float scalProd = position.x * dir.x + position.y * dir.y;
   float q = .9;
   gPos.x = position.x + q*amp*dir.x*cos(w*scalProd+time*animSpeed*w);
   gPos.y = position.y + q*amp*dir.y*cos(w*scalProd+time*animSpeed*w);
   gPos.z = amp*sin(scalProd*w + time * animSpeed * w);
   return gPos;
}