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

char *name = "Jmat.Real.erf";

double f_if(float x) {
        float r18905 = 1.0f;
        float r18906 = 0.3275911f;
        float r18907 = x;
        float r18908 = fabs(r18907);
        float r18909 = r18906 * r18908;
        float r18910 = r18905 + r18909;
        float r18911 = r18905 / r18910;
        float r18912 = 0.254829592f;
        float r18913 = -0.284496736f;
        float r18914 = 1.421413741f;
        float r18915 = -1.453152027f;
        float r18916 = 1.061405429f;
        float r18917 = r18911 * r18916;
        float r18918 = r18915 + r18917;
        float r18919 = r18911 * r18918;
        float r18920 = r18914 + r18919;
        float r18921 = r18911 * r18920;
        float r18922 = r18913 + r18921;
        float r18923 = r18911 * r18922;
        float r18924 = r18912 + r18923;
        float r18925 = r18911 * r18924;
        float r18926 = r18908 * r18908;
        float r18927 = -r18926;
        float r18928 = exp(r18927);
        float r18929 = r18925 * r18928;
        float r18930 = r18905 - r18929;
        return r18930;
}

double f_id(double x) {
        double r18931 = 1.0;
        double r18932 = 0.3275911;
        double r18933 = x;
        double r18934 = fabs(r18933);
        double r18935 = r18932 * r18934;
        double r18936 = r18931 + r18935;
        double r18937 = r18931 / r18936;
        double r18938 = 0.254829592;
        double r18939 = -0.284496736;
        double r18940 = 1.421413741;
        double r18941 = -1.453152027;
        double r18942 = 1.061405429;
        double r18943 = r18937 * r18942;
        double r18944 = r18941 + r18943;
        double r18945 = r18937 * r18944;
        double r18946 = r18940 + r18945;
        double r18947 = r18937 * r18946;
        double r18948 = r18939 + r18947;
        double r18949 = r18937 * r18948;
        double r18950 = r18938 + r18949;
        double r18951 = r18937 * r18950;
        double r18952 = r18934 * r18934;
        double r18953 = -r18952;
        double r18954 = exp(r18953);
        double r18955 = r18951 * r18954;
        double r18956 = r18931 - r18955;
        return r18956;
}


double f_of(float x) {
        float r18957 = 1.0f;
        float r18958 = 0.3275911f;
        float r18959 = x;
        float r18960 = fabs(r18959);
        float r18961 = r18958 * r18960;
        float r18962 = r18957 + r18961;
        float r18963 = r18957 / r18962;
        float r18964 = 0.254829592f;
        float r18965 = -0.284496736f;
        float r18966 = 1.421413741f;
        float r18967 = -1.453152027f;
        float r18968 = 1.061405429f;
        float r18969 = r18963 * r18968;
        float r18970 = r18967 + r18969;
        float r18971 = r18963 * r18970;
        float r18972 = r18966 + r18971;
        float r18973 = r18963 * r18972;
        float r18974 = r18965 + r18973;
        float r18975 = r18963 * r18974;
        float r18976 = cbrt(r18975);
        float r18977 = r18976 * (r18976 * r18976);
        float r18978 = r18964 + r18977;
        float r18979 = r18963 * r18978;
        float r18980 = r18960 * r18960;
        float r18981 = -r18980;
        float r18982 = exp(r18981);
        float r18983 = r18979 * r18982;
        float r18984 = r18957 - r18983;
        return r18984;
}

double f_od(double x) {
        double r18985 = 1.0;
        double r18986 = 0.3275911;
        double r18987 = x;
        double r18988 = fabs(r18987);
        double r18989 = r18986 * r18988;
        double r18990 = r18985 + r18989;
        double r18991 = r18985 / r18990;
        double r18992 = 0.254829592;
        double r18993 = -0.284496736;
        double r18994 = 1.421413741;
        double r18995 = -1.453152027;
        double r18996 = 1.061405429;
        double r18997 = r18991 * r18996;
        double r18998 = r18995 + r18997;
        double r18999 = r18991 * r18998;
        double r19000 = r18994 + r18999;
        double r19001 = r18991 * r19000;
        double r19002 = r18993 + r19001;
        double r19003 = r18991 * r19002;
        double r19004 = cbrt(r19003);
        double r19005 = r19004 * (r19004 * r19004);
        double r19006 = r18992 + r19005;
        double r19007 = r18991 * r19006;
        double r19008 = r18988 * r18988;
        double r19009 = -r19008;
        double r19010 = exp(r19009);
        double r19011 = r19007 * r19010;
        double r19012 = r18985 - r19011;
        return r19012;
}

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 r19013, r19014, r19015, r19016, r19017, r19018, r19019, r19020, r19021, r19022, r19023, r19024, r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033, r19034, r19035, r19036, r19037, r19038;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19013, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19014, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19015);
        mpfr_init(r19016);
        mpfr_init(r19017);
        mpfr_init(r19018);
        mpfr_init(r19019);
        mpfr_init_set_str(r19020, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19021, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19022, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19023, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19024, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19025);
        mpfr_init(r19026);
        mpfr_init(r19027);
        mpfr_init(r19028);
        mpfr_init(r19029);
        mpfr_init(r19030);
        mpfr_init(r19031);
        mpfr_init(r19032);
        mpfr_init(r19033);
        mpfr_init(r19034);
        mpfr_init(r19035);
        mpfr_init(r19036);
        mpfr_init(r19037);
        mpfr_init(r19038);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r19015, x, MPFR_RNDN);
        mpfr_abs(r19016, r19015, MPFR_RNDN);
        mpfr_mul(r19017, r19014, r19016, MPFR_RNDN);
        mpfr_add(r19018, r19013, r19017, MPFR_RNDN);
        mpfr_div(r19019, r19013, r19018, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r19025, r19019, r19024, MPFR_RNDN);
        mpfr_add(r19026, r19023, r19025, MPFR_RNDN);
        mpfr_mul(r19027, r19019, r19026, MPFR_RNDN);
        mpfr_add(r19028, r19022, r19027, MPFR_RNDN);
        mpfr_mul(r19029, r19019, r19028, MPFR_RNDN);
        mpfr_add(r19030, r19021, r19029, MPFR_RNDN);
        mpfr_mul(r19031, r19019, r19030, MPFR_RNDN);
        mpfr_add(r19032, r19020, r19031, MPFR_RNDN);
        mpfr_mul(r19033, r19019, r19032, MPFR_RNDN);
        mpfr_mul(r19034, r19016, r19016, MPFR_RNDN);
        mpfr_neg(r19035, r19034, MPFR_RNDN);
        mpfr_exp(r19036, r19035, MPFR_RNDN);
        mpfr_mul(r19037, r19033, r19036, MPFR_RNDN);
        mpfr_sub(r19038, r19013, r19037, MPFR_RNDN);
        return mpfr_get_d(r19038, MPFR_RNDN);
}

static mpfr_t r19039, r19040, r19041, r19042, r19043, r19044, r19045, r19046, r19047, r19048, r19049, r19050, r19051, r19052, r19053, r19054, r19055, r19056, r19057, r19058, r19059, r19060, r19061, r19062, r19063, r19064, r19065, r19066;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19039, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19040, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19041);
        mpfr_init(r19042);
        mpfr_init(r19043);
        mpfr_init(r19044);
        mpfr_init(r19045);
        mpfr_init_set_str(r19046, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19047, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19048, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19049, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19050, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19051);
        mpfr_init(r19052);
        mpfr_init(r19053);
        mpfr_init(r19054);
        mpfr_init(r19055);
        mpfr_init(r19056);
        mpfr_init(r19057);
        mpfr_init(r19058);
        mpfr_init(r19059);
        mpfr_init(r19060);
        mpfr_init(r19061);
        mpfr_init(r19062);
        mpfr_init(r19063);
        mpfr_init(r19064);
        mpfr_init(r19065);
        mpfr_init(r19066);
}

double f_fm(double x) {
        ;
        ;
        mpfr_set_d(r19041, x, MPFR_RNDN);
        mpfr_abs(r19042, r19041, MPFR_RNDN);
        mpfr_mul(r19043, r19040, r19042, MPFR_RNDN);
        mpfr_add(r19044, r19039, r19043, MPFR_RNDN);
        mpfr_div(r19045, r19039, r19044, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r19051, r19045, r19050, MPFR_RNDN);
        mpfr_add(r19052, r19049, r19051, MPFR_RNDN);
        mpfr_mul(r19053, r19045, r19052, MPFR_RNDN);
        mpfr_add(r19054, r19048, r19053, MPFR_RNDN);
        mpfr_mul(r19055, r19045, r19054, MPFR_RNDN);
        mpfr_add(r19056, r19047, r19055, MPFR_RNDN);
        mpfr_mul(r19057, r19045, r19056, MPFR_RNDN);
        mpfr_cbrt(r19058, r19057, MPFR_RNDN);
        mpfr_mul(r19059, r19058, r19058, MPFR_RNDN); mpfr_mul(r19059, r19059, r19058, MPFR_RNDN);
        mpfr_add(r19060, r19046, r19059, MPFR_RNDN);
        mpfr_mul(r19061, r19045, r19060, MPFR_RNDN);
        mpfr_mul(r19062, r19042, r19042, MPFR_RNDN);
        mpfr_neg(r19063, r19062, MPFR_RNDN);
        mpfr_exp(r19064, r19063, MPFR_RNDN);
        mpfr_mul(r19065, r19061, r19064, MPFR_RNDN);
        mpfr_sub(r19066, r19039, r19065, MPFR_RNDN);
        return mpfr_get_d(r19066, MPFR_RNDN);
}

static mpfr_t r19067, r19068, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086, r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19067, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19068, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19069);
        mpfr_init(r19070);
        mpfr_init(r19071);
        mpfr_init(r19072);
        mpfr_init(r19073);
        mpfr_init_set_str(r19074, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19075, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19076, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19077, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19078, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19079);
        mpfr_init(r19080);
        mpfr_init(r19081);
        mpfr_init(r19082);
        mpfr_init(r19083);
        mpfr_init(r19084);
        mpfr_init(r19085);
        mpfr_init(r19086);
        mpfr_init(r19087);
        mpfr_init(r19088);
        mpfr_init(r19089);
        mpfr_init(r19090);
        mpfr_init(r19091);
        mpfr_init(r19092);
        mpfr_init(r19093);
        mpfr_init(r19094);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r19069, x, MPFR_RNDN);
        mpfr_abs(r19070, r19069, MPFR_RNDN);
        mpfr_mul(r19071, r19068, r19070, MPFR_RNDN);
        mpfr_add(r19072, r19067, r19071, MPFR_RNDN);
        mpfr_div(r19073, r19067, r19072, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r19079, r19073, r19078, MPFR_RNDN);
        mpfr_add(r19080, r19077, r19079, MPFR_RNDN);
        mpfr_mul(r19081, r19073, r19080, MPFR_RNDN);
        mpfr_add(r19082, r19076, r19081, MPFR_RNDN);
        mpfr_mul(r19083, r19073, r19082, MPFR_RNDN);
        mpfr_add(r19084, r19075, r19083, MPFR_RNDN);
        mpfr_mul(r19085, r19073, r19084, MPFR_RNDN);
        mpfr_cbrt(r19086, r19085, MPFR_RNDN);
        mpfr_mul(r19087, r19086, r19086, MPFR_RNDN); mpfr_mul(r19087, r19087, r19086, MPFR_RNDN);
        mpfr_add(r19088, r19074, r19087, MPFR_RNDN);
        mpfr_mul(r19089, r19073, r19088, MPFR_RNDN);
        mpfr_mul(r19090, r19070, r19070, MPFR_RNDN);
        mpfr_neg(r19091, r19090, MPFR_RNDN);
        mpfr_exp(r19092, r19091, MPFR_RNDN);
        mpfr_mul(r19093, r19089, r19092, MPFR_RNDN);
        mpfr_sub(r19094, r19067, r19093, MPFR_RNDN);
        return mpfr_get_d(r19094, MPFR_RNDN);
}