#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 r18316 = b;
        float r18317 = -r18316;
        float r18318 = r18316 * r18316;
        float r18319 = 4.0f;
        float r18320 = a;
        float r18321 = r18319 * r18320;
        float r18322 = c;
        float r18323 = r18321 * r18322;
        float r18324 = r18318 - r18323;
        float r18325 = sqrt(r18324);
        float r18326 = r18317 + r18325;
        float r18327 = 2.0f;
        float r18328 = r18327 * r18320;
        float r18329 = r18326 / r18328;
        return r18329;
}

double f_id(double a, double b, double c) {
        double r18330 = b;
        double r18331 = -r18330;
        double r18332 = r18330 * r18330;
        double r18333 = 4.0;
        double r18334 = a;
        double r18335 = r18333 * r18334;
        double r18336 = c;
        double r18337 = r18335 * r18336;
        double r18338 = r18332 - r18337;
        double r18339 = sqrt(r18338);
        double r18340 = r18331 + r18339;
        double r18341 = 2.0;
        double r18342 = r18341 * r18334;
        double r18343 = r18340 / r18342;
        return r18343;
}


double f_of(float a, float b, float c) {
        float r18344 = b;
        float r18345 = -4.191830935951668e+150f;
        bool r18346 = r18344 <= r18345;
        float r18347 = c;
        float r18348 = r18347 / r18344;
        float r18349 = a;
        float r18350 = r18344 / r18349;
        float r18351 = r18348 - r18350;
        float r18352 = 7.444658085812687e-61f;
        bool r18353 = r18344 <= r18352;
        float r18354 = -r18344;
        float r18355 = r18344 * r18344;
        float r18356 = 4.0f;
        float r18357 = r18356 * r18349;
        float r18358 = r18357 * r18347;
        float r18359 = r18355 - r18358;
        float r18360 = sqrt(r18359);
        float r18361 = r18354 + r18360;
        float r18362 = 2.0f;
        float r18363 = r18362 * r18349;
        float r18364 = r18361 / r18363;
        float r18365 = r18356 / r18362;
        float r18366 = r18365 * r18347;
        float r18367 = r18349 * r18362;
        float r18368 = fma(r18367, r18348, r18354);
        float r18369 = r18368 - r18344;
        float r18370 = r18366 / r18369;
        float r18371 = r18353 ? r18364 : r18370;
        float r18372 = r18346 ? r18351 : r18371;
        return r18372;
}

double f_od(double a, double b, double c) {
        double r18373 = b;
        double r18374 = -4.191830935951668e+150;
        bool r18375 = r18373 <= r18374;
        double r18376 = c;
        double r18377 = r18376 / r18373;
        double r18378 = a;
        double r18379 = r18373 / r18378;
        double r18380 = r18377 - r18379;
        double r18381 = 7.444658085812687e-61;
        bool r18382 = r18373 <= r18381;
        double r18383 = -r18373;
        double r18384 = r18373 * r18373;
        double r18385 = 4.0;
        double r18386 = r18385 * r18378;
        double r18387 = r18386 * r18376;
        double r18388 = r18384 - r18387;
        double r18389 = sqrt(r18388);
        double r18390 = r18383 + r18389;
        double r18391 = 2.0;
        double r18392 = r18391 * r18378;
        double r18393 = r18390 / r18392;
        double r18394 = r18385 / r18391;
        double r18395 = r18394 * r18376;
        double r18396 = r18378 * r18391;
        double r18397 = fma(r18396, r18377, r18383);
        double r18398 = r18397 - r18373;
        double r18399 = r18395 / r18398;
        double r18400 = r18382 ? r18393 : r18399;
        double r18401 = r18375 ? r18380 : r18400;
        return r18401;
}

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 r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410, r18411, r18412, r18413, r18414, r18415;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18402);
        mpfr_init(r18403);
        mpfr_init(r18404);
        mpfr_init_set_str(r18405, "4", 10, MPFR_RNDN);
        mpfr_init(r18406);
        mpfr_init(r18407);
        mpfr_init(r18408);
        mpfr_init(r18409);
        mpfr_init(r18410);
        mpfr_init(r18411);
        mpfr_init(r18412);
        mpfr_init_set_str(r18413, "2", 10, MPFR_RNDN);
        mpfr_init(r18414);
        mpfr_init(r18415);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r18402, b, MPFR_RNDN);
        mpfr_neg(r18403, r18402, MPFR_RNDN);
        mpfr_sqr(r18404, r18402, MPFR_RNDN);
        ;
        mpfr_set_d(r18406, a, MPFR_RNDN);
        mpfr_mul(r18407, r18405, r18406, MPFR_RNDN);
        mpfr_set_d(r18408, c, MPFR_RNDN);
        mpfr_mul(r18409, r18407, r18408, MPFR_RNDN);
        mpfr_sub(r18410, r18404, r18409, MPFR_RNDN);
        mpfr_sqrt(r18411, r18410, MPFR_RNDN);
        mpfr_add(r18412, r18403, r18411, MPFR_RNDN);
        ;
        mpfr_mul(r18414, r18413, r18406, MPFR_RNDN);
        mpfr_div(r18415, r18412, r18414, MPFR_RNDN);
        return mpfr_get_d(r18415, MPFR_RNDN);
}

static mpfr_t r18416, r18417, r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432, r18433, r18434, r18435, r18436, r18437, r18438, r18439, r18440, r18441, r18442, r18443, r18444;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18416);
        mpfr_init_set_str(r18417, "-4.191830935951668e+150", 10, MPFR_RNDN);
        mpfr_init(r18418);
        mpfr_init(r18419);
        mpfr_init(r18420);
        mpfr_init(r18421);
        mpfr_init(r18422);
        mpfr_init(r18423);
        mpfr_init_set_str(r18424, "7.444658085812687e-61", 10, MPFR_RNDN);
        mpfr_init(r18425);
        mpfr_init(r18426);
        mpfr_init(r18427);
        mpfr_init_set_str(r18428, "4", 10, MPFR_RNDN);
        mpfr_init(r18429);
        mpfr_init(r18430);
        mpfr_init(r18431);
        mpfr_init(r18432);
        mpfr_init(r18433);
        mpfr_init_set_str(r18434, "2", 10, MPFR_RNDN);
        mpfr_init(r18435);
        mpfr_init(r18436);
        mpfr_init(r18437);
        mpfr_init(r18438);
        mpfr_init(r18439);
        mpfr_init(r18440);
        mpfr_init(r18441);
        mpfr_init(r18442);
        mpfr_init(r18443);
        mpfr_init(r18444);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r18416, b, MPFR_RNDN);
        ;
        mpfr_set_si(r18418, mpfr_cmp(r18416, r18417) <= 0, MPFR_RNDN);
        mpfr_set_d(r18419, c, MPFR_RNDN);
        mpfr_div(r18420, r18419, r18416, MPFR_RNDN);
        mpfr_set_d(r18421, a, MPFR_RNDN);
        mpfr_div(r18422, r18416, r18421, MPFR_RNDN);
        mpfr_sub(r18423, r18420, r18422, MPFR_RNDN);
        ;
        mpfr_set_si(r18425, mpfr_cmp(r18416, r18424) <= 0, MPFR_RNDN);
        mpfr_neg(r18426, r18416, MPFR_RNDN);
        mpfr_sqr(r18427, r18416, MPFR_RNDN);
        ;
        mpfr_mul(r18429, r18428, r18421, MPFR_RNDN);
        mpfr_mul(r18430, r18429, r18419, MPFR_RNDN);
        mpfr_sub(r18431, r18427, r18430, MPFR_RNDN);
        mpfr_sqrt(r18432, r18431, MPFR_RNDN);
        mpfr_add(r18433, r18426, r18432, MPFR_RNDN);
        ;
        mpfr_mul(r18435, r18434, r18421, MPFR_RNDN);
        mpfr_div(r18436, r18433, r18435, MPFR_RNDN);
        mpfr_div(r18437, r18428, r18434, MPFR_RNDN);
        mpfr_mul(r18438, r18437, r18419, MPFR_RNDN);
        mpfr_mul(r18439, r18421, r18434, MPFR_RNDN);
        mpfr_fma(r18440, r18439, r18420, r18426, MPFR_RNDN);
        mpfr_sub(r18441, r18440, r18416, MPFR_RNDN);
        mpfr_div(r18442, r18438, r18441, MPFR_RNDN);
        if (mpfr_get_si(r18425, MPFR_RNDN)) { mpfr_set(r18443, r18436, MPFR_RNDN); } else { mpfr_set(r18443, r18442, MPFR_RNDN); };
        if (mpfr_get_si(r18418, MPFR_RNDN)) { mpfr_set(r18444, r18423, MPFR_RNDN); } else { mpfr_set(r18444, r18443, MPFR_RNDN); };
        return mpfr_get_d(r18444, MPFR_RNDN);
}

static mpfr_t r18445, r18446, r18447, r18448, r18449, r18450, r18451, r18452, r18453, r18454, r18455, r18456, r18457, r18458, r18459, r18460, r18461, r18462, r18463, r18464, r18465, r18466, r18467, r18468, r18469, r18470, r18471, r18472, r18473;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18445);
        mpfr_init_set_str(r18446, "-4.191830935951668e+150", 10, MPFR_RNDN);
        mpfr_init(r18447);
        mpfr_init(r18448);
        mpfr_init(r18449);
        mpfr_init(r18450);
        mpfr_init(r18451);
        mpfr_init(r18452);
        mpfr_init_set_str(r18453, "7.444658085812687e-61", 10, MPFR_RNDN);
        mpfr_init(r18454);
        mpfr_init(r18455);
        mpfr_init(r18456);
        mpfr_init_set_str(r18457, "4", 10, MPFR_RNDN);
        mpfr_init(r18458);
        mpfr_init(r18459);
        mpfr_init(r18460);
        mpfr_init(r18461);
        mpfr_init(r18462);
        mpfr_init_set_str(r18463, "2", 10, MPFR_RNDN);
        mpfr_init(r18464);
        mpfr_init(r18465);
        mpfr_init(r18466);
        mpfr_init(r18467);
        mpfr_init(r18468);
        mpfr_init(r18469);
        mpfr_init(r18470);
        mpfr_init(r18471);
        mpfr_init(r18472);
        mpfr_init(r18473);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r18445, b, MPFR_RNDN);
        ;
        mpfr_set_si(r18447, mpfr_cmp(r18445, r18446) <= 0, MPFR_RNDN);
        mpfr_set_d(r18448, c, MPFR_RNDN);
        mpfr_div(r18449, r18448, r18445, MPFR_RNDN);
        mpfr_set_d(r18450, a, MPFR_RNDN);
        mpfr_div(r18451, r18445, r18450, MPFR_RNDN);
        mpfr_sub(r18452, r18449, r18451, MPFR_RNDN);
        ;
        mpfr_set_si(r18454, mpfr_cmp(r18445, r18453) <= 0, MPFR_RNDN);
        mpfr_neg(r18455, r18445, MPFR_RNDN);
        mpfr_sqr(r18456, r18445, MPFR_RNDN);
        ;
        mpfr_mul(r18458, r18457, r18450, MPFR_RNDN);
        mpfr_mul(r18459, r18458, r18448, MPFR_RNDN);
        mpfr_sub(r18460, r18456, r18459, MPFR_RNDN);
        mpfr_sqrt(r18461, r18460, MPFR_RNDN);
        mpfr_add(r18462, r18455, r18461, MPFR_RNDN);
        ;
        mpfr_mul(r18464, r18463, r18450, MPFR_RNDN);
        mpfr_div(r18465, r18462, r18464, MPFR_RNDN);
        mpfr_div(r18466, r18457, r18463, MPFR_RNDN);
        mpfr_mul(r18467, r18466, r18448, MPFR_RNDN);
        mpfr_mul(r18468, r18450, r18463, MPFR_RNDN);
        mpfr_fma(r18469, r18468, r18449, r18455, MPFR_RNDN);
        mpfr_sub(r18470, r18469, r18445, MPFR_RNDN);
        mpfr_div(r18471, r18467, r18470, MPFR_RNDN);
        if (mpfr_get_si(r18454, MPFR_RNDN)) { mpfr_set(r18472, r18465, MPFR_RNDN); } else { mpfr_set(r18472, r18471, MPFR_RNDN); };
        if (mpfr_get_si(r18447, MPFR_RNDN)) { mpfr_set(r18473, r18452, MPFR_RNDN); } else { mpfr_set(r18473, r18472, MPFR_RNDN); };
        return mpfr_get_d(r18473, MPFR_RNDN);
}