#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.284496736f;
        float r18979 = 0.3275911f;
        float r18980 = x;
        float r18981 = fabs(r18980);
        float r18982 = r18979 * r18981;
        float r18983 = r18977 + r18982;
        float r18984 = r18978 / r18983;
        float r18985 = 0.254829592f;
        float r18986 = r18984 + r18985;
        float r18987 = 1.421413741f;
        float r18988 = r18987 / r18983;
        float r18989 = r18988 / r18983;
        float r18990 = r18986 + r18989;
        float r18991 = -1.453152027f;
        float r18992 = 1.061405429f;
        float r18993 = r18992 / r18983;
        float r18994 = r18991 + r18993;
        float r18995 = r18983 * (r18983 * r18983);
        float r18996 = r18994 / r18995;
        float r18997 = r18990 + r18996;
        float r18998 = r18981 * r18981;
        float r18999 = exp(r18998);
        float r19000 = r18999 * r18983;
        float r19001 = r18997 / r19000;
        float r19002 = r18977 - r19001;
        return r19002;
}

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

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

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19029, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19030, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19031);
        mpfr_init(r19032);
        mpfr_init(r19033);
        mpfr_init(r19034);
        mpfr_init(r19035);
        mpfr_init_set_str(r19036, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r19037, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19038, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r19039, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19040, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19041);
        mpfr_init(r19042);
        mpfr_init(r19043);
        mpfr_init(r19044);
        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);
}

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

static mpfr_t r19055, r19056, r19057, r19058, r19059, r19060, r19061, r19062, r19063, r19064, r19065, r19066, r19067, r19068, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080;

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

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r19058, x, MPFR_RNDN);
        mpfr_abs(r19059, r19058, MPFR_RNDN);
        mpfr_mul(r19060, r19057, r19059, MPFR_RNDN);
        mpfr_add(r19061, r19055, r19060, MPFR_RNDN);
        mpfr_div(r19062, r19056, r19061, MPFR_RNDN);
        ;
        mpfr_add(r19064, r19062, r19063, MPFR_RNDN);
        ;
        mpfr_div(r19066, r19065, r19061, MPFR_RNDN);
        mpfr_div(r19067, r19066, r19061, MPFR_RNDN);
        mpfr_add(r19068, r19064, r19067, MPFR_RNDN);
        ;
        ;
        mpfr_div(r19071, r19070, r19061, MPFR_RNDN);
        mpfr_add(r19072, r19069, r19071, MPFR_RNDN);
        mpfr_mul(r19073, r19061, r19061, MPFR_RNDN); mpfr_mul(r19073, r19073, r19061, MPFR_RNDN);
        mpfr_div(r19074, r19072, r19073, MPFR_RNDN);
        mpfr_add(r19075, r19068, r19074, MPFR_RNDN);
        mpfr_sqr(r19076, r19059, MPFR_RNDN);
        mpfr_exp(r19077, r19076, MPFR_RNDN);
        mpfr_mul(r19078, r19077, r19061, MPFR_RNDN);
        mpfr_div(r19079, r19075, r19078, MPFR_RNDN);
        mpfr_sub(r19080, r19055, r19079, MPFR_RNDN);
        return mpfr_get_d(r19080, MPFR_RNDN);
}

static mpfr_t r19081, r19082, r19083, r19084, r19085, r19086, r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r19081, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r19082, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r19083, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r19084);
        mpfr_init(r19085);
        mpfr_init(r19086);
        mpfr_init(r19087);
        mpfr_init(r19088);
        mpfr_init_set_str(r19089, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r19090);
        mpfr_init_set_str(r19091, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r19092);
        mpfr_init(r19093);
        mpfr_init(r19094);
        mpfr_init_set_str(r19095, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r19096, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r19097);
        mpfr_init(r19098);
        mpfr_init(r19099);
        mpfr_init(r19100);
        mpfr_init(r19101);
        mpfr_init(r19102);
        mpfr_init(r19103);
        mpfr_init(r19104);
        mpfr_init(r19105);
        mpfr_init(r19106);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r19084, x, MPFR_RNDN);
        mpfr_abs(r19085, r19084, MPFR_RNDN);
        mpfr_mul(r19086, r19083, r19085, MPFR_RNDN);
        mpfr_add(r19087, r19081, r19086, MPFR_RNDN);
        mpfr_div(r19088, r19082, r19087, MPFR_RNDN);
        ;
        mpfr_add(r19090, r19088, r19089, MPFR_RNDN);
        ;
        mpfr_div(r19092, r19091, r19087, MPFR_RNDN);
        mpfr_div(r19093, r19092, r19087, MPFR_RNDN);
        mpfr_add(r19094, r19090, r19093, MPFR_RNDN);
        ;
        ;
        mpfr_div(r19097, r19096, r19087, MPFR_RNDN);
        mpfr_add(r19098, r19095, r19097, MPFR_RNDN);
        mpfr_mul(r19099, r19087, r19087, MPFR_RNDN); mpfr_mul(r19099, r19099, r19087, MPFR_RNDN);
        mpfr_div(r19100, r19098, r19099, MPFR_RNDN);
        mpfr_add(r19101, r19094, r19100, MPFR_RNDN);
        mpfr_sqr(r19102, r19085, MPFR_RNDN);
        mpfr_exp(r19103, r19102, MPFR_RNDN);
        mpfr_mul(r19104, r19103, r19087, MPFR_RNDN);
        mpfr_div(r19105, r19101, r19104, MPFR_RNDN);
        mpfr_sub(r19106, r19081, r19105, MPFR_RNDN);
        return mpfr_get_d(r19106, MPFR_RNDN);
}