#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 r26220 = b;
        float r26221 = -r26220;
        float r26222 = r26220 * r26220;
        float r26223 = 3;
        float r26224 = a;
        float r26225 = r26223 * r26224;
        float r26226 = c;
        float r26227 = r26225 * r26226;
        float r26228 = r26222 - r26227;
        float r26229 = sqrt(r26228);
        float r26230 = r26221 + r26229;
        float r26231 = r26230 / r26225;
        return r26231;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r26232 = b;
        double r26233 = -r26232;
        double r26234 = r26232 * r26232;
        double r26235 = 3;
        double r26236 = a;
        double r26237 = r26235 * r26236;
        double r26238 = c;
        double r26239 = r26237 * r26238;
        double r26240 = r26234 - r26239;
        double r26241 = sqrt(r26240);
        double r26242 = r26233 + r26241;
        double r26243 = r26242 / r26237;
        return r26243;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r26244 = b;
        float r26245 = -r26244;
        float r26246 = r26244 * r26244;
        float r26247 = 3;
        float r26248 = a;
        float r26249 = r26247 * r26248;
        float r26250 = c;
        float r26251 = r26249 * r26250;
        float r26252 = r26246 - r26251;
        float r26253 = sqrt(r26252);
        float r26254 = r26245 + r26253;
        float r26255 = r26254 / r26249;
        return r26255;
}

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

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 r26268, r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26268);
        mpfr_init(r26269);
        mpfr_init(r26270);
        mpfr_init_set_str(r26271, "3", 10, MPFR_RNDN);
        mpfr_init(r26272);
        mpfr_init(r26273);
        mpfr_init(r26274);
        mpfr_init(r26275);
        mpfr_init(r26276);
        mpfr_init(r26277);
        mpfr_init(r26278);
        mpfr_init(r26279);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26268, b, MPFR_RNDN);
        mpfr_neg(r26269, r26268, MPFR_RNDN);
        mpfr_mul(r26270, r26268, r26268, MPFR_RNDN);
        ;
        mpfr_set_d(r26272, a, MPFR_RNDN);
        mpfr_mul(r26273, r26271, r26272, MPFR_RNDN);
        mpfr_set_d(r26274, c, MPFR_RNDN);
        mpfr_mul(r26275, r26273, r26274, MPFR_RNDN);
        mpfr_sub(r26276, r26270, r26275, MPFR_RNDN);
        mpfr_sqrt(r26277, r26276, MPFR_RNDN);
        mpfr_add(r26278, r26269, r26277, MPFR_RNDN);
        mpfr_div(r26279, r26278, r26273, MPFR_RNDN);
        return mpfr_get_d(r26279, MPFR_RNDN);
}

static mpfr_t r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26280);
        mpfr_init(r26281);
        mpfr_init(r26282);
        mpfr_init_set_str(r26283, "3", 10, MPFR_RNDN);
        mpfr_init(r26284);
        mpfr_init(r26285);
        mpfr_init(r26286);
        mpfr_init(r26287);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
        mpfr_init(r26291);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26280, b, MPFR_RNDN);
        mpfr_neg(r26281, r26280, MPFR_RNDN);
        mpfr_mul(r26282, r26280, r26280, MPFR_RNDN);
        ;
        mpfr_set_d(r26284, a, MPFR_RNDN);
        mpfr_mul(r26285, r26283, r26284, MPFR_RNDN);
        mpfr_set_d(r26286, c, MPFR_RNDN);
        mpfr_mul(r26287, r26285, r26286, MPFR_RNDN);
        mpfr_sub(r26288, r26282, r26287, MPFR_RNDN);
        mpfr_sqrt(r26289, r26288, MPFR_RNDN);
        mpfr_add(r26290, r26281, r26289, MPFR_RNDN);
        mpfr_div(r26291, r26290, r26285, MPFR_RNDN);
        return mpfr_get_d(r26291, MPFR_RNDN);
}

static mpfr_t r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26292);
        mpfr_init(r26293);
        mpfr_init(r26294);
        mpfr_init_set_str(r26295, "3", 10, MPFR_RNDN);
        mpfr_init(r26296);
        mpfr_init(r26297);
        mpfr_init(r26298);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init(r26301);
        mpfr_init(r26302);
        mpfr_init(r26303);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26292, b, MPFR_RNDN);
        mpfr_neg(r26293, r26292, MPFR_RNDN);
        mpfr_mul(r26294, r26292, r26292, MPFR_RNDN);
        ;
        mpfr_set_d(r26296, a, MPFR_RNDN);
        mpfr_mul(r26297, r26295, r26296, MPFR_RNDN);
        mpfr_set_d(r26298, c, MPFR_RNDN);
        mpfr_mul(r26299, r26297, r26298, MPFR_RNDN);
        mpfr_sub(r26300, r26294, r26299, MPFR_RNDN);
        mpfr_sqrt(r26301, r26300, MPFR_RNDN);
        mpfr_add(r26302, r26293, r26301, MPFR_RNDN);
        mpfr_div(r26303, r26302, r26297, MPFR_RNDN);
        return mpfr_get_d(r26303, MPFR_RNDN);
}