Упрощение дробей

fract fract::reduce(fract ab)
{
    int tmp = 0;
    if ((prime(ab.numerator) == false) && (prime(ab.denominator) == false)) {
        int nod = NOD(ab.numerator, ab.denominator);
        ab.numerator /= nod;
        ab.denominator /= nod;
        return ab;
    }
    else
    {       
        ab.unit += ab.numerator / ab.denominator;
        while (ab.numerator>ab.denominator) 
        {
            ab.numerator -= ab.denominator;
        }
        return ab;
    }
}
 
bool fract::prime(int a)
{
    for (int i = 2; i<sqrt(a)+1; i++) 
    {
        if (a % i == 0) 
        {
            return false;
        }
    }
    return true;
}

fract fract::reduce(fract ab)
{
    while (ab.numerator > ab.denominator)
    {
        ab.numerator -= ab.denominator;
        ab.unit++;
    }
    int nod = NOD(ab.numerator, ab.denominator);
    ab.numerator /= nod;
    ab.denominator /= nod;
    ab.unit += ab.numerator / ab.denominator;
    return ab;
}

fract fract::reduce()
{
    while (numerator > denominator)
    {
        numerator -= denominator;
        unit++;
    }
    int nod = NOD(numerator, denominator);
    numerator /= nod;
    denominator /= nod;
    unit += numerator / denominator;
    return *this;
}