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

char *name = "quad2p (problem 3.2.1, positive)";

double f_if(float a, float b_2F2, float c) {
        float r18522 = b_2F2;
        float r18523 = -r18522;
        float r18524 = r18522 * r18522;
        float r18525 = a;
        float r18526 = c;
        float r18527 = r18525 * r18526;
        float r18528 = r18524 - r18527;
        float r18529 = sqrt(r18528);
        float r18530 = r18523 + r18529;
        float r18531 = r18530 / r18525;
        return r18531;
}

double f_id(double a, double b_2F2, double c) {
        double r18532 = b_2F2;
        double r18533 = -r18532;
        double r18534 = r18532 * r18532;
        double r18535 = a;
        double r18536 = c;
        double r18537 = r18535 * r18536;
        double r18538 = r18534 - r18537;
        double r18539 = sqrt(r18538);
        double r18540 = r18533 + r18539;
        double r18541 = r18540 / r18535;
        return r18541;
}


double f_of(float a, float b_2F2, float c) {
        float r18542 = b_2F2;
        float r18543 = -1.6573860688298954e+150f;
        bool r18544 = r18542 <= r18543;
        float r18545 = -2.0f;
        float r18546 = a;
        float r18547 = r18542 / r18546;
        float r18548 = r18545 * r18547;
        float r18549 = 1.828821921218125e-208f;
        bool r18550 = r18542 <= r18549;
        float r18551 = -r18542;
        float r18552 = r18542 * r18542;
        float r18553 = c;
        float r18554 = r18546 * r18553;
        float r18555 = r18552 - r18554;
        float r18556 = sqrt(r18555);
        float r18557 = r18551 + r18556;
        float r18558 = r18557 / r18546;
        float r18559 = 7.666985785588564e+91f;
        bool r18560 = r18542 <= r18559;
        float r18561 = r18553 * r18546;
        float r18562 = r18551 - r18556;
        float r18563 = r18561 / r18562;
        float r18564 = r18563 / r18546;
        float r18565 = r18553 / r18542;
        float r18566 = -0.5f;
        float r18567 = r18565 * r18566;
        float r18568 = r18560 ? r18564 : r18567;
        float r18569 = r18550 ? r18558 : r18568;
        float r18570 = r18544 ? r18548 : r18569;
        return r18570;
}

double f_od(double a, double b_2F2, double c) {
        double r18571 = b_2F2;
        double r18572 = -1.6573860688298954e+150;
        bool r18573 = r18571 <= r18572;
        double r18574 = -2.0;
        double r18575 = a;
        double r18576 = r18571 / r18575;
        double r18577 = r18574 * r18576;
        double r18578 = 1.828821921218125e-208;
        bool r18579 = r18571 <= r18578;
        double r18580 = -r18571;
        double r18581 = r18571 * r18571;
        double r18582 = c;
        double r18583 = r18575 * r18582;
        double r18584 = r18581 - r18583;
        double r18585 = sqrt(r18584);
        double r18586 = r18580 + r18585;
        double r18587 = r18586 / r18575;
        double r18588 = 7.666985785588564e+91;
        bool r18589 = r18571 <= r18588;
        double r18590 = r18582 * r18575;
        double r18591 = r18580 - r18585;
        double r18592 = r18590 / r18591;
        double r18593 = r18592 / r18575;
        double r18594 = r18582 / r18571;
        double r18595 = -0.5;
        double r18596 = r18594 * r18595;
        double r18597 = r18589 ? r18593 : r18596;
        double r18598 = r18579 ? r18587 : r18597;
        double r18599 = r18573 ? r18577 : r18598;
        return r18599;
}

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 r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2704);
        mpfr_init(r18600);
        mpfr_init(r18601);
        mpfr_init(r18602);
        mpfr_init(r18603);
        mpfr_init(r18604);
        mpfr_init(r18605);
        mpfr_init(r18606);
        mpfr_init(r18607);
        mpfr_init(r18608);
        mpfr_init(r18609);
}

double f_im(double a, double b_2F2, double c) {
        mpfr_set_d(r18600, b_2F2, MPFR_RNDN);
        mpfr_neg(r18601, r18600, MPFR_RNDN);
        mpfr_sqr(r18602, r18600, MPFR_RNDN);
        mpfr_set_d(r18603, a, MPFR_RNDN);
        mpfr_set_d(r18604, c, MPFR_RNDN);
        mpfr_mul(r18605, r18603, r18604, MPFR_RNDN);
        mpfr_sub(r18606, r18602, r18605, MPFR_RNDN);
        mpfr_sqrt(r18607, r18606, MPFR_RNDN);
        mpfr_add(r18608, r18601, r18607, MPFR_RNDN);
        mpfr_div(r18609, r18608, r18603, MPFR_RNDN);
        return mpfr_get_d(r18609, MPFR_RNDN);
}

static mpfr_t r18610, r18611, r18612, r18613, r18614, r18615, r18616, r18617, r18618, r18619, r18620, r18621, r18622, r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634, r18635, r18636, r18637, r18638;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2704);
        mpfr_init(r18610);
        mpfr_init_set_str(r18611, "-1.6573860688298954e+150", 10, MPFR_RNDN);
        mpfr_init(r18612);
        mpfr_init_set_str(r18613, "-2", 10, MPFR_RNDN);
        mpfr_init(r18614);
        mpfr_init(r18615);
        mpfr_init(r18616);
        mpfr_init_set_str(r18617, "1.828821921218125e-208", 10, MPFR_RNDN);
        mpfr_init(r18618);
        mpfr_init(r18619);
        mpfr_init(r18620);
        mpfr_init(r18621);
        mpfr_init(r18622);
        mpfr_init(r18623);
        mpfr_init(r18624);
        mpfr_init(r18625);
        mpfr_init(r18626);
        mpfr_init_set_str(r18627, "7.666985785588564e+91", 10, MPFR_RNDN);
        mpfr_init(r18628);
        mpfr_init(r18629);
        mpfr_init(r18630);
        mpfr_init(r18631);
        mpfr_init(r18632);
        mpfr_init(r18633);
        mpfr_init_set_str(r18634, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r18635);
        mpfr_init(r18636);
        mpfr_init(r18637);
        mpfr_init(r18638);
}

double f_fm(double a, double b_2F2, double c) {
        mpfr_set_d(r18610, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r18612, mpfr_cmp(r18610, r18611) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r18614, a, MPFR_RNDN);
        mpfr_div(r18615, r18610, r18614, MPFR_RNDN);
        mpfr_mul(r18616, r18613, r18615, MPFR_RNDN);
        ;
        mpfr_set_si(r18618, mpfr_cmp(r18610, r18617) <= 0, MPFR_RNDN);
        mpfr_neg(r18619, r18610, MPFR_RNDN);
        mpfr_sqr(r18620, r18610, MPFR_RNDN);
        mpfr_set_d(r18621, c, MPFR_RNDN);
        mpfr_mul(r18622, r18614, r18621, MPFR_RNDN);
        mpfr_sub(r18623, r18620, r18622, MPFR_RNDN);
        mpfr_sqrt(r18624, r18623, MPFR_RNDN);
        mpfr_add(r18625, r18619, r18624, MPFR_RNDN);
        mpfr_div(r18626, r18625, r18614, MPFR_RNDN);
        ;
        mpfr_set_si(r18628, mpfr_cmp(r18610, r18627) <= 0, MPFR_RNDN);
        mpfr_mul(r18629, r18621, r18614, MPFR_RNDN);
        mpfr_sub(r18630, r18619, r18624, MPFR_RNDN);
        mpfr_div(r18631, r18629, r18630, MPFR_RNDN);
        mpfr_div(r18632, r18631, r18614, MPFR_RNDN);
        mpfr_div(r18633, r18621, r18610, MPFR_RNDN);
        ;
        mpfr_mul(r18635, r18633, r18634, MPFR_RNDN);
        if (mpfr_get_si(r18628, MPFR_RNDN)) { mpfr_set(r18636, r18632, MPFR_RNDN); } else { mpfr_set(r18636, r18635, MPFR_RNDN); };
        if (mpfr_get_si(r18618, MPFR_RNDN)) { mpfr_set(r18637, r18626, MPFR_RNDN); } else { mpfr_set(r18637, r18636, MPFR_RNDN); };
        if (mpfr_get_si(r18612, MPFR_RNDN)) { mpfr_set(r18638, r18616, MPFR_RNDN); } else { mpfr_set(r18638, r18637, MPFR_RNDN); };
        return mpfr_get_d(r18638, MPFR_RNDN);
}

static mpfr_t r18639, r18640, r18641, r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2704);
        mpfr_init(r18639);
        mpfr_init_set_str(r18640, "-1.6573860688298954e+150", 10, MPFR_RNDN);
        mpfr_init(r18641);
        mpfr_init_set_str(r18642, "-2", 10, MPFR_RNDN);
        mpfr_init(r18643);
        mpfr_init(r18644);
        mpfr_init(r18645);
        mpfr_init_set_str(r18646, "1.828821921218125e-208", 10, MPFR_RNDN);
        mpfr_init(r18647);
        mpfr_init(r18648);
        mpfr_init(r18649);
        mpfr_init(r18650);
        mpfr_init(r18651);
        mpfr_init(r18652);
        mpfr_init(r18653);
        mpfr_init(r18654);
        mpfr_init(r18655);
        mpfr_init_set_str(r18656, "7.666985785588564e+91", 10, MPFR_RNDN);
        mpfr_init(r18657);
        mpfr_init(r18658);
        mpfr_init(r18659);
        mpfr_init(r18660);
        mpfr_init(r18661);
        mpfr_init(r18662);
        mpfr_init_set_str(r18663, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r18664);
        mpfr_init(r18665);
        mpfr_init(r18666);
        mpfr_init(r18667);
}

double f_dm(double a, double b_2F2, double c) {
        mpfr_set_d(r18639, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r18641, mpfr_cmp(r18639, r18640) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r18643, a, MPFR_RNDN);
        mpfr_div(r18644, r18639, r18643, MPFR_RNDN);
        mpfr_mul(r18645, r18642, r18644, MPFR_RNDN);
        ;
        mpfr_set_si(r18647, mpfr_cmp(r18639, r18646) <= 0, MPFR_RNDN);
        mpfr_neg(r18648, r18639, MPFR_RNDN);
        mpfr_sqr(r18649, r18639, MPFR_RNDN);
        mpfr_set_d(r18650, c, MPFR_RNDN);
        mpfr_mul(r18651, r18643, r18650, MPFR_RNDN);
        mpfr_sub(r18652, r18649, r18651, MPFR_RNDN);
        mpfr_sqrt(r18653, r18652, MPFR_RNDN);
        mpfr_add(r18654, r18648, r18653, MPFR_RNDN);
        mpfr_div(r18655, r18654, r18643, MPFR_RNDN);
        ;
        mpfr_set_si(r18657, mpfr_cmp(r18639, r18656) <= 0, MPFR_RNDN);
        mpfr_mul(r18658, r18650, r18643, MPFR_RNDN);
        mpfr_sub(r18659, r18648, r18653, MPFR_RNDN);
        mpfr_div(r18660, r18658, r18659, MPFR_RNDN);
        mpfr_div(r18661, r18660, r18643, MPFR_RNDN);
        mpfr_div(r18662, r18650, r18639, MPFR_RNDN);
        ;
        mpfr_mul(r18664, r18662, r18663, MPFR_RNDN);
        if (mpfr_get_si(r18657, MPFR_RNDN)) { mpfr_set(r18665, r18661, MPFR_RNDN); } else { mpfr_set(r18665, r18664, MPFR_RNDN); };
        if (mpfr_get_si(r18647, MPFR_RNDN)) { mpfr_set(r18666, r18655, MPFR_RNDN); } else { mpfr_set(r18666, r18665, MPFR_RNDN); };
        if (mpfr_get_si(r18641, MPFR_RNDN)) { mpfr_set(r18667, r18645, MPFR_RNDN); } else { mpfr_set(r18667, r18666, MPFR_RNDN); };
        return mpfr_get_d(r18667, MPFR_RNDN);
}