#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 r18925 = 1.0f;
        float r18926 = 0.3275911f;
        float r18927 = x;
        float r18928 = fabs(r18927);
        float r18929 = r18926 * r18928;
        float r18930 = r18925 + r18929;
        float r18931 = r18925 / r18930;
        float r18932 = 0.254829592f;
        float r18933 = -0.284496736f;
        float r18934 = 1.421413741f;
        float r18935 = -1.453152027f;
        float r18936 = 1.061405429f;
        float r18937 = r18931 * r18936;
        float r18938 = r18935 + r18937;
        float r18939 = r18931 * r18938;
        float r18940 = r18934 + r18939;
        float r18941 = r18931 * r18940;
        float r18942 = r18933 + r18941;
        float r18943 = r18931 * r18942;
        float r18944 = r18932 + r18943;
        float r18945 = r18931 * r18944;
        float r18946 = r18928 * r18928;
        float r18947 = -r18946;
        float r18948 = exp(r18947);
        float r18949 = r18945 * r18948;
        float r18950 = r18925 - r18949;
        return r18950;
}

double f_id(double x) {
        double r18951 = 1.0;
        double r18952 = 0.3275911;
        double r18953 = x;
        double r18954 = fabs(r18953);
        double r18955 = r18952 * r18954;
        double r18956 = r18951 + r18955;
        double r18957 = r18951 / r18956;
        double r18958 = 0.254829592;
        double r18959 = -0.284496736;
        double r18960 = 1.421413741;
        double r18961 = -1.453152027;
        double r18962 = 1.061405429;
        double r18963 = r18957 * r18962;
        double r18964 = r18961 + r18963;
        double r18965 = r18957 * r18964;
        double r18966 = r18960 + r18965;
        double r18967 = r18957 * r18966;
        double r18968 = r18959 + r18967;
        double r18969 = r18957 * r18968;
        double r18970 = r18958 + r18969;
        double r18971 = r18957 * r18970;
        double r18972 = r18954 * r18954;
        double r18973 = -r18972;
        double r18974 = exp(r18973);
        double r18975 = r18971 * r18974;
        double r18976 = r18951 - r18975;
        return r18976;
}


double f_of(float x) {
        float r18977 = 1.0f;
        float r18978 = 0.3275911f;
        float r18979 = x;
        float r18980 = fabs(r18979);
        float r18981 = r18978 * r18980;
        float r18982 = r18977 + r18981;
        float r18983 = r18977 / r18982;
        float r18984 = 0.254829592f;
        float r18985 = -0.284496736f;
        float r18986 = exp(r18983);
        float r18987 = log(r18986);
        float r18988 = 1.421413741f;
        float r18989 = -1.453152027f;
        float r18990 = 1.061405429f;
        float r18991 = r18983 * r18990;
        float r18992 = r18989 + r18991;
        float r18993 = r18983 * r18992;
        float r18994 = r18988 + r18993;
        float r18995 = r18987 * r18994;
        float r18996 = r18985 + r18995;
        float r18997 = r18983 * r18996;
        float r18998 = r18984 + r18997;
        float r18999 = r18983 * r18998;
        float r19000 = r18980 * r18980;
        float r19001 = -r19000;
        float r19002 = exp(r19001);
        float r19003 = r18999 * r19002;
        float r19004 = r18977 - r19003;
        return r19004;
}

double f_od(double x) {
        double r19005 = 1.0;
        double r19006 = 0.3275911;
        double r19007 = x;
        double r19008 = fabs(r19007);
        double r19009 = r19006 * r19008;
        double r19010 = r19005 + r19009;
        double r19011 = r19005 / r19010;
        double r19012 = 0.254829592;
        double r19013 = -0.284496736;
        double r19014 = exp(r19011);
        double r19015 = log(r19014);
        double r19016 = 1.421413741;
        double r19017 = -1.453152027;
        double r19018 = 1.061405429;
        double r19019 = r19011 * r19018;
        double r19020 = r19017 + r19019;
        double r19021 = r19011 * r19020;
        double r19022 = r19016 + r19021;
        double r19023 = r19015 * r19022;
        double r19024 = r19013 + r19023;
        double r19025 = r19011 * r19024;
        double r19026 = r19012 + r19025;
        double r19027 = r19011 * r19026;
        double r19028 = r19008 * r19008;
        double r19029 = -r19028;
        double r19030 = exp(r19029);
        double r19031 = r19027 * r19030;
        double r19032 = r19005 - r19031;
        return r19032;
}

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 r19033, r19034, r19035, r19036, r19037, r19038, r19039, r19040, r19041, r19042, r19043, r19044, r19045, r19046, r19047, r19048, r19049, r19050, r19051, r19052, r19053, r19054, r19055, r19056, r19057, r19058;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19033, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19034, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19035);
        mpfr_init(r19036);
        mpfr_init(r19037);
        mpfr_init(r19038);
        mpfr_init(r19039);
        mpfr_init_set_str(r19040, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19041, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19042, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19043, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19044, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19045);
        mpfr_init(r19046);
        mpfr_init(r19047);
        mpfr_init(r19048);
        mpfr_init(r19049);
        mpfr_init(r19050);
        mpfr_init(r19051);
        mpfr_init(r19052);
        mpfr_init(r19053);
        mpfr_init(r19054);
        mpfr_init(r19055);
        mpfr_init(r19056);
        mpfr_init(r19057);
        mpfr_init(r19058);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r19035, x, MPFR_RNDN);
        mpfr_abs(r19036, r19035, MPFR_RNDN);
        mpfr_mul(r19037, r19034, r19036, MPFR_RNDN);
        mpfr_add(r19038, r19033, r19037, MPFR_RNDN);
        mpfr_div(r19039, r19033, r19038, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r19045, r19039, r19044, MPFR_RNDN);
        mpfr_add(r19046, r19043, r19045, MPFR_RNDN);
        mpfr_mul(r19047, r19039, r19046, MPFR_RNDN);
        mpfr_add(r19048, r19042, r19047, MPFR_RNDN);
        mpfr_mul(r19049, r19039, r19048, MPFR_RNDN);
        mpfr_add(r19050, r19041, r19049, MPFR_RNDN);
        mpfr_mul(r19051, r19039, r19050, MPFR_RNDN);
        mpfr_add(r19052, r19040, r19051, MPFR_RNDN);
        mpfr_mul(r19053, r19039, r19052, MPFR_RNDN);
        mpfr_mul(r19054, r19036, r19036, MPFR_RNDN);
        mpfr_neg(r19055, r19054, MPFR_RNDN);
        mpfr_exp(r19056, r19055, MPFR_RNDN);
        mpfr_mul(r19057, r19053, r19056, MPFR_RNDN);
        mpfr_sub(r19058, r19033, r19057, MPFR_RNDN);
        return mpfr_get_d(r19058, MPFR_RNDN);
}

static mpfr_t r19059, r19060, r19061, r19062, r19063, r19064, r19065, r19066, r19067, r19068, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19059, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19060, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19061);
        mpfr_init(r19062);
        mpfr_init(r19063);
        mpfr_init(r19064);
        mpfr_init(r19065);
        mpfr_init_set_str(r19066, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19067, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r19068);
        mpfr_init(r19069);
        mpfr_init_set_str(r19070, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19071, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19072, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19073);
        mpfr_init(r19074);
        mpfr_init(r19075);
        mpfr_init(r19076);
        mpfr_init(r19077);
        mpfr_init(r19078);
        mpfr_init(r19079);
        mpfr_init(r19080);
        mpfr_init(r19081);
        mpfr_init(r19082);
        mpfr_init(r19083);
        mpfr_init(r19084);
        mpfr_init(r19085);
        mpfr_init(r19086);
}

double f_fm(double x) {
        ;
        ;
        mpfr_set_d(r19061, x, MPFR_RNDN);
        mpfr_abs(r19062, r19061, MPFR_RNDN);
        mpfr_mul(r19063, r19060, r19062, MPFR_RNDN);
        mpfr_add(r19064, r19059, r19063, MPFR_RNDN);
        mpfr_div(r19065, r19059, r19064, MPFR_RNDN);
        ;
        ;
        mpfr_exp(r19068, r19065, MPFR_RNDN);
        mpfr_log(r19069, r19068, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r19073, r19065, r19072, MPFR_RNDN);
        mpfr_add(r19074, r19071, r19073, MPFR_RNDN);
        mpfr_mul(r19075, r19065, r19074, MPFR_RNDN);
        mpfr_add(r19076, r19070, r19075, MPFR_RNDN);
        mpfr_mul(r19077, r19069, r19076, MPFR_RNDN);
        mpfr_add(r19078, r19067, r19077, MPFR_RNDN);
        mpfr_mul(r19079, r19065, r19078, MPFR_RNDN);
        mpfr_add(r19080, r19066, r19079, MPFR_RNDN);
        mpfr_mul(r19081, r19065, r19080, MPFR_RNDN);
        mpfr_mul(r19082, r19062, r19062, MPFR_RNDN);
        mpfr_neg(r19083, r19082, MPFR_RNDN);
        mpfr_exp(r19084, r19083, MPFR_RNDN);
        mpfr_mul(r19085, r19081, r19084, MPFR_RNDN);
        mpfr_sub(r19086, r19059, r19085, MPFR_RNDN);
        return mpfr_get_d(r19086, MPFR_RNDN);
}

static mpfr_t r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19087, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19088, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19089);
        mpfr_init(r19090);
        mpfr_init(r19091);
        mpfr_init(r19092);
        mpfr_init(r19093);
        mpfr_init_set_str(r19094, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19095, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r19096);
        mpfr_init(r19097);
        mpfr_init_set_str(r19098, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19099, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19100, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19101);
        mpfr_init(r19102);
        mpfr_init(r19103);
        mpfr_init(r19104);
        mpfr_init(r19105);
        mpfr_init(r19106);
        mpfr_init(r19107);
        mpfr_init(r19108);
        mpfr_init(r19109);
        mpfr_init(r19110);
        mpfr_init(r19111);
        mpfr_init(r19112);
        mpfr_init(r19113);
        mpfr_init(r19114);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r19089, x, MPFR_RNDN);
        mpfr_abs(r19090, r19089, MPFR_RNDN);
        mpfr_mul(r19091, r19088, r19090, MPFR_RNDN);
        mpfr_add(r19092, r19087, r19091, MPFR_RNDN);
        mpfr_div(r19093, r19087, r19092, MPFR_RNDN);
        ;
        ;
        mpfr_exp(r19096, r19093, MPFR_RNDN);
        mpfr_log(r19097, r19096, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r19101, r19093, r19100, MPFR_RNDN);
        mpfr_add(r19102, r19099, r19101, MPFR_RNDN);
        mpfr_mul(r19103, r19093, r19102, MPFR_RNDN);
        mpfr_add(r19104, r19098, r19103, MPFR_RNDN);
        mpfr_mul(r19105, r19097, r19104, MPFR_RNDN);
        mpfr_add(r19106, r19095, r19105, MPFR_RNDN);
        mpfr_mul(r19107, r19093, r19106, MPFR_RNDN);
        mpfr_add(r19108, r19094, r19107, MPFR_RNDN);
        mpfr_mul(r19109, r19093, r19108, MPFR_RNDN);
        mpfr_mul(r19110, r19090, r19090, MPFR_RNDN);
        mpfr_neg(r19111, r19110, MPFR_RNDN);
        mpfr_exp(r19112, r19111, MPFR_RNDN);
        mpfr_mul(r19113, r19109, r19112, MPFR_RNDN);
        mpfr_sub(r19114, r19087, r19113, MPFR_RNDN);
        return mpfr_get_d(r19114, MPFR_RNDN);
}