#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Henrywood and Agarwal, Equation (3)";

double f_if(float c0, float A, float V, float l) {
        float r9303 = c0;
        float r9304 = A;
        float r9305 = V;
        float r9306 = l;
        float r9307 = r9305 * r9306;
        float r9308 = r9304 / r9307;
        float r9309 = sqrt(r9308);
        float r9310 = r9303 * r9309;
        return r9310;
}

double f_id(double c0, double A, double V, double l) {
        double r9311 = c0;
        double r9312 = A;
        double r9313 = V;
        double r9314 = l;
        double r9315 = r9313 * r9314;
        double r9316 = r9312 / r9315;
        double r9317 = sqrt(r9316);
        double r9318 = r9311 * r9317;
        return r9318;
}


double f_of(float c0, float A, float V, float l) {
        float r9319 = A;
        float r9320 = V;
        float r9321 = r9319 / r9320;
        float r9322 = -2.004776637139463e+296f;
        bool r9323 = r9321 <= r9322;
        float r9324 = c0;
        float r9325 = l;
        float r9326 = r9325 * r9320;
        float r9327 = r9319 / r9326;
        float r9328 = sqrt(r9327);
        float r9329 = r9324 * r9328;
        float r9330 = -1.3501087918197612e-111f;
        bool r9331 = r9321 <= r9330;
        float r9332 = 1.0f;
        float r9333 = r9325 / r9321;
        float r9334 = r9332 / r9333;
        float r9335 = sqrt(r9334);
        float r9336 = r9324 * r9335;
        float r9337 = 9.043080568084718e-301f;
        bool r9338 = r9321 <= r9337;
        float r9339 = r9332 / r9325;
        float r9340 = sqrt(r9339);
        float r9341 = sqrt(r9321);
        float r9342 = r9340 * r9341;
        float r9343 = r9342 * r9324;
        float r9344 = r9338 ? r9329 : r9343;
        float r9345 = r9331 ? r9336 : r9344;
        float r9346 = r9323 ? r9329 : r9345;
        return r9346;
}

double f_od(double c0, double A, double V, double l) {
        double r9347 = A;
        double r9348 = V;
        double r9349 = r9347 / r9348;
        double r9350 = -2.004776637139463e+296;
        bool r9351 = r9349 <= r9350;
        double r9352 = c0;
        double r9353 = l;
        double r9354 = r9353 * r9348;
        double r9355 = r9347 / r9354;
        double r9356 = sqrt(r9355);
        double r9357 = r9352 * r9356;
        double r9358 = -1.3501087918197612e-111;
        bool r9359 = r9349 <= r9358;
        double r9360 = 1.0;
        double r9361 = r9353 / r9349;
        double r9362 = r9360 / r9361;
        double r9363 = sqrt(r9362);
        double r9364 = r9352 * r9363;
        double r9365 = 9.043080568084718e-301;
        bool r9366 = r9349 <= r9365;
        double r9367 = r9360 / r9353;
        double r9368 = sqrt(r9367);
        double r9369 = sqrt(r9349);
        double r9370 = r9368 * r9369;
        double r9371 = r9370 * r9352;
        double r9372 = r9366 ? r9357 : r9371;
        double r9373 = r9359 ? r9364 : r9372;
        double r9374 = r9351 ? r9357 : r9373;
        return r9374;
}

void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) {
        mpfr_fmod(r, n, d, rmd);
        if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd);
}


static mpfr_t r9375, r9376, r9377, r9378, r9379, r9380, r9381, r9382;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9375);
        mpfr_init(r9376);
        mpfr_init(r9377);
        mpfr_init(r9378);
        mpfr_init(r9379);
        mpfr_init(r9380);
        mpfr_init(r9381);
        mpfr_init(r9382);
}

double f_im(double c0, double A, double V, double l) {
        mpfr_set_d(r9375, c0, MPFR_RNDN);
        mpfr_set_d(r9376, A, MPFR_RNDN);
        mpfr_set_d(r9377, V, MPFR_RNDN);
        mpfr_set_d(r9378, l, MPFR_RNDN);
        mpfr_mul(r9379, r9377, r9378, MPFR_RNDN);
        mpfr_div(r9380, r9376, r9379, MPFR_RNDN);
        mpfr_sqrt(r9381, r9380, MPFR_RNDN);
        mpfr_mul(r9382, r9375, r9381, MPFR_RNDN);
        return mpfr_get_d(r9382, MPFR_RNDN);
}

static mpfr_t r9383, r9384, r9385, r9386, r9387, r9388, r9389, r9390, r9391, r9392, r9393, r9394, r9395, r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409, r9410;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9383);
        mpfr_init(r9384);
        mpfr_init(r9385);
        mpfr_init_set_str(r9386, "-2.004776637139463e+296", 10, MPFR_RNDN);
        mpfr_init(r9387);
        mpfr_init(r9388);
        mpfr_init(r9389);
        mpfr_init(r9390);
        mpfr_init(r9391);
        mpfr_init(r9392);
        mpfr_init(r9393);
        mpfr_init_set_str(r9394, "-1.3501087918197612e-111", 10, MPFR_RNDN);
        mpfr_init(r9395);
        mpfr_init_set_str(r9396, "1", 10, MPFR_RNDN);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init_set_str(r9401, "9.043080568084718e-301", 10, MPFR_RNDN);
        mpfr_init(r9402);
        mpfr_init(r9403);
        mpfr_init(r9404);
        mpfr_init(r9405);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
        mpfr_init(r9410);
}

double f_fm(double c0, double A, double V, double l) {
        mpfr_set_d(r9383, A, MPFR_RNDN);
        mpfr_set_d(r9384, V, MPFR_RNDN);
        mpfr_div(r9385, r9383, r9384, MPFR_RNDN);
        ;
        mpfr_set_si(r9387, mpfr_cmp(r9385, r9386) <= 0, MPFR_RNDN);
        mpfr_set_d(r9388, c0, MPFR_RNDN);
        mpfr_set_d(r9389, l, MPFR_RNDN);
        mpfr_mul(r9390, r9389, r9384, MPFR_RNDN);
        mpfr_div(r9391, r9383, r9390, MPFR_RNDN);
        mpfr_sqrt(r9392, r9391, MPFR_RNDN);
        mpfr_mul(r9393, r9388, r9392, MPFR_RNDN);
        ;
        mpfr_set_si(r9395, mpfr_cmp(r9385, r9394) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9397, r9389, r9385, MPFR_RNDN);
        mpfr_div(r9398, r9396, r9397, MPFR_RNDN);
        mpfr_sqrt(r9399, r9398, MPFR_RNDN);
        mpfr_mul(r9400, r9388, r9399, MPFR_RNDN);
        ;
        mpfr_set_si(r9402, mpfr_cmp(r9385, r9401) <= 0, MPFR_RNDN);
        mpfr_div(r9403, r9396, r9389, MPFR_RNDN);
        mpfr_sqrt(r9404, r9403, MPFR_RNDN);
        mpfr_sqrt(r9405, r9385, MPFR_RNDN);
        mpfr_mul(r9406, r9404, r9405, MPFR_RNDN);
        mpfr_mul(r9407, r9406, r9388, MPFR_RNDN);
        if (mpfr_get_si(r9402, MPFR_RNDN)) { mpfr_set(r9408, r9393, MPFR_RNDN); } else { mpfr_set(r9408, r9407, MPFR_RNDN); };
        if (mpfr_get_si(r9395, MPFR_RNDN)) { mpfr_set(r9409, r9400, MPFR_RNDN); } else { mpfr_set(r9409, r9408, MPFR_RNDN); };
        if (mpfr_get_si(r9387, MPFR_RNDN)) { mpfr_set(r9410, r9393, MPFR_RNDN); } else { mpfr_set(r9410, r9409, MPFR_RNDN); };
        return mpfr_get_d(r9410, MPFR_RNDN);
}

static mpfr_t r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423, r9424, r9425, r9426, r9427, r9428, r9429, r9430, r9431, r9432, r9433, r9434, r9435, r9436, r9437, r9438;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9411);
        mpfr_init(r9412);
        mpfr_init(r9413);
        mpfr_init_set_str(r9414, "-2.004776637139463e+296", 10, MPFR_RNDN);
        mpfr_init(r9415);
        mpfr_init(r9416);
        mpfr_init(r9417);
        mpfr_init(r9418);
        mpfr_init(r9419);
        mpfr_init(r9420);
        mpfr_init(r9421);
        mpfr_init_set_str(r9422, "-1.3501087918197612e-111", 10, MPFR_RNDN);
        mpfr_init(r9423);
        mpfr_init_set_str(r9424, "1", 10, MPFR_RNDN);
        mpfr_init(r9425);
        mpfr_init(r9426);
        mpfr_init(r9427);
        mpfr_init(r9428);
        mpfr_init_set_str(r9429, "9.043080568084718e-301", 10, MPFR_RNDN);
        mpfr_init(r9430);
        mpfr_init(r9431);
        mpfr_init(r9432);
        mpfr_init(r9433);
        mpfr_init(r9434);
        mpfr_init(r9435);
        mpfr_init(r9436);
        mpfr_init(r9437);
        mpfr_init(r9438);
}

double f_dm(double c0, double A, double V, double l) {
        mpfr_set_d(r9411, A, MPFR_RNDN);
        mpfr_set_d(r9412, V, MPFR_RNDN);
        mpfr_div(r9413, r9411, r9412, MPFR_RNDN);
        ;
        mpfr_set_si(r9415, mpfr_cmp(r9413, r9414) <= 0, MPFR_RNDN);
        mpfr_set_d(r9416, c0, MPFR_RNDN);
        mpfr_set_d(r9417, l, MPFR_RNDN);
        mpfr_mul(r9418, r9417, r9412, MPFR_RNDN);
        mpfr_div(r9419, r9411, r9418, MPFR_RNDN);
        mpfr_sqrt(r9420, r9419, MPFR_RNDN);
        mpfr_mul(r9421, r9416, r9420, MPFR_RNDN);
        ;
        mpfr_set_si(r9423, mpfr_cmp(r9413, r9422) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9425, r9417, r9413, MPFR_RNDN);
        mpfr_div(r9426, r9424, r9425, MPFR_RNDN);
        mpfr_sqrt(r9427, r9426, MPFR_RNDN);
        mpfr_mul(r9428, r9416, r9427, MPFR_RNDN);
        ;
        mpfr_set_si(r9430, mpfr_cmp(r9413, r9429) <= 0, MPFR_RNDN);
        mpfr_div(r9431, r9424, r9417, MPFR_RNDN);
        mpfr_sqrt(r9432, r9431, MPFR_RNDN);
        mpfr_sqrt(r9433, r9413, MPFR_RNDN);
        mpfr_mul(r9434, r9432, r9433, MPFR_RNDN);
        mpfr_mul(r9435, r9434, r9416, MPFR_RNDN);
        if (mpfr_get_si(r9430, MPFR_RNDN)) { mpfr_set(r9436, r9421, MPFR_RNDN); } else { mpfr_set(r9436, r9435, MPFR_RNDN); };
        if (mpfr_get_si(r9423, MPFR_RNDN)) { mpfr_set(r9437, r9428, MPFR_RNDN); } else { mpfr_set(r9437, r9436, MPFR_RNDN); };
        if (mpfr_get_si(r9415, MPFR_RNDN)) { mpfr_set(r9438, r9421, MPFR_RNDN); } else { mpfr_set(r9438, r9437, MPFR_RNDN); };
        return mpfr_get_d(r9438, MPFR_RNDN);
}