#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 r25090 = b;
        float r25091 = -r25090;
        float r25092 = r25090 * r25090;
        float r25093 = 4;
        float r25094 = a;
        float r25095 = r25093 * r25094;
        float r25096 = c;
        float r25097 = r25095 * r25096;
        float r25098 = r25092 - r25097;
        float r25099 = sqrt(r25098);
        float r25100 = r25091 + r25099;
        float r25101 = 2;
        float r25102 = r25101 * r25094;
        float r25103 = r25100 / r25102;
        return r25103;
}

double f_id(double a, double b, double c) {
        double r25104 = b;
        double r25105 = -r25104;
        double r25106 = r25104 * r25104;
        double r25107 = 4;
        double r25108 = a;
        double r25109 = r25107 * r25108;
        double r25110 = c;
        double r25111 = r25109 * r25110;
        double r25112 = r25106 - r25111;
        double r25113 = sqrt(r25112);
        double r25114 = r25105 + r25113;
        double r25115 = 2;
        double r25116 = r25115 * r25108;
        double r25117 = r25114 / r25116;
        return r25117;
}


double f_of(float a, float b, float c) {
        float r25118 = b;
        float r25119 = -1.3563250571528177e+154;
        bool r25120 = r25118 <= r25119;
        float r25121 = -1/2;
        float r25122 = a;
        float r25123 = r25118 / r25122;
        float r25124 = r25121 * r25123;
        float r25125 = 2.14399979241458e-120;
        bool r25126 = r25118 <= r25125;
        float r25127 = 4;
        float r25128 = r25127 * r25122;
        float r25129 = c;
        float r25130 = -r25129;
        float r25131 = r25118 * r25118;
        float r25132 = fma(r25128, r25130, r25131);
        float r25133 = sqrt(r25132);
        float r25134 = 2;
        float r25135 = r25134 * r25122;
        float r25136 = r25133 / r25135;
        float r25137 = r25118 / r25135;
        float r25138 = r25136 - r25137;
        float r25139 = r25122 * r25127;
        float r25140 = r25139 * r25130;
        float r25141 = r25140 / r25135;
        float r25142 = r25133 + r25118;
        float r25143 = r25141 / r25142;
        float r25144 = r25126 ? r25138 : r25143;
        float r25145 = r25120 ? r25124 : r25144;
        return r25145;
}

double f_od(double a, double b, double c) {
        double r25146 = b;
        double r25147 = -1.3563250571528177e+154;
        bool r25148 = r25146 <= r25147;
        double r25149 = -1/2;
        double r25150 = a;
        double r25151 = r25146 / r25150;
        double r25152 = r25149 * r25151;
        double r25153 = 2.14399979241458e-120;
        bool r25154 = r25146 <= r25153;
        double r25155 = 4;
        double r25156 = r25155 * r25150;
        double r25157 = c;
        double r25158 = -r25157;
        double r25159 = r25146 * r25146;
        double r25160 = fma(r25156, r25158, r25159);
        double r25161 = sqrt(r25160);
        double r25162 = 2;
        double r25163 = r25162 * r25150;
        double r25164 = r25161 / r25163;
        double r25165 = r25146 / r25163;
        double r25166 = r25164 - r25165;
        double r25167 = r25150 * r25155;
        double r25168 = r25167 * r25158;
        double r25169 = r25168 / r25163;
        double r25170 = r25161 + r25146;
        double r25171 = r25169 / r25170;
        double r25172 = r25154 ? r25166 : r25171;
        double r25173 = r25148 ? r25152 : r25172;
        return r25173;
}

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 r25174, r25175, r25176, r25177, r25178, r25179, r25180, r25181, r25182, r25183, r25184, r25185, r25186, r25187;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3152);
        mpfr_init(r25174);
        mpfr_init(r25175);
        mpfr_init(r25176);
        mpfr_init_set_str(r25177, "4", 10, MPFR_RNDN);
        mpfr_init(r25178);
        mpfr_init(r25179);
        mpfr_init(r25180);
        mpfr_init(r25181);
        mpfr_init(r25182);
        mpfr_init(r25183);
        mpfr_init(r25184);
        mpfr_init_set_str(r25185, "2", 10, MPFR_RNDN);
        mpfr_init(r25186);
        mpfr_init(r25187);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r25174, b, MPFR_RNDN);
        mpfr_neg(r25175, r25174, MPFR_RNDN);
        mpfr_mul(r25176, r25174, r25174, MPFR_RNDN);
        ;
        mpfr_set_d(r25178, a, MPFR_RNDN);
        mpfr_mul(r25179, r25177, r25178, MPFR_RNDN);
        mpfr_set_d(r25180, c, MPFR_RNDN);
        mpfr_mul(r25181, r25179, r25180, MPFR_RNDN);
        mpfr_sub(r25182, r25176, r25181, MPFR_RNDN);
        mpfr_sqrt(r25183, r25182, MPFR_RNDN);
        mpfr_add(r25184, r25175, r25183, MPFR_RNDN);
        ;
        mpfr_mul(r25186, r25185, r25178, MPFR_RNDN);
        mpfr_div(r25187, r25184, r25186, MPFR_RNDN);
        return mpfr_get_d(r25187, MPFR_RNDN);
}

static mpfr_t r25188, r25189, r25190, r25191, r25192, r25193, r25194, r25195, r25196, r25197, r25198, r25199, r25200, r25201, r25202, r25203, r25204, r25205, r25206, r25207, r25208, r25209, r25210, r25211, r25212, r25213, r25214, r25215;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r25188);
        mpfr_init_set_str(r25189, "-1.3563250571528177e+154", 10, MPFR_RNDN);
        mpfr_init(r25190);
        mpfr_init_set_str(r25191, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r25192);
        mpfr_init(r25193);
        mpfr_init(r25194);
        mpfr_init_set_str(r25195, "2.14399979241458e-120", 10, MPFR_RNDN);
        mpfr_init(r25196);
        mpfr_init_set_str(r25197, "4", 10, MPFR_RNDN);
        mpfr_init(r25198);
        mpfr_init(r25199);
        mpfr_init(r25200);
        mpfr_init(r25201);
        mpfr_init(r25202);
        mpfr_init(r25203);
        mpfr_init_set_str(r25204, "2", 10, MPFR_RNDN);
        mpfr_init(r25205);
        mpfr_init(r25206);
        mpfr_init(r25207);
        mpfr_init(r25208);
        mpfr_init(r25209);
        mpfr_init(r25210);
        mpfr_init(r25211);
        mpfr_init(r25212);
        mpfr_init(r25213);
        mpfr_init(r25214);
        mpfr_init(r25215);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r25188, b, MPFR_RNDN);
        ;
        mpfr_set_si(r25190, mpfr_cmp(r25188, r25189) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r25192, a, MPFR_RNDN);
        mpfr_div(r25193, r25188, r25192, MPFR_RNDN);
        mpfr_mul(r25194, r25191, r25193, MPFR_RNDN);
        ;
        mpfr_set_si(r25196, mpfr_cmp(r25188, r25195) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r25198, r25197, r25192, MPFR_RNDN);
        mpfr_set_d(r25199, c, MPFR_RNDN);
        mpfr_neg(r25200, r25199, MPFR_RNDN);
        mpfr_mul(r25201, r25188, r25188, MPFR_RNDN);
        mpfr_fma(r25202, r25198, r25200, r25201, MPFR_RNDN);
        mpfr_sqrt(r25203, r25202, MPFR_RNDN);
        ;
        mpfr_mul(r25205, r25204, r25192, MPFR_RNDN);
        mpfr_div(r25206, r25203, r25205, MPFR_RNDN);
        mpfr_div(r25207, r25188, r25205, MPFR_RNDN);
        mpfr_sub(r25208, r25206, r25207, MPFR_RNDN);
        mpfr_mul(r25209, r25192, r25197, MPFR_RNDN);
        mpfr_mul(r25210, r25209, r25200, MPFR_RNDN);
        mpfr_div(r25211, r25210, r25205, MPFR_RNDN);
        mpfr_add(r25212, r25203, r25188, MPFR_RNDN);
        mpfr_div(r25213, r25211, r25212, MPFR_RNDN);
        if (mpfr_get_si(r25196, MPFR_RNDN)) { mpfr_set(r25214, r25208, MPFR_RNDN); } else { mpfr_set(r25214, r25213, MPFR_RNDN); };
        if (mpfr_get_si(r25190, MPFR_RNDN)) { mpfr_set(r25215, r25194, MPFR_RNDN); } else { mpfr_set(r25215, r25214, MPFR_RNDN); };
        return mpfr_get_d(r25215, MPFR_RNDN);
}

static mpfr_t r25216, r25217, r25218, r25219, r25220, r25221, r25222, r25223, r25224, r25225, r25226, r25227, r25228, r25229, r25230, r25231, r25232, r25233, r25234, r25235, r25236, r25237, r25238, r25239, r25240, r25241, r25242, r25243;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r25216);
        mpfr_init_set_str(r25217, "-1.3563250571528177e+154", 10, MPFR_RNDN);
        mpfr_init(r25218);
        mpfr_init_set_str(r25219, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r25220);
        mpfr_init(r25221);
        mpfr_init(r25222);
        mpfr_init_set_str(r25223, "2.14399979241458e-120", 10, MPFR_RNDN);
        mpfr_init(r25224);
        mpfr_init_set_str(r25225, "4", 10, MPFR_RNDN);
        mpfr_init(r25226);
        mpfr_init(r25227);
        mpfr_init(r25228);
        mpfr_init(r25229);
        mpfr_init(r25230);
        mpfr_init(r25231);
        mpfr_init_set_str(r25232, "2", 10, MPFR_RNDN);
        mpfr_init(r25233);
        mpfr_init(r25234);
        mpfr_init(r25235);
        mpfr_init(r25236);
        mpfr_init(r25237);
        mpfr_init(r25238);
        mpfr_init(r25239);
        mpfr_init(r25240);
        mpfr_init(r25241);
        mpfr_init(r25242);
        mpfr_init(r25243);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r25216, b, MPFR_RNDN);
        ;
        mpfr_set_si(r25218, mpfr_cmp(r25216, r25217) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r25220, a, MPFR_RNDN);
        mpfr_div(r25221, r25216, r25220, MPFR_RNDN);
        mpfr_mul(r25222, r25219, r25221, MPFR_RNDN);
        ;
        mpfr_set_si(r25224, mpfr_cmp(r25216, r25223) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r25226, r25225, r25220, MPFR_RNDN);
        mpfr_set_d(r25227, c, MPFR_RNDN);
        mpfr_neg(r25228, r25227, MPFR_RNDN);
        mpfr_mul(r25229, r25216, r25216, MPFR_RNDN);
        mpfr_fma(r25230, r25226, r25228, r25229, MPFR_RNDN);
        mpfr_sqrt(r25231, r25230, MPFR_RNDN);
        ;
        mpfr_mul(r25233, r25232, r25220, MPFR_RNDN);
        mpfr_div(r25234, r25231, r25233, MPFR_RNDN);
        mpfr_div(r25235, r25216, r25233, MPFR_RNDN);
        mpfr_sub(r25236, r25234, r25235, MPFR_RNDN);
        mpfr_mul(r25237, r25220, r25225, MPFR_RNDN);
        mpfr_mul(r25238, r25237, r25228, MPFR_RNDN);
        mpfr_div(r25239, r25238, r25233, MPFR_RNDN);
        mpfr_add(r25240, r25231, r25216, MPFR_RNDN);
        mpfr_div(r25241, r25239, r25240, MPFR_RNDN);
        if (mpfr_get_si(r25224, MPFR_RNDN)) { mpfr_set(r25242, r25236, MPFR_RNDN); } else { mpfr_set(r25242, r25241, MPFR_RNDN); };
        if (mpfr_get_si(r25218, MPFR_RNDN)) { mpfr_set(r25243, r25222, MPFR_RNDN); } else { mpfr_set(r25243, r25242, MPFR_RNDN); };
        return mpfr_get_d(r25243, MPFR_RNDN);
}