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

char *name = "quad2m (problem 3.2.1, negative)";

double f_if(float a, float b_2F2, float c) {
        float r20798 = b_2F2;
        float r20799 = -r20798;
        float r20800 = r20798 * r20798;
        float r20801 = a;
        float r20802 = c;
        float r20803 = r20801 * r20802;
        float r20804 = r20800 - r20803;
        float r20805 = sqrt(r20804);
        float r20806 = r20799 - r20805;
        float r20807 = r20806 / r20801;
        return r20807;
}

double f_id(double a, double b_2F2, double c) {
        double r20808 = b_2F2;
        double r20809 = -r20808;
        double r20810 = r20808 * r20808;
        double r20811 = a;
        double r20812 = c;
        double r20813 = r20811 * r20812;
        double r20814 = r20810 - r20813;
        double r20815 = sqrt(r20814);
        double r20816 = r20809 - r20815;
        double r20817 = r20816 / r20811;
        return r20817;
}


double f_of(float a, float b_2F2, float c) {
        float r20818 = b_2F2;
        float r20819 = -4.005537006288397e-85;
        bool r20820 = r20818 <= r20819;
        float r20821 = c;
        float r20822 = r20821 / r20818;
        float r20823 = -1/2;
        float r20824 = r20822 * r20823;
        float r20825 = 1.026219862232508e+81;
        bool r20826 = r20818 <= r20825;
        float r20827 = 1;
        float r20828 = a;
        float r20829 = -r20818;
        float r20830 = r20818 * r20818;
        float r20831 = r20828 * r20821;
        float r20832 = r20830 - r20831;
        float r20833 = sqrt(r20832);
        float r20834 = r20829 - r20833;
        float r20835 = r20828 / r20834;
        float r20836 = r20827 / r20835;
        float r20837 = 1/2;
        float r20838 = r20837 * r20822;
        float r20839 = r20829 - r20818;
        float r20840 = r20839 / r20828;
        float r20841 = r20838 + r20840;
        float r20842 = r20826 ? r20836 : r20841;
        float r20843 = r20820 ? r20824 : r20842;
        return r20843;
}

double f_od(double a, double b_2F2, double c) {
        double r20844 = b_2F2;
        double r20845 = -4.005537006288397e-85;
        bool r20846 = r20844 <= r20845;
        double r20847 = c;
        double r20848 = r20847 / r20844;
        double r20849 = -1/2;
        double r20850 = r20848 * r20849;
        double r20851 = 1.026219862232508e+81;
        bool r20852 = r20844 <= r20851;
        double r20853 = 1;
        double r20854 = a;
        double r20855 = -r20844;
        double r20856 = r20844 * r20844;
        double r20857 = r20854 * r20847;
        double r20858 = r20856 - r20857;
        double r20859 = sqrt(r20858);
        double r20860 = r20855 - r20859;
        double r20861 = r20854 / r20860;
        double r20862 = r20853 / r20861;
        double r20863 = 1/2;
        double r20864 = r20863 * r20848;
        double r20865 = r20855 - r20844;
        double r20866 = r20865 / r20854;
        double r20867 = r20864 + r20866;
        double r20868 = r20852 ? r20862 : r20867;
        double r20869 = r20846 ? r20850 : r20868;
        return r20869;
}

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 r20870, r20871, r20872, r20873, r20874, r20875, r20876, r20877, r20878, r20879;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3472);
        mpfr_init(r20870);
        mpfr_init(r20871);
        mpfr_init(r20872);
        mpfr_init(r20873);
        mpfr_init(r20874);
        mpfr_init(r20875);
        mpfr_init(r20876);
        mpfr_init(r20877);
        mpfr_init(r20878);
        mpfr_init(r20879);
}

double f_im(double a, double b_2F2, double c) {
        mpfr_set_d(r20870, b_2F2, MPFR_RNDN);
        mpfr_neg(r20871, r20870, MPFR_RNDN);
        mpfr_mul(r20872, r20870, r20870, MPFR_RNDN);
        mpfr_set_d(r20873, a, MPFR_RNDN);
        mpfr_set_d(r20874, c, MPFR_RNDN);
        mpfr_mul(r20875, r20873, r20874, MPFR_RNDN);
        mpfr_sub(r20876, r20872, r20875, MPFR_RNDN);
        mpfr_sqrt(r20877, r20876, MPFR_RNDN);
        mpfr_sub(r20878, r20871, r20877, MPFR_RNDN);
        mpfr_div(r20879, r20878, r20873, MPFR_RNDN);
        return mpfr_get_d(r20879, MPFR_RNDN);
}

static mpfr_t r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903, r20904, r20905;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r20880);
        mpfr_init_set_str(r20881, "-4.005537006288397e-85", 10, MPFR_RNDN);
        mpfr_init(r20882);
        mpfr_init(r20883);
        mpfr_init(r20884);
        mpfr_init_set_str(r20885, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r20886);
        mpfr_init_set_str(r20887, "1.026219862232508e+81", 10, MPFR_RNDN);
        mpfr_init(r20888);
        mpfr_init_set_str(r20889, "1", 10, MPFR_RNDN);
        mpfr_init(r20890);
        mpfr_init(r20891);
        mpfr_init(r20892);
        mpfr_init(r20893);
        mpfr_init(r20894);
        mpfr_init(r20895);
        mpfr_init(r20896);
        mpfr_init(r20897);
        mpfr_init(r20898);
        mpfr_init_set_str(r20899, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20900);
        mpfr_init(r20901);
        mpfr_init(r20902);
        mpfr_init(r20903);
        mpfr_init(r20904);
        mpfr_init(r20905);
}

double f_fm(double a, double b_2F2, double c) {
        mpfr_set_d(r20880, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r20882, mpfr_cmp(r20880, r20881) <= 0, MPFR_RNDN);
        mpfr_set_d(r20883, c, MPFR_RNDN);
        mpfr_div(r20884, r20883, r20880, MPFR_RNDN);
        ;
        mpfr_mul(r20886, r20884, r20885, MPFR_RNDN);
        ;
        mpfr_set_si(r20888, mpfr_cmp(r20880, r20887) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r20890, a, MPFR_RNDN);
        mpfr_neg(r20891, r20880, MPFR_RNDN);
        mpfr_mul(r20892, r20880, r20880, MPFR_RNDN);
        mpfr_mul(r20893, r20890, r20883, MPFR_RNDN);
        mpfr_sub(r20894, r20892, r20893, MPFR_RNDN);
        mpfr_sqrt(r20895, r20894, MPFR_RNDN);
        mpfr_sub(r20896, r20891, r20895, MPFR_RNDN);
        mpfr_div(r20897, r20890, r20896, MPFR_RNDN);
        mpfr_div(r20898, r20889, r20897, MPFR_RNDN);
        ;
        mpfr_mul(r20900, r20899, r20884, MPFR_RNDN);
        mpfr_sub(r20901, r20891, r20880, MPFR_RNDN);
        mpfr_div(r20902, r20901, r20890, MPFR_RNDN);
        mpfr_add(r20903, r20900, r20902, MPFR_RNDN);
        if (mpfr_get_si(r20888, MPFR_RNDN)) { mpfr_set(r20904, r20898, MPFR_RNDN); } else { mpfr_set(r20904, r20903, MPFR_RNDN); };
        if (mpfr_get_si(r20882, MPFR_RNDN)) { mpfr_set(r20905, r20886, MPFR_RNDN); } else { mpfr_set(r20905, r20904, MPFR_RNDN); };
        return mpfr_get_d(r20905, MPFR_RNDN);
}

static mpfr_t r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930, r20931;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3472);
        mpfr_init(r20906);
        mpfr_init_set_str(r20907, "-4.005537006288397e-85", 10, MPFR_RNDN);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init(r20910);
        mpfr_init_set_str(r20911, "-1/2", 10, MPFR_RNDN);
        mpfr_init(r20912);
        mpfr_init_set_str(r20913, "1.026219862232508e+81", 10, MPFR_RNDN);
        mpfr_init(r20914);
        mpfr_init_set_str(r20915, "1", 10, MPFR_RNDN);
        mpfr_init(r20916);
        mpfr_init(r20917);
        mpfr_init(r20918);
        mpfr_init(r20919);
        mpfr_init(r20920);
        mpfr_init(r20921);
        mpfr_init(r20922);
        mpfr_init(r20923);
        mpfr_init(r20924);
        mpfr_init_set_str(r20925, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20926);
        mpfr_init(r20927);
        mpfr_init(r20928);
        mpfr_init(r20929);
        mpfr_init(r20930);
        mpfr_init(r20931);
}

double f_dm(double a, double b_2F2, double c) {
        mpfr_set_d(r20906, b_2F2, MPFR_RNDN);
        ;
        mpfr_set_si(r20908, mpfr_cmp(r20906, r20907) <= 0, MPFR_RNDN);
        mpfr_set_d(r20909, c, MPFR_RNDN);
        mpfr_div(r20910, r20909, r20906, MPFR_RNDN);
        ;
        mpfr_mul(r20912, r20910, r20911, MPFR_RNDN);
        ;
        mpfr_set_si(r20914, mpfr_cmp(r20906, r20913) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r20916, a, MPFR_RNDN);
        mpfr_neg(r20917, r20906, MPFR_RNDN);
        mpfr_mul(r20918, r20906, r20906, MPFR_RNDN);
        mpfr_mul(r20919, r20916, r20909, MPFR_RNDN);
        mpfr_sub(r20920, r20918, r20919, MPFR_RNDN);
        mpfr_sqrt(r20921, r20920, MPFR_RNDN);
        mpfr_sub(r20922, r20917, r20921, MPFR_RNDN);
        mpfr_div(r20923, r20916, r20922, MPFR_RNDN);
        mpfr_div(r20924, r20915, r20923, MPFR_RNDN);
        ;
        mpfr_mul(r20926, r20925, r20910, MPFR_RNDN);
        mpfr_sub(r20927, r20917, r20906, MPFR_RNDN);
        mpfr_div(r20928, r20927, r20916, MPFR_RNDN);
        mpfr_add(r20929, r20926, r20928, MPFR_RNDN);
        if (mpfr_get_si(r20914, MPFR_RNDN)) { mpfr_set(r20930, r20924, MPFR_RNDN); } else { mpfr_set(r20930, r20929, MPFR_RNDN); };
        if (mpfr_get_si(r20908, MPFR_RNDN)) { mpfr_set(r20931, r20912, MPFR_RNDN); } else { mpfr_set(r20931, r20930, MPFR_RNDN); };
        return mpfr_get_d(r20931, MPFR_RNDN);
}