#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 r18397 = 1.0f;
        float r18398 = 0.3275911f;
        float r18399 = x;
        float r18400 = fabs(r18399);
        float r18401 = r18398 * r18400;
        float r18402 = r18397 + r18401;
        float r18403 = r18397 / r18402;
        float r18404 = 0.254829592f;
        float r18405 = -0.284496736f;
        float r18406 = 1.421413741f;
        float r18407 = -1.453152027f;
        float r18408 = 1.061405429f;
        float r18409 = r18403 * r18408;
        float r18410 = r18407 + r18409;
        float r18411 = r18403 * r18410;
        float r18412 = r18406 + r18411;
        float r18413 = r18403 * r18412;
        float r18414 = r18405 + r18413;
        float r18415 = r18403 * r18414;
        float r18416 = r18404 + r18415;
        float r18417 = r18403 * r18416;
        float r18418 = r18400 * r18400;
        float r18419 = -r18418;
        float r18420 = exp(r18419);
        float r18421 = r18417 * r18420;
        float r18422 = r18397 - r18421;
        return r18422;
}

double f_id(double x) {
        double r18423 = 1.0;
        double r18424 = 0.3275911;
        double r18425 = x;
        double r18426 = fabs(r18425);
        double r18427 = r18424 * r18426;
        double r18428 = r18423 + r18427;
        double r18429 = r18423 / r18428;
        double r18430 = 0.254829592;
        double r18431 = -0.284496736;
        double r18432 = 1.421413741;
        double r18433 = -1.453152027;
        double r18434 = 1.061405429;
        double r18435 = r18429 * r18434;
        double r18436 = r18433 + r18435;
        double r18437 = r18429 * r18436;
        double r18438 = r18432 + r18437;
        double r18439 = r18429 * r18438;
        double r18440 = r18431 + r18439;
        double r18441 = r18429 * r18440;
        double r18442 = r18430 + r18441;
        double r18443 = r18429 * r18442;
        double r18444 = r18426 * r18426;
        double r18445 = -r18444;
        double r18446 = exp(r18445);
        double r18447 = r18443 * r18446;
        double r18448 = r18423 - r18447;
        return r18448;
}


double f_of(float x) {
        float r18449 = 1.0f;
        float r18450 = 1.061405429f;
        float r18451 = 0.3275911f;
        float r18452 = x;
        float r18453 = fabs(r18452);
        float r18454 = fma(r18451, r18453, r18449);
        float r18455 = r18450 / r18454;
        float r18456 = -1.453152027f;
        float r18457 = r18455 + r18456;
        float r18458 = r18449 / r18454;
        float r18459 = r18458 / r18454;
        float r18460 = 1.421413741f;
        float r18461 = r18460 / r18454;
        float r18462 = -0.284496736f;
        float r18463 = r18461 + r18462;
        float r18464 = fma(r18457, r18459, r18463);
        float r18465 = log1p(r18464);
        float r18466 = expm1(r18465);
        float r18467 = 0.254829592f;
        float r18468 = r18467 / r18454;
        float r18469 = fma(r18466, r18459, r18468);
        float r18470 = r18453 * r18453;
        float r18471 = exp(r18470);
        float r18472 = r18469 / r18471;
        float r18473 = r18449 - r18472;
        float r18474 = exp(r18473);
        float r18475 = log(r18474);
        return r18475;
}

double f_od(double x) {
        double r18476 = 1.0;
        double r18477 = 1.061405429;
        double r18478 = 0.3275911;
        double r18479 = x;
        double r18480 = fabs(r18479);
        double r18481 = fma(r18478, r18480, r18476);
        double r18482 = r18477 / r18481;
        double r18483 = -1.453152027;
        double r18484 = r18482 + r18483;
        double r18485 = r18476 / r18481;
        double r18486 = r18485 / r18481;
        double r18487 = 1.421413741;
        double r18488 = r18487 / r18481;
        double r18489 = -0.284496736;
        double r18490 = r18488 + r18489;
        double r18491 = fma(r18484, r18486, r18490);
        double r18492 = log1p(r18491);
        double r18493 = expm1(r18492);
        double r18494 = 0.254829592;
        double r18495 = r18494 / r18481;
        double r18496 = fma(r18493, r18486, r18495);
        double r18497 = r18480 * r18480;
        double r18498 = exp(r18497);
        double r18499 = r18496 / r18498;
        double r18500 = r18476 - r18499;
        double r18501 = exp(r18500);
        double r18502 = log(r18501);
        return r18502;
}

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 r18503, r18504, r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513, r18514, r18515, r18516, r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525, r18526, r18527, r18528;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18503, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18504, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18505);
        mpfr_init(r18506);
        mpfr_init(r18507);
        mpfr_init(r18508);
        mpfr_init(r18509);
        mpfr_init_set_str(r18510, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18511, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18512, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18513, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18514, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18515);
        mpfr_init(r18516);
        mpfr_init(r18517);
        mpfr_init(r18518);
        mpfr_init(r18519);
        mpfr_init(r18520);
        mpfr_init(r18521);
        mpfr_init(r18522);
        mpfr_init(r18523);
        mpfr_init(r18524);
        mpfr_init(r18525);
        mpfr_init(r18526);
        mpfr_init(r18527);
        mpfr_init(r18528);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18505, x, MPFR_RNDN);
        mpfr_abs(r18506, r18505, MPFR_RNDN);
        mpfr_mul(r18507, r18504, r18506, MPFR_RNDN);
        mpfr_add(r18508, r18503, r18507, MPFR_RNDN);
        mpfr_div(r18509, r18503, r18508, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18515, r18509, r18514, MPFR_RNDN);
        mpfr_add(r18516, r18513, r18515, MPFR_RNDN);
        mpfr_mul(r18517, r18509, r18516, MPFR_RNDN);
        mpfr_add(r18518, r18512, r18517, MPFR_RNDN);
        mpfr_mul(r18519, r18509, r18518, MPFR_RNDN);
        mpfr_add(r18520, r18511, r18519, MPFR_RNDN);
        mpfr_mul(r18521, r18509, r18520, MPFR_RNDN);
        mpfr_add(r18522, r18510, r18521, MPFR_RNDN);
        mpfr_mul(r18523, r18509, r18522, MPFR_RNDN);
        mpfr_mul(r18524, r18506, r18506, MPFR_RNDN);
        mpfr_neg(r18525, r18524, MPFR_RNDN);
        mpfr_exp(r18526, r18525, MPFR_RNDN);
        mpfr_mul(r18527, r18523, r18526, MPFR_RNDN);
        mpfr_sub(r18528, r18503, r18527, MPFR_RNDN);
        return mpfr_get_d(r18528, MPFR_RNDN);
}

static mpfr_t r18529, r18530, r18531, r18532, r18533, r18534, r18535, r18536, r18537, r18538, r18539, r18540, r18541, r18542, r18543, r18544, r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18529, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18530, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18531, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18532);
        mpfr_init(r18533);
        mpfr_init(r18534);
        mpfr_init(r18535);
        mpfr_init_set_str(r18536, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18537);
        mpfr_init(r18538);
        mpfr_init(r18539);
        mpfr_init_set_str(r18540, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18541);
        mpfr_init_set_str(r18542, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18543);
        mpfr_init(r18544);
        mpfr_init(r18545);
        mpfr_init(r18546);
        mpfr_init_set_str(r18547, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18548);
        mpfr_init(r18549);
        mpfr_init(r18550);
        mpfr_init(r18551);
        mpfr_init(r18552);
        mpfr_init(r18553);
        mpfr_init(r18554);
        mpfr_init(r18555);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18532, x, MPFR_RNDN);
        mpfr_abs(r18533, r18532, MPFR_RNDN);
        mpfr_fma(r18534, r18531, r18533, r18529, MPFR_RNDN);
        mpfr_div(r18535, r18530, r18534, MPFR_RNDN);
        ;
        mpfr_add(r18537, r18535, r18536, MPFR_RNDN);
        mpfr_div(r18538, r18529, r18534, MPFR_RNDN);
        mpfr_div(r18539, r18538, r18534, MPFR_RNDN);
        ;
        mpfr_div(r18541, r18540, r18534, MPFR_RNDN);
        ;
        mpfr_add(r18543, r18541, r18542, MPFR_RNDN);
        mpfr_fma(r18544, r18537, r18539, r18543, MPFR_RNDN);
        mpfr_log1p(r18545, r18544, MPFR_RNDN);
        mpfr_expm1(r18546, r18545, MPFR_RNDN);
        ;
        mpfr_div(r18548, r18547, r18534, MPFR_RNDN);
        mpfr_fma(r18549, r18546, r18539, r18548, MPFR_RNDN);
        mpfr_mul(r18550, r18533, r18533, MPFR_RNDN);
        mpfr_exp(r18551, r18550, MPFR_RNDN);
        mpfr_div(r18552, r18549, r18551, MPFR_RNDN);
        mpfr_sub(r18553, r18529, r18552, MPFR_RNDN);
        mpfr_exp(r18554, r18553, MPFR_RNDN);
        mpfr_log(r18555, r18554, MPFR_RNDN);
        return mpfr_get_d(r18555, MPFR_RNDN);
}

static mpfr_t r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567, r18568, r18569, r18570, r18571, r18572, r18573, r18574, r18575, r18576, r18577, r18578, r18579, r18580, r18581, r18582;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18556, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18557, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18558, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18559);
        mpfr_init(r18560);
        mpfr_init(r18561);
        mpfr_init(r18562);
        mpfr_init_set_str(r18563, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18564);
        mpfr_init(r18565);
        mpfr_init(r18566);
        mpfr_init_set_str(r18567, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18568);
        mpfr_init_set_str(r18569, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18570);
        mpfr_init(r18571);
        mpfr_init(r18572);
        mpfr_init(r18573);
        mpfr_init_set_str(r18574, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18575);
        mpfr_init(r18576);
        mpfr_init(r18577);
        mpfr_init(r18578);
        mpfr_init(r18579);
        mpfr_init(r18580);
        mpfr_init(r18581);
        mpfr_init(r18582);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18559, x, MPFR_RNDN);
        mpfr_abs(r18560, r18559, MPFR_RNDN);
        mpfr_fma(r18561, r18558, r18560, r18556, MPFR_RNDN);
        mpfr_div(r18562, r18557, r18561, MPFR_RNDN);
        ;
        mpfr_add(r18564, r18562, r18563, MPFR_RNDN);
        mpfr_div(r18565, r18556, r18561, MPFR_RNDN);
        mpfr_div(r18566, r18565, r18561, MPFR_RNDN);
        ;
        mpfr_div(r18568, r18567, r18561, MPFR_RNDN);
        ;
        mpfr_add(r18570, r18568, r18569, MPFR_RNDN);
        mpfr_fma(r18571, r18564, r18566, r18570, MPFR_RNDN);
        mpfr_log1p(r18572, r18571, MPFR_RNDN);
        mpfr_expm1(r18573, r18572, MPFR_RNDN);
        ;
        mpfr_div(r18575, r18574, r18561, MPFR_RNDN);
        mpfr_fma(r18576, r18573, r18566, r18575, MPFR_RNDN);
        mpfr_mul(r18577, r18560, r18560, MPFR_RNDN);
        mpfr_exp(r18578, r18577, MPFR_RNDN);
        mpfr_div(r18579, r18576, r18578, MPFR_RNDN);
        mpfr_sub(r18580, r18556, r18579, MPFR_RNDN);
        mpfr_exp(r18581, r18580, MPFR_RNDN);
        mpfr_log(r18582, r18581, MPFR_RNDN);
        return mpfr_get_d(r18582, MPFR_RNDN);
}