#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 r26246 = b;
        float r26247 = -r26246;
        float r26248 = r26246 * r26246;
        float r26249 = 3;
        float r26250 = a;
        float r26251 = r26249 * r26250;
        float r26252 = c;
        float r26253 = r26251 * r26252;
        float r26254 = r26248 - r26253;
        float r26255 = sqrt(r26254);
        float r26256 = r26247 + r26255;
        float r26257 = r26256 / r26251;
        return r26257;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r26258 = b;
        double r26259 = -r26258;
        double r26260 = r26258 * r26258;
        double r26261 = 3;
        double r26262 = a;
        double r26263 = r26261 * r26262;
        double r26264 = c;
        double r26265 = r26263 * r26264;
        double r26266 = r26260 - r26265;
        double r26267 = sqrt(r26266);
        double r26268 = r26259 + r26267;
        double r26269 = r26268 / r26263;
        return r26269;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r26270 = c;
        float r26271 = 3;
        float r26272 = a;
        float r26273 = r26271 * r26272;
        float r26274 = r26270 * r26273;
        float r26275 = -r26270;
        float r26276 = b;
        float r26277 = r26276 * r26276;
        float r26278 = fma(r26275, r26273, r26277);
        float r26279 = sqrt(r26278);
        float r26280 = r26279 + r26276;
        float r26281 = r26274 / r26280;
        float r26282 = -r26273;
        float r26283 = r26281 / r26282;
        return r26283;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r26284 = c;
        double r26285 = 3;
        double r26286 = a;
        double r26287 = r26285 * r26286;
        double r26288 = r26284 * r26287;
        double r26289 = -r26284;
        double r26290 = b;
        double r26291 = r26290 * r26290;
        double r26292 = fma(r26289, r26287, r26291);
        double r26293 = sqrt(r26292);
        double r26294 = r26293 + r26290;
        double r26295 = r26288 / r26294;
        double r26296 = -r26287;
        double r26297 = r26295 / r26296;
        return r26297;
}

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 r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26298);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init_set_str(r26301, "3", 10, MPFR_RNDN);
        mpfr_init(r26302);
        mpfr_init(r26303);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26298, b, MPFR_RNDN);
        mpfr_neg(r26299, r26298, MPFR_RNDN);
        mpfr_mul(r26300, r26298, r26298, MPFR_RNDN);
        ;
        mpfr_set_d(r26302, a, MPFR_RNDN);
        mpfr_mul(r26303, r26301, r26302, MPFR_RNDN);
        mpfr_set_d(r26304, c, MPFR_RNDN);
        mpfr_mul(r26305, r26303, r26304, MPFR_RNDN);
        mpfr_sub(r26306, r26300, r26305, MPFR_RNDN);
        mpfr_sqrt(r26307, r26306, MPFR_RNDN);
        mpfr_add(r26308, r26299, r26307, MPFR_RNDN);
        mpfr_div(r26309, r26308, r26303, MPFR_RNDN);
        return mpfr_get_d(r26309, MPFR_RNDN);
}

static mpfr_t r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26310);
        mpfr_init_set_str(r26311, "3", 10, MPFR_RNDN);
        mpfr_init(r26312);
        mpfr_init(r26313);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init(r26317);
        mpfr_init(r26318);
        mpfr_init(r26319);
        mpfr_init(r26320);
        mpfr_init(r26321);
        mpfr_init(r26322);
        mpfr_init(r26323);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26310, c, MPFR_RNDN);
        ;
        mpfr_set_d(r26312, a, MPFR_RNDN);
        mpfr_mul(r26313, r26311, r26312, MPFR_RNDN);
        mpfr_mul(r26314, r26310, r26313, MPFR_RNDN);
        mpfr_neg(r26315, r26310, MPFR_RNDN);
        mpfr_set_d(r26316, b, MPFR_RNDN);
        mpfr_mul(r26317, r26316, r26316, MPFR_RNDN);
        mpfr_fma(r26318, r26315, r26313, r26317, MPFR_RNDN);
        mpfr_sqrt(r26319, r26318, MPFR_RNDN);
        mpfr_add(r26320, r26319, r26316, MPFR_RNDN);
        mpfr_div(r26321, r26314, r26320, MPFR_RNDN);
        mpfr_neg(r26322, r26313, MPFR_RNDN);
        mpfr_div(r26323, r26321, r26322, MPFR_RNDN);
        return mpfr_get_d(r26323, MPFR_RNDN);
}

static mpfr_t r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26324);
        mpfr_init_set_str(r26325, "3", 10, MPFR_RNDN);
        mpfr_init(r26326);
        mpfr_init(r26327);
        mpfr_init(r26328);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26324, c, MPFR_RNDN);
        ;
        mpfr_set_d(r26326, a, MPFR_RNDN);
        mpfr_mul(r26327, r26325, r26326, MPFR_RNDN);
        mpfr_mul(r26328, r26324, r26327, MPFR_RNDN);
        mpfr_neg(r26329, r26324, MPFR_RNDN);
        mpfr_set_d(r26330, b, MPFR_RNDN);
        mpfr_mul(r26331, r26330, r26330, MPFR_RNDN);
        mpfr_fma(r26332, r26329, r26327, r26331, MPFR_RNDN);
        mpfr_sqrt(r26333, r26332, MPFR_RNDN);
        mpfr_add(r26334, r26333, r26330, MPFR_RNDN);
        mpfr_div(r26335, r26328, r26334, MPFR_RNDN);
        mpfr_neg(r26336, r26327, MPFR_RNDN);
        mpfr_div(r26337, r26335, r26336, MPFR_RNDN);
        return mpfr_get_d(r26337, MPFR_RNDN);
}