int_power.c

Go to the documentation of this file.

 
#include "grib2_int.h"
 
double
int_power(double x, g2int y)
{
    double value;
 
    if (y < 0)
    {
        y = -y;
        x = 1.0 / x;
    }
    value = 1.0;
 
    while (y)
    {
        if (y & 1)
        {
            value *= x;
        }
        x = x * x;
        y >>= 1;
    }
    return value;
}