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

char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5";

double f_if(float x) {
        float r21883 = 1.0f;
        float r21884 = atan2(1.0, 0.0);
        float r21885 = sqrt(r21884);
        float r21886 = r21883 / r21885;
        float r21887 = 2.0f;
        float r21888 = x;
        float r21889 = fabs(r21888);
        float r21890 = r21887 * r21889;
        float r21891 = 3.0f;
        float r21892 = r21887 / r21891;
        float r21893 = r21889 * r21889;
        float r21894 = r21893 * r21889;
        float r21895 = r21892 * r21894;
        float r21896 = r21890 + r21895;
        float r21897 = 5.0f;
        float r21898 = r21883 / r21897;
        float r21899 = r21894 * r21889;
        float r21900 = r21899 * r21889;
        float r21901 = r21898 * r21900;
        float r21902 = r21896 + r21901;
        float r21903 = 21.0f;
        float r21904 = r21883 / r21903;
        float r21905 = r21900 * r21889;
        float r21906 = r21905 * r21889;
        float r21907 = r21904 * r21906;
        float r21908 = r21902 + r21907;
        float r21909 = r21886 * r21908;
        float r21910 = fabs(r21909);
        return r21910;
}

double f_id(double x) {
        double r21911 = 1.0;
        double r21912 = atan2(1.0, 0.0);
        double r21913 = sqrt(r21912);
        double r21914 = r21911 / r21913;
        double r21915 = 2.0;
        double r21916 = x;
        double r21917 = fabs(r21916);
        double r21918 = r21915 * r21917;
        double r21919 = 3.0;
        double r21920 = r21915 / r21919;
        double r21921 = r21917 * r21917;
        double r21922 = r21921 * r21917;
        double r21923 = r21920 * r21922;
        double r21924 = r21918 + r21923;
        double r21925 = 5.0;
        double r21926 = r21911 / r21925;
        double r21927 = r21922 * r21917;
        double r21928 = r21927 * r21917;
        double r21929 = r21926 * r21928;
        double r21930 = r21924 + r21929;
        double r21931 = 21.0;
        double r21932 = r21911 / r21931;
        double r21933 = r21928 * r21917;
        double r21934 = r21933 * r21917;
        double r21935 = r21932 * r21934;
        double r21936 = r21930 + r21935;
        double r21937 = r21914 * r21936;
        double r21938 = fabs(r21937);
        return r21938;
}


double f_of(float x) {
        float r21939 = 0.6666666666666666f;
        float r21940 = x;
        float r21941 = fabs(r21940);
        float r21942 = 3.0f;
        float r21943 = pow(r21941, r21942);
        float r21944 = r21939 * r21943;
        float r21945 = 2.0f;
        float r21946 = r21945 * r21941;
        float r21947 = 0.047619047619047616f;
        float r21948 = 7.0f;
        float r21949 = pow(r21941, r21948);
        float r21950 = r21947 * r21949;
        float r21951 = 0.2f;
        float r21952 = 5.0f;
        float r21953 = pow(r21941, r21952);
        float r21954 = r21951 * r21953;
        float r21955 = r21950 + r21954;
        float r21956 = r21946 + r21955;
        float r21957 = r21944 + r21956;
        float r21958 = 1.0f;
        float r21959 = atan2(1.0, 0.0);
        float r21960 = r21958 / r21959;
        float r21961 = sqrt(r21960);
        float r21962 = r21957 * r21961;
        float r21963 = fabs(r21962);
        return r21963;
}

double f_od(double x) {
        double r21964 = 0.6666666666666666;
        double r21965 = x;
        double r21966 = fabs(r21965);
        double r21967 = 3.0;
        double r21968 = pow(r21966, r21967);
        double r21969 = r21964 * r21968;
        double r21970 = 2.0;
        double r21971 = r21970 * r21966;
        double r21972 = 0.047619047619047616;
        double r21973 = 7.0;
        double r21974 = pow(r21966, r21973);
        double r21975 = r21972 * r21974;
        double r21976 = 0.2;
        double r21977 = 5.0;
        double r21978 = pow(r21966, r21977);
        double r21979 = r21976 * r21978;
        double r21980 = r21975 + r21979;
        double r21981 = r21971 + r21980;
        double r21982 = r21969 + r21981;
        double r21983 = 1.0;
        double r21984 = atan2(1.0, 0.0);
        double r21985 = r21983 / r21984;
        double r21986 = sqrt(r21985);
        double r21987 = r21982 * r21986;
        double r21988 = fabs(r21987);
        return r21988;
}

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 r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000, r22001, r22002, r22003, r22004, r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r21989, "1", 10, MPFR_RNDN);
        mpfr_init(r21990);
        mpfr_init(r21991);
        mpfr_init(r21992);
        mpfr_init_set_str(r21993, "2", 10, MPFR_RNDN);
        mpfr_init(r21994);
        mpfr_init(r21995);
        mpfr_init(r21996);
        mpfr_init_set_str(r21997, "3", 10, MPFR_RNDN);
        mpfr_init(r21998);
        mpfr_init(r21999);
        mpfr_init(r22000);
        mpfr_init(r22001);
        mpfr_init(r22002);
        mpfr_init_set_str(r22003, "5", 10, MPFR_RNDN);
        mpfr_init(r22004);
        mpfr_init(r22005);
        mpfr_init(r22006);
        mpfr_init(r22007);
        mpfr_init(r22008);
        mpfr_init_set_str(r22009, "21", 10, MPFR_RNDN);
        mpfr_init(r22010);
        mpfr_init(r22011);
        mpfr_init(r22012);
        mpfr_init(r22013);
        mpfr_init(r22014);
        mpfr_init(r22015);
        mpfr_init(r22016);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r21990, MPFR_RNDN);
        mpfr_sqrt(r21991, r21990, MPFR_RNDN);
        mpfr_div(r21992, r21989, r21991, MPFR_RNDN);
        ;
        mpfr_set_d(r21994, x, MPFR_RNDN);
        mpfr_abs(r21995, r21994, MPFR_RNDN);
        mpfr_mul(r21996, r21993, r21995, MPFR_RNDN);
        ;
        mpfr_div(r21998, r21993, r21997, MPFR_RNDN);
        mpfr_mul(r21999, r21995, r21995, MPFR_RNDN);
        mpfr_mul(r22000, r21999, r21995, MPFR_RNDN);
        mpfr_mul(r22001, r21998, r22000, MPFR_RNDN);
        mpfr_add(r22002, r21996, r22001, MPFR_RNDN);
        ;
        mpfr_div(r22004, r21989, r22003, MPFR_RNDN);
        mpfr_mul(r22005, r22000, r21995, MPFR_RNDN);
        mpfr_mul(r22006, r22005, r21995, MPFR_RNDN);
        mpfr_mul(r22007, r22004, r22006, MPFR_RNDN);
        mpfr_add(r22008, r22002, r22007, MPFR_RNDN);
        ;
        mpfr_div(r22010, r21989, r22009, MPFR_RNDN);
        mpfr_mul(r22011, r22006, r21995, MPFR_RNDN);
        mpfr_mul(r22012, r22011, r21995, MPFR_RNDN);
        mpfr_mul(r22013, r22010, r22012, MPFR_RNDN);
        mpfr_add(r22014, r22008, r22013, MPFR_RNDN);
        mpfr_mul(r22015, r21992, r22014, MPFR_RNDN);
        mpfr_abs(r22016, r22015, MPFR_RNDN);
        return mpfr_get_d(r22016, MPFR_RNDN);
}

static mpfr_t r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039, r22040, r22041;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22017, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22018);
        mpfr_init(r22019);
        mpfr_init_set_str(r22020, "3", 10, MPFR_RNDN);
        mpfr_init(r22021);
        mpfr_init(r22022);
        mpfr_init_set_str(r22023, "2", 10, MPFR_RNDN);
        mpfr_init(r22024);
        mpfr_init_set_str(r22025, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r22026, "7", 10, MPFR_RNDN);
        mpfr_init(r22027);
        mpfr_init(r22028);
        mpfr_init_set_str(r22029, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r22030, "5", 10, MPFR_RNDN);
        mpfr_init(r22031);
        mpfr_init(r22032);
        mpfr_init(r22033);
        mpfr_init(r22034);
        mpfr_init(r22035);
        mpfr_init_set_str(r22036, "1", 10, MPFR_RNDN);
        mpfr_init(r22037);
        mpfr_init(r22038);
        mpfr_init(r22039);
        mpfr_init(r22040);
        mpfr_init(r22041);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r22018, x, MPFR_RNDN);
        mpfr_abs(r22019, r22018, MPFR_RNDN);
        ;
        mpfr_pow(r22021, r22019, r22020, MPFR_RNDN);
        mpfr_mul(r22022, r22017, r22021, MPFR_RNDN);
        ;
        mpfr_mul(r22024, r22023, r22019, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22027, r22019, r22026, MPFR_RNDN);
        mpfr_mul(r22028, r22025, r22027, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22031, r22019, r22030, MPFR_RNDN);
        mpfr_mul(r22032, r22029, r22031, MPFR_RNDN);
        mpfr_add(r22033, r22028, r22032, MPFR_RNDN);
        mpfr_add(r22034, r22024, r22033, MPFR_RNDN);
        mpfr_add(r22035, r22022, r22034, MPFR_RNDN);
        ;
        mpfr_const_pi(r22037, MPFR_RNDN);
        mpfr_div(r22038, r22036, r22037, MPFR_RNDN);
        mpfr_sqrt(r22039, r22038, MPFR_RNDN);
        mpfr_mul(r22040, r22035, r22039, MPFR_RNDN);
        mpfr_abs(r22041, r22040, MPFR_RNDN);
        return mpfr_get_d(r22041, MPFR_RNDN);
}

static mpfr_t r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22042, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22043);
        mpfr_init(r22044);
        mpfr_init_set_str(r22045, "3", 10, MPFR_RNDN);
        mpfr_init(r22046);
        mpfr_init(r22047);
        mpfr_init_set_str(r22048, "2", 10, MPFR_RNDN);
        mpfr_init(r22049);
        mpfr_init_set_str(r22050, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r22051, "7", 10, MPFR_RNDN);
        mpfr_init(r22052);
        mpfr_init(r22053);
        mpfr_init_set_str(r22054, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r22055, "5", 10, MPFR_RNDN);
        mpfr_init(r22056);
        mpfr_init(r22057);
        mpfr_init(r22058);
        mpfr_init(r22059);
        mpfr_init(r22060);
        mpfr_init_set_str(r22061, "1", 10, MPFR_RNDN);
        mpfr_init(r22062);
        mpfr_init(r22063);
        mpfr_init(r22064);
        mpfr_init(r22065);
        mpfr_init(r22066);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r22043, x, MPFR_RNDN);
        mpfr_abs(r22044, r22043, MPFR_RNDN);
        ;
        mpfr_pow(r22046, r22044, r22045, MPFR_RNDN);
        mpfr_mul(r22047, r22042, r22046, MPFR_RNDN);
        ;
        mpfr_mul(r22049, r22048, r22044, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22052, r22044, r22051, MPFR_RNDN);
        mpfr_mul(r22053, r22050, r22052, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22056, r22044, r22055, MPFR_RNDN);
        mpfr_mul(r22057, r22054, r22056, MPFR_RNDN);
        mpfr_add(r22058, r22053, r22057, MPFR_RNDN);
        mpfr_add(r22059, r22049, r22058, MPFR_RNDN);
        mpfr_add(r22060, r22047, r22059, MPFR_RNDN);
        ;
        mpfr_const_pi(r22062, MPFR_RNDN);
        mpfr_div(r22063, r22061, r22062, MPFR_RNDN);
        mpfr_sqrt(r22064, r22063, MPFR_RNDN);
        mpfr_mul(r22065, r22060, r22064, MPFR_RNDN);
        mpfr_abs(r22066, r22065, MPFR_RNDN);
        return mpfr_get_d(r22066, MPFR_RNDN);
}