Function powInt

src/CPnumerics.cpp:40–69

    }
}
double powInt(double x, int y) {
    // Raise a double to an integer power
    // Overload not provided in math.h
    int i = 0;
    double product = 1.0;
    double x_in = NAN;
    int y_in = 0;

    if (y == 0) {
        return 1.0;
    }

    if (y < 0) {
        x_in = 1 / x;
        y_in = -y;
    } else {
        x_in = x;
        y_in = y;
    }

    if (y_in == 1) {
        return x_in;
    }

    product = x_in;
    for (i = 1; i < y_in; i++) {
        product = product * x_in;
    }
    return product;
}

void MatInv_2(double A[2][2], double B[2][2]) {
    double Det = NAN;