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

char *name = "The quadratic formula (r1)";

double f_if(float a, float b, float c) {
        float r24528 = b;
        float r24529 = -r24528;
        float r24530 = r24528 * r24528;
        float r24531 = 4;
        float r24532 = a;
        float r24533 = r24531 * r24532;
        float r24534 = c;
        float r24535 = r24533 * r24534;
        float r24536 = r24530 - r24535;
        float r24537 = sqrt(r24536);
        float r24538 = r24529 + r24537;
        float r24539 = 2;
        float r24540 = r24539 * r24532;
        float r24541 = r24538 / r24540;
        return r24541;
}

double f_id(double a, double b, double c) {
        double r24542 = b;
        double r24543 = -r24542;
        double r24544 = r24542 * r24542;
        double r24545 = 4;
        double r24546 = a;
        double r24547 = r24545 * r24546;
        double r24548 = c;
        double r24549 = r24547 * r24548;
        double r24550 = r24544 - r24549;
        double r24551 = sqrt(r24550);
        double r24552 = r24543 + r24551;
        double r24553 = 2;
        double r24554 = r24553 * r24546;
        double r24555 = r24552 / r24554;
        return r24555;
}


double f_of(float a, float b, float c) {
        float r24556 = b;
        float r24557 = -1.5023809489540946e+154;
        bool r24558 = r24556 <= r24557;
        float r24559 = -1/2;
        float r24560 = a;
        float r24561 = r24556 / r24560;
        float r24562 = r24559 * r24561;
        float r24563 = 5.761969426420715e-61;
        bool r24564 = r24556 <= r24563;
        float r24565 = 4;
        float r24566 = r24565 * r24560;
        float r24567 = c;
        float r24568 = -r24567;
        float r24569 = r24556 * r24556;
        float r24570 = fma(r24566, r24568, r24569);
        float r24571 = sqrt(r24570);
        float r24572 = sqrt(r24571);
        float r24573 = -r24556;
        float r24574 = fma(r24572, r24572, r24573);
        float r24575 = 2;
        float r24576 = r24575 * r24560;
        float r24577 = r24574 / r24576;
        float r24578 = r24567 / r24573;
        float r24579 = r24564 ? r24577 : r24578;
        float r24580 = r24558 ? r24562 : r24579;
        return r24580;
}

double f_od(double a, double b, double c) {
        double r24581 = b;
        double r24582 = -1.5023809489540946e+154;
        bool r24583 = r24581 <= r24582;
        double r24584 = -1/2;
        double r24585 = a;
        double r24586 = r24581 / r24585;
        double r24587 = r24584 * r24586;
        double r24588 = 5.761969426420715e-61;
        bool r24589 = r24581 <= r24588;
        double r24590 = 4;
        double r24591 = r24590 * r24585;
        double r24592 = c;
        double r24593 = -r24592;
        double r24594 = r24581 * r24581;
        double r24595 = fma(r24591, r24593, r24594);
        double r24596 = sqrt(r24595);
        double r24597 = sqrt(r24596);
        double r24598 = -r24581;
        double r24599 = fma(r24597, r24597, r24598);
        double r24600 = 2;
        double r24601 = r24600 * r24585;
        double r24602 = r24599 / r24601;
        double r24603 = r24592 / r24598;
        double r24604 = r24589 ? r24602 : r24603;
        double r24605 = r24583 ? r24587 : r24604;
        return r24605;
}

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 r24606, r24607, r24608, r24609, r24610, r24611, r24612, r24613, r24614, r24615, r24616, r24617, r24618, r24619;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3472);
        mpfr_init(r24606);
        mpfr_init(r24607);
        mpfr_init(r24608);
        mpfr_init_set_str(r24609, "4", 10, MPFR_RNDN);
        mpfr_init(r24610);
        mpfr_init(r24611);
        mpfr_init(r24612);
        mpfr_init(r24613);
        mpfr_init(r24614);
        mpfr_init(r24615);
        mpfr_init(r24616);
        mpfr_init_set_str(r24617, "2", 10, MPFR_RNDN);
        mpfr_init(r24618);
        mpfr_init(r24619);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r24606, b, MPFR_RNDN);
        mpfr_neg(r24607, r24606, MPFR_RNDN);
        mpfr_mul(r24608, r24606, r24606, MPFR_RNDN);
        ;
        mpfr_set_d(r24610, a, MPFR_RNDN);
        mpfr_mul(r24611, r24609, r24610, MPFR_RNDN);
        mpfr_set_d(r24612, c, MPFR_RNDN);
        mpfr_mul(r24613, r24611, r24612, MPFR_RNDN);
        mpfr_sub(r24614, r24608, r24613, MPFR_RNDN);
        mpfr_sqrt(r24615, r24614, MPFR_RNDN);
        mpfr_add(r24616, r24607, r24615, MPFR_RNDN);
        ;
        mpfr_mul(r24618, r24617, r24610, MPFR_RNDN);
        mpfr_div(r24619, r24616, r24618, MPFR_RNDN);
        return mpfr_get_d(r24619, MPFR_RNDN);
}

static mpfr_t r24620, r24621, r24622, r24623, r24624, r24625, r24626, r24627, r24628, r24629, r24630, r24631, r24632, r24633, r24634, r24635, r24636, r24637, r24638, r24639, r24640, r24641, r24642, r24643, r24644;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r24620);
        mpfr_init_set_str(r24621, "-1.5023809489540946e+154", 10, MPFR_RNDN);
        mpfr_init(r24622);
        mpfr_init_set_str(r24623, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r24624);
        mpfr_init(r24625);
        mpfr_init(r24626);
        mpfr_init_set_str(r24627, "5.761969426420715e-61", 10, MPFR_RNDN);
        mpfr_init(r24628);
        mpfr_init_set_str(r24629, "4", 10, MPFR_RNDN);
        mpfr_init(r24630);
        mpfr_init(r24631);
        mpfr_init(r24632);
        mpfr_init(r24633);
        mpfr_init(r24634);
        mpfr_init(r24635);
        mpfr_init(r24636);
        mpfr_init(r24637);
        mpfr_init(r24638);
        mpfr_init_set_str(r24639, "2", 10, MPFR_RNDN);
        mpfr_init(r24640);
        mpfr_init(r24641);
        mpfr_init(r24642);
        mpfr_init(r24643);
        mpfr_init(r24644);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r24620, b, MPFR_RNDN);
        ;
        mpfr_set_si(r24622, mpfr_cmp(r24620, r24621) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r24624, a, MPFR_RNDN);
        mpfr_div(r24625, r24620, r24624, MPFR_RNDN);
        mpfr_mul(r24626, r24623, r24625, MPFR_RNDN);
        ;
        mpfr_set_si(r24628, mpfr_cmp(r24620, r24627) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r24630, r24629, r24624, MPFR_RNDN);
        mpfr_set_d(r24631, c, MPFR_RNDN);
        mpfr_neg(r24632, r24631, MPFR_RNDN);
        mpfr_mul(r24633, r24620, r24620, MPFR_RNDN);
        mpfr_fma(r24634, r24630, r24632, r24633, MPFR_RNDN);
        mpfr_sqrt(r24635, r24634, MPFR_RNDN);
        mpfr_sqrt(r24636, r24635, MPFR_RNDN);
        mpfr_neg(r24637, r24620, MPFR_RNDN);
        mpfr_fma(r24638, r24636, r24636, r24637, MPFR_RNDN);
        ;
        mpfr_mul(r24640, r24639, r24624, MPFR_RNDN);
        mpfr_div(r24641, r24638, r24640, MPFR_RNDN);
        mpfr_div(r24642, r24631, r24637, MPFR_RNDN);
        if (mpfr_get_si(r24628, MPFR_RNDN)) { mpfr_set(r24643, r24641, MPFR_RNDN); } else { mpfr_set(r24643, r24642, MPFR_RNDN); };
        if (mpfr_get_si(r24622, MPFR_RNDN)) { mpfr_set(r24644, r24626, MPFR_RNDN); } else { mpfr_set(r24644, r24643, MPFR_RNDN); };
        return mpfr_get_d(r24644, MPFR_RNDN);
}

static mpfr_t r24645, r24646, r24647, r24648, r24649, r24650, r24651, r24652, r24653, r24654, r24655, r24656, r24657, r24658, r24659, r24660, r24661, r24662, r24663, r24664, r24665, r24666, r24667, r24668, r24669;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r24645);
        mpfr_init_set_str(r24646, "-1.5023809489540946e+154", 10, MPFR_RNDN);
        mpfr_init(r24647);
        mpfr_init_set_str(r24648, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r24649);
        mpfr_init(r24650);
        mpfr_init(r24651);
        mpfr_init_set_str(r24652, "5.761969426420715e-61", 10, MPFR_RNDN);
        mpfr_init(r24653);
        mpfr_init_set_str(r24654, "4", 10, MPFR_RNDN);
        mpfr_init(r24655);
        mpfr_init(r24656);
        mpfr_init(r24657);
        mpfr_init(r24658);
        mpfr_init(r24659);
        mpfr_init(r24660);
        mpfr_init(r24661);
        mpfr_init(r24662);
        mpfr_init(r24663);
        mpfr_init_set_str(r24664, "2", 10, MPFR_RNDN);
        mpfr_init(r24665);
        mpfr_init(r24666);
        mpfr_init(r24667);
        mpfr_init(r24668);
        mpfr_init(r24669);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r24645, b, MPFR_RNDN);
        ;
        mpfr_set_si(r24647, mpfr_cmp(r24645, r24646) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r24649, a, MPFR_RNDN);
        mpfr_div(r24650, r24645, r24649, MPFR_RNDN);
        mpfr_mul(r24651, r24648, r24650, MPFR_RNDN);
        ;
        mpfr_set_si(r24653, mpfr_cmp(r24645, r24652) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r24655, r24654, r24649, MPFR_RNDN);
        mpfr_set_d(r24656, c, MPFR_RNDN);
        mpfr_neg(r24657, r24656, MPFR_RNDN);
        mpfr_mul(r24658, r24645, r24645, MPFR_RNDN);
        mpfr_fma(r24659, r24655, r24657, r24658, MPFR_RNDN);
        mpfr_sqrt(r24660, r24659, MPFR_RNDN);
        mpfr_sqrt(r24661, r24660, MPFR_RNDN);
        mpfr_neg(r24662, r24645, MPFR_RNDN);
        mpfr_fma(r24663, r24661, r24661, r24662, MPFR_RNDN);
        ;
        mpfr_mul(r24665, r24664, r24649, MPFR_RNDN);
        mpfr_div(r24666, r24663, r24665, MPFR_RNDN);
        mpfr_div(r24667, r24656, r24662, MPFR_RNDN);
        if (mpfr_get_si(r24653, MPFR_RNDN)) { mpfr_set(r24668, r24666, MPFR_RNDN); } else { mpfr_set(r24668, r24667, MPFR_RNDN); };
        if (mpfr_get_si(r24647, MPFR_RNDN)) { mpfr_set(r24669, r24651, MPFR_RNDN); } else { mpfr_set(r24669, r24668, MPFR_RNDN); };
        return mpfr_get_d(r24669, MPFR_RNDN);
}