#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 r17448 = b;
        float r17449 = -r17448;
        float r17450 = r17448 * r17448;
        float r17451 = 4.0f;
        float r17452 = a;
        float r17453 = r17451 * r17452;
        float r17454 = c;
        float r17455 = r17453 * r17454;
        float r17456 = r17450 - r17455;
        float r17457 = sqrt(r17456);
        float r17458 = r17449 + r17457;
        float r17459 = 2.0f;
        float r17460 = r17459 * r17452;
        float r17461 = r17458 / r17460;
        return r17461;
}

double f_id(double a, double b, double c) {
        double r17462 = b;
        double r17463 = -r17462;
        double r17464 = r17462 * r17462;
        double r17465 = 4.0;
        double r17466 = a;
        double r17467 = r17465 * r17466;
        double r17468 = c;
        double r17469 = r17467 * r17468;
        double r17470 = r17464 - r17469;
        double r17471 = sqrt(r17470);
        double r17472 = r17463 + r17471;
        double r17473 = 2.0;
        double r17474 = r17473 * r17466;
        double r17475 = r17472 / r17474;
        return r17475;
}


double f_of(float a, float b, float c) {
        float r17476 = b;
        float r17477 = -r17476;
        float r17478 = r17476 * r17476;
        float r17479 = 4.0f;
        float r17480 = a;
        float r17481 = r17479 * r17480;
        float r17482 = c;
        float r17483 = r17481 * r17482;
        float r17484 = r17478 - r17483;
        float r17485 = sqrt(r17484);
        float r17486 = r17477 + r17485;
        float r17487 = 2.0f;
        float r17488 = r17487 * r17480;
        float r17489 = r17486 / r17488;
        return r17489;
}

double f_od(double a, double b, double c) {
        double r17490 = b;
        double r17491 = -r17490;
        double r17492 = r17490 * r17490;
        double r17493 = 4.0;
        double r17494 = a;
        double r17495 = r17493 * r17494;
        double r17496 = c;
        double r17497 = r17495 * r17496;
        double r17498 = r17492 - r17497;
        double r17499 = sqrt(r17498);
        double r17500 = r17491 + r17499;
        double r17501 = 2.0;
        double r17502 = r17501 * r17494;
        double r17503 = r17500 / r17502;
        return r17503;
}

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 r17504, r17505, r17506, r17507, r17508, r17509, r17510, r17511, r17512, r17513, r17514, r17515, r17516, r17517;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17504);
        mpfr_init(r17505);
        mpfr_init(r17506);
        mpfr_init_set_str(r17507, "4", 10, MPFR_RNDN);
        mpfr_init(r17508);
        mpfr_init(r17509);
        mpfr_init(r17510);
        mpfr_init(r17511);
        mpfr_init(r17512);
        mpfr_init(r17513);
        mpfr_init(r17514);
        mpfr_init_set_str(r17515, "2", 10, MPFR_RNDN);
        mpfr_init(r17516);
        mpfr_init(r17517);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r17504, b, MPFR_RNDN);
        mpfr_neg(r17505, r17504, MPFR_RNDN);
        mpfr_sqr(r17506, r17504, MPFR_RNDN);
        ;
        mpfr_set_d(r17508, a, MPFR_RNDN);
        mpfr_mul(r17509, r17507, r17508, MPFR_RNDN);
        mpfr_set_d(r17510, c, MPFR_RNDN);
        mpfr_mul(r17511, r17509, r17510, MPFR_RNDN);
        mpfr_sub(r17512, r17506, r17511, MPFR_RNDN);
        mpfr_sqrt(r17513, r17512, MPFR_RNDN);
        mpfr_add(r17514, r17505, r17513, MPFR_RNDN);
        ;
        mpfr_mul(r17516, r17515, r17508, MPFR_RNDN);
        mpfr_div(r17517, r17514, r17516, MPFR_RNDN);
        return mpfr_get_d(r17517, MPFR_RNDN);
}

static mpfr_t r17518, r17519, r17520, r17521, r17522, r17523, r17524, r17525, r17526, r17527, r17528, r17529, r17530, r17531;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17518);
        mpfr_init(r17519);
        mpfr_init(r17520);
        mpfr_init_set_str(r17521, "4", 10, MPFR_RNDN);
        mpfr_init(r17522);
        mpfr_init(r17523);
        mpfr_init(r17524);
        mpfr_init(r17525);
        mpfr_init(r17526);
        mpfr_init(r17527);
        mpfr_init(r17528);
        mpfr_init_set_str(r17529, "2", 10, MPFR_RNDN);
        mpfr_init(r17530);
        mpfr_init(r17531);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r17518, b, MPFR_RNDN);
        mpfr_neg(r17519, r17518, MPFR_RNDN);
        mpfr_sqr(r17520, r17518, MPFR_RNDN);
        ;
        mpfr_set_d(r17522, a, MPFR_RNDN);
        mpfr_mul(r17523, r17521, r17522, MPFR_RNDN);
        mpfr_set_d(r17524, c, MPFR_RNDN);
        mpfr_mul(r17525, r17523, r17524, MPFR_RNDN);
        mpfr_sub(r17526, r17520, r17525, MPFR_RNDN);
        mpfr_sqrt(r17527, r17526, MPFR_RNDN);
        mpfr_add(r17528, r17519, r17527, MPFR_RNDN);
        ;
        mpfr_mul(r17530, r17529, r17522, MPFR_RNDN);
        mpfr_div(r17531, r17528, r17530, MPFR_RNDN);
        return mpfr_get_d(r17531, MPFR_RNDN);
}

static mpfr_t r17532, r17533, r17534, r17535, r17536, r17537, r17538, r17539, r17540, r17541, r17542, r17543, r17544, r17545;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17532);
        mpfr_init(r17533);
        mpfr_init(r17534);
        mpfr_init_set_str(r17535, "4", 10, MPFR_RNDN);
        mpfr_init(r17536);
        mpfr_init(r17537);
        mpfr_init(r17538);
        mpfr_init(r17539);
        mpfr_init(r17540);
        mpfr_init(r17541);
        mpfr_init(r17542);
        mpfr_init_set_str(r17543, "2", 10, MPFR_RNDN);
        mpfr_init(r17544);
        mpfr_init(r17545);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r17532, b, MPFR_RNDN);
        mpfr_neg(r17533, r17532, MPFR_RNDN);
        mpfr_sqr(r17534, r17532, MPFR_RNDN);
        ;
        mpfr_set_d(r17536, a, MPFR_RNDN);
        mpfr_mul(r17537, r17535, r17536, MPFR_RNDN);
        mpfr_set_d(r17538, c, MPFR_RNDN);
        mpfr_mul(r17539, r17537, r17538, MPFR_RNDN);
        mpfr_sub(r17540, r17534, r17539, MPFR_RNDN);
        mpfr_sqrt(r17541, r17540, MPFR_RNDN);
        mpfr_add(r17542, r17533, r17541, MPFR_RNDN);
        ;
        mpfr_mul(r17544, r17543, r17536, MPFR_RNDN);
        mpfr_div(r17545, r17542, r17544, MPFR_RNDN);
        return mpfr_get_d(r17545, MPFR_RNDN);
}