#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 r18185 = 1.0f;
        float r18186 = 0.3275911f;
        float r18187 = x;
        float r18188 = fabs(r18187);
        float r18189 = r18186 * r18188;
        float r18190 = r18185 + r18189;
        float r18191 = r18185 / r18190;
        float r18192 = 0.254829592f;
        float r18193 = -0.284496736f;
        float r18194 = 1.421413741f;
        float r18195 = -1.453152027f;
        float r18196 = 1.061405429f;
        float r18197 = r18191 * r18196;
        float r18198 = r18195 + r18197;
        float r18199 = r18191 * r18198;
        float r18200 = r18194 + r18199;
        float r18201 = r18191 * r18200;
        float r18202 = r18193 + r18201;
        float r18203 = r18191 * r18202;
        float r18204 = r18192 + r18203;
        float r18205 = r18191 * r18204;
        float r18206 = r18188 * r18188;
        float r18207 = -r18206;
        float r18208 = exp(r18207);
        float r18209 = r18205 * r18208;
        float r18210 = r18185 - r18209;
        return r18210;
}

double f_id(double x) {
        double r18211 = 1.0;
        double r18212 = 0.3275911;
        double r18213 = x;
        double r18214 = fabs(r18213);
        double r18215 = r18212 * r18214;
        double r18216 = r18211 + r18215;
        double r18217 = r18211 / r18216;
        double r18218 = 0.254829592;
        double r18219 = -0.284496736;
        double r18220 = 1.421413741;
        double r18221 = -1.453152027;
        double r18222 = 1.061405429;
        double r18223 = r18217 * r18222;
        double r18224 = r18221 + r18223;
        double r18225 = r18217 * r18224;
        double r18226 = r18220 + r18225;
        double r18227 = r18217 * r18226;
        double r18228 = r18219 + r18227;
        double r18229 = r18217 * r18228;
        double r18230 = r18218 + r18229;
        double r18231 = r18217 * r18230;
        double r18232 = r18214 * r18214;
        double r18233 = -r18232;
        double r18234 = exp(r18233);
        double r18235 = r18231 * r18234;
        double r18236 = r18211 - r18235;
        return r18236;
}


double f_of(float x) {
        float r18237 = 1.0f;
        float r18238 = 1.061405429f;
        float r18239 = 0.3275911f;
        float r18240 = x;
        float r18241 = fabs(r18240);
        float r18242 = fma(r18239, r18241, r18237);
        float r18243 = r18237 / r18242;
        float r18244 = r18238 * r18243;
        float r18245 = 1.453152027f;
        float r18246 = r18244 - r18245;
        float r18247 = r18242 * r18242;
        float r18248 = r18237 / r18247;
        float r18249 = 1.421413741f;
        float r18250 = r18249 * r18243;
        float r18251 = 0.284496736f;
        float r18252 = r18250 - r18251;
        float r18253 = fma(r18246, r18248, r18252);
        float r18254 = 0.254829592f;
        float r18255 = r18254 / r18242;
        float r18256 = fma(r18253, r18248, r18255);
        float r18257 = exp(r18256);
        float r18258 = log(r18257);
        float r18259 = r18241 * r18241;
        float r18260 = exp(r18259);
        float r18261 = r18258 / r18260;
        float r18262 = r18237 - r18261;
        float r18263 = log(r18262);
        float r18264 = exp(r18263);
        return r18264;
}

double f_od(double x) {
        double r18265 = 1.0;
        double r18266 = 1.061405429;
        double r18267 = 0.3275911;
        double r18268 = x;
        double r18269 = fabs(r18268);
        double r18270 = fma(r18267, r18269, r18265);
        double r18271 = r18265 / r18270;
        double r18272 = r18266 * r18271;
        double r18273 = 1.453152027;
        double r18274 = r18272 - r18273;
        double r18275 = r18270 * r18270;
        double r18276 = r18265 / r18275;
        double r18277 = 1.421413741;
        double r18278 = r18277 * r18271;
        double r18279 = 0.284496736;
        double r18280 = r18278 - r18279;
        double r18281 = fma(r18274, r18276, r18280);
        double r18282 = 0.254829592;
        double r18283 = r18282 / r18270;
        double r18284 = fma(r18281, r18276, r18283);
        double r18285 = exp(r18284);
        double r18286 = log(r18285);
        double r18287 = r18269 * r18269;
        double r18288 = exp(r18287);
        double r18289 = r18286 / r18288;
        double r18290 = r18265 - r18289;
        double r18291 = log(r18290);
        double r18292 = exp(r18291);
        return r18292;
}

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 r18293, r18294, r18295, r18296, r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309, r18310, r18311, r18312, r18313, r18314, r18315, r18316, r18317, r18318;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18293, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18294, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18295);
        mpfr_init(r18296);
        mpfr_init(r18297);
        mpfr_init(r18298);
        mpfr_init(r18299);
        mpfr_init_set_str(r18300, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18301, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18302, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18303, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18304, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18305);
        mpfr_init(r18306);
        mpfr_init(r18307);
        mpfr_init(r18308);
        mpfr_init(r18309);
        mpfr_init(r18310);
        mpfr_init(r18311);
        mpfr_init(r18312);
        mpfr_init(r18313);
        mpfr_init(r18314);
        mpfr_init(r18315);
        mpfr_init(r18316);
        mpfr_init(r18317);
        mpfr_init(r18318);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18295, x, MPFR_RNDN);
        mpfr_abs(r18296, r18295, MPFR_RNDN);
        mpfr_mul(r18297, r18294, r18296, MPFR_RNDN);
        mpfr_add(r18298, r18293, r18297, MPFR_RNDN);
        mpfr_div(r18299, r18293, r18298, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18305, r18299, r18304, MPFR_RNDN);
        mpfr_add(r18306, r18303, r18305, MPFR_RNDN);
        mpfr_mul(r18307, r18299, r18306, MPFR_RNDN);
        mpfr_add(r18308, r18302, r18307, MPFR_RNDN);
        mpfr_mul(r18309, r18299, r18308, MPFR_RNDN);
        mpfr_add(r18310, r18301, r18309, MPFR_RNDN);
        mpfr_mul(r18311, r18299, r18310, MPFR_RNDN);
        mpfr_add(r18312, r18300, r18311, MPFR_RNDN);
        mpfr_mul(r18313, r18299, r18312, MPFR_RNDN);
        mpfr_mul(r18314, r18296, r18296, MPFR_RNDN);
        mpfr_neg(r18315, r18314, MPFR_RNDN);
        mpfr_exp(r18316, r18315, MPFR_RNDN);
        mpfr_mul(r18317, r18313, r18316, MPFR_RNDN);
        mpfr_sub(r18318, r18293, r18317, MPFR_RNDN);
        return mpfr_get_d(r18318, MPFR_RNDN);
}

static mpfr_t r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326, r18327, r18328, r18329, r18330, r18331, r18332, r18333, r18334, r18335, r18336, r18337, r18338, r18339, r18340, r18341, r18342, r18343, r18344, r18345, r18346;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18319, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18320, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18321, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18322);
        mpfr_init(r18323);
        mpfr_init(r18324);
        mpfr_init(r18325);
        mpfr_init(r18326);
        mpfr_init_set_str(r18327, "1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18328);
        mpfr_init(r18329);
        mpfr_init(r18330);
        mpfr_init_set_str(r18331, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18332);
        mpfr_init_set_str(r18333, "0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18334);
        mpfr_init(r18335);
        mpfr_init_set_str(r18336, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18337);
        mpfr_init(r18338);
        mpfr_init(r18339);
        mpfr_init(r18340);
        mpfr_init(r18341);
        mpfr_init(r18342);
        mpfr_init(r18343);
        mpfr_init(r18344);
        mpfr_init(r18345);
        mpfr_init(r18346);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18322, x, MPFR_RNDN);
        mpfr_abs(r18323, r18322, MPFR_RNDN);
        mpfr_fma(r18324, r18321, r18323, r18319, MPFR_RNDN);
        mpfr_div(r18325, r18319, r18324, MPFR_RNDN);
        mpfr_mul(r18326, r18320, r18325, MPFR_RNDN);
        ;
        mpfr_sub(r18328, r18326, r18327, MPFR_RNDN);
        mpfr_sqr(r18329, r18324, MPFR_RNDN);
        mpfr_div(r18330, r18319, r18329, MPFR_RNDN);
        ;
        mpfr_mul(r18332, r18331, r18325, MPFR_RNDN);
        ;
        mpfr_sub(r18334, r18332, r18333, MPFR_RNDN);
        mpfr_fma(r18335, r18328, r18330, r18334, MPFR_RNDN);
        ;
        mpfr_div(r18337, r18336, r18324, MPFR_RNDN);
        mpfr_fma(r18338, r18335, r18330, r18337, MPFR_RNDN);
        mpfr_exp(r18339, r18338, MPFR_RNDN);
        mpfr_log(r18340, r18339, MPFR_RNDN);
        mpfr_sqr(r18341, r18323, MPFR_RNDN);
        mpfr_exp(r18342, r18341, MPFR_RNDN);
        mpfr_div(r18343, r18340, r18342, MPFR_RNDN);
        mpfr_sub(r18344, r18319, r18343, MPFR_RNDN);
        mpfr_log(r18345, r18344, MPFR_RNDN);
        mpfr_exp(r18346, r18345, MPFR_RNDN);
        return mpfr_get_d(r18346, MPFR_RNDN);
}

static mpfr_t r18347, r18348, r18349, r18350, r18351, r18352, r18353, r18354, r18355, r18356, r18357, r18358, r18359, r18360, r18361, r18362, r18363, r18364, r18365, r18366, r18367, r18368, r18369, r18370, r18371, r18372, r18373, r18374;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18347, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18348, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18349, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18350);
        mpfr_init(r18351);
        mpfr_init(r18352);
        mpfr_init(r18353);
        mpfr_init(r18354);
        mpfr_init_set_str(r18355, "1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18356);
        mpfr_init(r18357);
        mpfr_init(r18358);
        mpfr_init_set_str(r18359, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18360);
        mpfr_init_set_str(r18361, "0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18362);
        mpfr_init(r18363);
        mpfr_init_set_str(r18364, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18365);
        mpfr_init(r18366);
        mpfr_init(r18367);
        mpfr_init(r18368);
        mpfr_init(r18369);
        mpfr_init(r18370);
        mpfr_init(r18371);
        mpfr_init(r18372);
        mpfr_init(r18373);
        mpfr_init(r18374);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18350, x, MPFR_RNDN);
        mpfr_abs(r18351, r18350, MPFR_RNDN);
        mpfr_fma(r18352, r18349, r18351, r18347, MPFR_RNDN);
        mpfr_div(r18353, r18347, r18352, MPFR_RNDN);
        mpfr_mul(r18354, r18348, r18353, MPFR_RNDN);
        ;
        mpfr_sub(r18356, r18354, r18355, MPFR_RNDN);
        mpfr_sqr(r18357, r18352, MPFR_RNDN);
        mpfr_div(r18358, r18347, r18357, MPFR_RNDN);
        ;
        mpfr_mul(r18360, r18359, r18353, MPFR_RNDN);
        ;
        mpfr_sub(r18362, r18360, r18361, MPFR_RNDN);
        mpfr_fma(r18363, r18356, r18358, r18362, MPFR_RNDN);
        ;
        mpfr_div(r18365, r18364, r18352, MPFR_RNDN);
        mpfr_fma(r18366, r18363, r18358, r18365, MPFR_RNDN);
        mpfr_exp(r18367, r18366, MPFR_RNDN);
        mpfr_log(r18368, r18367, MPFR_RNDN);
        mpfr_sqr(r18369, r18351, MPFR_RNDN);
        mpfr_exp(r18370, r18369, MPFR_RNDN);
        mpfr_div(r18371, r18368, r18370, MPFR_RNDN);
        mpfr_sub(r18372, r18347, r18371, MPFR_RNDN);
        mpfr_log(r18373, r18372, MPFR_RNDN);
        mpfr_exp(r18374, r18373, MPFR_RNDN);
        return mpfr_get_d(r18374, MPFR_RNDN);
}