#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 r9010 = b;
        float r9011 = -r9010;
        float r9012 = r9010 * r9010;
        float r9013 = 4.0f;
        float r9014 = a;
        float r9015 = r9013 * r9014;
        float r9016 = c;
        float r9017 = r9015 * r9016;
        float r9018 = r9012 - r9017;
        float r9019 = sqrt(r9018);
        float r9020 = r9011 + r9019;
        float r9021 = 2.0f;
        float r9022 = r9021 * r9014;
        float r9023 = r9020 / r9022;
        return r9023;
}

double f_id(double a, double b, double c) {
        double r9024 = b;
        double r9025 = -r9024;
        double r9026 = r9024 * r9024;
        double r9027 = 4.0;
        double r9028 = a;
        double r9029 = r9027 * r9028;
        double r9030 = c;
        double r9031 = r9029 * r9030;
        double r9032 = r9026 - r9031;
        double r9033 = sqrt(r9032);
        double r9034 = r9025 + r9033;
        double r9035 = 2.0;
        double r9036 = r9035 * r9028;
        double r9037 = r9034 / r9036;
        return r9037;
}


double f_of(float a, float b, float c) {
        float r9038 = b;
        float r9039 = -1.338815475246526e+154f;
        bool r9040 = r9038 <= r9039;
        float r9041 = -0.5f;
        float r9042 = a;
        float r9043 = r9038 / r9042;
        float r9044 = r9041 * r9043;
        float r9045 = 3.690239526200231e-95f;
        bool r9046 = r9038 <= r9045;
        float r9047 = 4.0f;
        float r9048 = r9042 * r9047;
        float r9049 = c;
        float r9050 = -r9049;
        float r9051 = r9038 * r9038;
        float r9052 = fma(r9048, r9050, r9051);
        float r9053 = sqrt(r9052);
        float r9054 = 2.0f;
        float r9055 = r9054 * r9042;
        float r9056 = r9053 / r9055;
        float r9057 = r9038 / r9055;
        float r9058 = r9056 - r9057;
        float r9059 = r9050 * r9048;
        float r9060 = r9053 + r9038;
        float r9061 = r9055 * r9060;
        float r9062 = r9059 / r9061;
        float r9063 = r9046 ? r9058 : r9062;
        float r9064 = r9040 ? r9044 : r9063;
        return r9064;
}

double f_od(double a, double b, double c) {
        double r9065 = b;
        double r9066 = -1.338815475246526e+154;
        bool r9067 = r9065 <= r9066;
        double r9068 = -0.5;
        double r9069 = a;
        double r9070 = r9065 / r9069;
        double r9071 = r9068 * r9070;
        double r9072 = 3.690239526200231e-95;
        bool r9073 = r9065 <= r9072;
        double r9074 = 4.0;
        double r9075 = r9069 * r9074;
        double r9076 = c;
        double r9077 = -r9076;
        double r9078 = r9065 * r9065;
        double r9079 = fma(r9075, r9077, r9078);
        double r9080 = sqrt(r9079);
        double r9081 = 2.0;
        double r9082 = r9081 * r9069;
        double r9083 = r9080 / r9082;
        double r9084 = r9065 / r9082;
        double r9085 = r9083 - r9084;
        double r9086 = r9077 * r9075;
        double r9087 = r9080 + r9065;
        double r9088 = r9082 * r9087;
        double r9089 = r9086 / r9088;
        double r9090 = r9073 ? r9085 : r9089;
        double r9091 = r9067 ? r9071 : r9090;
        return r9091;
}

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 r9092, r9093, r9094, r9095, r9096, r9097, r9098, r9099, r9100, r9101, r9102, r9103, r9104, r9105;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3152);
        mpfr_init(r9092);
        mpfr_init(r9093);
        mpfr_init(r9094);
        mpfr_init_set_str(r9095, "4", 10, MPFR_RNDN);
        mpfr_init(r9096);
        mpfr_init(r9097);
        mpfr_init(r9098);
        mpfr_init(r9099);
        mpfr_init(r9100);
        mpfr_init(r9101);
        mpfr_init(r9102);
        mpfr_init_set_str(r9103, "2", 10, MPFR_RNDN);
        mpfr_init(r9104);
        mpfr_init(r9105);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r9092, b, MPFR_RNDN);
        mpfr_neg(r9093, r9092, MPFR_RNDN);
        mpfr_mul(r9094, r9092, r9092, MPFR_RNDN);
        ;
        mpfr_set_d(r9096, a, MPFR_RNDN);
        mpfr_mul(r9097, r9095, r9096, MPFR_RNDN);
        mpfr_set_d(r9098, c, MPFR_RNDN);
        mpfr_mul(r9099, r9097, r9098, MPFR_RNDN);
        mpfr_sub(r9100, r9094, r9099, MPFR_RNDN);
        mpfr_sqrt(r9101, r9100, MPFR_RNDN);
        mpfr_add(r9102, r9093, r9101, MPFR_RNDN);
        ;
        mpfr_mul(r9104, r9103, r9096, MPFR_RNDN);
        mpfr_div(r9105, r9102, r9104, MPFR_RNDN);
        return mpfr_get_d(r9105, MPFR_RNDN);
}

static mpfr_t r9106, r9107, r9108, r9109, r9110, r9111, r9112, r9113, r9114, r9115, r9116, r9117, r9118, r9119, r9120, r9121, r9122, r9123, r9124, r9125, r9126, r9127, r9128, r9129, r9130, r9131, r9132;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r9106);
        mpfr_init_set_str(r9107, "-1.338815475246526e+154", 10, MPFR_RNDN);
        mpfr_init(r9108);
        mpfr_init_set_str(r9109, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r9110);
        mpfr_init(r9111);
        mpfr_init(r9112);
        mpfr_init_set_str(r9113, "3.690239526200231e-95", 10, MPFR_RNDN);
        mpfr_init(r9114);
        mpfr_init_set_str(r9115, "4", 10, MPFR_RNDN);
        mpfr_init(r9116);
        mpfr_init(r9117);
        mpfr_init(r9118);
        mpfr_init(r9119);
        mpfr_init(r9120);
        mpfr_init(r9121);
        mpfr_init_set_str(r9122, "2", 10, MPFR_RNDN);
        mpfr_init(r9123);
        mpfr_init(r9124);
        mpfr_init(r9125);
        mpfr_init(r9126);
        mpfr_init(r9127);
        mpfr_init(r9128);
        mpfr_init(r9129);
        mpfr_init(r9130);
        mpfr_init(r9131);
        mpfr_init(r9132);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r9106, b, MPFR_RNDN);
        ;
        mpfr_set_si(r9108, mpfr_cmp(r9106, r9107) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r9110, a, MPFR_RNDN);
        mpfr_div(r9111, r9106, r9110, MPFR_RNDN);
        mpfr_mul(r9112, r9109, r9111, MPFR_RNDN);
        ;
        mpfr_set_si(r9114, mpfr_cmp(r9106, r9113) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r9116, r9110, r9115, MPFR_RNDN);
        mpfr_set_d(r9117, c, MPFR_RNDN);
        mpfr_neg(r9118, r9117, MPFR_RNDN);
        mpfr_mul(r9119, r9106, r9106, MPFR_RNDN);
        mpfr_fma(r9120, r9116, r9118, r9119, MPFR_RNDN);
        mpfr_sqrt(r9121, r9120, MPFR_RNDN);
        ;
        mpfr_mul(r9123, r9122, r9110, MPFR_RNDN);
        mpfr_div(r9124, r9121, r9123, MPFR_RNDN);
        mpfr_div(r9125, r9106, r9123, MPFR_RNDN);
        mpfr_sub(r9126, r9124, r9125, MPFR_RNDN);
        mpfr_mul(r9127, r9118, r9116, MPFR_RNDN);
        mpfr_add(r9128, r9121, r9106, MPFR_RNDN);
        mpfr_mul(r9129, r9123, r9128, MPFR_RNDN);
        mpfr_div(r9130, r9127, r9129, MPFR_RNDN);
        if (mpfr_get_si(r9114, MPFR_RNDN)) { mpfr_set(r9131, r9126, MPFR_RNDN); } else { mpfr_set(r9131, r9130, MPFR_RNDN); };
        if (mpfr_get_si(r9108, MPFR_RNDN)) { mpfr_set(r9132, r9112, MPFR_RNDN); } else { mpfr_set(r9132, r9131, MPFR_RNDN); };
        return mpfr_get_d(r9132, MPFR_RNDN);
}

static mpfr_t r9133, r9134, r9135, r9136, r9137, r9138, r9139, r9140, r9141, r9142, r9143, r9144, r9145, r9146, r9147, r9148, r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3152);
        mpfr_init(r9133);
        mpfr_init_set_str(r9134, "-1.338815475246526e+154", 10, MPFR_RNDN);
        mpfr_init(r9135);
        mpfr_init_set_str(r9136, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r9137);
        mpfr_init(r9138);
        mpfr_init(r9139);
        mpfr_init_set_str(r9140, "3.690239526200231e-95", 10, MPFR_RNDN);
        mpfr_init(r9141);
        mpfr_init_set_str(r9142, "4", 10, MPFR_RNDN);
        mpfr_init(r9143);
        mpfr_init(r9144);
        mpfr_init(r9145);
        mpfr_init(r9146);
        mpfr_init(r9147);
        mpfr_init(r9148);
        mpfr_init_set_str(r9149, "2", 10, MPFR_RNDN);
        mpfr_init(r9150);
        mpfr_init(r9151);
        mpfr_init(r9152);
        mpfr_init(r9153);
        mpfr_init(r9154);
        mpfr_init(r9155);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
        mpfr_init(r9159);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r9133, b, MPFR_RNDN);
        ;
        mpfr_set_si(r9135, mpfr_cmp(r9133, r9134) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r9137, a, MPFR_RNDN);
        mpfr_div(r9138, r9133, r9137, MPFR_RNDN);
        mpfr_mul(r9139, r9136, r9138, MPFR_RNDN);
        ;
        mpfr_set_si(r9141, mpfr_cmp(r9133, r9140) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r9143, r9137, r9142, MPFR_RNDN);
        mpfr_set_d(r9144, c, MPFR_RNDN);
        mpfr_neg(r9145, r9144, MPFR_RNDN);
        mpfr_mul(r9146, r9133, r9133, MPFR_RNDN);
        mpfr_fma(r9147, r9143, r9145, r9146, MPFR_RNDN);
        mpfr_sqrt(r9148, r9147, MPFR_RNDN);
        ;
        mpfr_mul(r9150, r9149, r9137, MPFR_RNDN);
        mpfr_div(r9151, r9148, r9150, MPFR_RNDN);
        mpfr_div(r9152, r9133, r9150, MPFR_RNDN);
        mpfr_sub(r9153, r9151, r9152, MPFR_RNDN);
        mpfr_mul(r9154, r9145, r9143, MPFR_RNDN);
        mpfr_add(r9155, r9148, r9133, MPFR_RNDN);
        mpfr_mul(r9156, r9150, r9155, MPFR_RNDN);
        mpfr_div(r9157, r9154, r9156, MPFR_RNDN);
        if (mpfr_get_si(r9141, MPFR_RNDN)) { mpfr_set(r9158, r9153, MPFR_RNDN); } else { mpfr_set(r9158, r9157, MPFR_RNDN); };
        if (mpfr_get_si(r9135, MPFR_RNDN)) { mpfr_set(r9159, r9139, MPFR_RNDN); } else { mpfr_set(r9159, r9158, MPFR_RNDN); };
        return mpfr_get_d(r9159, MPFR_RNDN);
}