#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 r25928 = b;
        float r25929 = -r25928;
        float r25930 = r25928 * r25928;
        float r25931 = 3;
        float r25932 = a;
        float r25933 = r25931 * r25932;
        float r25934 = c;
        float r25935 = r25933 * r25934;
        float r25936 = r25930 - r25935;
        float r25937 = sqrt(r25936);
        float r25938 = r25929 + r25937;
        float r25939 = r25938 / r25933;
        return r25939;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r25940 = b;
        double r25941 = -r25940;
        double r25942 = r25940 * r25940;
        double r25943 = 3;
        double r25944 = a;
        double r25945 = r25943 * r25944;
        double r25946 = c;
        double r25947 = r25945 * r25946;
        double r25948 = r25942 - r25947;
        double r25949 = sqrt(r25948);
        double r25950 = r25941 + r25949;
        double r25951 = r25950 / r25945;
        return r25951;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r25952 = c;
        float r25953 = a;
        float r25954 = r25952 * r25953;
        float r25955 = -r25954;
        float r25956 = -r25952;
        float r25957 = 3;
        float r25958 = r25957 * r25953;
        float r25959 = b;
        float r25960 = r25959 * r25959;
        float r25961 = fma(r25956, r25958, r25960);
        float r25962 = sqrt(r25961);
        float r25963 = r25962 + r25959;
        float r25964 = r25955 / r25963;
        float r25965 = r25964 / r25953;
        return r25965;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r25966 = c;
        double r25967 = a;
        double r25968 = r25966 * r25967;
        double r25969 = -r25968;
        double r25970 = -r25966;
        double r25971 = 3;
        double r25972 = r25971 * r25967;
        double r25973 = b;
        double r25974 = r25973 * r25973;
        double r25975 = fma(r25970, r25972, r25974);
        double r25976 = sqrt(r25975);
        double r25977 = r25976 + r25973;
        double r25978 = r25969 / r25977;
        double r25979 = r25978 / r25967;
        return r25979;
}

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 r25980, r25981, r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r25980);
        mpfr_init(r25981);
        mpfr_init(r25982);
        mpfr_init_set_str(r25983, "3", 10, MPFR_RNDN);
        mpfr_init(r25984);
        mpfr_init(r25985);
        mpfr_init(r25986);
        mpfr_init(r25987);
        mpfr_init(r25988);
        mpfr_init(r25989);
        mpfr_init(r25990);
        mpfr_init(r25991);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r25980, b, MPFR_RNDN);
        mpfr_neg(r25981, r25980, MPFR_RNDN);
        mpfr_mul(r25982, r25980, r25980, MPFR_RNDN);
        ;
        mpfr_set_d(r25984, a, MPFR_RNDN);
        mpfr_mul(r25985, r25983, r25984, MPFR_RNDN);
        mpfr_set_d(r25986, c, MPFR_RNDN);
        mpfr_mul(r25987, r25985, r25986, MPFR_RNDN);
        mpfr_sub(r25988, r25982, r25987, MPFR_RNDN);
        mpfr_sqrt(r25989, r25988, MPFR_RNDN);
        mpfr_add(r25990, r25981, r25989, MPFR_RNDN);
        mpfr_div(r25991, r25990, r25985, MPFR_RNDN);
        return mpfr_get_d(r25991, MPFR_RNDN);
}

static mpfr_t r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000, r26001, r26002, r26003, r26004, r26005;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r25992);
        mpfr_init(r25993);
        mpfr_init(r25994);
        mpfr_init(r25995);
        mpfr_init(r25996);
        mpfr_init_set_str(r25997, "3", 10, MPFR_RNDN);
        mpfr_init(r25998);
        mpfr_init(r25999);
        mpfr_init(r26000);
        mpfr_init(r26001);
        mpfr_init(r26002);
        mpfr_init(r26003);
        mpfr_init(r26004);
        mpfr_init(r26005);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r25992, c, MPFR_RNDN);
        mpfr_set_d(r25993, a, MPFR_RNDN);
        mpfr_mul(r25994, r25992, r25993, MPFR_RNDN);
        mpfr_neg(r25995, r25994, MPFR_RNDN);
        mpfr_neg(r25996, r25992, MPFR_RNDN);
        ;
        mpfr_mul(r25998, r25997, r25993, MPFR_RNDN);
        mpfr_set_d(r25999, b, MPFR_RNDN);
        mpfr_mul(r26000, r25999, r25999, MPFR_RNDN);
        mpfr_fma(r26001, r25996, r25998, r26000, MPFR_RNDN);
        mpfr_sqrt(r26002, r26001, MPFR_RNDN);
        mpfr_add(r26003, r26002, r25999, MPFR_RNDN);
        mpfr_div(r26004, r25995, r26003, MPFR_RNDN);
        mpfr_div(r26005, r26004, r25993, MPFR_RNDN);
        return mpfr_get_d(r26005, MPFR_RNDN);
}

static mpfr_t r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017, r26018, r26019;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26006);
        mpfr_init(r26007);
        mpfr_init(r26008);
        mpfr_init(r26009);
        mpfr_init(r26010);
        mpfr_init_set_str(r26011, "3", 10, MPFR_RNDN);
        mpfr_init(r26012);
        mpfr_init(r26013);
        mpfr_init(r26014);
        mpfr_init(r26015);
        mpfr_init(r26016);
        mpfr_init(r26017);
        mpfr_init(r26018);
        mpfr_init(r26019);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26006, c, MPFR_RNDN);
        mpfr_set_d(r26007, a, MPFR_RNDN);
        mpfr_mul(r26008, r26006, r26007, MPFR_RNDN);
        mpfr_neg(r26009, r26008, MPFR_RNDN);
        mpfr_neg(r26010, r26006, MPFR_RNDN);
        ;
        mpfr_mul(r26012, r26011, r26007, MPFR_RNDN);
        mpfr_set_d(r26013, b, MPFR_RNDN);
        mpfr_mul(r26014, r26013, r26013, MPFR_RNDN);
        mpfr_fma(r26015, r26010, r26012, r26014, MPFR_RNDN);
        mpfr_sqrt(r26016, r26015, MPFR_RNDN);
        mpfr_add(r26017, r26016, r26013, MPFR_RNDN);
        mpfr_div(r26018, r26009, r26017, MPFR_RNDN);
        mpfr_div(r26019, r26018, r26007, MPFR_RNDN);
        return mpfr_get_d(r26019, MPFR_RNDN);
}