#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 = -7.833327672762858e+152;
        bool r19673 = r19671 <= r19672;
        float r19674 = -2;
        float r19675 = a;
        float r19676 = r19671 / r19675;
        float r19677 = r19674 * r19676;
        float r19678 = 7.698366662298276e-138;
        bool r19679 = r19671 <= r19678;
        float r19680 = -r19671;
        float r19681 = r19671 * r19671;
        float r19682 = c;
        float r19683 = r19675 * r19682;
        float r19684 = r19681 - r19683;
        float r19685 = sqrt(r19684);
        float r19686 = r19680 + r19685;
        float r19687 = 1;
        float r19688 = r19687 / r19675;
        float r19689 = r19686 * r19688;
        float r19690 = r19682 / r19671;
        float r19691 = -1/2;
        float r19692 = r19690 * r19691;
        float r19693 = r19679 ? r19689 : r19692;
        float r19694 = r19673 ? r19677 : r19693;
        return r19694;
}

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

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

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2960);
        mpfr_init(r19719);
        mpfr_init(r19720);
        mpfr_init(r19721);
        mpfr_init(r19722);
        mpfr_init(r19723);
        mpfr_init(r19724);
        mpfr_init(r19725);
        mpfr_init(r19726);
        mpfr_init(r19727);
        mpfr_init(r19728);
}

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

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

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2960);
        mpfr_init(r19729);
        mpfr_init_set_str(r19730, "-7.833327672762858e+152", 10, MPFR_RNDN);
        mpfr_init(r19731);
        mpfr_init_set_str(r19732, "-2", 10, MPFR_RNDN);
        mpfr_init(r19733);
        mpfr_init(r19734);
        mpfr_init(r19735);
        mpfr_init_set_str(r19736, "7.698366662298276e-138", 10, MPFR_RNDN);
        mpfr_init(r19737);
        mpfr_init(r19738);
        mpfr_init(r19739);
        mpfr_init(r19740);
        mpfr_init(r19741);
        mpfr_init(r19742);
        mpfr_init(r19743);
        mpfr_init(r19744);
        mpfr_init_set_str(r19745, "1", 10, MPFR_RNDN);
        mpfr_init(r19746);
        mpfr_init(r19747);
        mpfr_init(r19748);
        mpfr_init_set_str(r19749, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r19750);
        mpfr_init(r19751);
        mpfr_init(r19752);
}

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

static mpfr_t r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2960);
        mpfr_init(r19753);
        mpfr_init_set_str(r19754, "-7.833327672762858e+152", 10, MPFR_RNDN);
        mpfr_init(r19755);
        mpfr_init_set_str(r19756, "-2", 10, MPFR_RNDN);
        mpfr_init(r19757);
        mpfr_init(r19758);
        mpfr_init(r19759);
        mpfr_init_set_str(r19760, "7.698366662298276e-138", 10, MPFR_RNDN);
        mpfr_init(r19761);
        mpfr_init(r19762);
        mpfr_init(r19763);
        mpfr_init(r19764);
        mpfr_init(r19765);
        mpfr_init(r19766);
        mpfr_init(r19767);
        mpfr_init(r19768);
        mpfr_init_set_str(r19769, "1", 10, MPFR_RNDN);
        mpfr_init(r19770);
        mpfr_init(r19771);
        mpfr_init(r19772);
        mpfr_init_set_str(r19773, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r19774);
        mpfr_init(r19775);
        mpfr_init(r19776);
}

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