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

char *name = "Cubic critical, medium range";

double f_if(float a, float b, float c, float __attribute__((unused)) d) {
        float r23535 = b;
        float r23536 = -r23535;
        float r23537 = r23535 * r23535;
        float r23538 = 3;
        float r23539 = a;
        float r23540 = r23538 * r23539;
        float r23541 = c;
        float r23542 = r23540 * r23541;
        float r23543 = r23537 - r23542;
        float r23544 = sqrt(r23543);
        float r23545 = r23536 + r23544;
        float r23546 = r23545 / r23540;
        return r23546;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r23547 = b;
        double r23548 = -r23547;
        double r23549 = r23547 * r23547;
        double r23550 = 3;
        double r23551 = a;
        double r23552 = r23550 * r23551;
        double r23553 = c;
        double r23554 = r23552 * r23553;
        double r23555 = r23549 - r23554;
        double r23556 = sqrt(r23555);
        double r23557 = r23548 + r23556;
        double r23558 = r23557 / r23552;
        return r23558;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r23559 = b;
        float r23560 = r23559 * r23559;
        float r23561 = a;
        float r23562 = c;
        float r23563 = 3;
        float r23564 = r23562 * r23563;
        float r23565 = r23561 * r23564;
        float r23566 = r23560 - r23565;
        float r23567 = r23565 + r23560;
        float r23568 = r23566 * r23567;
        float r23569 = sqrt(r23568);
        float r23570 = sqrt(r23567);
        float r23571 = r23569 / r23570;
        float r23572 = r23571 - r23559;
        float r23573 = r23563 * r23561;
        float r23574 = r23572 / r23573;
        float r23575 = -1.850810227248647e-21;
        bool r23576 = r23574 <= r23575;
        float r23577 = -r23559;
        float r23578 = r23577 - r23559;
        float r23579 = r23578 / r23573;
        float r23580 = r23576 ? r23574 : r23579;
        return r23580;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r23581 = b;
        double r23582 = r23581 * r23581;
        double r23583 = a;
        double r23584 = c;
        double r23585 = 3;
        double r23586 = r23584 * r23585;
        double r23587 = r23583 * r23586;
        double r23588 = r23582 - r23587;
        double r23589 = r23587 + r23582;
        double r23590 = r23588 * r23589;
        double r23591 = sqrt(r23590);
        double r23592 = sqrt(r23589);
        double r23593 = r23591 / r23592;
        double r23594 = r23593 - r23581;
        double r23595 = r23585 * r23583;
        double r23596 = r23594 / r23595;
        double r23597 = -1.850810227248647e-21;
        bool r23598 = r23596 <= r23597;
        double r23599 = -r23581;
        double r23600 = r23599 - r23581;
        double r23601 = r23600 / r23595;
        double r23602 = r23598 ? r23596 : r23601;
        return r23602;
}

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 r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23603);
        mpfr_init(r23604);
        mpfr_init(r23605);
        mpfr_init_set_str(r23606, "3", 10, MPFR_RNDN);
        mpfr_init(r23607);
        mpfr_init(r23608);
        mpfr_init(r23609);
        mpfr_init(r23610);
        mpfr_init(r23611);
        mpfr_init(r23612);
        mpfr_init(r23613);
        mpfr_init(r23614);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r23603, b, MPFR_RNDN);
        mpfr_neg(r23604, r23603, MPFR_RNDN);
        mpfr_mul(r23605, r23603, r23603, MPFR_RNDN);
        ;
        mpfr_set_d(r23607, a, MPFR_RNDN);
        mpfr_mul(r23608, r23606, r23607, MPFR_RNDN);
        mpfr_set_d(r23609, c, MPFR_RNDN);
        mpfr_mul(r23610, r23608, r23609, MPFR_RNDN);
        mpfr_sub(r23611, r23605, r23610, MPFR_RNDN);
        mpfr_sqrt(r23612, r23611, MPFR_RNDN);
        mpfr_add(r23613, r23604, r23612, MPFR_RNDN);
        mpfr_div(r23614, r23613, r23608, MPFR_RNDN);
        return mpfr_get_d(r23614, MPFR_RNDN);
}

static mpfr_t r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23615);
        mpfr_init(r23616);
        mpfr_init(r23617);
        mpfr_init(r23618);
        mpfr_init_set_str(r23619, "3", 10, MPFR_RNDN);
        mpfr_init(r23620);
        mpfr_init(r23621);
        mpfr_init(r23622);
        mpfr_init(r23623);
        mpfr_init(r23624);
        mpfr_init(r23625);
        mpfr_init(r23626);
        mpfr_init(r23627);
        mpfr_init(r23628);
        mpfr_init(r23629);
        mpfr_init(r23630);
        mpfr_init_set_str(r23631, "-1.850810227248647e-21", 10, MPFR_RNDN);
        mpfr_init(r23632);
        mpfr_init(r23633);
        mpfr_init(r23634);
        mpfr_init(r23635);
        mpfr_init(r23636);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r23615, b, MPFR_RNDN);
        mpfr_mul(r23616, r23615, r23615, MPFR_RNDN);
        mpfr_set_d(r23617, a, MPFR_RNDN);
        mpfr_set_d(r23618, c, MPFR_RNDN);
        ;
        mpfr_mul(r23620, r23618, r23619, MPFR_RNDN);
        mpfr_mul(r23621, r23617, r23620, MPFR_RNDN);
        mpfr_sub(r23622, r23616, r23621, MPFR_RNDN);
        mpfr_add(r23623, r23621, r23616, MPFR_RNDN);
        mpfr_mul(r23624, r23622, r23623, MPFR_RNDN);
        mpfr_sqrt(r23625, r23624, MPFR_RNDN);
        mpfr_sqrt(r23626, r23623, MPFR_RNDN);
        mpfr_div(r23627, r23625, r23626, MPFR_RNDN);
        mpfr_sub(r23628, r23627, r23615, MPFR_RNDN);
        mpfr_mul(r23629, r23619, r23617, MPFR_RNDN);
        mpfr_div(r23630, r23628, r23629, MPFR_RNDN);
        ;
        mpfr_set_si(r23632, mpfr_cmp(r23630, r23631) <= 0, MPFR_RNDN);
        mpfr_neg(r23633, r23615, MPFR_RNDN);
        mpfr_sub(r23634, r23633, r23615, MPFR_RNDN);
        mpfr_div(r23635, r23634, r23629, MPFR_RNDN);
        if (mpfr_get_si(r23632, MPFR_RNDN)) { mpfr_set(r23636, r23630, MPFR_RNDN); } else { mpfr_set(r23636, r23635, MPFR_RNDN); };
        return mpfr_get_d(r23636, MPFR_RNDN);
}

static mpfr_t r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648, r23649, r23650, r23651, r23652, r23653, r23654, r23655, r23656, r23657, r23658;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23637);
        mpfr_init(r23638);
        mpfr_init(r23639);
        mpfr_init(r23640);
        mpfr_init_set_str(r23641, "3", 10, MPFR_RNDN);
        mpfr_init(r23642);
        mpfr_init(r23643);
        mpfr_init(r23644);
        mpfr_init(r23645);
        mpfr_init(r23646);
        mpfr_init(r23647);
        mpfr_init(r23648);
        mpfr_init(r23649);
        mpfr_init(r23650);
        mpfr_init(r23651);
        mpfr_init(r23652);
        mpfr_init_set_str(r23653, "-1.850810227248647e-21", 10, MPFR_RNDN);
        mpfr_init(r23654);
        mpfr_init(r23655);
        mpfr_init(r23656);
        mpfr_init(r23657);
        mpfr_init(r23658);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r23637, b, MPFR_RNDN);
        mpfr_mul(r23638, r23637, r23637, MPFR_RNDN);
        mpfr_set_d(r23639, a, MPFR_RNDN);
        mpfr_set_d(r23640, c, MPFR_RNDN);
        ;
        mpfr_mul(r23642, r23640, r23641, MPFR_RNDN);
        mpfr_mul(r23643, r23639, r23642, MPFR_RNDN);
        mpfr_sub(r23644, r23638, r23643, MPFR_RNDN);
        mpfr_add(r23645, r23643, r23638, MPFR_RNDN);
        mpfr_mul(r23646, r23644, r23645, MPFR_RNDN);
        mpfr_sqrt(r23647, r23646, MPFR_RNDN);
        mpfr_sqrt(r23648, r23645, MPFR_RNDN);
        mpfr_div(r23649, r23647, r23648, MPFR_RNDN);
        mpfr_sub(r23650, r23649, r23637, MPFR_RNDN);
        mpfr_mul(r23651, r23641, r23639, MPFR_RNDN);
        mpfr_div(r23652, r23650, r23651, MPFR_RNDN);
        ;
        mpfr_set_si(r23654, mpfr_cmp(r23652, r23653) <= 0, MPFR_RNDN);
        mpfr_neg(r23655, r23637, MPFR_RNDN);
        mpfr_sub(r23656, r23655, r23637, MPFR_RNDN);
        mpfr_div(r23657, r23656, r23651, MPFR_RNDN);
        if (mpfr_get_si(r23654, MPFR_RNDN)) { mpfr_set(r23658, r23652, MPFR_RNDN); } else { mpfr_set(r23658, r23657, MPFR_RNDN); };
        return mpfr_get_d(r23658, MPFR_RNDN);
}