#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 r19741 = b_2F2;
        float r19742 = -r19741;
        float r19743 = r19741 * r19741;
        float r19744 = a;
        float r19745 = c;
        float r19746 = r19744 * r19745;
        float r19747 = r19743 - r19746;
        float r19748 = sqrt(r19747);
        float r19749 = r19742 + r19748;
        float r19750 = r19749 / r19744;
        return r19750;
}

double f_id(double a, double b_2F2, double c) {
        double r19751 = b_2F2;
        double r19752 = -r19751;
        double r19753 = r19751 * r19751;
        double r19754 = a;
        double r19755 = c;
        double r19756 = r19754 * r19755;
        double r19757 = r19753 - r19756;
        double r19758 = sqrt(r19757);
        double r19759 = r19752 + r19758;
        double r19760 = r19759 / r19754;
        return r19760;
}


double f_of(float a, float b_2F2, float c) {
        float r19761 = b_2F2;
        float r19762 = -5.190855692504989e+126;
        bool r19763 = r19761 <= r19762;
        float r19764 = -2;
        float r19765 = a;
        float r19766 = r19761 / r19765;
        float r19767 = r19764 * r19766;
        float r19768 = 9.165338891167788e-118;
        bool r19769 = r19761 <= r19768;
        float r19770 = 1;
        float r19771 = -r19761;
        float r19772 = r19761 * r19761;
        float r19773 = c;
        float r19774 = r19765 * r19773;
        float r19775 = r19772 - r19774;
        float r19776 = sqrt(r19775);
        float r19777 = r19771 + r19776;
        float r19778 = r19765 / r19777;
        float r19779 = r19770 / r19778;
        float r19780 = r19771 - r19761;
        float r19781 = 1/2;
        float r19782 = r19781 * r19773;
        float r19783 = r19782 / r19766;
        float r19784 = r19780 + r19783;
        float r19785 = r19773 / r19784;
        float r19786 = r19769 ? r19779 : r19785;
        float r19787 = r19763 ? r19767 : r19786;
        return r19787;
}

double f_od(double a, double b_2F2, double c) {
        double r19788 = b_2F2;
        double r19789 = -5.190855692504989e+126;
        bool r19790 = r19788 <= r19789;
        double r19791 = -2;
        double r19792 = a;
        double r19793 = r19788 / r19792;
        double r19794 = r19791 * r19793;
        double r19795 = 9.165338891167788e-118;
        bool r19796 = r19788 <= r19795;
        double r19797 = 1;
        double r19798 = -r19788;
        double r19799 = r19788 * r19788;
        double r19800 = c;
        double r19801 = r19792 * r19800;
        double r19802 = r19799 - r19801;
        double r19803 = sqrt(r19802);
        double r19804 = r19798 + r19803;
        double r19805 = r19792 / r19804;
        double r19806 = r19797 / r19805;
        double r19807 = r19798 - r19788;
        double r19808 = 1/2;
        double r19809 = r19808 * r19800;
        double r19810 = r19809 / r19793;
        double r19811 = r19807 + r19810;
        double r19812 = r19800 / r19811;
        double r19813 = r19796 ? r19806 : r19812;
        double r19814 = r19790 ? r19794 : r19813;
        return r19814;
}

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 r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19815);
        mpfr_init(r19816);
        mpfr_init(r19817);
        mpfr_init(r19818);
        mpfr_init(r19819);
        mpfr_init(r19820);
        mpfr_init(r19821);
        mpfr_init(r19822);
        mpfr_init(r19823);
        mpfr_init(r19824);
}

double f_im(double a, double b_2F2, double c) {
        mpfr_set_d(r19815, b_2F2, MPFR_RNDN);
        mpfr_neg(r19816, r19815, MPFR_RNDN);
        mpfr_mul(r19817, r19815, r19815, MPFR_RNDN);
        mpfr_set_d(r19818, a, MPFR_RNDN);
        mpfr_set_d(r19819, c, MPFR_RNDN);
        mpfr_mul(r19820, r19818, r19819, MPFR_RNDN);
        mpfr_sub(r19821, r19817, r19820, MPFR_RNDN);
        mpfr_sqrt(r19822, r19821, MPFR_RNDN);
        mpfr_add(r19823, r19816, r19822, MPFR_RNDN);
        mpfr_div(r19824, r19823, r19818, MPFR_RNDN);
        return mpfr_get_d(r19824, MPFR_RNDN);
}

static mpfr_t r19825, r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19825);
        mpfr_init_set_str(r19826, "-5.190855692504989e+126", 10, MPFR_RNDN);
        mpfr_init(r19827);
        mpfr_init_set_str(r19828, "-2", 10, MPFR_RNDN);
        mpfr_init(r19829);
        mpfr_init(r19830);
        mpfr_init(r19831);
        mpfr_init_set_str(r19832, "9.165338891167788e-118", 10, MPFR_RNDN);
        mpfr_init(r19833);
        mpfr_init_set_str(r19834, "1", 10, MPFR_RNDN);
        mpfr_init(r19835);
        mpfr_init(r19836);
        mpfr_init(r19837);
        mpfr_init(r19838);
        mpfr_init(r19839);
        mpfr_init(r19840);
        mpfr_init(r19841);
        mpfr_init(r19842);
        mpfr_init(r19843);
        mpfr_init(r19844);
        mpfr_init_set_str(r19845, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19846);
        mpfr_init(r19847);
        mpfr_init(r19848);
        mpfr_init(r19849);
        mpfr_init(r19850);
        mpfr_init(r19851);
}

double f_fm(double a, double b_2F2, double c) {
        mpfr_set_d(r19825, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r19827, mpfr_cmp(r19825, r19826) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19829, a, MPFR_RNDN);
        mpfr_div(r19830, r19825, r19829, MPFR_RNDN);
        mpfr_mul(r19831, r19828, r19830, MPFR_RNDN);
        ;
        mpfr_set_si(r19833, mpfr_cmp(r19825, r19832) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r19835, r19825, MPFR_RNDN);
        mpfr_mul(r19836, r19825, r19825, MPFR_RNDN);
        mpfr_set_d(r19837, c, MPFR_RNDN);
        mpfr_mul(r19838, r19829, r19837, MPFR_RNDN);
        mpfr_sub(r19839, r19836, r19838, MPFR_RNDN);
        mpfr_sqrt(r19840, r19839, MPFR_RNDN);
        mpfr_add(r19841, r19835, r19840, MPFR_RNDN);
        mpfr_div(r19842, r19829, r19841, MPFR_RNDN);
        mpfr_div(r19843, r19834, r19842, MPFR_RNDN);
        mpfr_sub(r19844, r19835, r19825, MPFR_RNDN);
        ;
        mpfr_mul(r19846, r19845, r19837, MPFR_RNDN);
        mpfr_div(r19847, r19846, r19830, MPFR_RNDN);
        mpfr_add(r19848, r19844, r19847, MPFR_RNDN);
        mpfr_div(r19849, r19837, r19848, MPFR_RNDN);
        if (mpfr_get_si(r19833, MPFR_RNDN)) { mpfr_set(r19850, r19843, MPFR_RNDN); } else { mpfr_set(r19850, r19849, MPFR_RNDN); };
        if (mpfr_get_si(r19827, MPFR_RNDN)) { mpfr_set(r19851, r19831, MPFR_RNDN); } else { mpfr_set(r19851, r19850, MPFR_RNDN); };
        return mpfr_get_d(r19851, MPFR_RNDN);
}

static mpfr_t r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866, r19867, r19868, r19869, r19870, r19871, r19872, r19873, r19874, r19875, r19876, r19877, r19878;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19852);
        mpfr_init_set_str(r19853, "-5.190855692504989e+126", 10, MPFR_RNDN);
        mpfr_init(r19854);
        mpfr_init_set_str(r19855, "-2", 10, MPFR_RNDN);
        mpfr_init(r19856);
        mpfr_init(r19857);
        mpfr_init(r19858);
        mpfr_init_set_str(r19859, "9.165338891167788e-118", 10, MPFR_RNDN);
        mpfr_init(r19860);
        mpfr_init_set_str(r19861, "1", 10, MPFR_RNDN);
        mpfr_init(r19862);
        mpfr_init(r19863);
        mpfr_init(r19864);
        mpfr_init(r19865);
        mpfr_init(r19866);
        mpfr_init(r19867);
        mpfr_init(r19868);
        mpfr_init(r19869);
        mpfr_init(r19870);
        mpfr_init(r19871);
        mpfr_init_set_str(r19872, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19873);
        mpfr_init(r19874);
        mpfr_init(r19875);
        mpfr_init(r19876);
        mpfr_init(r19877);
        mpfr_init(r19878);
}

double f_dm(double a, double b_2F2, double c) {
        mpfr_set_d(r19852, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r19854, mpfr_cmp(r19852, r19853) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r19856, a, MPFR_RNDN);
        mpfr_div(r19857, r19852, r19856, MPFR_RNDN);
        mpfr_mul(r19858, r19855, r19857, MPFR_RNDN);
        ;
        mpfr_set_si(r19860, mpfr_cmp(r19852, r19859) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r19862, r19852, MPFR_RNDN);
        mpfr_mul(r19863, r19852, r19852, MPFR_RNDN);
        mpfr_set_d(r19864, c, MPFR_RNDN);
        mpfr_mul(r19865, r19856, r19864, MPFR_RNDN);
        mpfr_sub(r19866, r19863, r19865, MPFR_RNDN);
        mpfr_sqrt(r19867, r19866, MPFR_RNDN);
        mpfr_add(r19868, r19862, r19867, MPFR_RNDN);
        mpfr_div(r19869, r19856, r19868, MPFR_RNDN);
        mpfr_div(r19870, r19861, r19869, MPFR_RNDN);
        mpfr_sub(r19871, r19862, r19852, MPFR_RNDN);
        ;
        mpfr_mul(r19873, r19872, r19864, MPFR_RNDN);
        mpfr_div(r19874, r19873, r19857, MPFR_RNDN);
        mpfr_add(r19875, r19871, r19874, MPFR_RNDN);
        mpfr_div(r19876, r19864, r19875, MPFR_RNDN);
        if (mpfr_get_si(r19860, MPFR_RNDN)) { mpfr_set(r19877, r19870, MPFR_RNDN); } else { mpfr_set(r19877, r19876, MPFR_RNDN); };
        if (mpfr_get_si(r19854, MPFR_RNDN)) { mpfr_set(r19878, r19858, MPFR_RNDN); } else { mpfr_set(r19878, r19877, MPFR_RNDN); };
        return mpfr_get_d(r19878, MPFR_RNDN);
}