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

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

double f_if(float a, float b, float c, float __attribute__((unused)) d) {
        float r9602 = b;
        float r9603 = -r9602;
        float r9604 = r9602 * r9602;
        float r9605 = 3.0f;
        float r9606 = a;
        float r9607 = r9605 * r9606;
        float r9608 = c;
        float r9609 = r9607 * r9608;
        float r9610 = r9604 - r9609;
        float r9611 = sqrt(r9610);
        float r9612 = r9603 + r9611;
        float r9613 = r9612 / r9607;
        return r9613;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r9614 = b;
        double r9615 = -r9614;
        double r9616 = r9614 * r9614;
        double r9617 = 3.0;
        double r9618 = a;
        double r9619 = r9617 * r9618;
        double r9620 = c;
        double r9621 = r9619 * r9620;
        double r9622 = r9616 - r9621;
        double r9623 = sqrt(r9622);
        double r9624 = r9615 + r9623;
        double r9625 = r9624 / r9619;
        return r9625;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r9626 = c;
        float r9627 = -r9626;
        float r9628 = b;
        float r9629 = a;
        float r9630 = 3.0f;
        float r9631 = r9629 * r9630;
        float r9632 = r9628 * r9628;
        float r9633 = fma(r9631, r9627, r9632);
        float r9634 = sqrt(r9633);
        float r9635 = r9628 + r9634;
        float r9636 = r9627 / r9635;
        return r9636;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r9637 = c;
        double r9638 = -r9637;
        double r9639 = b;
        double r9640 = a;
        double r9641 = 3.0;
        double r9642 = r9640 * r9641;
        double r9643 = r9639 * r9639;
        double r9644 = fma(r9642, r9638, r9643);
        double r9645 = sqrt(r9644);
        double r9646 = r9639 + r9645;
        double r9647 = r9638 / r9646;
        return r9647;
}

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 r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658, r9659;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init_set_str(r9651, "3", 10, MPFR_RNDN);
        mpfr_init(r9652);
        mpfr_init(r9653);
        mpfr_init(r9654);
        mpfr_init(r9655);
        mpfr_init(r9656);
        mpfr_init(r9657);
        mpfr_init(r9658);
        mpfr_init(r9659);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r9648, b, MPFR_RNDN);
        mpfr_neg(r9649, r9648, MPFR_RNDN);
        mpfr_mul(r9650, r9648, r9648, MPFR_RNDN);
        ;
        mpfr_set_d(r9652, a, MPFR_RNDN);
        mpfr_mul(r9653, r9651, r9652, MPFR_RNDN);
        mpfr_set_d(r9654, c, MPFR_RNDN);
        mpfr_mul(r9655, r9653, r9654, MPFR_RNDN);
        mpfr_sub(r9656, r9650, r9655, MPFR_RNDN);
        mpfr_sqrt(r9657, r9656, MPFR_RNDN);
        mpfr_add(r9658, r9649, r9657, MPFR_RNDN);
        mpfr_div(r9659, r9658, r9653, MPFR_RNDN);
        return mpfr_get_d(r9659, MPFR_RNDN);
}

static mpfr_t r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init(r9660);
        mpfr_init(r9661);
        mpfr_init(r9662);
        mpfr_init(r9663);
        mpfr_init_set_str(r9664, "3", 10, MPFR_RNDN);
        mpfr_init(r9665);
        mpfr_init(r9666);
        mpfr_init(r9667);
        mpfr_init(r9668);
        mpfr_init(r9669);
        mpfr_init(r9670);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r9660, c, MPFR_RNDN);
        mpfr_neg(r9661, r9660, MPFR_RNDN);
        mpfr_set_d(r9662, b, MPFR_RNDN);
        mpfr_set_d(r9663, a, MPFR_RNDN);
        ;
        mpfr_mul(r9665, r9663, r9664, MPFR_RNDN);
        mpfr_mul(r9666, r9662, r9662, MPFR_RNDN);
        mpfr_fma(r9667, r9665, r9661, r9666, MPFR_RNDN);
        mpfr_sqrt(r9668, r9667, MPFR_RNDN);
        mpfr_add(r9669, r9662, r9668, MPFR_RNDN);
        mpfr_div(r9670, r9661, r9669, MPFR_RNDN);
        return mpfr_get_d(r9670, MPFR_RNDN);
}

static mpfr_t r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680, r9681;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init(r9671);
        mpfr_init(r9672);
        mpfr_init(r9673);
        mpfr_init(r9674);
        mpfr_init_set_str(r9675, "3", 10, MPFR_RNDN);
        mpfr_init(r9676);
        mpfr_init(r9677);
        mpfr_init(r9678);
        mpfr_init(r9679);
        mpfr_init(r9680);
        mpfr_init(r9681);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r9671, c, MPFR_RNDN);
        mpfr_neg(r9672, r9671, MPFR_RNDN);
        mpfr_set_d(r9673, b, MPFR_RNDN);
        mpfr_set_d(r9674, a, MPFR_RNDN);
        ;
        mpfr_mul(r9676, r9674, r9675, MPFR_RNDN);
        mpfr_mul(r9677, r9673, r9673, MPFR_RNDN);
        mpfr_fma(r9678, r9676, r9672, r9677, MPFR_RNDN);
        mpfr_sqrt(r9679, r9678, MPFR_RNDN);
        mpfr_add(r9680, r9673, r9679, MPFR_RNDN);
        mpfr_div(r9681, r9672, r9680, MPFR_RNDN);
        return mpfr_get_d(r9681, MPFR_RNDN);
}