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

char *name = "quadp (p42, positive)";

double f_if(float a, float b, float c) {
        float r9565 = b;
        float r9566 = -r9565;
        float r9567 = r9565 * r9565;
        float r9568 = 4.0f;
        float r9569 = a;
        float r9570 = c;
        float r9571 = r9569 * r9570;
        float r9572 = r9568 * r9571;
        float r9573 = r9567 - r9572;
        float r9574 = sqrt(r9573);
        float r9575 = r9566 + r9574;
        float r9576 = 2.0f;
        float r9577 = r9576 * r9569;
        float r9578 = r9575 / r9577;
        return r9578;
}

double f_id(double a, double b, double c) {
        double r9579 = b;
        double r9580 = -r9579;
        double r9581 = r9579 * r9579;
        double r9582 = 4.0;
        double r9583 = a;
        double r9584 = c;
        double r9585 = r9583 * r9584;
        double r9586 = r9582 * r9585;
        double r9587 = r9581 - r9586;
        double r9588 = sqrt(r9587);
        double r9589 = r9580 + r9588;
        double r9590 = 2.0;
        double r9591 = r9590 * r9583;
        double r9592 = r9589 / r9591;
        return r9592;
}


double f_of(float a, float b, float c) {
        float r9593 = b;
        float r9594 = -1.3395204081552134e+154f;
        bool r9595 = r9593 <= r9594;
        float r9596 = -0.5f;
        float r9597 = a;
        float r9598 = r9593 / r9597;
        float r9599 = r9596 * r9598;
        float r9600 = -1.96292444197926e-264f;
        bool r9601 = r9593 <= r9600;
        float r9602 = 4.0f;
        float r9603 = -r9602;
        float r9604 = c;
        float r9605 = r9604 * r9597;
        float r9606 = r9593 * r9593;
        float r9607 = fma(r9603, r9605, r9606);
        float r9608 = sqrt(r9607);
        float r9609 = r9608 - r9593;
        float r9610 = 1.0f;
        float r9611 = 2.0f;
        float r9612 = r9597 * r9611;
        float r9613 = r9610 / r9612;
        float r9614 = r9609 * r9613;
        float r9615 = 1.473326279994868e+84f;
        bool r9616 = r9593 <= r9615;
        float r9617 = -2.0f;
        float r9618 = r9604 * r9617;
        float r9619 = pow(r9593, r9611);
        float r9620 = r9602 * r9605;
        float r9621 = r9619 - r9620;
        float r9622 = sqrt(r9621);
        float r9623 = r9622 + r9593;
        float r9624 = r9618 / r9623;
        float r9625 = r9593 + r9593;
        float r9626 = r9618 / r9625;
        float r9627 = r9616 ? r9624 : r9626;
        float r9628 = r9601 ? r9614 : r9627;
        float r9629 = r9595 ? r9599 : r9628;
        return r9629;
}

double f_od(double a, double b, double c) {
        double r9630 = b;
        double r9631 = -1.3395204081552134e+154;
        bool r9632 = r9630 <= r9631;
        double r9633 = -0.5;
        double r9634 = a;
        double r9635 = r9630 / r9634;
        double r9636 = r9633 * r9635;
        double r9637 = -1.96292444197926e-264;
        bool r9638 = r9630 <= r9637;
        double r9639 = 4.0;
        double r9640 = -r9639;
        double r9641 = c;
        double r9642 = r9641 * r9634;
        double r9643 = r9630 * r9630;
        double r9644 = fma(r9640, r9642, r9643);
        double r9645 = sqrt(r9644);
        double r9646 = r9645 - r9630;
        double r9647 = 1.0;
        double r9648 = 2.0;
        double r9649 = r9634 * r9648;
        double r9650 = r9647 / r9649;
        double r9651 = r9646 * r9650;
        double r9652 = 1.473326279994868e+84;
        bool r9653 = r9630 <= r9652;
        double r9654 = -2.0;
        double r9655 = r9641 * r9654;
        double r9656 = pow(r9630, r9648);
        double r9657 = r9639 * r9642;
        double r9658 = r9656 - r9657;
        double r9659 = sqrt(r9658);
        double r9660 = r9659 + r9630;
        double r9661 = r9655 / r9660;
        double r9662 = r9630 + r9630;
        double r9663 = r9655 / r9662;
        double r9664 = r9653 ? r9661 : r9663;
        double r9665 = r9638 ? r9651 : r9664;
        double r9666 = r9632 ? r9636 : r9665;
        return r9666;
}

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 r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r9667);
        mpfr_init(r9668);
        mpfr_init(r9669);
        mpfr_init_set_str(r9670, "4", 10, MPFR_RNDN);
        mpfr_init(r9671);
        mpfr_init(r9672);
        mpfr_init(r9673);
        mpfr_init(r9674);
        mpfr_init(r9675);
        mpfr_init(r9676);
        mpfr_init(r9677);
        mpfr_init_set_str(r9678, "2", 10, MPFR_RNDN);
        mpfr_init(r9679);
        mpfr_init(r9680);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r9667, b, MPFR_RNDN);
        mpfr_neg(r9668, r9667, MPFR_RNDN);
        mpfr_mul(r9669, r9667, r9667, MPFR_RNDN);
        ;
        mpfr_set_d(r9671, a, MPFR_RNDN);
        mpfr_set_d(r9672, c, MPFR_RNDN);
        mpfr_mul(r9673, r9671, r9672, MPFR_RNDN);
        mpfr_mul(r9674, r9670, r9673, MPFR_RNDN);
        mpfr_sub(r9675, r9669, r9674, MPFR_RNDN);
        mpfr_sqrt(r9676, r9675, MPFR_RNDN);
        mpfr_add(r9677, r9668, r9676, MPFR_RNDN);
        ;
        mpfr_mul(r9679, r9678, r9671, MPFR_RNDN);
        mpfr_div(r9680, r9677, r9679, MPFR_RNDN);
        return mpfr_get_d(r9680, MPFR_RNDN);
}

static mpfr_t r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689, r9690, r9691, r9692, r9693, r9694, r9695, r9696, r9697, r9698, r9699, r9700, r9701, r9702, r9703, r9704, r9705, r9706, r9707, r9708, r9709, r9710, r9711, r9712, r9713, r9714, r9715, r9716, r9717;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r9681);
        mpfr_init_set_str(r9682, "-1.3395204081552134e+154", 10, MPFR_RNDN);
        mpfr_init(r9683);
        mpfr_init_set_str(r9684, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r9685);
        mpfr_init(r9686);
        mpfr_init(r9687);
        mpfr_init_set_str(r9688, "-1.96292444197926e-264", 10, MPFR_RNDN);
        mpfr_init(r9689);
        mpfr_init_set_str(r9690, "4", 10, MPFR_RNDN);
        mpfr_init(r9691);
        mpfr_init(r9692);
        mpfr_init(r9693);
        mpfr_init(r9694);
        mpfr_init(r9695);
        mpfr_init(r9696);
        mpfr_init(r9697);
        mpfr_init_set_str(r9698, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r9699, "2", 10, MPFR_RNDN);
        mpfr_init(r9700);
        mpfr_init(r9701);
        mpfr_init(r9702);
        mpfr_init_set_str(r9703, "1.473326279994868e+84", 10, MPFR_RNDN);
        mpfr_init(r9704);
        mpfr_init_set_str(r9705, "-2", 10, MPFR_RNDN);
        mpfr_init(r9706);
        mpfr_init(r9707);
        mpfr_init(r9708);
        mpfr_init(r9709);
        mpfr_init(r9710);
        mpfr_init(r9711);
        mpfr_init(r9712);
        mpfr_init(r9713);
        mpfr_init(r9714);
        mpfr_init(r9715);
        mpfr_init(r9716);
        mpfr_init(r9717);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r9681, b, MPFR_RNDN);
        ;
        mpfr_set_si(r9683, mpfr_cmp(r9681, r9682) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r9685, a, MPFR_RNDN);
        mpfr_div(r9686, r9681, r9685, MPFR_RNDN);
        mpfr_mul(r9687, r9684, r9686, MPFR_RNDN);
        ;
        mpfr_set_si(r9689, mpfr_cmp(r9681, r9688) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r9691, r9690, MPFR_RNDN);
        mpfr_set_d(r9692, c, MPFR_RNDN);
        mpfr_mul(r9693, r9692, r9685, MPFR_RNDN);
        mpfr_mul(r9694, r9681, r9681, MPFR_RNDN);
        mpfr_fma(r9695, r9691, r9693, r9694, MPFR_RNDN);
        mpfr_sqrt(r9696, r9695, MPFR_RNDN);
        mpfr_sub(r9697, r9696, r9681, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r9700, r9685, r9699, MPFR_RNDN);
        mpfr_div(r9701, r9698, r9700, MPFR_RNDN);
        mpfr_mul(r9702, r9697, r9701, MPFR_RNDN);
        ;
        mpfr_set_si(r9704, mpfr_cmp(r9681, r9703) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r9706, r9692, r9705, MPFR_RNDN);
        mpfr_pow(r9707, r9681, r9699, MPFR_RNDN);
        mpfr_mul(r9708, r9690, r9693, MPFR_RNDN);
        mpfr_sub(r9709, r9707, r9708, MPFR_RNDN);
        mpfr_sqrt(r9710, r9709, MPFR_RNDN);
        mpfr_add(r9711, r9710, r9681, MPFR_RNDN);
        mpfr_div(r9712, r9706, r9711, MPFR_RNDN);
        mpfr_add(r9713, r9681, r9681, MPFR_RNDN);
        mpfr_div(r9714, r9706, r9713, MPFR_RNDN);
        if (mpfr_get_si(r9704, MPFR_RNDN)) { mpfr_set(r9715, r9712, MPFR_RNDN); } else { mpfr_set(r9715, r9714, MPFR_RNDN); };
        if (mpfr_get_si(r9689, MPFR_RNDN)) { mpfr_set(r9716, r9702, MPFR_RNDN); } else { mpfr_set(r9716, r9715, MPFR_RNDN); };
        if (mpfr_get_si(r9683, MPFR_RNDN)) { mpfr_set(r9717, r9687, MPFR_RNDN); } else { mpfr_set(r9717, r9716, MPFR_RNDN); };
        return mpfr_get_d(r9717, MPFR_RNDN);
}

static mpfr_t r9718, r9719, r9720, r9721, r9722, r9723, r9724, r9725, r9726, r9727, r9728, r9729, r9730, r9731, r9732, r9733, r9734, r9735, r9736, r9737, r9738, r9739, r9740, r9741, r9742, r9743, r9744, r9745, r9746, r9747, r9748, r9749, r9750, r9751, r9752, r9753, r9754;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r9718);
        mpfr_init_set_str(r9719, "-1.3395204081552134e+154", 10, MPFR_RNDN);
        mpfr_init(r9720);
        mpfr_init_set_str(r9721, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r9722);
        mpfr_init(r9723);
        mpfr_init(r9724);
        mpfr_init_set_str(r9725, "-1.96292444197926e-264", 10, MPFR_RNDN);
        mpfr_init(r9726);
        mpfr_init_set_str(r9727, "4", 10, MPFR_RNDN);
        mpfr_init(r9728);
        mpfr_init(r9729);
        mpfr_init(r9730);
        mpfr_init(r9731);
        mpfr_init(r9732);
        mpfr_init(r9733);
        mpfr_init(r9734);
        mpfr_init_set_str(r9735, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r9736, "2", 10, MPFR_RNDN);
        mpfr_init(r9737);
        mpfr_init(r9738);
        mpfr_init(r9739);
        mpfr_init_set_str(r9740, "1.473326279994868e+84", 10, MPFR_RNDN);
        mpfr_init(r9741);
        mpfr_init_set_str(r9742, "-2", 10, MPFR_RNDN);
        mpfr_init(r9743);
        mpfr_init(r9744);
        mpfr_init(r9745);
        mpfr_init(r9746);
        mpfr_init(r9747);
        mpfr_init(r9748);
        mpfr_init(r9749);
        mpfr_init(r9750);
        mpfr_init(r9751);
        mpfr_init(r9752);
        mpfr_init(r9753);
        mpfr_init(r9754);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r9718, b, MPFR_RNDN);
        ;
        mpfr_set_si(r9720, mpfr_cmp(r9718, r9719) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r9722, a, MPFR_RNDN);
        mpfr_div(r9723, r9718, r9722, MPFR_RNDN);
        mpfr_mul(r9724, r9721, r9723, MPFR_RNDN);
        ;
        mpfr_set_si(r9726, mpfr_cmp(r9718, r9725) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r9728, r9727, MPFR_RNDN);
        mpfr_set_d(r9729, c, MPFR_RNDN);
        mpfr_mul(r9730, r9729, r9722, MPFR_RNDN);
        mpfr_mul(r9731, r9718, r9718, MPFR_RNDN);
        mpfr_fma(r9732, r9728, r9730, r9731, MPFR_RNDN);
        mpfr_sqrt(r9733, r9732, MPFR_RNDN);
        mpfr_sub(r9734, r9733, r9718, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r9737, r9722, r9736, MPFR_RNDN);
        mpfr_div(r9738, r9735, r9737, MPFR_RNDN);
        mpfr_mul(r9739, r9734, r9738, MPFR_RNDN);
        ;
        mpfr_set_si(r9741, mpfr_cmp(r9718, r9740) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r9743, r9729, r9742, MPFR_RNDN);
        mpfr_pow(r9744, r9718, r9736, MPFR_RNDN);
        mpfr_mul(r9745, r9727, r9730, MPFR_RNDN);
        mpfr_sub(r9746, r9744, r9745, MPFR_RNDN);
        mpfr_sqrt(r9747, r9746, MPFR_RNDN);
        mpfr_add(r9748, r9747, r9718, MPFR_RNDN);
        mpfr_div(r9749, r9743, r9748, MPFR_RNDN);
        mpfr_add(r9750, r9718, r9718, MPFR_RNDN);
        mpfr_div(r9751, r9743, r9750, MPFR_RNDN);
        if (mpfr_get_si(r9741, MPFR_RNDN)) { mpfr_set(r9752, r9749, MPFR_RNDN); } else { mpfr_set(r9752, r9751, MPFR_RNDN); };
        if (mpfr_get_si(r9726, MPFR_RNDN)) { mpfr_set(r9753, r9739, MPFR_RNDN); } else { mpfr_set(r9753, r9752, MPFR_RNDN); };
        if (mpfr_get_si(r9720, MPFR_RNDN)) { mpfr_set(r9754, r9724, MPFR_RNDN); } else { mpfr_set(r9754, r9753, MPFR_RNDN); };
        return mpfr_get_d(r9754, MPFR_RNDN);
}