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

char *name = "quad2p (problem 3.2.1, positive)";

double f_if(float a, float b_2F2, float c) {
        float r19651 = b_2F2;
        float r19652 = -r19651;
        float r19653 = r19651 * r19651;
        float r19654 = a;
        float r19655 = c;
        float r19656 = r19654 * r19655;
        float r19657 = r19653 - r19656;
        float r19658 = sqrt(r19657);
        float r19659 = r19652 + r19658;
        float r19660 = r19659 / r19654;
        return r19660;
}

double f_id(double a, double b_2F2, double c) {
        double r19661 = b_2F2;
        double r19662 = -r19661;
        double r19663 = r19661 * r19661;
        double r19664 = a;
        double r19665 = c;
        double r19666 = r19664 * r19665;
        double r19667 = r19663 - r19666;
        double r19668 = sqrt(r19667);
        double r19669 = r19662 + r19668;
        double r19670 = r19669 / r19664;
        return r19670;
}


double f_of(float a, float b_2F2, float c) {
        float r19671 = b_2F2;
        float r19672 = -2.788370916558726e+153;
        bool r19673 = r19671 <= r19672;
        float r19674 = 1/2;
        float r19675 = c;
        float r19676 = r19675 / r19671;
        float r19677 = r19674 * r19676;
        float r19678 = 2;
        float r19679 = a;
        float r19680 = r19671 / r19679;
        float r19681 = r19678 * r19680;
        float r19682 = r19677 - r19681;
        float r19683 = 7.580412743766101e-138;
        bool r19684 = r19671 <= r19683;
        float r19685 = -r19671;
        float r19686 = r19671 * r19671;
        float r19687 = r19679 * r19675;
        float r19688 = r19686 - r19687;
        float r19689 = sqrt(r19688);
        float r19690 = r19685 + r19689;
        float r19691 = r19690 / r19679;
        float r19692 = -1/2;
        float r19693 = r19676 * r19692;
        float r19694 = r19684 ? r19691 : r19693;
        float r19695 = r19673 ? r19682 : r19694;
        return r19695;
}

double f_od(double a, double b_2F2, double c) {
        double r19696 = b_2F2;
        double r19697 = -2.788370916558726e+153;
        bool r19698 = r19696 <= r19697;
        double r19699 = 1/2;
        double r19700 = c;
        double r19701 = r19700 / r19696;
        double r19702 = r19699 * r19701;
        double r19703 = 2;
        double r19704 = a;
        double r19705 = r19696 / r19704;
        double r19706 = r19703 * r19705;
        double r19707 = r19702 - r19706;
        double r19708 = 7.580412743766101e-138;
        bool r19709 = r19696 <= r19708;
        double r19710 = -r19696;
        double r19711 = r19696 * r19696;
        double r19712 = r19704 * r19700;
        double r19713 = r19711 - r19712;
        double r19714 = sqrt(r19713);
        double r19715 = r19710 + r19714;
        double r19716 = r19715 / r19704;
        double r19717 = -1/2;
        double r19718 = r19701 * r19717;
        double r19719 = r19709 ? r19716 : r19718;
        double r19720 = r19698 ? r19707 : r19719;
        return r19720;
}

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 r19721, r19722, r19723, r19724, r19725, r19726, r19727, r19728, r19729, r19730;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3472);
        mpfr_init(r19721);
        mpfr_init(r19722);
        mpfr_init(r19723);
        mpfr_init(r19724);
        mpfr_init(r19725);
        mpfr_init(r19726);
        mpfr_init(r19727);
        mpfr_init(r19728);
        mpfr_init(r19729);
        mpfr_init(r19730);
}

double f_im(double a, double b_2F2, double c) {
        mpfr_set_d(r19721, b_2F2, MPFR_RNDN);
        mpfr_neg(r19722, r19721, MPFR_RNDN);
        mpfr_mul(r19723, r19721, r19721, MPFR_RNDN);
        mpfr_set_d(r19724, a, MPFR_RNDN);
        mpfr_set_d(r19725, c, MPFR_RNDN);
        mpfr_mul(r19726, r19724, r19725, MPFR_RNDN);
        mpfr_sub(r19727, r19723, r19726, MPFR_RNDN);
        mpfr_sqrt(r19728, r19727, MPFR_RNDN);
        mpfr_add(r19729, r19722, r19728, MPFR_RNDN);
        mpfr_div(r19730, r19729, r19724, MPFR_RNDN);
        return mpfr_get_d(r19730, MPFR_RNDN);
}

static mpfr_t r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r19731);
        mpfr_init_set_str(r19732, "-2.788370916558726e+153", 10, MPFR_RNDN);
        mpfr_init(r19733);
        mpfr_init_set_str(r19734, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19735);
        mpfr_init(r19736);
        mpfr_init(r19737);
        mpfr_init_set_str(r19738, "2", 10, MPFR_RNDN);
        mpfr_init(r19739);
        mpfr_init(r19740);
        mpfr_init(r19741);
        mpfr_init(r19742);
        mpfr_init_set_str(r19743, "7.580412743766101e-138", 10, MPFR_RNDN);
        mpfr_init(r19744);
        mpfr_init(r19745);
        mpfr_init(r19746);
        mpfr_init(r19747);
        mpfr_init(r19748);
        mpfr_init(r19749);
        mpfr_init(r19750);
        mpfr_init(r19751);
        mpfr_init_set_str(r19752, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r19753);
        mpfr_init(r19754);
        mpfr_init(r19755);
}

double f_fm(double a, double b_2F2, double c) {
        mpfr_set_d(r19731, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r19733, mpfr_cmp(r19731, r19732) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19735, c, MPFR_RNDN);
        mpfr_div(r19736, r19735, r19731, MPFR_RNDN);
        mpfr_mul(r19737, r19734, r19736, MPFR_RNDN);
        ;
        mpfr_set_d(r19739, a, MPFR_RNDN);
        mpfr_div(r19740, r19731, r19739, MPFR_RNDN);
        mpfr_mul(r19741, r19738, r19740, MPFR_RNDN);
        mpfr_sub(r19742, r19737, r19741, MPFR_RNDN);
        ;
        mpfr_set_si(r19744, mpfr_cmp(r19731, r19743) <= 0, MPFR_RNDN);
        mpfr_neg(r19745, r19731, MPFR_RNDN);
        mpfr_mul(r19746, r19731, r19731, MPFR_RNDN);
        mpfr_mul(r19747, r19739, r19735, MPFR_RNDN);
        mpfr_sub(r19748, r19746, r19747, MPFR_RNDN);
        mpfr_sqrt(r19749, r19748, MPFR_RNDN);
        mpfr_add(r19750, r19745, r19749, MPFR_RNDN);
        mpfr_div(r19751, r19750, r19739, MPFR_RNDN);
        ;
        mpfr_mul(r19753, r19736, r19752, MPFR_RNDN);
        if (mpfr_get_si(r19744, MPFR_RNDN)) { mpfr_set(r19754, r19751, MPFR_RNDN); } else { mpfr_set(r19754, r19753, MPFR_RNDN); };
        if (mpfr_get_si(r19733, MPFR_RNDN)) { mpfr_set(r19755, r19742, MPFR_RNDN); } else { mpfr_set(r19755, r19754, MPFR_RNDN); };
        return mpfr_get_d(r19755, MPFR_RNDN);
}

static mpfr_t r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r19756);
        mpfr_init_set_str(r19757, "-2.788370916558726e+153", 10, MPFR_RNDN);
        mpfr_init(r19758);
        mpfr_init_set_str(r19759, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19760);
        mpfr_init(r19761);
        mpfr_init(r19762);
        mpfr_init_set_str(r19763, "2", 10, MPFR_RNDN);
        mpfr_init(r19764);
        mpfr_init(r19765);
        mpfr_init(r19766);
        mpfr_init(r19767);
        mpfr_init_set_str(r19768, "7.580412743766101e-138", 10, MPFR_RNDN);
        mpfr_init(r19769);
        mpfr_init(r19770);
        mpfr_init(r19771);
        mpfr_init(r19772);
        mpfr_init(r19773);
        mpfr_init(r19774);
        mpfr_init(r19775);
        mpfr_init(r19776);
        mpfr_init_set_str(r19777, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r19778);
        mpfr_init(r19779);
        mpfr_init(r19780);
}

double f_dm(double a, double b_2F2, double c) {
        mpfr_set_d(r19756, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r19758, mpfr_cmp(r19756, r19757) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19760, c, MPFR_RNDN);
        mpfr_div(r19761, r19760, r19756, MPFR_RNDN);
        mpfr_mul(r19762, r19759, r19761, MPFR_RNDN);
        ;
        mpfr_set_d(r19764, a, MPFR_RNDN);
        mpfr_div(r19765, r19756, r19764, MPFR_RNDN);
        mpfr_mul(r19766, r19763, r19765, MPFR_RNDN);
        mpfr_sub(r19767, r19762, r19766, MPFR_RNDN);
        ;
        mpfr_set_si(r19769, mpfr_cmp(r19756, r19768) <= 0, MPFR_RNDN);
        mpfr_neg(r19770, r19756, MPFR_RNDN);
        mpfr_mul(r19771, r19756, r19756, MPFR_RNDN);
        mpfr_mul(r19772, r19764, r19760, MPFR_RNDN);
        mpfr_sub(r19773, r19771, r19772, MPFR_RNDN);
        mpfr_sqrt(r19774, r19773, MPFR_RNDN);
        mpfr_add(r19775, r19770, r19774, MPFR_RNDN);
        mpfr_div(r19776, r19775, r19764, MPFR_RNDN);
        ;
        mpfr_mul(r19778, r19761, r19777, MPFR_RNDN);
        if (mpfr_get_si(r19769, MPFR_RNDN)) { mpfr_set(r19779, r19776, MPFR_RNDN); } else { mpfr_set(r19779, r19778, MPFR_RNDN); };
        if (mpfr_get_si(r19758, MPFR_RNDN)) { mpfr_set(r19780, r19767, MPFR_RNDN); } else { mpfr_set(r19780, r19779, MPFR_RNDN); };
        return mpfr_get_d(r19780, MPFR_RNDN);
}