PinkGenerator.cpp

Go to the documentation of this file.

#include "PinkGenerator.h"


double PinkGenerator::av[3][5] = {{4.3606e-3, 5.9961e-3, 8.3586e-3, 0.0, 0.0},
                                  {0.0100, 0.00763, 0.00832, 0.01036, 0.0},
                                  {0.00224, 0.00821, 0.00798, 0.00694, 0.00992}};

double PinkGenerator::pv[3][5] = {{3.1878e-1, 7.7686e-1, 9.7785e-1, 0.0, 0.0},
                                  {0.3030, 0.74170, 0.91680, 0.97820, 0.0},
                                  {0.15571, 0.30194, 0.74115, 0.93003, 0.98035}};

//ThreadSpecificRandomDistribution< UniformDistribution > PinkGenerator::white_noise;

UniformDistribution PinkGenerator::white_noise;
MersenneTwister19937 PinkGenerator::mt;

PinkGenerator::PinkGenerator(int numGen)
{
    setNumGen(numGen);
}


void PinkGenerator::setNumGen(int numGen)
{
    if (numGen>=3 && numGen <= 5)
        nGen=numGen;
    else
        nGen=3;

    reset();
}

PinkGenerator::~PinkGenerator()
{
}


int PinkGenerator::adjust(double dt)
{
//    white_noise.set(UniformDistribution(0.0, 1.0));
    return 0;
}


int PinkGenerator::reset()
{
    for (int i=0; i<nGen; i++)
//        randreg[i]=av[nGen-3][i]*(2*white_noise()-0.5);               
        randreg[i]=av[nGen-3][i]*(2*white_noise(mt)-0.5);               

    return 0;
}


void PinkGenerator::advance()
{
//    double rv=2*white_noise();
    double rv=2*white_noise(mt);
    value=0.0;

    for (int i=0; i<nGen; i++)
    {
        if (rv > pv[nGen-3][i])
//            randreg[i] = av[nGen-3][i]*(2*white_noise()-0.5);
            randreg[i] = av[nGen-3][i]*(2*white_noise(mt)-0.5);

        value += randreg[i];
    }
}

---