#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 r21875 = 1.0f;
        float r21876 = atan2(1.0, 0.0);
        float r21877 = sqrt(r21876);
        float r21878 = r21875 / r21877;
        float r21879 = 2.0f;
        float r21880 = x;
        float r21881 = fabs(r21880);
        float r21882 = r21879 * r21881;
        float r21883 = 3.0f;
        float r21884 = r21879 / r21883;
        float r21885 = r21881 * r21881;
        float r21886 = r21885 * r21881;
        float r21887 = r21884 * r21886;
        float r21888 = r21882 + r21887;
        float r21889 = 5.0f;
        float r21890 = r21875 / r21889;
        float r21891 = r21886 * r21881;
        float r21892 = r21891 * r21881;
        float r21893 = r21890 * r21892;
        float r21894 = r21888 + r21893;
        float r21895 = 21.0f;
        float r21896 = r21875 / r21895;
        float r21897 = r21892 * r21881;
        float r21898 = r21897 * r21881;
        float r21899 = r21896 * r21898;
        float r21900 = r21894 + r21899;
        float r21901 = r21878 * r21900;
        float r21902 = fabs(r21901);
        return r21902;
}

double f_id(double x) {
        double r21903 = 1.0;
        double r21904 = atan2(1.0, 0.0);
        double r21905 = sqrt(r21904);
        double r21906 = r21903 / r21905;
        double r21907 = 2.0;
        double r21908 = x;
        double r21909 = fabs(r21908);
        double r21910 = r21907 * r21909;
        double r21911 = 3.0;
        double r21912 = r21907 / r21911;
        double r21913 = r21909 * r21909;
        double r21914 = r21913 * r21909;
        double r21915 = r21912 * r21914;
        double r21916 = r21910 + r21915;
        double r21917 = 5.0;
        double r21918 = r21903 / r21917;
        double r21919 = r21914 * r21909;
        double r21920 = r21919 * r21909;
        double r21921 = r21918 * r21920;
        double r21922 = r21916 + r21921;
        double r21923 = 21.0;
        double r21924 = r21903 / r21923;
        double r21925 = r21920 * r21909;
        double r21926 = r21925 * r21909;
        double r21927 = r21924 * r21926;
        double r21928 = r21922 + r21927;
        double r21929 = r21906 * r21928;
        double r21930 = fabs(r21929);
        return r21930;
}


double f_of(float x) {
        float r21931 = 0.6666666666666666f;
        float r21932 = x;
        float r21933 = fabs(r21932);
        float r21934 = 3.0f;
        float r21935 = pow(r21933, r21934);
        float r21936 = r21931 * r21935;
        float r21937 = 2.0f;
        float r21938 = r21937 * r21933;
        float r21939 = 0.047619047619047616f;
        float r21940 = 7.0f;
        float r21941 = pow(r21933, r21940);
        float r21942 = r21939 * r21941;
        float r21943 = 0.2f;
        float r21944 = 5.0f;
        float r21945 = pow(r21933, r21944);
        float r21946 = r21943 * r21945;
        float r21947 = r21942 + r21946;
        float r21948 = r21938 + r21947;
        float r21949 = r21936 + r21948;
        float r21950 = 1.0f;
        float r21951 = atan2(1.0, 0.0);
        float r21952 = r21950 / r21951;
        float r21953 = sqrt(r21952);
        float r21954 = r21949 * r21953;
        float r21955 = fabs(r21954);
        return r21955;
}

double f_od(double x) {
        double r21956 = 0.6666666666666666;
        double r21957 = x;
        double r21958 = fabs(r21957);
        double r21959 = 3.0;
        double r21960 = pow(r21958, r21959);
        double r21961 = r21956 * r21960;
        double r21962 = 2.0;
        double r21963 = r21962 * r21958;
        double r21964 = 0.047619047619047616;
        double r21965 = 7.0;
        double r21966 = pow(r21958, r21965);
        double r21967 = r21964 * r21966;
        double r21968 = 0.2;
        double r21969 = 5.0;
        double r21970 = pow(r21958, r21969);
        double r21971 = r21968 * r21970;
        double r21972 = r21967 + r21971;
        double r21973 = r21963 + r21972;
        double r21974 = r21961 + r21973;
        double r21975 = 1.0;
        double r21976 = atan2(1.0, 0.0);
        double r21977 = r21975 / r21976;
        double r21978 = sqrt(r21977);
        double r21979 = r21974 * r21978;
        double r21980 = fabs(r21979);
        return r21980;
}

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 r21981, r21982, r21983, r21984, r21985, r21986, r21987, r21988, r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000, r22001, r22002, r22003, r22004, r22005, r22006, r22007, r22008;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r21981, "1", 10, MPFR_RNDN);
        mpfr_init(r21982);
        mpfr_init(r21983);
        mpfr_init(r21984);
        mpfr_init_set_str(r21985, "2", 10, MPFR_RNDN);
        mpfr_init(r21986);
        mpfr_init(r21987);
        mpfr_init(r21988);
        mpfr_init_set_str(r21989, "3", 10, MPFR_RNDN);
        mpfr_init(r21990);
        mpfr_init(r21991);
        mpfr_init(r21992);
        mpfr_init(r21993);
        mpfr_init(r21994);
        mpfr_init_set_str(r21995, "5", 10, MPFR_RNDN);
        mpfr_init(r21996);
        mpfr_init(r21997);
        mpfr_init(r21998);
        mpfr_init(r21999);
        mpfr_init(r22000);
        mpfr_init_set_str(r22001, "21", 10, MPFR_RNDN);
        mpfr_init(r22002);
        mpfr_init(r22003);
        mpfr_init(r22004);
        mpfr_init(r22005);
        mpfr_init(r22006);
        mpfr_init(r22007);
        mpfr_init(r22008);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r21982, MPFR_RNDN);
        mpfr_sqrt(r21983, r21982, MPFR_RNDN);
        mpfr_div(r21984, r21981, r21983, MPFR_RNDN);
        ;
        mpfr_set_d(r21986, x, MPFR_RNDN);
        mpfr_abs(r21987, r21986, MPFR_RNDN);
        mpfr_mul(r21988, r21985, r21987, MPFR_RNDN);
        ;
        mpfr_div(r21990, r21985, r21989, MPFR_RNDN);
        mpfr_mul(r21991, r21987, r21987, MPFR_RNDN);
        mpfr_mul(r21992, r21991, r21987, MPFR_RNDN);
        mpfr_mul(r21993, r21990, r21992, MPFR_RNDN);
        mpfr_add(r21994, r21988, r21993, MPFR_RNDN);
        ;
        mpfr_div(r21996, r21981, r21995, MPFR_RNDN);
        mpfr_mul(r21997, r21992, r21987, MPFR_RNDN);
        mpfr_mul(r21998, r21997, r21987, MPFR_RNDN);
        mpfr_mul(r21999, r21996, r21998, MPFR_RNDN);
        mpfr_add(r22000, r21994, r21999, MPFR_RNDN);
        ;
        mpfr_div(r22002, r21981, r22001, MPFR_RNDN);
        mpfr_mul(r22003, r21998, r21987, MPFR_RNDN);
        mpfr_mul(r22004, r22003, r21987, MPFR_RNDN);
        mpfr_mul(r22005, r22002, r22004, MPFR_RNDN);
        mpfr_add(r22006, r22000, r22005, MPFR_RNDN);
        mpfr_mul(r22007, r21984, r22006, MPFR_RNDN);
        mpfr_abs(r22008, r22007, MPFR_RNDN);
        return mpfr_get_d(r22008, MPFR_RNDN);
}

static mpfr_t r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22009, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22010);
        mpfr_init(r22011);
        mpfr_init_set_str(r22012, "3", 10, MPFR_RNDN);
        mpfr_init(r22013);
        mpfr_init(r22014);
        mpfr_init_set_str(r22015, "2", 10, MPFR_RNDN);
        mpfr_init(r22016);
        mpfr_init_set_str(r22017, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r22018, "7", 10, MPFR_RNDN);
        mpfr_init(r22019);
        mpfr_init(r22020);
        mpfr_init_set_str(r22021, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r22022, "5", 10, MPFR_RNDN);
        mpfr_init(r22023);
        mpfr_init(r22024);
        mpfr_init(r22025);
        mpfr_init(r22026);
        mpfr_init(r22027);
        mpfr_init_set_str(r22028, "1", 10, MPFR_RNDN);
        mpfr_init(r22029);
        mpfr_init(r22030);
        mpfr_init(r22031);
        mpfr_init(r22032);
        mpfr_init(r22033);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r22010, x, MPFR_RNDN);
        mpfr_abs(r22011, r22010, MPFR_RNDN);
        ;
        mpfr_pow(r22013, r22011, r22012, MPFR_RNDN);
        mpfr_mul(r22014, r22009, r22013, MPFR_RNDN);
        ;
        mpfr_mul(r22016, r22015, r22011, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22019, r22011, r22018, MPFR_RNDN);
        mpfr_mul(r22020, r22017, r22019, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22023, r22011, r22022, MPFR_RNDN);
        mpfr_mul(r22024, r22021, r22023, MPFR_RNDN);
        mpfr_add(r22025, r22020, r22024, MPFR_RNDN);
        mpfr_add(r22026, r22016, r22025, MPFR_RNDN);
        mpfr_add(r22027, r22014, r22026, MPFR_RNDN);
        ;
        mpfr_const_pi(r22029, MPFR_RNDN);
        mpfr_div(r22030, r22028, r22029, MPFR_RNDN);
        mpfr_sqrt(r22031, r22030, MPFR_RNDN);
        mpfr_mul(r22032, r22027, r22031, MPFR_RNDN);
        mpfr_abs(r22033, r22032, MPFR_RNDN);
        return mpfr_get_d(r22033, MPFR_RNDN);
}

static mpfr_t r22034, r22035, r22036, r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22034, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22035);
        mpfr_init(r22036);
        mpfr_init_set_str(r22037, "3", 10, MPFR_RNDN);
        mpfr_init(r22038);
        mpfr_init(r22039);
        mpfr_init_set_str(r22040, "2", 10, MPFR_RNDN);
        mpfr_init(r22041);
        mpfr_init_set_str(r22042, "1/21", 10, MPFR_RNDN);
        mpfr_init_set_str(r22043, "7", 10, MPFR_RNDN);
        mpfr_init(r22044);
        mpfr_init(r22045);
        mpfr_init_set_str(r22046, "1/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r22047, "5", 10, MPFR_RNDN);
        mpfr_init(r22048);
        mpfr_init(r22049);
        mpfr_init(r22050);
        mpfr_init(r22051);
        mpfr_init(r22052);
        mpfr_init_set_str(r22053, "1", 10, MPFR_RNDN);
        mpfr_init(r22054);
        mpfr_init(r22055);
        mpfr_init(r22056);
        mpfr_init(r22057);
        mpfr_init(r22058);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r22035, x, MPFR_RNDN);
        mpfr_abs(r22036, r22035, MPFR_RNDN);
        ;
        mpfr_pow(r22038, r22036, r22037, MPFR_RNDN);
        mpfr_mul(r22039, r22034, r22038, MPFR_RNDN);
        ;
        mpfr_mul(r22041, r22040, r22036, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22044, r22036, r22043, MPFR_RNDN);
        mpfr_mul(r22045, r22042, r22044, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22048, r22036, r22047, MPFR_RNDN);
        mpfr_mul(r22049, r22046, r22048, MPFR_RNDN);
        mpfr_add(r22050, r22045, r22049, MPFR_RNDN);
        mpfr_add(r22051, r22041, r22050, MPFR_RNDN);
        mpfr_add(r22052, r22039, r22051, MPFR_RNDN);
        ;
        mpfr_const_pi(r22054, MPFR_RNDN);
        mpfr_div(r22055, r22053, r22054, MPFR_RNDN);
        mpfr_sqrt(r22056, r22055, MPFR_RNDN);
        mpfr_mul(r22057, r22052, r22056, MPFR_RNDN);
        mpfr_abs(r22058, r22057, MPFR_RNDN);
        return mpfr_get_d(r22058, MPFR_RNDN);
}