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

char *name = "Cubic critical";

double f_if(float a, float b, float c, float __attribute__((unused)) d) {
        float r26206 = b;
        float r26207 = -r26206;
        float r26208 = r26206 * r26206;
        float r26209 = 3;
        float r26210 = a;
        float r26211 = r26209 * r26210;
        float r26212 = c;
        float r26213 = r26211 * r26212;
        float r26214 = r26208 - r26213;
        float r26215 = sqrt(r26214);
        float r26216 = r26207 + r26215;
        float r26217 = r26216 / r26211;
        return r26217;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r26218 = b;
        double r26219 = -r26218;
        double r26220 = r26218 * r26218;
        double r26221 = 3;
        double r26222 = a;
        double r26223 = r26221 * r26222;
        double r26224 = c;
        double r26225 = r26223 * r26224;
        double r26226 = r26220 - r26225;
        double r26227 = sqrt(r26226);
        double r26228 = r26219 + r26227;
        double r26229 = r26228 / r26223;
        return r26229;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r26230 = b;
        float r26231 = -3.273058773164801e+121;
        bool r26232 = r26230 <= r26231;
        float r26233 = 3/2;
        float r26234 = a;
        float r26235 = r26233 * r26234;
        float r26236 = c;
        float r26237 = r26230 / r26236;
        float r26238 = r26235 / r26237;
        float r26239 = r26230 + r26230;
        float r26240 = r26238 - r26239;
        float r26241 = 3;
        float r26242 = r26241 * r26234;
        float r26243 = r26240 / r26242;
        float r26244 = 1.0696681173689353e-137;
        bool r26245 = r26230 <= r26244;
        float r26246 = -r26236;
        float r26247 = r26230 * r26230;
        float r26248 = fma(r26246, r26242, r26247);
        float r26249 = sqrt(r26248);
        float r26250 = r26249 - r26230;
        float r26251 = r26250 / r26241;
        float r26252 = r26251 / r26234;
        float r26253 = r26236 / r26230;
        float r26254 = -3/2;
        float r26255 = r26254 / r26241;
        float r26256 = r26253 * r26255;
        float r26257 = r26245 ? r26252 : r26256;
        float r26258 = r26232 ? r26243 : r26257;
        return r26258;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r26259 = b;
        double r26260 = -3.273058773164801e+121;
        bool r26261 = r26259 <= r26260;
        double r26262 = 3/2;
        double r26263 = a;
        double r26264 = r26262 * r26263;
        double r26265 = c;
        double r26266 = r26259 / r26265;
        double r26267 = r26264 / r26266;
        double r26268 = r26259 + r26259;
        double r26269 = r26267 - r26268;
        double r26270 = 3;
        double r26271 = r26270 * r26263;
        double r26272 = r26269 / r26271;
        double r26273 = 1.0696681173689353e-137;
        bool r26274 = r26259 <= r26273;
        double r26275 = -r26265;
        double r26276 = r26259 * r26259;
        double r26277 = fma(r26275, r26271, r26276);
        double r26278 = sqrt(r26277);
        double r26279 = r26278 - r26259;
        double r26280 = r26279 / r26270;
        double r26281 = r26280 / r26263;
        double r26282 = r26265 / r26259;
        double r26283 = -3/2;
        double r26284 = r26283 / r26270;
        double r26285 = r26282 * r26284;
        double r26286 = r26274 ? r26281 : r26285;
        double r26287 = r26261 ? r26272 : r26286;
        return r26287;
}

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 r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297, r26298, r26299;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
        mpfr_init_set_str(r26291, "3", 10, MPFR_RNDN);
        mpfr_init(r26292);
        mpfr_init(r26293);
        mpfr_init(r26294);
        mpfr_init(r26295);
        mpfr_init(r26296);
        mpfr_init(r26297);
        mpfr_init(r26298);
        mpfr_init(r26299);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26288, b, MPFR_RNDN);
        mpfr_neg(r26289, r26288, MPFR_RNDN);
        mpfr_mul(r26290, r26288, r26288, MPFR_RNDN);
        ;
        mpfr_set_d(r26292, a, MPFR_RNDN);
        mpfr_mul(r26293, r26291, r26292, MPFR_RNDN);
        mpfr_set_d(r26294, c, MPFR_RNDN);
        mpfr_mul(r26295, r26293, r26294, MPFR_RNDN);
        mpfr_sub(r26296, r26290, r26295, MPFR_RNDN);
        mpfr_sqrt(r26297, r26296, MPFR_RNDN);
        mpfr_add(r26298, r26289, r26297, MPFR_RNDN);
        mpfr_div(r26299, r26298, r26293, MPFR_RNDN);
        return mpfr_get_d(r26299, MPFR_RNDN);
}

static mpfr_t r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26300);
        mpfr_init_set_str(r26301, "-3.273058773164801e+121", 10, MPFR_RNDN);
        mpfr_init(r26302);
        mpfr_init_set_str(r26303, "3/2", 10, MPFR_RNDN);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
        mpfr_init(r26310);
        mpfr_init_set_str(r26311, "3", 10, MPFR_RNDN);
        mpfr_init(r26312);
        mpfr_init(r26313);
        mpfr_init_set_str(r26314, "1.0696681173689353e-137", 10, MPFR_RNDN);
        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);
        mpfr_init_set_str(r26324, "-3/2", 10, MPFR_RNDN);
        mpfr_init(r26325);
        mpfr_init(r26326);
        mpfr_init(r26327);
        mpfr_init(r26328);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26300, b, MPFR_RNDN);
        ;
        mpfr_set_si(r26302, mpfr_cmp(r26300, r26301) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r26304, a, MPFR_RNDN);
        mpfr_mul(r26305, r26303, r26304, MPFR_RNDN);
        mpfr_set_d(r26306, c, MPFR_RNDN);
        mpfr_div(r26307, r26300, r26306, MPFR_RNDN);
        mpfr_div(r26308, r26305, r26307, MPFR_RNDN);
        mpfr_add(r26309, r26300, r26300, MPFR_RNDN);
        mpfr_sub(r26310, r26308, r26309, MPFR_RNDN);
        ;
        mpfr_mul(r26312, r26311, r26304, MPFR_RNDN);
        mpfr_div(r26313, r26310, r26312, MPFR_RNDN);
        ;
        mpfr_set_si(r26315, mpfr_cmp(r26300, r26314) <= 0, MPFR_RNDN);
        mpfr_neg(r26316, r26306, MPFR_RNDN);
        mpfr_mul(r26317, r26300, r26300, MPFR_RNDN);
        mpfr_fma(r26318, r26316, r26312, r26317, MPFR_RNDN);
        mpfr_sqrt(r26319, r26318, MPFR_RNDN);
        mpfr_sub(r26320, r26319, r26300, MPFR_RNDN);
        mpfr_div(r26321, r26320, r26311, MPFR_RNDN);
        mpfr_div(r26322, r26321, r26304, MPFR_RNDN);
        mpfr_div(r26323, r26306, r26300, MPFR_RNDN);
        ;
        mpfr_div(r26325, r26324, r26311, MPFR_RNDN);
        mpfr_mul(r26326, r26323, r26325, MPFR_RNDN);
        if (mpfr_get_si(r26315, MPFR_RNDN)) { mpfr_set(r26327, r26322, MPFR_RNDN); } else { mpfr_set(r26327, r26326, MPFR_RNDN); };
        if (mpfr_get_si(r26302, MPFR_RNDN)) { mpfr_set(r26328, r26313, MPFR_RNDN); } else { mpfr_set(r26328, r26327, MPFR_RNDN); };
        return mpfr_get_d(r26328, MPFR_RNDN);
}

static mpfr_t r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351, r26352, r26353, r26354, r26355, r26356, r26357;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26329);
        mpfr_init_set_str(r26330, "-3.273058773164801e+121", 10, MPFR_RNDN);
        mpfr_init(r26331);
        mpfr_init_set_str(r26332, "3/2", 10, MPFR_RNDN);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init(r26339);
        mpfr_init_set_str(r26340, "3", 10, MPFR_RNDN);
        mpfr_init(r26341);
        mpfr_init(r26342);
        mpfr_init_set_str(r26343, "1.0696681173689353e-137", 10, MPFR_RNDN);
        mpfr_init(r26344);
        mpfr_init(r26345);
        mpfr_init(r26346);
        mpfr_init(r26347);
        mpfr_init(r26348);
        mpfr_init(r26349);
        mpfr_init(r26350);
        mpfr_init(r26351);
        mpfr_init(r26352);
        mpfr_init_set_str(r26353, "-3/2", 10, MPFR_RNDN);
        mpfr_init(r26354);
        mpfr_init(r26355);
        mpfr_init(r26356);
        mpfr_init(r26357);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26329, b, MPFR_RNDN);
        ;
        mpfr_set_si(r26331, mpfr_cmp(r26329, r26330) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r26333, a, MPFR_RNDN);
        mpfr_mul(r26334, r26332, r26333, MPFR_RNDN);
        mpfr_set_d(r26335, c, MPFR_RNDN);
        mpfr_div(r26336, r26329, r26335, MPFR_RNDN);
        mpfr_div(r26337, r26334, r26336, MPFR_RNDN);
        mpfr_add(r26338, r26329, r26329, MPFR_RNDN);
        mpfr_sub(r26339, r26337, r26338, MPFR_RNDN);
        ;
        mpfr_mul(r26341, r26340, r26333, MPFR_RNDN);
        mpfr_div(r26342, r26339, r26341, MPFR_RNDN);
        ;
        mpfr_set_si(r26344, mpfr_cmp(r26329, r26343) <= 0, MPFR_RNDN);
        mpfr_neg(r26345, r26335, MPFR_RNDN);
        mpfr_mul(r26346, r26329, r26329, MPFR_RNDN);
        mpfr_fma(r26347, r26345, r26341, r26346, MPFR_RNDN);
        mpfr_sqrt(r26348, r26347, MPFR_RNDN);
        mpfr_sub(r26349, r26348, r26329, MPFR_RNDN);
        mpfr_div(r26350, r26349, r26340, MPFR_RNDN);
        mpfr_div(r26351, r26350, r26333, MPFR_RNDN);
        mpfr_div(r26352, r26335, r26329, MPFR_RNDN);
        ;
        mpfr_div(r26354, r26353, r26340, MPFR_RNDN);
        mpfr_mul(r26355, r26352, r26354, MPFR_RNDN);
        if (mpfr_get_si(r26344, MPFR_RNDN)) { mpfr_set(r26356, r26351, MPFR_RNDN); } else { mpfr_set(r26356, r26355, MPFR_RNDN); };
        if (mpfr_get_si(r26331, MPFR_RNDN)) { mpfr_set(r26357, r26342, MPFR_RNDN); } else { mpfr_set(r26357, r26356, MPFR_RNDN); };
        return mpfr_get_d(r26357, MPFR_RNDN);
}