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

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

double f_if(float a, float b, float c, float __attribute__((unused)) d) {
        float r25993 = b;
        float r25994 = -r25993;
        float r25995 = r25993 * r25993;
        float r25996 = 3;
        float r25997 = a;
        float r25998 = r25996 * r25997;
        float r25999 = c;
        float r26000 = r25998 * r25999;
        float r26001 = r25995 - r26000;
        float r26002 = sqrt(r26001);
        float r26003 = r25994 + r26002;
        float r26004 = r26003 / r25998;
        return r26004;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r26005 = b;
        double r26006 = -r26005;
        double r26007 = r26005 * r26005;
        double r26008 = 3;
        double r26009 = a;
        double r26010 = r26008 * r26009;
        double r26011 = c;
        double r26012 = r26010 * r26011;
        double r26013 = r26007 - r26012;
        double r26014 = sqrt(r26013);
        double r26015 = r26006 + r26014;
        double r26016 = r26015 / r26010;
        return r26016;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r26017 = b;
        float r26018 = -r26017;
        float r26019 = r26017 * r26017;
        float r26020 = 3;
        float r26021 = a;
        float r26022 = r26020 * r26021;
        float r26023 = c;
        float r26024 = r26022 * r26023;
        float r26025 = r26019 - r26024;
        float r26026 = sqrt(r26025);
        float r26027 = r26018 + r26026;
        float r26028 = r26027 / r26022;
        return r26028;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r26029 = b;
        double r26030 = -r26029;
        double r26031 = r26029 * r26029;
        double r26032 = 3;
        double r26033 = a;
        double r26034 = r26032 * r26033;
        double r26035 = c;
        double r26036 = r26034 * r26035;
        double r26037 = r26031 - r26036;
        double r26038 = sqrt(r26037);
        double r26039 = r26030 + r26038;
        double r26040 = r26039 / r26034;
        return r26040;
}

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 r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26041);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init_set_str(r26044, "3", 10, MPFR_RNDN);
        mpfr_init(r26045);
        mpfr_init(r26046);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init(r26049);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26041, b, MPFR_RNDN);
        mpfr_neg(r26042, r26041, MPFR_RNDN);
        mpfr_mul(r26043, r26041, r26041, MPFR_RNDN);
        ;
        mpfr_set_d(r26045, a, MPFR_RNDN);
        mpfr_mul(r26046, r26044, r26045, MPFR_RNDN);
        mpfr_set_d(r26047, c, MPFR_RNDN);
        mpfr_mul(r26048, r26046, r26047, MPFR_RNDN);
        mpfr_sub(r26049, r26043, r26048, MPFR_RNDN);
        mpfr_sqrt(r26050, r26049, MPFR_RNDN);
        mpfr_add(r26051, r26042, r26050, MPFR_RNDN);
        mpfr_div(r26052, r26051, r26046, MPFR_RNDN);
        return mpfr_get_d(r26052, MPFR_RNDN);
}

static mpfr_t r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
        mpfr_init_set_str(r26056, "3", 10, MPFR_RNDN);
        mpfr_init(r26057);
        mpfr_init(r26058);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init(r26061);
        mpfr_init(r26062);
        mpfr_init(r26063);
        mpfr_init(r26064);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26053, b, MPFR_RNDN);
        mpfr_neg(r26054, r26053, MPFR_RNDN);
        mpfr_mul(r26055, r26053, r26053, MPFR_RNDN);
        ;
        mpfr_set_d(r26057, a, MPFR_RNDN);
        mpfr_mul(r26058, r26056, r26057, MPFR_RNDN);
        mpfr_set_d(r26059, c, MPFR_RNDN);
        mpfr_mul(r26060, r26058, r26059, MPFR_RNDN);
        mpfr_sub(r26061, r26055, r26060, MPFR_RNDN);
        mpfr_sqrt(r26062, r26061, MPFR_RNDN);
        mpfr_add(r26063, r26054, r26062, MPFR_RNDN);
        mpfr_div(r26064, r26063, r26058, MPFR_RNDN);
        return mpfr_get_d(r26064, MPFR_RNDN);
}

static mpfr_t r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26065);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init_set_str(r26068, "3", 10, MPFR_RNDN);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
        mpfr_init(r26073);
        mpfr_init(r26074);
        mpfr_init(r26075);
        mpfr_init(r26076);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26065, b, MPFR_RNDN);
        mpfr_neg(r26066, r26065, MPFR_RNDN);
        mpfr_mul(r26067, r26065, r26065, MPFR_RNDN);
        ;
        mpfr_set_d(r26069, a, MPFR_RNDN);
        mpfr_mul(r26070, r26068, r26069, MPFR_RNDN);
        mpfr_set_d(r26071, c, MPFR_RNDN);
        mpfr_mul(r26072, r26070, r26071, MPFR_RNDN);
        mpfr_sub(r26073, r26067, r26072, MPFR_RNDN);
        mpfr_sqrt(r26074, r26073, MPFR_RNDN);
        mpfr_add(r26075, r26066, r26074, MPFR_RNDN);
        mpfr_div(r26076, r26075, r26070, MPFR_RNDN);
        return mpfr_get_d(r26076, MPFR_RNDN);
}