#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 = r18964 + r18975;
        float r18977 = r18976 * (r18976 * r18976);
        float r18978 = cbrt(r18977);
        float r18979 = r18963 * r18978;
        float r18980 = r18960 * r18960;
        float r18981 = -r18980;
        float r18982 = exp(r18981);
        float r18983 = r18979 * r18982;
        float r18984 = r18957 - r18983;
        float r18985 = exp(r18984);
        float r18986 = log(r18985);
        return r18986;
}

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

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

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19017, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19018, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19019);
        mpfr_init(r19020);
        mpfr_init(r19021);
        mpfr_init(r19022);
        mpfr_init(r19023);
        mpfr_init_set_str(r19024, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19025, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19026, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19027, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19028, "1.061405429", 10, MPFR_RNDN);
        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);
        mpfr_init(r19039);
        mpfr_init(r19040);
        mpfr_init(r19041);
        mpfr_init(r19042);
}

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

static mpfr_t r19043, r19044, r19045, r19046, r19047, r19048, r19049, r19050, r19051, r19052, r19053, r19054, r19055, r19056, r19057, r19058, r19059, r19060, r19061, r19062, r19063, r19064, r19065, r19066, r19067, r19068, r19069, r19070, r19071, r19072;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19043, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19044, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19045);
        mpfr_init(r19046);
        mpfr_init(r19047);
        mpfr_init(r19048);
        mpfr_init(r19049);
        mpfr_init_set_str(r19050, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19051, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19052, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19053, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19054, "1.061405429", 10, MPFR_RNDN);
        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);
        mpfr_init(r19067);
        mpfr_init(r19068);
        mpfr_init(r19069);
        mpfr_init(r19070);
        mpfr_init(r19071);
        mpfr_init(r19072);
}

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

static mpfr_t r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086, r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19073, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19074, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19075);
        mpfr_init(r19076);
        mpfr_init(r19077);
        mpfr_init(r19078);
        mpfr_init(r19079);
        mpfr_init_set_str(r19080, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19081, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19082, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19083, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19084, "1.061405429", 10, MPFR_RNDN);
        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);
        mpfr_init(r19095);
        mpfr_init(r19096);
        mpfr_init(r19097);
        mpfr_init(r19098);
        mpfr_init(r19099);
        mpfr_init(r19100);
        mpfr_init(r19101);
        mpfr_init(r19102);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r19075, x, MPFR_RNDN);
        mpfr_abs(r19076, r19075, MPFR_RNDN);
        mpfr_mul(r19077, r19074, r19076, MPFR_RNDN);
        mpfr_add(r19078, r19073, r19077, MPFR_RNDN);
        mpfr_div(r19079, r19073, r19078, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r19085, r19079, r19084, MPFR_RNDN);
        mpfr_add(r19086, r19083, r19085, MPFR_RNDN);
        mpfr_mul(r19087, r19079, r19086, MPFR_RNDN);
        mpfr_add(r19088, r19082, r19087, MPFR_RNDN);
        mpfr_mul(r19089, r19079, r19088, MPFR_RNDN);
        mpfr_add(r19090, r19081, r19089, MPFR_RNDN);
        mpfr_mul(r19091, r19079, r19090, MPFR_RNDN);
        mpfr_add(r19092, r19080, r19091, MPFR_RNDN);
        mpfr_mul(r19093, r19092, r19092, MPFR_RNDN); mpfr_mul(r19093, r19093, r19092, MPFR_RNDN);
        mpfr_cbrt(r19094, r19093, MPFR_RNDN);
        mpfr_mul(r19095, r19079, r19094, MPFR_RNDN);
        mpfr_mul(r19096, r19076, r19076, MPFR_RNDN);
        mpfr_neg(r19097, r19096, MPFR_RNDN);
        mpfr_exp(r19098, r19097, MPFR_RNDN);
        mpfr_mul(r19099, r19095, r19098, MPFR_RNDN);
        mpfr_sub(r19100, r19073, r19099, MPFR_RNDN);
        mpfr_exp(r19101, r19100, MPFR_RNDN);
        mpfr_log(r19102, r19101, MPFR_RNDN);
        return mpfr_get_d(r19102, MPFR_RNDN);
}