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

char *name = "quadp (p42, positive)";

double f_if(float a, float b, float c) {
        float r25974 = b;
        float r25975 = -r25974;
        float r25976 = r25974 * r25974;
        float r25977 = 4;
        float r25978 = a;
        float r25979 = c;
        float r25980 = r25978 * r25979;
        float r25981 = r25977 * r25980;
        float r25982 = r25976 - r25981;
        float r25983 = sqrt(r25982);
        float r25984 = r25975 + r25983;
        float r25985 = 2;
        float r25986 = r25985 * r25978;
        float r25987 = r25984 / r25986;
        return r25987;
}

double f_id(double a, double b, double c) {
        double r25988 = b;
        double r25989 = -r25988;
        double r25990 = r25988 * r25988;
        double r25991 = 4;
        double r25992 = a;
        double r25993 = c;
        double r25994 = r25992 * r25993;
        double r25995 = r25991 * r25994;
        double r25996 = r25990 - r25995;
        double r25997 = sqrt(r25996);
        double r25998 = r25989 + r25997;
        double r25999 = 2;
        double r26000 = r25999 * r25992;
        double r26001 = r25998 / r26000;
        return r26001;
}


double f_of(float a, float b, float c) {
        float r26002 = b;
        float r26003 = -r26002;
        float r26004 = -5.7583530328199956e-65;
        bool r26005 = r26003 <= r26004;
        float r26006 = c;
        float r26007 = r26006 / r26002;
        float r26008 = -r26007;
        float r26009 = 7.455938382452604e+91;
        bool r26010 = r26003 <= r26009;
        float r26011 = r26002 * r26002;
        float r26012 = 4;
        float r26013 = a;
        float r26014 = r26006 * r26013;
        float r26015 = r26012 * r26014;
        float r26016 = r26011 - r26015;
        float r26017 = sqrt(r26016);
        float r26018 = r26017 + r26003;
        float r26019 = 2;
        float r26020 = r26019 * r26013;
        float r26021 = r26018 / r26020;
        float r26022 = r26002 + r26002;
        float r26023 = r26022 / r26020;
        float r26024 = r26007 - r26023;
        float r26025 = r26010 ? r26021 : r26024;
        float r26026 = r26005 ? r26008 : r26025;
        return r26026;
}

double f_od(double a, double b, double c) {
        double r26027 = b;
        double r26028 = -r26027;
        double r26029 = -5.7583530328199956e-65;
        bool r26030 = r26028 <= r26029;
        double r26031 = c;
        double r26032 = r26031 / r26027;
        double r26033 = -r26032;
        double r26034 = 7.455938382452604e+91;
        bool r26035 = r26028 <= r26034;
        double r26036 = r26027 * r26027;
        double r26037 = 4;
        double r26038 = a;
        double r26039 = r26031 * r26038;
        double r26040 = r26037 * r26039;
        double r26041 = r26036 - r26040;
        double r26042 = sqrt(r26041);
        double r26043 = r26042 + r26028;
        double r26044 = 2;
        double r26045 = r26044 * r26038;
        double r26046 = r26043 / r26045;
        double r26047 = r26027 + r26027;
        double r26048 = r26047 / r26045;
        double r26049 = r26032 - r26048;
        double r26050 = r26035 ? r26046 : r26049;
        double r26051 = r26030 ? r26033 : r26050;
        return r26051;
}

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 r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init_set_str(r26055, "4", 10, MPFR_RNDN);
        mpfr_init(r26056);
        mpfr_init(r26057);
        mpfr_init(r26058);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init(r26061);
        mpfr_init(r26062);
        mpfr_init_set_str(r26063, "2", 10, MPFR_RNDN);
        mpfr_init(r26064);
        mpfr_init(r26065);
}

double f_im(double a, double b, double c) {
        mpfr_set_d(r26052, b, MPFR_RNDN);
        mpfr_neg(r26053, r26052, MPFR_RNDN);
        mpfr_mul(r26054, r26052, r26052, MPFR_RNDN);
        ;
        mpfr_set_d(r26056, a, MPFR_RNDN);
        mpfr_set_d(r26057, c, MPFR_RNDN);
        mpfr_mul(r26058, r26056, r26057, MPFR_RNDN);
        mpfr_mul(r26059, r26055, r26058, MPFR_RNDN);
        mpfr_sub(r26060, r26054, r26059, MPFR_RNDN);
        mpfr_sqrt(r26061, r26060, MPFR_RNDN);
        mpfr_add(r26062, r26053, r26061, MPFR_RNDN);
        ;
        mpfr_mul(r26064, r26063, r26056, MPFR_RNDN);
        mpfr_div(r26065, r26062, r26064, MPFR_RNDN);
        return mpfr_get_d(r26065, MPFR_RNDN);
}

static mpfr_t r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26066);
        mpfr_init(r26067);
        mpfr_init_set_str(r26068, "-5.7583530328199956e-65", 10, MPFR_RNDN);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
        mpfr_init_set_str(r26073, "7.455938382452604e+91", 10, MPFR_RNDN);
        mpfr_init(r26074);
        mpfr_init(r26075);
        mpfr_init_set_str(r26076, "4", 10, MPFR_RNDN);
        mpfr_init(r26077);
        mpfr_init(r26078);
        mpfr_init(r26079);
        mpfr_init(r26080);
        mpfr_init(r26081);
        mpfr_init(r26082);
        mpfr_init_set_str(r26083, "2", 10, MPFR_RNDN);
        mpfr_init(r26084);
        mpfr_init(r26085);
        mpfr_init(r26086);
        mpfr_init(r26087);
        mpfr_init(r26088);
        mpfr_init(r26089);
        mpfr_init(r26090);
}

double f_fm(double a, double b, double c) {
        mpfr_set_d(r26066, b, MPFR_RNDN);
        mpfr_neg(r26067, r26066, MPFR_RNDN);
        ;
        mpfr_set_si(r26069, mpfr_cmp(r26067, r26068) <= 0, MPFR_RNDN);
        mpfr_set_d(r26070, c, MPFR_RNDN);
        mpfr_div(r26071, r26070, r26066, MPFR_RNDN);
        mpfr_neg(r26072, r26071, MPFR_RNDN);
        ;
        mpfr_set_si(r26074, mpfr_cmp(r26067, r26073) <= 0, MPFR_RNDN);
        mpfr_mul(r26075, r26066, r26066, MPFR_RNDN);
        ;
        mpfr_set_d(r26077, a, MPFR_RNDN);
        mpfr_mul(r26078, r26070, r26077, MPFR_RNDN);
        mpfr_mul(r26079, r26076, r26078, MPFR_RNDN);
        mpfr_sub(r26080, r26075, r26079, MPFR_RNDN);
        mpfr_sqrt(r26081, r26080, MPFR_RNDN);
        mpfr_add(r26082, r26081, r26067, MPFR_RNDN);
        ;
        mpfr_mul(r26084, r26083, r26077, MPFR_RNDN);
        mpfr_div(r26085, r26082, r26084, MPFR_RNDN);
        mpfr_add(r26086, r26066, r26066, MPFR_RNDN);
        mpfr_div(r26087, r26086, r26084, MPFR_RNDN);
        mpfr_sub(r26088, r26071, r26087, MPFR_RNDN);
        if (mpfr_get_si(r26074, MPFR_RNDN)) { mpfr_set(r26089, r26085, MPFR_RNDN); } else { mpfr_set(r26089, r26088, MPFR_RNDN); };
        if (mpfr_get_si(r26069, MPFR_RNDN)) { mpfr_set(r26090, r26072, MPFR_RNDN); } else { mpfr_set(r26090, r26089, MPFR_RNDN); };
        return mpfr_get_d(r26090, MPFR_RNDN);
}

static mpfr_t r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3408);
        mpfr_init(r26091);
        mpfr_init(r26092);
        mpfr_init_set_str(r26093, "-5.7583530328199956e-65", 10, MPFR_RNDN);
        mpfr_init(r26094);
        mpfr_init(r26095);
        mpfr_init(r26096);
        mpfr_init(r26097);
        mpfr_init_set_str(r26098, "7.455938382452604e+91", 10, MPFR_RNDN);
        mpfr_init(r26099);
        mpfr_init(r26100);
        mpfr_init_set_str(r26101, "4", 10, MPFR_RNDN);
        mpfr_init(r26102);
        mpfr_init(r26103);
        mpfr_init(r26104);
        mpfr_init(r26105);
        mpfr_init(r26106);
        mpfr_init(r26107);
        mpfr_init_set_str(r26108, "2", 10, MPFR_RNDN);
        mpfr_init(r26109);
        mpfr_init(r26110);
        mpfr_init(r26111);
        mpfr_init(r26112);
        mpfr_init(r26113);
        mpfr_init(r26114);
        mpfr_init(r26115);
}

double f_dm(double a, double b, double c) {
        mpfr_set_d(r26091, b, MPFR_RNDN);
        mpfr_neg(r26092, r26091, MPFR_RNDN);
        ;
        mpfr_set_si(r26094, mpfr_cmp(r26092, r26093) <= 0, MPFR_RNDN);
        mpfr_set_d(r26095, c, MPFR_RNDN);
        mpfr_div(r26096, r26095, r26091, MPFR_RNDN);
        mpfr_neg(r26097, r26096, MPFR_RNDN);
        ;
        mpfr_set_si(r26099, mpfr_cmp(r26092, r26098) <= 0, MPFR_RNDN);
        mpfr_mul(r26100, r26091, r26091, MPFR_RNDN);
        ;
        mpfr_set_d(r26102, a, MPFR_RNDN);
        mpfr_mul(r26103, r26095, r26102, MPFR_RNDN);
        mpfr_mul(r26104, r26101, r26103, MPFR_RNDN);
        mpfr_sub(r26105, r26100, r26104, MPFR_RNDN);
        mpfr_sqrt(r26106, r26105, MPFR_RNDN);
        mpfr_add(r26107, r26106, r26092, MPFR_RNDN);
        ;
        mpfr_mul(r26109, r26108, r26102, MPFR_RNDN);
        mpfr_div(r26110, r26107, r26109, MPFR_RNDN);
        mpfr_add(r26111, r26091, r26091, MPFR_RNDN);
        mpfr_div(r26112, r26111, r26109, MPFR_RNDN);
        mpfr_sub(r26113, r26096, r26112, MPFR_RNDN);
        if (mpfr_get_si(r26099, MPFR_RNDN)) { mpfr_set(r26114, r26110, MPFR_RNDN); } else { mpfr_set(r26114, r26113, MPFR_RNDN); };
        if (mpfr_get_si(r26094, MPFR_RNDN)) { mpfr_set(r26115, r26097, MPFR_RNDN); } else { mpfr_set(r26115, r26114, MPFR_RNDN); };
        return mpfr_get_d(r26115, MPFR_RNDN);
}