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

char *name = "NMSE Section 6.1 mentioned, A";

double f_if(float x, float eps) {
        float r20893 = 1;
        float r20894 = eps;
        float r20895 = r20893 / r20894;
        float r20896 = r20893 + r20895;
        float r20897 = r20893 - r20894;
        float r20898 = x;
        float r20899 = r20897 * r20898;
        float r20900 = -r20899;
        float r20901 = exp(r20900);
        float r20902 = r20896 * r20901;
        float r20903 = r20895 - r20893;
        float r20904 = r20893 + r20894;
        float r20905 = r20904 * r20898;
        float r20906 = -r20905;
        float r20907 = exp(r20906);
        float r20908 = r20903 * r20907;
        float r20909 = r20902 - r20908;
        float r20910 = 2;
        float r20911 = r20909 / r20910;
        return r20911;
}

double f_id(double x, double eps) {
        double r20912 = 1;
        double r20913 = eps;
        double r20914 = r20912 / r20913;
        double r20915 = r20912 + r20914;
        double r20916 = r20912 - r20913;
        double r20917 = x;
        double r20918 = r20916 * r20917;
        double r20919 = -r20918;
        double r20920 = exp(r20919);
        double r20921 = r20915 * r20920;
        double r20922 = r20914 - r20912;
        double r20923 = r20912 + r20913;
        double r20924 = r20923 * r20917;
        double r20925 = -r20924;
        double r20926 = exp(r20925);
        double r20927 = r20922 * r20926;
        double r20928 = r20921 - r20927;
        double r20929 = 2;
        double r20930 = r20928 / r20929;
        return r20930;
}


double f_of(float x, float eps) {
        float r20931 = x;
        float r20932 = 35.794330612085844;
        bool r20933 = r20931 <= r20932;
        float r20934 = 2/3;
        float r20935 = r20931 * r20934;
        float r20936 = r20931 * r20931;
        float r20937 = 2;
        float r20938 = r20937 - r20936;
        float r20939 = fma(r20935, r20936, r20938);
        float r20940 = r20939 / r20937;
        float r20941 = 1;
        float r20942 = eps;
        float r20943 = r20941 / r20942;
        float r20944 = r20943 + r20941;
        float r20945 = exp(r20931);
        float r20946 = r20941 - r20942;
        float r20947 = pow(r20945, r20946);
        float r20948 = r20944 / r20947;
        float r20949 = r20943 - r20941;
        float r20950 = fma(r20942, r20931, r20931);
        float r20951 = exp(r20950);
        float r20952 = r20949 / r20951;
        float r20953 = r20948 - r20952;
        float r20954 = r20953 / r20937;
        float r20955 = r20933 ? r20940 : r20954;
        return r20955;
}

double f_od(double x, double eps) {
        double r20956 = x;
        double r20957 = 35.794330612085844;
        bool r20958 = r20956 <= r20957;
        double r20959 = 2/3;
        double r20960 = r20956 * r20959;
        double r20961 = r20956 * r20956;
        double r20962 = 2;
        double r20963 = r20962 - r20961;
        double r20964 = fma(r20960, r20961, r20963);
        double r20965 = r20964 / r20962;
        double r20966 = 1;
        double r20967 = eps;
        double r20968 = r20966 / r20967;
        double r20969 = r20968 + r20966;
        double r20970 = exp(r20956);
        double r20971 = r20966 - r20967;
        double r20972 = pow(r20970, r20971);
        double r20973 = r20969 / r20972;
        double r20974 = r20968 - r20966;
        double r20975 = fma(r20967, r20956, r20956);
        double r20976 = exp(r20975);
        double r20977 = r20974 / r20976;
        double r20978 = r20973 - r20977;
        double r20979 = r20978 / r20962;
        double r20980 = r20958 ? r20965 : r20979;
        return r20980;
}

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 r20981, r20982, r20983, r20984, r20985, r20986, r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994, r20995, r20996, r20997, r20998, r20999;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20981, "1", 10, MPFR_RNDN);
        mpfr_init(r20982);
        mpfr_init(r20983);
        mpfr_init(r20984);
        mpfr_init(r20985);
        mpfr_init(r20986);
        mpfr_init(r20987);
        mpfr_init(r20988);
        mpfr_init(r20989);
        mpfr_init(r20990);
        mpfr_init(r20991);
        mpfr_init(r20992);
        mpfr_init(r20993);
        mpfr_init(r20994);
        mpfr_init(r20995);
        mpfr_init(r20996);
        mpfr_init(r20997);
        mpfr_init_set_str(r20998, "2", 10, MPFR_RNDN);
        mpfr_init(r20999);
}

double f_im(double x, double eps) {
        ;
        mpfr_set_d(r20982, eps, MPFR_RNDN);
        mpfr_div(r20983, r20981, r20982, MPFR_RNDN);
        mpfr_add(r20984, r20981, r20983, MPFR_RNDN);
        mpfr_sub(r20985, r20981, r20982, MPFR_RNDN);
        mpfr_set_d(r20986, x, MPFR_RNDN);
        mpfr_mul(r20987, r20985, r20986, MPFR_RNDN);
        mpfr_neg(r20988, r20987, MPFR_RNDN);
        mpfr_exp(r20989, r20988, MPFR_RNDN);
        mpfr_mul(r20990, r20984, r20989, MPFR_RNDN);
        mpfr_sub(r20991, r20983, r20981, MPFR_RNDN);
        mpfr_add(r20992, r20981, r20982, MPFR_RNDN);
        mpfr_mul(r20993, r20992, r20986, MPFR_RNDN);
        mpfr_neg(r20994, r20993, MPFR_RNDN);
        mpfr_exp(r20995, r20994, MPFR_RNDN);
        mpfr_mul(r20996, r20991, r20995, MPFR_RNDN);
        mpfr_sub(r20997, r20990, r20996, MPFR_RNDN);
        ;
        mpfr_div(r20999, r20997, r20998, MPFR_RNDN);
        return mpfr_get_d(r20999, MPFR_RNDN);
}

static mpfr_t r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020, r21021, r21022, r21023, r21024;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r21000);
        mpfr_init_set_str(r21001, "35.794330612085844", 10, MPFR_RNDN);
        mpfr_init(r21002);
        mpfr_init_set_str(r21003, "2/3", 10, MPFR_RNDN);
        mpfr_init(r21004);
        mpfr_init(r21005);
        mpfr_init_set_str(r21006, "2", 10, MPFR_RNDN);
        mpfr_init(r21007);
        mpfr_init(r21008);
        mpfr_init(r21009);
        mpfr_init_set_str(r21010, "1", 10, MPFR_RNDN);
        mpfr_init(r21011);
        mpfr_init(r21012);
        mpfr_init(r21013);
        mpfr_init(r21014);
        mpfr_init(r21015);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
        mpfr_init(r21019);
        mpfr_init(r21020);
        mpfr_init(r21021);
        mpfr_init(r21022);
        mpfr_init(r21023);
        mpfr_init(r21024);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r21000, x, MPFR_RNDN);
        ;
        mpfr_set_si(r21002, mpfr_cmp(r21000, r21001) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r21004, r21000, r21003, MPFR_RNDN);
        mpfr_mul(r21005, r21000, r21000, MPFR_RNDN);
        ;
        mpfr_sub(r21007, r21006, r21005, MPFR_RNDN);
        mpfr_fma(r21008, r21004, r21005, r21007, MPFR_RNDN);
        mpfr_div(r21009, r21008, r21006, MPFR_RNDN);
        ;
        mpfr_set_d(r21011, eps, MPFR_RNDN);
        mpfr_div(r21012, r21010, r21011, MPFR_RNDN);
        mpfr_add(r21013, r21012, r21010, MPFR_RNDN);
        mpfr_exp(r21014, r21000, MPFR_RNDN);
        mpfr_sub(r21015, r21010, r21011, MPFR_RNDN);
        mpfr_pow(r21016, r21014, r21015, MPFR_RNDN);
        mpfr_div(r21017, r21013, r21016, MPFR_RNDN);
        mpfr_sub(r21018, r21012, r21010, MPFR_RNDN);
        mpfr_fma(r21019, r21011, r21000, r21000, MPFR_RNDN);
        mpfr_exp(r21020, r21019, MPFR_RNDN);
        mpfr_div(r21021, r21018, r21020, MPFR_RNDN);
        mpfr_sub(r21022, r21017, r21021, MPFR_RNDN);
        mpfr_div(r21023, r21022, r21006, MPFR_RNDN);
        if (mpfr_get_si(r21002, MPFR_RNDN)) { mpfr_set(r21024, r21009, MPFR_RNDN); } else { mpfr_set(r21024, r21023, MPFR_RNDN); };
        return mpfr_get_d(r21024, MPFR_RNDN);
}

static mpfr_t r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r21025);
        mpfr_init_set_str(r21026, "35.794330612085844", 10, MPFR_RNDN);
        mpfr_init(r21027);
        mpfr_init_set_str(r21028, "2/3", 10, MPFR_RNDN);
        mpfr_init(r21029);
        mpfr_init(r21030);
        mpfr_init_set_str(r21031, "2", 10, MPFR_RNDN);
        mpfr_init(r21032);
        mpfr_init(r21033);
        mpfr_init(r21034);
        mpfr_init_set_str(r21035, "1", 10, MPFR_RNDN);
        mpfr_init(r21036);
        mpfr_init(r21037);
        mpfr_init(r21038);
        mpfr_init(r21039);
        mpfr_init(r21040);
        mpfr_init(r21041);
        mpfr_init(r21042);
        mpfr_init(r21043);
        mpfr_init(r21044);
        mpfr_init(r21045);
        mpfr_init(r21046);
        mpfr_init(r21047);
        mpfr_init(r21048);
        mpfr_init(r21049);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r21025, x, MPFR_RNDN);
        ;
        mpfr_set_si(r21027, mpfr_cmp(r21025, r21026) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r21029, r21025, r21028, MPFR_RNDN);
        mpfr_mul(r21030, r21025, r21025, MPFR_RNDN);
        ;
        mpfr_sub(r21032, r21031, r21030, MPFR_RNDN);
        mpfr_fma(r21033, r21029, r21030, r21032, MPFR_RNDN);
        mpfr_div(r21034, r21033, r21031, MPFR_RNDN);
        ;
        mpfr_set_d(r21036, eps, MPFR_RNDN);
        mpfr_div(r21037, r21035, r21036, MPFR_RNDN);
        mpfr_add(r21038, r21037, r21035, MPFR_RNDN);
        mpfr_exp(r21039, r21025, MPFR_RNDN);
        mpfr_sub(r21040, r21035, r21036, MPFR_RNDN);
        mpfr_pow(r21041, r21039, r21040, MPFR_RNDN);
        mpfr_div(r21042, r21038, r21041, MPFR_RNDN);
        mpfr_sub(r21043, r21037, r21035, MPFR_RNDN);
        mpfr_fma(r21044, r21036, r21025, r21025, MPFR_RNDN);
        mpfr_exp(r21045, r21044, MPFR_RNDN);
        mpfr_div(r21046, r21043, r21045, MPFR_RNDN);
        mpfr_sub(r21047, r21042, r21046, MPFR_RNDN);
        mpfr_div(r21048, r21047, r21031, MPFR_RNDN);
        if (mpfr_get_si(r21027, MPFR_RNDN)) { mpfr_set(r21049, r21034, MPFR_RNDN); } else { mpfr_set(r21049, r21048, MPFR_RNDN); };
        return mpfr_get_d(r21049, MPFR_RNDN);
}