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

char *name = "jeff quadratic root 2";

double f_if(float a, float b, float c) {
        float r8303 = b;
        float r8304 = 0;
        bool r8305 = r8303 >= r8304;
        float r8306 = 2;
        float r8307 = c;
        float r8308 = r8306 * r8307;
        float r8309 = -r8303;
        float r8310 = r8303 * r8303;
        float r8311 = 4;
        float r8312 = a;
        float r8313 = r8311 * r8312;
        float r8314 = r8313 * r8307;
        float r8315 = r8310 - r8314;
        float r8316 = sqrt(r8315);
        float r8317 = r8309 - r8316;
        float r8318 = r8308 / r8317;
        float r8319 = r8309 + r8316;
        float r8320 = r8306 * r8312;
        float r8321 = r8319 / r8320;
        float r8322 = r8305 ? r8318 : r8321;
        return r8322;
}

double f_id(double a, double b, double c) {
        double r8323 = b;
        double r8324 = 0;
        bool r8325 = r8323 >= r8324;
        double r8326 = 2;
        double r8327 = c;
        double r8328 = r8326 * r8327;
        double r8329 = -r8323;
        double r8330 = r8323 * r8323;
        double r8331 = 4;
        double r8332 = a;
        double r8333 = r8331 * r8332;
        double r8334 = r8333 * r8327;
        double r8335 = r8330 - r8334;
        double r8336 = sqrt(r8335);
        double r8337 = r8329 - r8336;
        double r8338 = r8328 / r8337;
        double r8339 = r8329 + r8336;
        double r8340 = r8326 * r8332;
        double r8341 = r8339 / r8340;
        double r8342 = r8325 ? r8338 : r8341;
        return r8342;
}


double f_of(float a, float b, float c) {
        float r8343 = b;
        float r8344 = -1.3172728653882788e+154;
        bool r8345 = r8343 <= r8344;
        float r8346 = 0;
        bool r8347 = r8343 >= r8346;
        float r8348 = c;
        float r8349 = 2;
        float r8350 = r8348 * r8349;
        float r8351 = -r8343;
        float r8352 = r8343 * r8343;
        float r8353 = 4;
        float r8354 = a;
        float r8355 = r8353 * r8354;
        float r8356 = r8348 * r8355;
        float r8357 = r8352 - r8356;
        float r8358 = sqrt(r8357);
        float r8359 = r8351 - r8358;
        float r8360 = r8350 / r8359;
        float r8361 = r8348 / r8343;
        float r8362 = r8343 / r8354;
        float r8363 = r8361 - r8362;
        float r8364 = r8347 ? r8360 : r8363;
        float r8365 = 8.022931310201715e+36;
        bool r8366 = r8343 <= r8365;
        float r8367 = sqrt(r8358);
        float r8368 = r8367 * r8367;
        float r8369 = r8351 - r8368;
        float r8370 = r8350 / r8369;
        float r8371 = r8358 + r8351;
        float r8372 = r8354 * r8349;
        float r8373 = r8371 / r8372;
        float r8374 = r8347 ? r8370 : r8373;
        float r8375 = r8354 / r8343;
        float r8376 = r8350 * r8375;
        float r8377 = r8343 + r8343;
        float r8378 = r8376 - r8377;
        float r8379 = r8350 / r8378;
        float r8380 = r8353 * r8348;
        float r8381 = r8380 * r8354;
        float r8382 = r8352 - r8381;
        float r8383 = sqrt(r8382);
        float r8384 = r8383 - r8343;
        float r8385 = r8384 / r8372;
        float r8386 = r8347 ? r8379 : r8385;
        float r8387 = r8366 ? r8374 : r8386;
        float r8388 = r8345 ? r8364 : r8387;
        return r8388;
}

double f_od(double a, double b, double c) {
        double r8389 = b;
        double r8390 = -1.3172728653882788e+154;
        bool r8391 = r8389 <= r8390;
        double r8392 = 0;
        bool r8393 = r8389 >= r8392;
        double r8394 = c;
        double r8395 = 2;
        double r8396 = r8394 * r8395;
        double r8397 = -r8389;
        double r8398 = r8389 * r8389;
        double r8399 = 4;
        double r8400 = a;
        double r8401 = r8399 * r8400;
        double r8402 = r8394 * r8401;
        double r8403 = r8398 - r8402;
        double r8404 = sqrt(r8403);
        double r8405 = r8397 - r8404;
        double r8406 = r8396 / r8405;
        double r8407 = r8394 / r8389;
        double r8408 = r8389 / r8400;
        double r8409 = r8407 - r8408;
        double r8410 = r8393 ? r8406 : r8409;
        double r8411 = 8.022931310201715e+36;
        bool r8412 = r8389 <= r8411;
        double r8413 = sqrt(r8404);
        double r8414 = r8413 * r8413;
        double r8415 = r8397 - r8414;
        double r8416 = r8396 / r8415;
        double r8417 = r8404 + r8397;
        double r8418 = r8400 * r8395;
        double r8419 = r8417 / r8418;
        double r8420 = r8393 ? r8416 : r8419;
        double r8421 = r8400 / r8389;
        double r8422 = r8396 * r8421;
        double r8423 = r8389 + r8389;
        double r8424 = r8422 - r8423;
        double r8425 = r8396 / r8424;
        double r8426 = r8399 * r8394;
        double r8427 = r8426 * r8400;
        double r8428 = r8398 - r8427;
        double r8429 = sqrt(r8428);
        double r8430 = r8429 - r8389;
        double r8431 = r8430 / r8418;
        double r8432 = r8393 ? r8425 : r8431;
        double r8433 = r8412 ? r8420 : r8432;
        double r8434 = r8391 ? r8410 : r8433;
        return r8434;
}

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 r8435, r8436, r8437, r8438, r8439, r8440, r8441, r8442, r8443, r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451, r8452, r8453, r8454;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r8435);
        mpfr_init_set_str(r8436, "0", 10, MPFR_RNDN);
        mpfr_init(r8437);
        mpfr_init_set_str(r8438, "2", 10, MPFR_RNDN);
        mpfr_init(r8439);
        mpfr_init(r8440);
        mpfr_init(r8441);
        mpfr_init(r8442);
        mpfr_init_set_str(r8443, "4", 10, MPFR_RNDN);
        mpfr_init(r8444);
        mpfr_init(r8445);
        mpfr_init(r8446);
        mpfr_init(r8447);
        mpfr_init(r8448);
        mpfr_init(r8449);
        mpfr_init(r8450);
        mpfr_init(r8451);
        mpfr_init(r8452);
        mpfr_init(r8453);
        mpfr_init(r8454);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r8435, b, MPFR_RNDN);
        ;
        mpfr_set_si(r8437, mpfr_cmp(r8435, r8436) >= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r8439, c, MPFR_RNDN);
        mpfr_mul(r8440, r8438, r8439, MPFR_RNDN);
        mpfr_neg(r8441, r8435, MPFR_RNDN);
        mpfr_mul(r8442, r8435, r8435, MPFR_RNDN);
        ;
        mpfr_set_d(r8444, a, MPFR_RNDN);
        mpfr_mul(r8445, r8443, r8444, MPFR_RNDN);
        mpfr_mul(r8446, r8445, r8439, MPFR_RNDN);
        mpfr_sub(r8447, r8442, r8446, MPFR_RNDN);
        mpfr_sqrt(r8448, r8447, MPFR_RNDN);
        mpfr_sub(r8449, r8441, r8448, MPFR_RNDN);
        mpfr_div(r8450, r8440, r8449, MPFR_RNDN);
        mpfr_add(r8451, r8441, r8448, MPFR_RNDN);
        mpfr_mul(r8452, r8438, r8444, MPFR_RNDN);
        mpfr_div(r8453, r8451, r8452, MPFR_RNDN);
        if (mpfr_get_si(r8437, MPFR_RNDN)) { mpfr_set(r8454, r8450, MPFR_RNDN); } else { mpfr_set(r8454, r8453, MPFR_RNDN); };
        return mpfr_get_d(r8454, MPFR_RNDN);
}

static mpfr_t r8455, r8456, r8457, r8458, r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466, r8467, r8468, r8469, r8470, r8471, r8472, r8473, r8474, r8475, r8476, r8477, r8478, r8479, r8480, r8481, r8482, r8483, r8484, r8485, r8486, r8487, r8488, r8489, r8490, r8491, r8492, r8493, r8494, r8495, r8496, r8497, r8498, r8499, r8500;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8455);
        mpfr_init_set_str(r8456, "-1.3172728653882788e+154", 10, MPFR_RNDN);
        mpfr_init(r8457);
        mpfr_init_set_str(r8458, "0", 10, MPFR_RNDN);
        mpfr_init(r8459);
        mpfr_init(r8460);
        mpfr_init_set_str(r8461, "2", 10, MPFR_RNDN);
        mpfr_init(r8462);
        mpfr_init(r8463);
        mpfr_init(r8464);
        mpfr_init_set_str(r8465, "4", 10, MPFR_RNDN);
        mpfr_init(r8466);
        mpfr_init(r8467);
        mpfr_init(r8468);
        mpfr_init(r8469);
        mpfr_init(r8470);
        mpfr_init(r8471);
        mpfr_init(r8472);
        mpfr_init(r8473);
        mpfr_init(r8474);
        mpfr_init(r8475);
        mpfr_init(r8476);
        mpfr_init_set_str(r8477, "8.022931310201715e+36", 10, MPFR_RNDN);
        mpfr_init(r8478);
        mpfr_init(r8479);
        mpfr_init(r8480);
        mpfr_init(r8481);
        mpfr_init(r8482);
        mpfr_init(r8483);
        mpfr_init(r8484);
        mpfr_init(r8485);
        mpfr_init(r8486);
        mpfr_init(r8487);
        mpfr_init(r8488);
        mpfr_init(r8489);
        mpfr_init(r8490);
        mpfr_init(r8491);
        mpfr_init(r8492);
        mpfr_init(r8493);
        mpfr_init(r8494);
        mpfr_init(r8495);
        mpfr_init(r8496);
        mpfr_init(r8497);
        mpfr_init(r8498);
        mpfr_init(r8499);
        mpfr_init(r8500);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r8455, b, MPFR_RNDN);
        ;
        mpfr_set_si(r8457, mpfr_cmp(r8455, r8456) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r8459, mpfr_cmp(r8455, r8458) >= 0, MPFR_RNDN);
        mpfr_set_d(r8460, c, MPFR_RNDN);
        ;
        mpfr_mul(r8462, r8460, r8461, MPFR_RNDN);
        mpfr_neg(r8463, r8455, MPFR_RNDN);
        mpfr_mul(r8464, r8455, r8455, MPFR_RNDN);
        ;
        mpfr_set_d(r8466, a, MPFR_RNDN);
        mpfr_mul(r8467, r8465, r8466, MPFR_RNDN);
        mpfr_mul(r8468, r8460, r8467, MPFR_RNDN);
        mpfr_sub(r8469, r8464, r8468, MPFR_RNDN);
        mpfr_sqrt(r8470, r8469, MPFR_RNDN);
        mpfr_sub(r8471, r8463, r8470, MPFR_RNDN);
        mpfr_div(r8472, r8462, r8471, MPFR_RNDN);
        mpfr_div(r8473, r8460, r8455, MPFR_RNDN);
        mpfr_div(r8474, r8455, r8466, MPFR_RNDN);
        mpfr_sub(r8475, r8473, r8474, MPFR_RNDN);
        if (mpfr_get_si(r8459, MPFR_RNDN)) { mpfr_set(r8476, r8472, MPFR_RNDN); } else { mpfr_set(r8476, r8475, MPFR_RNDN); };
        ;
        mpfr_set_si(r8478, mpfr_cmp(r8455, r8477) <= 0, MPFR_RNDN);
        mpfr_sqrt(r8479, r8470, MPFR_RNDN);
        mpfr_mul(r8480, r8479, r8479, MPFR_RNDN);
        mpfr_sub(r8481, r8463, r8480, MPFR_RNDN);
        mpfr_div(r8482, r8462, r8481, MPFR_RNDN);
        mpfr_add(r8483, r8470, r8463, MPFR_RNDN);
        mpfr_mul(r8484, r8466, r8461, MPFR_RNDN);
        mpfr_div(r8485, r8483, r8484, MPFR_RNDN);
        if (mpfr_get_si(r8459, MPFR_RNDN)) { mpfr_set(r8486, r8482, MPFR_RNDN); } else { mpfr_set(r8486, r8485, MPFR_RNDN); };
        mpfr_div(r8487, r8466, r8455, MPFR_RNDN);
        mpfr_mul(r8488, r8462, r8487, MPFR_RNDN);
        mpfr_add(r8489, r8455, r8455, MPFR_RNDN);
        mpfr_sub(r8490, r8488, r8489, MPFR_RNDN);
        mpfr_div(r8491, r8462, r8490, MPFR_RNDN);
        mpfr_mul(r8492, r8465, r8460, MPFR_RNDN);
        mpfr_mul(r8493, r8492, r8466, MPFR_RNDN);
        mpfr_sub(r8494, r8464, r8493, MPFR_RNDN);
        mpfr_sqrt(r8495, r8494, MPFR_RNDN);
        mpfr_sub(r8496, r8495, r8455, MPFR_RNDN);
        mpfr_div(r8497, r8496, r8484, MPFR_RNDN);
        if (mpfr_get_si(r8459, MPFR_RNDN)) { mpfr_set(r8498, r8491, MPFR_RNDN); } else { mpfr_set(r8498, r8497, MPFR_RNDN); };
        if (mpfr_get_si(r8478, MPFR_RNDN)) { mpfr_set(r8499, r8486, MPFR_RNDN); } else { mpfr_set(r8499, r8498, MPFR_RNDN); };
        if (mpfr_get_si(r8457, MPFR_RNDN)) { mpfr_set(r8500, r8476, MPFR_RNDN); } else { mpfr_set(r8500, r8499, MPFR_RNDN); };
        return mpfr_get_d(r8500, MPFR_RNDN);
}

static mpfr_t r8501, r8502, r8503, r8504, r8505, r8506, r8507, r8508, r8509, r8510, r8511, r8512, r8513, r8514, r8515, r8516, r8517, r8518, r8519, r8520, r8521, r8522, r8523, r8524, r8525, r8526, r8527, r8528, r8529, r8530, r8531, r8532, r8533, r8534, r8535, r8536, r8537, r8538, r8539, r8540, r8541, r8542, r8543, r8544, r8545, r8546;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r8501);
        mpfr_init_set_str(r8502, "-1.3172728653882788e+154", 10, MPFR_RNDN);
        mpfr_init(r8503);
        mpfr_init_set_str(r8504, "0", 10, MPFR_RNDN);
        mpfr_init(r8505);
        mpfr_init(r8506);
        mpfr_init_set_str(r8507, "2", 10, MPFR_RNDN);
        mpfr_init(r8508);
        mpfr_init(r8509);
        mpfr_init(r8510);
        mpfr_init_set_str(r8511, "4", 10, MPFR_RNDN);
        mpfr_init(r8512);
        mpfr_init(r8513);
        mpfr_init(r8514);
        mpfr_init(r8515);
        mpfr_init(r8516);
        mpfr_init(r8517);
        mpfr_init(r8518);
        mpfr_init(r8519);
        mpfr_init(r8520);
        mpfr_init(r8521);
        mpfr_init(r8522);
        mpfr_init_set_str(r8523, "8.022931310201715e+36", 10, MPFR_RNDN);
        mpfr_init(r8524);
        mpfr_init(r8525);
        mpfr_init(r8526);
        mpfr_init(r8527);
        mpfr_init(r8528);
        mpfr_init(r8529);
        mpfr_init(r8530);
        mpfr_init(r8531);
        mpfr_init(r8532);
        mpfr_init(r8533);
        mpfr_init(r8534);
        mpfr_init(r8535);
        mpfr_init(r8536);
        mpfr_init(r8537);
        mpfr_init(r8538);
        mpfr_init(r8539);
        mpfr_init(r8540);
        mpfr_init(r8541);
        mpfr_init(r8542);
        mpfr_init(r8543);
        mpfr_init(r8544);
        mpfr_init(r8545);
        mpfr_init(r8546);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r8501, b, MPFR_RNDN);
        ;
        mpfr_set_si(r8503, mpfr_cmp(r8501, r8502) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r8505, mpfr_cmp(r8501, r8504) >= 0, MPFR_RNDN);
        mpfr_set_d(r8506, c, MPFR_RNDN);
        ;
        mpfr_mul(r8508, r8506, r8507, MPFR_RNDN);
        mpfr_neg(r8509, r8501, MPFR_RNDN);
        mpfr_mul(r8510, r8501, r8501, MPFR_RNDN);
        ;
        mpfr_set_d(r8512, a, MPFR_RNDN);
        mpfr_mul(r8513, r8511, r8512, MPFR_RNDN);
        mpfr_mul(r8514, r8506, r8513, MPFR_RNDN);
        mpfr_sub(r8515, r8510, r8514, MPFR_RNDN);
        mpfr_sqrt(r8516, r8515, MPFR_RNDN);
        mpfr_sub(r8517, r8509, r8516, MPFR_RNDN);
        mpfr_div(r8518, r8508, r8517, MPFR_RNDN);
        mpfr_div(r8519, r8506, r8501, MPFR_RNDN);
        mpfr_div(r8520, r8501, r8512, MPFR_RNDN);
        mpfr_sub(r8521, r8519, r8520, MPFR_RNDN);
        if (mpfr_get_si(r8505, MPFR_RNDN)) { mpfr_set(r8522, r8518, MPFR_RNDN); } else { mpfr_set(r8522, r8521, MPFR_RNDN); };
        ;
        mpfr_set_si(r8524, mpfr_cmp(r8501, r8523) <= 0, MPFR_RNDN);
        mpfr_sqrt(r8525, r8516, MPFR_RNDN);
        mpfr_mul(r8526, r8525, r8525, MPFR_RNDN);
        mpfr_sub(r8527, r8509, r8526, MPFR_RNDN);
        mpfr_div(r8528, r8508, r8527, MPFR_RNDN);
        mpfr_add(r8529, r8516, r8509, MPFR_RNDN);
        mpfr_mul(r8530, r8512, r8507, MPFR_RNDN);
        mpfr_div(r8531, r8529, r8530, MPFR_RNDN);
        if (mpfr_get_si(r8505, MPFR_RNDN)) { mpfr_set(r8532, r8528, MPFR_RNDN); } else { mpfr_set(r8532, r8531, MPFR_RNDN); };
        mpfr_div(r8533, r8512, r8501, MPFR_RNDN);
        mpfr_mul(r8534, r8508, r8533, MPFR_RNDN);
        mpfr_add(r8535, r8501, r8501, MPFR_RNDN);
        mpfr_sub(r8536, r8534, r8535, MPFR_RNDN);
        mpfr_div(r8537, r8508, r8536, MPFR_RNDN);
        mpfr_mul(r8538, r8511, r8506, MPFR_RNDN);
        mpfr_mul(r8539, r8538, r8512, MPFR_RNDN);
        mpfr_sub(r8540, r8510, r8539, MPFR_RNDN);
        mpfr_sqrt(r8541, r8540, MPFR_RNDN);
        mpfr_sub(r8542, r8541, r8501, MPFR_RNDN);
        mpfr_div(r8543, r8542, r8530, MPFR_RNDN);
        if (mpfr_get_si(r8505, MPFR_RNDN)) { mpfr_set(r8544, r8537, MPFR_RNDN); } else { mpfr_set(r8544, r8543, MPFR_RNDN); };
        if (mpfr_get_si(r8524, MPFR_RNDN)) { mpfr_set(r8545, r8532, MPFR_RNDN); } else { mpfr_set(r8545, r8544, MPFR_RNDN); };
        if (mpfr_get_si(r8503, MPFR_RNDN)) { mpfr_set(r8546, r8522, MPFR_RNDN); } else { mpfr_set(r8546, r8545, MPFR_RNDN); };
        return mpfr_get_d(r8546, MPFR_RNDN);
}