bubbleRadius.h

Go to the documentation of this file.

bool isNaN(double var);
double bubbleRadius (const double m0, const double m1);
double nodeRadius(const double &v);

bool isNaN(double var)
{
    volatile double d = var;
    return d != d;
}
double bubbleRadius (const double m0, const double m1)
{
    double R;
    R   = pow((3.0*m1/(4.0*M_PI*m0)), 1.0/3.0);
    if (isNaN(R)) {
        R=init_size;
    }
    return R;
}

double nodeRadius(const double &v)
{
    return (pow((6.0*v/M_PI), 1/3.0));
}