ModelChannelPowerlaw.cpp

#include "Headers.h"
#include "ModelChannelPowerLaw.h"
#include "Rand.h"

CModelChannelPowerLaw::CModelChannelPowerLaw(
    unsigned int    numChannels,
    double          durationMin,
    double          durationMax,
    double          alpha
    ) : CModelChannel(numChannels)
{
    assert(durationMin < durationMax);

    this->durationMin = durationMin;
    this->durationMax = durationMax;
    this->alpha = alpha;
}

void CModelChannelPowerLaw::GenerateUptime(
    unsigned int&   channel,
    __time&         duration
    )
{
    channel = CRand::Generate(this->numChannels);
    assert(channel < this->numChannels);
    double rnd = CRand::Generate();
    
    duration = pow(rnd * (pow(this->durationMax,1-this->alpha) - pow(this->durationMin,1-this->alpha))
            + pow(this->durationMin, 1-this->alpha), 1/(1-this->alpha));
}

void CModelChannelPowerLaw::GenerateUptime(
    unsigned int    oldChannel,
    __time          oldDuration,
    unsigned int&   newChannel,
    __time&         newDuration
    )
{
    newChannel = CRand::Generate(this->numChannels-1);
    if(newChannel >= oldChannel) newChannel++;
    assert(newChannel < this->numChannels);
    double rnd = CRand::Generate();

    newDuration = pow(rnd * (pow(this->durationMax,1-this->alpha) - pow(this->durationMin,1-this->alpha))
            + pow(this->durationMin, 1-this->alpha), 1/(1-this->alpha));
}