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

char *name = "The quadratic formula (r1)";

double f_if(float a, float b, float c) {
        float r27274 = b;
        float r27275 = -r27274;
        float r27276 = r27274 * r27274;
        float r27277 = 4;
        float r27278 = a;
        float r27279 = r27277 * r27278;
        float r27280 = c;
        float r27281 = r27279 * r27280;
        float r27282 = r27276 - r27281;
        float r27283 = sqrt(r27282);
        float r27284 = r27275 + r27283;
        float r27285 = 2;
        float r27286 = r27285 * r27278;
        float r27287 = r27284 / r27286;
        return r27287;
}

double f_id(double a, double b, double c) {
        double r27288 = b;
        double r27289 = -r27288;
        double r27290 = r27288 * r27288;
        double r27291 = 4;
        double r27292 = a;
        double r27293 = r27291 * r27292;
        double r27294 = c;
        double r27295 = r27293 * r27294;
        double r27296 = r27290 - r27295;
        double r27297 = sqrt(r27296);
        double r27298 = r27289 + r27297;
        double r27299 = 2;
        double r27300 = r27299 * r27292;
        double r27301 = r27298 / r27300;
        return r27301;
}


double f_of(float a, float b, float c) {
        float r27302 = b;
        float r27303 = -r27302;
        float r27304 = -5.7583530328199956e-65;
        bool r27305 = r27303 <= r27304;
        float r27306 = c;
        float r27307 = -r27306;
        float r27308 = r27307 / r27302;
        float r27309 = 5.9344413950308096e+91;
        bool r27310 = r27303 <= r27309;
        float r27311 = 1;
        float r27312 = a;
        float r27313 = 2;
        float r27314 = r27312 * r27313;
        float r27315 = 4;
        float r27316 = -r27315;
        float r27317 = r27312 * r27306;
        float r27318 = r27302 * r27302;
        float r27319 = fma(r27316, r27317, r27318);
        float r27320 = sqrt(r27319);
        float r27321 = r27320 - r27302;
        float r27322 = r27314 / r27321;
        float r27323 = r27311 / r27322;
        float r27324 = r27306 / r27302;
        float r27325 = r27302 + r27302;
        float r27326 = r27325 / r27314;
        float r27327 = r27324 - r27326;
        float r27328 = r27310 ? r27323 : r27327;
        float r27329 = r27305 ? r27308 : r27328;
        return r27329;
}

double f_od(double a, double b, double c) {
        double r27330 = b;
        double r27331 = -r27330;
        double r27332 = -5.7583530328199956e-65;
        bool r27333 = r27331 <= r27332;
        double r27334 = c;
        double r27335 = -r27334;
        double r27336 = r27335 / r27330;
        double r27337 = 5.9344413950308096e+91;
        bool r27338 = r27331 <= r27337;
        double r27339 = 1;
        double r27340 = a;
        double r27341 = 2;
        double r27342 = r27340 * r27341;
        double r27343 = 4;
        double r27344 = -r27343;
        double r27345 = r27340 * r27334;
        double r27346 = r27330 * r27330;
        double r27347 = fma(r27344, r27345, r27346);
        double r27348 = sqrt(r27347);
        double r27349 = r27348 - r27330;
        double r27350 = r27342 / r27349;
        double r27351 = r27339 / r27350;
        double r27352 = r27334 / r27330;
        double r27353 = r27330 + r27330;
        double r27354 = r27353 / r27342;
        double r27355 = r27352 - r27354;
        double r27356 = r27338 ? r27351 : r27355;
        double r27357 = r27333 ? r27336 : r27356;
        return r27357;
}

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 r27358, r27359, r27360, r27361, r27362, r27363, r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r27358);
        mpfr_init(r27359);
        mpfr_init(r27360);
        mpfr_init_set_str(r27361, "4", 10, MPFR_RNDN);
        mpfr_init(r27362);
        mpfr_init(r27363);
        mpfr_init(r27364);
        mpfr_init(r27365);
        mpfr_init(r27366);
        mpfr_init(r27367);
        mpfr_init(r27368);
        mpfr_init_set_str(r27369, "2", 10, MPFR_RNDN);
        mpfr_init(r27370);
        mpfr_init(r27371);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r27358, b, MPFR_RNDN);
        mpfr_neg(r27359, r27358, MPFR_RNDN);
        mpfr_mul(r27360, r27358, r27358, MPFR_RNDN);
        ;
        mpfr_set_d(r27362, a, MPFR_RNDN);
        mpfr_mul(r27363, r27361, r27362, MPFR_RNDN);
        mpfr_set_d(r27364, c, MPFR_RNDN);
        mpfr_mul(r27365, r27363, r27364, MPFR_RNDN);
        mpfr_sub(r27366, r27360, r27365, MPFR_RNDN);
        mpfr_sqrt(r27367, r27366, MPFR_RNDN);
        mpfr_add(r27368, r27359, r27367, MPFR_RNDN);
        ;
        mpfr_mul(r27370, r27369, r27362, MPFR_RNDN);
        mpfr_div(r27371, r27368, r27370, MPFR_RNDN);
        return mpfr_get_d(r27371, MPFR_RNDN);
}

static mpfr_t r27372, r27373, r27374, r27375, r27376, r27377, r27378, r27379, r27380, r27381, r27382, r27383, r27384, r27385, r27386, r27387, r27388, r27389, r27390, r27391, r27392, r27393, r27394, r27395, r27396, r27397, r27398, r27399;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r27372);
        mpfr_init(r27373);
        mpfr_init_set_str(r27374, "-5.7583530328199956e-65", 10, MPFR_RNDN);
        mpfr_init(r27375);
        mpfr_init(r27376);
        mpfr_init(r27377);
        mpfr_init(r27378);
        mpfr_init_set_str(r27379, "5.9344413950308096e+91", 10, MPFR_RNDN);
        mpfr_init(r27380);
        mpfr_init_set_str(r27381, "1", 10, MPFR_RNDN);
        mpfr_init(r27382);
        mpfr_init_set_str(r27383, "2", 10, MPFR_RNDN);
        mpfr_init(r27384);
        mpfr_init_set_str(r27385, "4", 10, MPFR_RNDN);
        mpfr_init(r27386);
        mpfr_init(r27387);
        mpfr_init(r27388);
        mpfr_init(r27389);
        mpfr_init(r27390);
        mpfr_init(r27391);
        mpfr_init(r27392);
        mpfr_init(r27393);
        mpfr_init(r27394);
        mpfr_init(r27395);
        mpfr_init(r27396);
        mpfr_init(r27397);
        mpfr_init(r27398);
        mpfr_init(r27399);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r27372, b, MPFR_RNDN);
        mpfr_neg(r27373, r27372, MPFR_RNDN);
        ;
        mpfr_set_si(r27375, mpfr_cmp(r27373, r27374) <= 0, MPFR_RNDN);
        mpfr_set_d(r27376, c, MPFR_RNDN);
        mpfr_neg(r27377, r27376, MPFR_RNDN);
        mpfr_div(r27378, r27377, r27372, MPFR_RNDN);
        ;
        mpfr_set_si(r27380, mpfr_cmp(r27373, r27379) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r27382, a, MPFR_RNDN);
        ;
        mpfr_mul(r27384, r27382, r27383, MPFR_RNDN);
        ;
        mpfr_neg(r27386, r27385, MPFR_RNDN);
        mpfr_mul(r27387, r27382, r27376, MPFR_RNDN);
        mpfr_mul(r27388, r27372, r27372, MPFR_RNDN);
        mpfr_fma(r27389, r27386, r27387, r27388, MPFR_RNDN);
        mpfr_sqrt(r27390, r27389, MPFR_RNDN);
        mpfr_sub(r27391, r27390, r27372, MPFR_RNDN);
        mpfr_div(r27392, r27384, r27391, MPFR_RNDN);
        mpfr_div(r27393, r27381, r27392, MPFR_RNDN);
        mpfr_div(r27394, r27376, r27372, MPFR_RNDN);
        mpfr_add(r27395, r27372, r27372, MPFR_RNDN);
        mpfr_div(r27396, r27395, r27384, MPFR_RNDN);
        mpfr_sub(r27397, r27394, r27396, MPFR_RNDN);
        if (mpfr_get_si(r27380, MPFR_RNDN)) { mpfr_set(r27398, r27393, MPFR_RNDN); } else { mpfr_set(r27398, r27397, MPFR_RNDN); };
        if (mpfr_get_si(r27375, MPFR_RNDN)) { mpfr_set(r27399, r27378, MPFR_RNDN); } else { mpfr_set(r27399, r27398, MPFR_RNDN); };
        return mpfr_get_d(r27399, MPFR_RNDN);
}

static mpfr_t r27400, r27401, r27402, r27403, r27404, r27405, r27406, r27407, r27408, r27409, r27410, r27411, r27412, r27413, r27414, r27415, r27416, r27417, r27418, r27419, r27420, r27421, r27422, r27423, r27424, r27425, r27426, r27427;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r27400);
        mpfr_init(r27401);
        mpfr_init_set_str(r27402, "-5.7583530328199956e-65", 10, MPFR_RNDN);
        mpfr_init(r27403);
        mpfr_init(r27404);
        mpfr_init(r27405);
        mpfr_init(r27406);
        mpfr_init_set_str(r27407, "5.9344413950308096e+91", 10, MPFR_RNDN);
        mpfr_init(r27408);
        mpfr_init_set_str(r27409, "1", 10, MPFR_RNDN);
        mpfr_init(r27410);
        mpfr_init_set_str(r27411, "2", 10, MPFR_RNDN);
        mpfr_init(r27412);
        mpfr_init_set_str(r27413, "4", 10, MPFR_RNDN);
        mpfr_init(r27414);
        mpfr_init(r27415);
        mpfr_init(r27416);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init(r27420);
        mpfr_init(r27421);
        mpfr_init(r27422);
        mpfr_init(r27423);
        mpfr_init(r27424);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init(r27427);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r27400, b, MPFR_RNDN);
        mpfr_neg(r27401, r27400, MPFR_RNDN);
        ;
        mpfr_set_si(r27403, mpfr_cmp(r27401, r27402) <= 0, MPFR_RNDN);
        mpfr_set_d(r27404, c, MPFR_RNDN);
        mpfr_neg(r27405, r27404, MPFR_RNDN);
        mpfr_div(r27406, r27405, r27400, MPFR_RNDN);
        ;
        mpfr_set_si(r27408, mpfr_cmp(r27401, r27407) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r27410, a, MPFR_RNDN);
        ;
        mpfr_mul(r27412, r27410, r27411, MPFR_RNDN);
        ;
        mpfr_neg(r27414, r27413, MPFR_RNDN);
        mpfr_mul(r27415, r27410, r27404, MPFR_RNDN);
        mpfr_mul(r27416, r27400, r27400, MPFR_RNDN);
        mpfr_fma(r27417, r27414, r27415, r27416, MPFR_RNDN);
        mpfr_sqrt(r27418, r27417, MPFR_RNDN);
        mpfr_sub(r27419, r27418, r27400, MPFR_RNDN);
        mpfr_div(r27420, r27412, r27419, MPFR_RNDN);
        mpfr_div(r27421, r27409, r27420, MPFR_RNDN);
        mpfr_div(r27422, r27404, r27400, MPFR_RNDN);
        mpfr_add(r27423, r27400, r27400, MPFR_RNDN);
        mpfr_div(r27424, r27423, r27412, MPFR_RNDN);
        mpfr_sub(r27425, r27422, r27424, MPFR_RNDN);
        if (mpfr_get_si(r27408, MPFR_RNDN)) { mpfr_set(r27426, r27421, MPFR_RNDN); } else { mpfr_set(r27426, r27425, MPFR_RNDN); };
        if (mpfr_get_si(r27403, MPFR_RNDN)) { mpfr_set(r27427, r27406, MPFR_RNDN); } else { mpfr_set(r27427, r27426, MPFR_RNDN); };
        return mpfr_get_d(r27427, MPFR_RNDN);
}