#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 r26396 = b;
        float r26397 = -r26396;
        float r26398 = r26396 * r26396;
        float r26399 = 3;
        float r26400 = a;
        float r26401 = r26399 * r26400;
        float r26402 = c;
        float r26403 = r26401 * r26402;
        float r26404 = r26398 - r26403;
        float r26405 = sqrt(r26404);
        float r26406 = r26397 + r26405;
        float r26407 = r26406 / r26401;
        return r26407;
}

double f_id(double a, double b, double c, double __attribute__((unused)) d) {
        double r26408 = b;
        double r26409 = -r26408;
        double r26410 = r26408 * r26408;
        double r26411 = 3;
        double r26412 = a;
        double r26413 = r26411 * r26412;
        double r26414 = c;
        double r26415 = r26413 * r26414;
        double r26416 = r26410 - r26415;
        double r26417 = sqrt(r26416);
        double r26418 = r26409 + r26417;
        double r26419 = r26418 / r26413;
        return r26419;
}


double f_of(float a, float b, float c, float __attribute__((unused)) d) {
        float r26420 = b;
        float r26421 = -r26420;
        float r26422 = r26420 * r26420;
        float r26423 = 3;
        float r26424 = a;
        float r26425 = r26423 * r26424;
        float r26426 = c;
        float r26427 = r26425 * r26426;
        float r26428 = r26422 - r26427;
        float r26429 = sqrt(r26428);
        float r26430 = r26421 + r26429;
        float r26431 = r26430 / r26425;
        return r26431;
}

double f_od(double a, double b, double c, double __attribute__((unused)) d) {
        double r26432 = b;
        double r26433 = -r26432;
        double r26434 = r26432 * r26432;
        double r26435 = 3;
        double r26436 = a;
        double r26437 = r26435 * r26436;
        double r26438 = c;
        double r26439 = r26437 * r26438;
        double r26440 = r26434 - r26439;
        double r26441 = sqrt(r26440);
        double r26442 = r26433 + r26441;
        double r26443 = r26442 / r26437;
        return r26443;
}

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 r26444, r26445, r26446, r26447, r26448, r26449, r26450, r26451, r26452, r26453, r26454, r26455;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26444);
        mpfr_init(r26445);
        mpfr_init(r26446);
        mpfr_init_set_str(r26447, "3", 10, MPFR_RNDN);
        mpfr_init(r26448);
        mpfr_init(r26449);
        mpfr_init(r26450);
        mpfr_init(r26451);
        mpfr_init(r26452);
        mpfr_init(r26453);
        mpfr_init(r26454);
        mpfr_init(r26455);
}

double f_im(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26444, b, MPFR_RNDN);
        mpfr_neg(r26445, r26444, MPFR_RNDN);
        mpfr_mul(r26446, r26444, r26444, MPFR_RNDN);
        ;
        mpfr_set_d(r26448, a, MPFR_RNDN);
        mpfr_mul(r26449, r26447, r26448, MPFR_RNDN);
        mpfr_set_d(r26450, c, MPFR_RNDN);
        mpfr_mul(r26451, r26449, r26450, MPFR_RNDN);
        mpfr_sub(r26452, r26446, r26451, MPFR_RNDN);
        mpfr_sqrt(r26453, r26452, MPFR_RNDN);
        mpfr_add(r26454, r26445, r26453, MPFR_RNDN);
        mpfr_div(r26455, r26454, r26449, MPFR_RNDN);
        return mpfr_get_d(r26455, MPFR_RNDN);
}

static mpfr_t r26456, r26457, r26458, r26459, r26460, r26461, r26462, r26463, r26464, r26465, r26466, r26467;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26456);
        mpfr_init(r26457);
        mpfr_init(r26458);
        mpfr_init_set_str(r26459, "3", 10, MPFR_RNDN);
        mpfr_init(r26460);
        mpfr_init(r26461);
        mpfr_init(r26462);
        mpfr_init(r26463);
        mpfr_init(r26464);
        mpfr_init(r26465);
        mpfr_init(r26466);
        mpfr_init(r26467);
}

double f_fm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26456, b, MPFR_RNDN);
        mpfr_neg(r26457, r26456, MPFR_RNDN);
        mpfr_mul(r26458, r26456, r26456, MPFR_RNDN);
        ;
        mpfr_set_d(r26460, a, MPFR_RNDN);
        mpfr_mul(r26461, r26459, r26460, MPFR_RNDN);
        mpfr_set_d(r26462, c, MPFR_RNDN);
        mpfr_mul(r26463, r26461, r26462, MPFR_RNDN);
        mpfr_sub(r26464, r26458, r26463, MPFR_RNDN);
        mpfr_sqrt(r26465, r26464, MPFR_RNDN);
        mpfr_add(r26466, r26457, r26465, MPFR_RNDN);
        mpfr_div(r26467, r26466, r26461, MPFR_RNDN);
        return mpfr_get_d(r26467, MPFR_RNDN);
}

static mpfr_t r26468, r26469, r26470, r26471, r26472, r26473, r26474, r26475, r26476, r26477, r26478, r26479;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26468);
        mpfr_init(r26469);
        mpfr_init(r26470);
        mpfr_init_set_str(r26471, "3", 10, MPFR_RNDN);
        mpfr_init(r26472);
        mpfr_init(r26473);
        mpfr_init(r26474);
        mpfr_init(r26475);
        mpfr_init(r26476);
        mpfr_init(r26477);
        mpfr_init(r26478);
        mpfr_init(r26479);
}

double f_dm(double a, double b, double c, double __attribute__((unused)) d) {
        mpfr_set_d(r26468, b, MPFR_RNDN);
        mpfr_neg(r26469, r26468, MPFR_RNDN);
        mpfr_mul(r26470, r26468, r26468, MPFR_RNDN);
        ;
        mpfr_set_d(r26472, a, MPFR_RNDN);
        mpfr_mul(r26473, r26471, r26472, MPFR_RNDN);
        mpfr_set_d(r26474, c, MPFR_RNDN);
        mpfr_mul(r26475, r26473, r26474, MPFR_RNDN);
        mpfr_sub(r26476, r26470, r26475, MPFR_RNDN);
        mpfr_sqrt(r26477, r26476, MPFR_RNDN);
        mpfr_add(r26478, r26469, r26477, MPFR_RNDN);
        mpfr_div(r26479, r26478, r26473, MPFR_RNDN);
        return mpfr_get_d(r26479, MPFR_RNDN);
}