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

char *name = "Quadratic roots, medium range";

double f_if(float a, float b, float c) {
        float r8365 = b;
        float r8366 = -r8365;
        float r8367 = r8365 * r8365;
        float r8368 = 4;
        float r8369 = a;
        float r8370 = r8368 * r8369;
        float r8371 = c;
        float r8372 = r8370 * r8371;
        float r8373 = r8367 - r8372;
        float r8374 = sqrt(r8373);
        float r8375 = r8366 + r8374;
        float r8376 = 2;
        float r8377 = r8376 * r8369;
        float r8378 = r8375 / r8377;
        return r8378;
}

double f_id(double a, double b, double c) {
        double r8379 = b;
        double r8380 = -r8379;
        double r8381 = r8379 * r8379;
        double r8382 = 4;
        double r8383 = a;
        double r8384 = r8382 * r8383;
        double r8385 = c;
        double r8386 = r8384 * r8385;
        double r8387 = r8381 - r8386;
        double r8388 = sqrt(r8387);
        double r8389 = r8380 + r8388;
        double r8390 = 2;
        double r8391 = r8390 * r8383;
        double r8392 = r8389 / r8391;
        return r8392;
}


double f_of(float a, float b, float c) {
        float r8393 = 4;
        float r8394 = 2;
        float r8395 = r8393 / r8394;
        float r8396 = c;
        float r8397 = b;
        float r8398 = -r8397;
        float r8399 = r8397 * r8397;
        float r8400 = r8396 * r8393;
        float r8401 = a;
        float r8402 = r8400 * r8401;
        float r8403 = r8399 - r8402;
        float r8404 = sqrt(r8403);
        float r8405 = r8398 - r8404;
        float r8406 = r8396 / r8405;
        float r8407 = r8395 * r8406;
        return r8407;
}

double f_od(double a, double b, double c) {
        double r8408 = 4;
        double r8409 = 2;
        double r8410 = r8408 / r8409;
        double r8411 = c;
        double r8412 = b;
        double r8413 = -r8412;
        double r8414 = r8412 * r8412;
        double r8415 = r8411 * r8408;
        double r8416 = a;
        double r8417 = r8415 * r8416;
        double r8418 = r8414 - r8417;
        double r8419 = sqrt(r8418);
        double r8420 = r8413 - r8419;
        double r8421 = r8411 / r8420;
        double r8422 = r8410 * r8421;
        return r8422;
}

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 r8423, r8424, r8425, r8426, r8427, r8428, r8429, r8430, r8431, r8432, r8433, r8434, r8435, r8436;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8423);
        mpfr_init(r8424);
        mpfr_init(r8425);
        mpfr_init_set_str(r8426, "4", 10, MPFR_RNDN);
        mpfr_init(r8427);
        mpfr_init(r8428);
        mpfr_init(r8429);
        mpfr_init(r8430);
        mpfr_init(r8431);
        mpfr_init(r8432);
        mpfr_init(r8433);
        mpfr_init_set_str(r8434, "2", 10, MPFR_RNDN);
        mpfr_init(r8435);
        mpfr_init(r8436);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r8423, b, MPFR_RNDN);
        mpfr_neg(r8424, r8423, MPFR_RNDN);
        mpfr_mul(r8425, r8423, r8423, MPFR_RNDN);
        ;
        mpfr_set_d(r8427, a, MPFR_RNDN);
        mpfr_mul(r8428, r8426, r8427, MPFR_RNDN);
        mpfr_set_d(r8429, c, MPFR_RNDN);
        mpfr_mul(r8430, r8428, r8429, MPFR_RNDN);
        mpfr_sub(r8431, r8425, r8430, MPFR_RNDN);
        mpfr_sqrt(r8432, r8431, MPFR_RNDN);
        mpfr_add(r8433, r8424, r8432, MPFR_RNDN);
        ;
        mpfr_mul(r8435, r8434, r8427, MPFR_RNDN);
        mpfr_div(r8436, r8433, r8435, MPFR_RNDN);
        return mpfr_get_d(r8436, MPFR_RNDN);
}

static mpfr_t r8437, r8438, r8439, r8440, r8441, r8442, r8443, r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8437, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r8438, "2", 10, MPFR_RNDN);
        mpfr_init(r8439);
        mpfr_init(r8440);
        mpfr_init(r8441);
        mpfr_init(r8442);
        mpfr_init(r8443);
        mpfr_init(r8444);
        mpfr_init(r8445);
        mpfr_init(r8446);
        mpfr_init(r8447);
        mpfr_init(r8448);
        mpfr_init(r8449);
        mpfr_init(r8450);
        mpfr_init(r8451);
}

double f_fm(double a, double b, double c) {
        ;
        ;
        mpfr_div(r8439, r8437, r8438, MPFR_RNDN);
        mpfr_set_d(r8440, c, MPFR_RNDN);
        mpfr_set_d(r8441, b, MPFR_RNDN);
        mpfr_neg(r8442, r8441, MPFR_RNDN);
        mpfr_mul(r8443, r8441, r8441, MPFR_RNDN);
        mpfr_mul(r8444, r8440, r8437, MPFR_RNDN);
        mpfr_set_d(r8445, a, MPFR_RNDN);
        mpfr_mul(r8446, r8444, r8445, MPFR_RNDN);
        mpfr_sub(r8447, r8443, r8446, MPFR_RNDN);
        mpfr_sqrt(r8448, r8447, MPFR_RNDN);
        mpfr_sub(r8449, r8442, r8448, MPFR_RNDN);
        mpfr_div(r8450, r8440, r8449, MPFR_RNDN);
        mpfr_mul(r8451, r8439, r8450, MPFR_RNDN);
        return mpfr_get_d(r8451, MPFR_RNDN);
}

static mpfr_t r8452, r8453, r8454, r8455, r8456, r8457, r8458, r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8452, "4", 10, MPFR_RNDN);
        mpfr_init_set_str(r8453, "2", 10, MPFR_RNDN);
        mpfr_init(r8454);
        mpfr_init(r8455);
        mpfr_init(r8456);
        mpfr_init(r8457);
        mpfr_init(r8458);
        mpfr_init(r8459);
        mpfr_init(r8460);
        mpfr_init(r8461);
        mpfr_init(r8462);
        mpfr_init(r8463);
        mpfr_init(r8464);
        mpfr_init(r8465);
        mpfr_init(r8466);
}

double f_dm(double a, double b, double c) {
        ;
        ;
        mpfr_div(r8454, r8452, r8453, MPFR_RNDN);
        mpfr_set_d(r8455, c, MPFR_RNDN);
        mpfr_set_d(r8456, b, MPFR_RNDN);
        mpfr_neg(r8457, r8456, MPFR_RNDN);
        mpfr_mul(r8458, r8456, r8456, MPFR_RNDN);
        mpfr_mul(r8459, r8455, r8452, MPFR_RNDN);
        mpfr_set_d(r8460, a, MPFR_RNDN);
        mpfr_mul(r8461, r8459, r8460, MPFR_RNDN);
        mpfr_sub(r8462, r8458, r8461, MPFR_RNDN);
        mpfr_sqrt(r8463, r8462, MPFR_RNDN);
        mpfr_sub(r8464, r8457, r8463, MPFR_RNDN);
        mpfr_div(r8465, r8455, r8464, MPFR_RNDN);
        mpfr_mul(r8466, r8454, r8465, MPFR_RNDN);
        return mpfr_get_d(r8466, MPFR_RNDN);
}