#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 r18041 = 1.0f;
        float r18042 = 0.3275911f;
        float r18043 = x;
        float r18044 = fabs(r18043);
        float r18045 = r18042 * r18044;
        float r18046 = r18041 + r18045;
        float r18047 = r18041 / r18046;
        float r18048 = 0.254829592f;
        float r18049 = -0.284496736f;
        float r18050 = 1.421413741f;
        float r18051 = -1.453152027f;
        float r18052 = 1.061405429f;
        float r18053 = r18047 * r18052;
        float r18054 = r18051 + r18053;
        float r18055 = r18047 * r18054;
        float r18056 = r18050 + r18055;
        float r18057 = r18047 * r18056;
        float r18058 = r18049 + r18057;
        float r18059 = r18047 * r18058;
        float r18060 = r18048 + r18059;
        float r18061 = r18047 * r18060;
        float r18062 = r18044 * r18044;
        float r18063 = -r18062;
        float r18064 = exp(r18063);
        float r18065 = r18061 * r18064;
        float r18066 = r18041 - r18065;
        return r18066;
}

double f_id(double x) {
        double r18067 = 1.0;
        double r18068 = 0.3275911;
        double r18069 = x;
        double r18070 = fabs(r18069);
        double r18071 = r18068 * r18070;
        double r18072 = r18067 + r18071;
        double r18073 = r18067 / r18072;
        double r18074 = 0.254829592;
        double r18075 = -0.284496736;
        double r18076 = 1.421413741;
        double r18077 = -1.453152027;
        double r18078 = 1.061405429;
        double r18079 = r18073 * r18078;
        double r18080 = r18077 + r18079;
        double r18081 = r18073 * r18080;
        double r18082 = r18076 + r18081;
        double r18083 = r18073 * r18082;
        double r18084 = r18075 + r18083;
        double r18085 = r18073 * r18084;
        double r18086 = r18074 + r18085;
        double r18087 = r18073 * r18086;
        double r18088 = r18070 * r18070;
        double r18089 = -r18088;
        double r18090 = exp(r18089);
        double r18091 = r18087 * r18090;
        double r18092 = r18067 - r18091;
        return r18092;
}


double f_of(float x) {
        float r18093 = 1.0f;
        float r18094 = 1.061405429f;
        float r18095 = 0.3275911f;
        float r18096 = x;
        float r18097 = fabs(r18096);
        float r18098 = fma(r18095, r18097, r18093);
        float r18099 = r18094 / r18098;
        float r18100 = -1.453152027f;
        float r18101 = r18099 + r18100;
        float r18102 = r18093 / r18098;
        float r18103 = r18102 / r18098;
        float r18104 = 1.421413741f;
        float r18105 = r18104 / r18098;
        float r18106 = -0.284496736f;
        float r18107 = r18105 + r18106;
        float r18108 = fma(r18101, r18103, r18107);
        float r18109 = cbrt(r18102);
        float r18110 = exp(r18109);
        float r18111 = log(r18110);
        float r18112 = r18111 * (r18111 * r18111);
        float r18113 = r18112 / r18098;
        float r18114 = 0.254829592f;
        float r18115 = r18114 / r18098;
        float r18116 = fma(r18108, r18113, r18115);
        float r18117 = r18097 * r18097;
        float r18118 = exp(r18117);
        float r18119 = r18116 / r18118;
        float r18120 = r18093 - r18119;
        return r18120;
}

double f_od(double x) {
        double r18121 = 1.0;
        double r18122 = 1.061405429;
        double r18123 = 0.3275911;
        double r18124 = x;
        double r18125 = fabs(r18124);
        double r18126 = fma(r18123, r18125, r18121);
        double r18127 = r18122 / r18126;
        double r18128 = -1.453152027;
        double r18129 = r18127 + r18128;
        double r18130 = r18121 / r18126;
        double r18131 = r18130 / r18126;
        double r18132 = 1.421413741;
        double r18133 = r18132 / r18126;
        double r18134 = -0.284496736;
        double r18135 = r18133 + r18134;
        double r18136 = fma(r18129, r18131, r18135);
        double r18137 = cbrt(r18130);
        double r18138 = exp(r18137);
        double r18139 = log(r18138);
        double r18140 = r18139 * (r18139 * r18139);
        double r18141 = r18140 / r18126;
        double r18142 = 0.254829592;
        double r18143 = r18142 / r18126;
        double r18144 = fma(r18136, r18141, r18143);
        double r18145 = r18125 * r18125;
        double r18146 = exp(r18145);
        double r18147 = r18144 / r18146;
        double r18148 = r18121 - r18147;
        return r18148;
}

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 r18149, r18150, r18151, r18152, r18153, r18154, r18155, r18156, r18157, r18158, r18159, r18160, r18161, r18162, r18163, r18164, r18165, r18166, r18167, r18168, r18169, r18170, r18171, r18172, r18173, r18174;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18149, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18150, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18151);
        mpfr_init(r18152);
        mpfr_init(r18153);
        mpfr_init(r18154);
        mpfr_init(r18155);
        mpfr_init_set_str(r18156, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18157, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18158, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18159, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18160, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18161);
        mpfr_init(r18162);
        mpfr_init(r18163);
        mpfr_init(r18164);
        mpfr_init(r18165);
        mpfr_init(r18166);
        mpfr_init(r18167);
        mpfr_init(r18168);
        mpfr_init(r18169);
        mpfr_init(r18170);
        mpfr_init(r18171);
        mpfr_init(r18172);
        mpfr_init(r18173);
        mpfr_init(r18174);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18151, x, MPFR_RNDN);
        mpfr_abs(r18152, r18151, MPFR_RNDN);
        mpfr_mul(r18153, r18150, r18152, MPFR_RNDN);
        mpfr_add(r18154, r18149, r18153, MPFR_RNDN);
        mpfr_div(r18155, r18149, r18154, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18161, r18155, r18160, MPFR_RNDN);
        mpfr_add(r18162, r18159, r18161, MPFR_RNDN);
        mpfr_mul(r18163, r18155, r18162, MPFR_RNDN);
        mpfr_add(r18164, r18158, r18163, MPFR_RNDN);
        mpfr_mul(r18165, r18155, r18164, MPFR_RNDN);
        mpfr_add(r18166, r18157, r18165, MPFR_RNDN);
        mpfr_mul(r18167, r18155, r18166, MPFR_RNDN);
        mpfr_add(r18168, r18156, r18167, MPFR_RNDN);
        mpfr_mul(r18169, r18155, r18168, MPFR_RNDN);
        mpfr_mul(r18170, r18152, r18152, MPFR_RNDN);
        mpfr_neg(r18171, r18170, MPFR_RNDN);
        mpfr_exp(r18172, r18171, MPFR_RNDN);
        mpfr_mul(r18173, r18169, r18172, MPFR_RNDN);
        mpfr_sub(r18174, r18149, r18173, MPFR_RNDN);
        return mpfr_get_d(r18174, MPFR_RNDN);
}

static mpfr_t r18175, r18176, r18177, r18178, r18179, r18180, r18181, r18182, r18183, r18184, r18185, r18186, r18187, r18188, r18189, r18190, r18191, r18192, r18193, r18194, r18195, r18196, r18197, r18198, r18199, r18200, r18201, r18202;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18175, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18176, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18177, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18178);
        mpfr_init(r18179);
        mpfr_init(r18180);
        mpfr_init(r18181);
        mpfr_init_set_str(r18182, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18183);
        mpfr_init(r18184);
        mpfr_init(r18185);
        mpfr_init_set_str(r18186, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18187);
        mpfr_init_set_str(r18188, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18189);
        mpfr_init(r18190);
        mpfr_init(r18191);
        mpfr_init(r18192);
        mpfr_init(r18193);
        mpfr_init(r18194);
        mpfr_init(r18195);
        mpfr_init_set_str(r18196, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18197);
        mpfr_init(r18198);
        mpfr_init(r18199);
        mpfr_init(r18200);
        mpfr_init(r18201);
        mpfr_init(r18202);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18178, x, MPFR_RNDN);
        mpfr_abs(r18179, r18178, MPFR_RNDN);
        mpfr_fma(r18180, r18177, r18179, r18175, MPFR_RNDN);
        mpfr_div(r18181, r18176, r18180, MPFR_RNDN);
        ;
        mpfr_add(r18183, r18181, r18182, MPFR_RNDN);
        mpfr_div(r18184, r18175, r18180, MPFR_RNDN);
        mpfr_div(r18185, r18184, r18180, MPFR_RNDN);
        ;
        mpfr_div(r18187, r18186, r18180, MPFR_RNDN);
        ;
        mpfr_add(r18189, r18187, r18188, MPFR_RNDN);
        mpfr_fma(r18190, r18183, r18185, r18189, MPFR_RNDN);
        mpfr_cbrt(r18191, r18184, MPFR_RNDN);
        mpfr_exp(r18192, r18191, MPFR_RNDN);
        mpfr_log(r18193, r18192, MPFR_RNDN);
        mpfr_mul(r18194, r18193, r18193, MPFR_RNDN); mpfr_mul(r18194, r18194, r18193, MPFR_RNDN);
        mpfr_div(r18195, r18194, r18180, MPFR_RNDN);
        ;
        mpfr_div(r18197, r18196, r18180, MPFR_RNDN);
        mpfr_fma(r18198, r18190, r18195, r18197, MPFR_RNDN);
        mpfr_mul(r18199, r18179, r18179, MPFR_RNDN);
        mpfr_exp(r18200, r18199, MPFR_RNDN);
        mpfr_div(r18201, r18198, r18200, MPFR_RNDN);
        mpfr_sub(r18202, r18175, r18201, MPFR_RNDN);
        return mpfr_get_d(r18202, MPFR_RNDN);
}

static mpfr_t r18203, r18204, r18205, r18206, r18207, r18208, r18209, r18210, r18211, r18212, r18213, r18214, r18215, r18216, r18217, r18218, r18219, r18220, r18221, r18222, r18223, r18224, r18225, r18226, r18227, r18228, r18229, r18230;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18203, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18204, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18205, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18206);
        mpfr_init(r18207);
        mpfr_init(r18208);
        mpfr_init(r18209);
        mpfr_init_set_str(r18210, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18211);
        mpfr_init(r18212);
        mpfr_init(r18213);
        mpfr_init_set_str(r18214, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18215);
        mpfr_init_set_str(r18216, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18217);
        mpfr_init(r18218);
        mpfr_init(r18219);
        mpfr_init(r18220);
        mpfr_init(r18221);
        mpfr_init(r18222);
        mpfr_init(r18223);
        mpfr_init_set_str(r18224, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18225);
        mpfr_init(r18226);
        mpfr_init(r18227);
        mpfr_init(r18228);
        mpfr_init(r18229);
        mpfr_init(r18230);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18206, x, MPFR_RNDN);
        mpfr_abs(r18207, r18206, MPFR_RNDN);
        mpfr_fma(r18208, r18205, r18207, r18203, MPFR_RNDN);
        mpfr_div(r18209, r18204, r18208, MPFR_RNDN);
        ;
        mpfr_add(r18211, r18209, r18210, MPFR_RNDN);
        mpfr_div(r18212, r18203, r18208, MPFR_RNDN);
        mpfr_div(r18213, r18212, r18208, MPFR_RNDN);
        ;
        mpfr_div(r18215, r18214, r18208, MPFR_RNDN);
        ;
        mpfr_add(r18217, r18215, r18216, MPFR_RNDN);
        mpfr_fma(r18218, r18211, r18213, r18217, MPFR_RNDN);
        mpfr_cbrt(r18219, r18212, MPFR_RNDN);
        mpfr_exp(r18220, r18219, MPFR_RNDN);
        mpfr_log(r18221, r18220, MPFR_RNDN);
        mpfr_mul(r18222, r18221, r18221, MPFR_RNDN); mpfr_mul(r18222, r18222, r18221, MPFR_RNDN);
        mpfr_div(r18223, r18222, r18208, MPFR_RNDN);
        ;
        mpfr_div(r18225, r18224, r18208, MPFR_RNDN);
        mpfr_fma(r18226, r18218, r18223, r18225, MPFR_RNDN);
        mpfr_mul(r18227, r18207, r18207, MPFR_RNDN);
        mpfr_exp(r18228, r18227, MPFR_RNDN);
        mpfr_div(r18229, r18226, r18228, MPFR_RNDN);
        mpfr_sub(r18230, r18203, r18229, MPFR_RNDN);
        return mpfr_get_d(r18230, MPFR_RNDN);
}