LISASMDCalcMethod1.H

Go to the documentation of this file.

{
    // calculate the new diameter with a Rosin Rammler distribution
 
    scalar minValue = min(d, dD/10.0);
    scalar maxValue = dD;
 
    if (maxValue - minValue < small)
    {
        minValue = d/10.0;
    }
 
    scalar range = maxValue - minValue;
 
    scalar y = 0;
 
    bool success = false;
 
    while(!success)
    {
 
        x = minValue + range*rndGen.sample01<scalar>();
        y = rndGen.sample01<scalar>();
        scalar p = 0.0;
        scalar nExp = 1;
        scalar xx = pow(x/dD, nExp);
 
        p = xx*exp(-xx);
        if (y<p)
        {
            success = true;
        }
    }
}