#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 r18409 = 1.0f;
        float r18410 = 0.3275911f;
        float r18411 = x;
        float r18412 = fabs(r18411);
        float r18413 = r18410 * r18412;
        float r18414 = r18409 + r18413;
        float r18415 = r18409 / r18414;
        float r18416 = 0.254829592f;
        float r18417 = -0.284496736f;
        float r18418 = 1.421413741f;
        float r18419 = -1.453152027f;
        float r18420 = 1.061405429f;
        float r18421 = r18415 * r18420;
        float r18422 = r18419 + r18421;
        float r18423 = r18415 * r18422;
        float r18424 = r18418 + r18423;
        float r18425 = r18415 * r18424;
        float r18426 = r18417 + r18425;
        float r18427 = r18415 * r18426;
        float r18428 = r18416 + r18427;
        float r18429 = r18415 * r18428;
        float r18430 = r18412 * r18412;
        float r18431 = -r18430;
        float r18432 = exp(r18431);
        float r18433 = r18429 * r18432;
        float r18434 = r18409 - r18433;
        return r18434;
}

double f_id(double x) {
        double r18435 = 1.0;
        double r18436 = 0.3275911;
        double r18437 = x;
        double r18438 = fabs(r18437);
        double r18439 = r18436 * r18438;
        double r18440 = r18435 + r18439;
        double r18441 = r18435 / r18440;
        double r18442 = 0.254829592;
        double r18443 = -0.284496736;
        double r18444 = 1.421413741;
        double r18445 = -1.453152027;
        double r18446 = 1.061405429;
        double r18447 = r18441 * r18446;
        double r18448 = r18445 + r18447;
        double r18449 = r18441 * r18448;
        double r18450 = r18444 + r18449;
        double r18451 = r18441 * r18450;
        double r18452 = r18443 + r18451;
        double r18453 = r18441 * r18452;
        double r18454 = r18442 + r18453;
        double r18455 = r18441 * r18454;
        double r18456 = r18438 * r18438;
        double r18457 = -r18456;
        double r18458 = exp(r18457);
        double r18459 = r18455 * r18458;
        double r18460 = r18435 - r18459;
        return r18460;
}


double f_of(float x) {
        float r18461 = 1.0f;
        float r18462 = 1.061405429f;
        float r18463 = 0.3275911f;
        float r18464 = x;
        float r18465 = fabs(r18464);
        float r18466 = fma(r18463, r18465, r18461);
        float r18467 = r18462 / r18466;
        float r18468 = -1.453152027f;
        float r18469 = r18467 + r18468;
        float r18470 = r18461 / r18466;
        float r18471 = r18470 / r18466;
        float r18472 = 1.421413741f;
        float r18473 = r18472 / r18466;
        float r18474 = -0.284496736f;
        float r18475 = r18473 + r18474;
        float r18476 = fma(r18469, r18471, r18475);
        float r18477 = cbrt(r18470);
        float r18478 = exp(r18477);
        float r18479 = log(r18478);
        float r18480 = r18479 * (r18479 * r18479);
        float r18481 = r18480 / r18466;
        float r18482 = 0.254829592f;
        float r18483 = r18482 / r18466;
        float r18484 = fma(r18476, r18481, r18483);
        float r18485 = r18465 * r18465;
        float r18486 = exp(r18485);
        float r18487 = r18484 / r18486;
        float r18488 = r18461 - r18487;
        return r18488;
}

double f_od(double x) {
        double r18489 = 1.0;
        double r18490 = 1.061405429;
        double r18491 = 0.3275911;
        double r18492 = x;
        double r18493 = fabs(r18492);
        double r18494 = fma(r18491, r18493, r18489);
        double r18495 = r18490 / r18494;
        double r18496 = -1.453152027;
        double r18497 = r18495 + r18496;
        double r18498 = r18489 / r18494;
        double r18499 = r18498 / r18494;
        double r18500 = 1.421413741;
        double r18501 = r18500 / r18494;
        double r18502 = -0.284496736;
        double r18503 = r18501 + r18502;
        double r18504 = fma(r18497, r18499, r18503);
        double r18505 = cbrt(r18498);
        double r18506 = exp(r18505);
        double r18507 = log(r18506);
        double r18508 = r18507 * (r18507 * r18507);
        double r18509 = r18508 / r18494;
        double r18510 = 0.254829592;
        double r18511 = r18510 / r18494;
        double r18512 = fma(r18504, r18509, r18511);
        double r18513 = r18493 * r18493;
        double r18514 = exp(r18513);
        double r18515 = r18512 / r18514;
        double r18516 = r18489 - r18515;
        return r18516;
}

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 r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525, r18526, r18527, r18528, r18529, r18530, r18531, r18532, r18533, r18534, r18535, r18536, r18537, r18538, r18539, r18540, r18541, r18542;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18517, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18518, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18519);
        mpfr_init(r18520);
        mpfr_init(r18521);
        mpfr_init(r18522);
        mpfr_init(r18523);
        mpfr_init_set_str(r18524, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18525, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18526, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18527, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18528, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18529);
        mpfr_init(r18530);
        mpfr_init(r18531);
        mpfr_init(r18532);
        mpfr_init(r18533);
        mpfr_init(r18534);
        mpfr_init(r18535);
        mpfr_init(r18536);
        mpfr_init(r18537);
        mpfr_init(r18538);
        mpfr_init(r18539);
        mpfr_init(r18540);
        mpfr_init(r18541);
        mpfr_init(r18542);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18519, x, MPFR_RNDN);
        mpfr_abs(r18520, r18519, MPFR_RNDN);
        mpfr_mul(r18521, r18518, r18520, MPFR_RNDN);
        mpfr_add(r18522, r18517, r18521, MPFR_RNDN);
        mpfr_div(r18523, r18517, r18522, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18529, r18523, r18528, MPFR_RNDN);
        mpfr_add(r18530, r18527, r18529, MPFR_RNDN);
        mpfr_mul(r18531, r18523, r18530, MPFR_RNDN);
        mpfr_add(r18532, r18526, r18531, MPFR_RNDN);
        mpfr_mul(r18533, r18523, r18532, MPFR_RNDN);
        mpfr_add(r18534, r18525, r18533, MPFR_RNDN);
        mpfr_mul(r18535, r18523, r18534, MPFR_RNDN);
        mpfr_add(r18536, r18524, r18535, MPFR_RNDN);
        mpfr_mul(r18537, r18523, r18536, MPFR_RNDN);
        mpfr_mul(r18538, r18520, r18520, MPFR_RNDN);
        mpfr_neg(r18539, r18538, MPFR_RNDN);
        mpfr_exp(r18540, r18539, MPFR_RNDN);
        mpfr_mul(r18541, r18537, r18540, MPFR_RNDN);
        mpfr_sub(r18542, r18517, r18541, MPFR_RNDN);
        return mpfr_get_d(r18542, MPFR_RNDN);
}

static mpfr_t r18543, r18544, r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567, r18568, r18569, r18570;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18543, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18544, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18545, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18546);
        mpfr_init(r18547);
        mpfr_init(r18548);
        mpfr_init(r18549);
        mpfr_init_set_str(r18550, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18551);
        mpfr_init(r18552);
        mpfr_init(r18553);
        mpfr_init_set_str(r18554, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18555);
        mpfr_init_set_str(r18556, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18557);
        mpfr_init(r18558);
        mpfr_init(r18559);
        mpfr_init(r18560);
        mpfr_init(r18561);
        mpfr_init(r18562);
        mpfr_init(r18563);
        mpfr_init_set_str(r18564, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18565);
        mpfr_init(r18566);
        mpfr_init(r18567);
        mpfr_init(r18568);
        mpfr_init(r18569);
        mpfr_init(r18570);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18546, x, MPFR_RNDN);
        mpfr_abs(r18547, r18546, MPFR_RNDN);
        mpfr_fma(r18548, r18545, r18547, r18543, MPFR_RNDN);
        mpfr_div(r18549, r18544, r18548, MPFR_RNDN);
        ;
        mpfr_add(r18551, r18549, r18550, MPFR_RNDN);
        mpfr_div(r18552, r18543, r18548, MPFR_RNDN);
        mpfr_div(r18553, r18552, r18548, MPFR_RNDN);
        ;
        mpfr_div(r18555, r18554, r18548, MPFR_RNDN);
        ;
        mpfr_add(r18557, r18555, r18556, MPFR_RNDN);
        mpfr_fma(r18558, r18551, r18553, r18557, MPFR_RNDN);
        mpfr_cbrt(r18559, r18552, MPFR_RNDN);
        mpfr_exp(r18560, r18559, MPFR_RNDN);
        mpfr_log(r18561, r18560, MPFR_RNDN);
        mpfr_mul(r18562, r18561, r18561, MPFR_RNDN); mpfr_mul(r18562, r18562, r18561, MPFR_RNDN);
        mpfr_div(r18563, r18562, r18548, MPFR_RNDN);
        ;
        mpfr_div(r18565, r18564, r18548, MPFR_RNDN);
        mpfr_fma(r18566, r18558, r18563, r18565, MPFR_RNDN);
        mpfr_mul(r18567, r18547, r18547, MPFR_RNDN);
        mpfr_exp(r18568, r18567, MPFR_RNDN);
        mpfr_div(r18569, r18566, r18568, MPFR_RNDN);
        mpfr_sub(r18570, r18543, r18569, MPFR_RNDN);
        return mpfr_get_d(r18570, MPFR_RNDN);
}

static mpfr_t r18571, r18572, r18573, r18574, r18575, r18576, r18577, r18578, r18579, r18580, r18581, r18582, r18583, r18584, r18585, r18586, r18587, r18588, r18589, r18590, r18591, r18592, r18593, r18594, r18595, r18596, r18597, r18598;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18571, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18572, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18573, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18574);
        mpfr_init(r18575);
        mpfr_init(r18576);
        mpfr_init(r18577);
        mpfr_init_set_str(r18578, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18579);
        mpfr_init(r18580);
        mpfr_init(r18581);
        mpfr_init_set_str(r18582, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18583);
        mpfr_init_set_str(r18584, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18585);
        mpfr_init(r18586);
        mpfr_init(r18587);
        mpfr_init(r18588);
        mpfr_init(r18589);
        mpfr_init(r18590);
        mpfr_init(r18591);
        mpfr_init_set_str(r18592, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18593);
        mpfr_init(r18594);
        mpfr_init(r18595);
        mpfr_init(r18596);
        mpfr_init(r18597);
        mpfr_init(r18598);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18574, x, MPFR_RNDN);
        mpfr_abs(r18575, r18574, MPFR_RNDN);
        mpfr_fma(r18576, r18573, r18575, r18571, MPFR_RNDN);
        mpfr_div(r18577, r18572, r18576, MPFR_RNDN);
        ;
        mpfr_add(r18579, r18577, r18578, MPFR_RNDN);
        mpfr_div(r18580, r18571, r18576, MPFR_RNDN);
        mpfr_div(r18581, r18580, r18576, MPFR_RNDN);
        ;
        mpfr_div(r18583, r18582, r18576, MPFR_RNDN);
        ;
        mpfr_add(r18585, r18583, r18584, MPFR_RNDN);
        mpfr_fma(r18586, r18579, r18581, r18585, MPFR_RNDN);
        mpfr_cbrt(r18587, r18580, MPFR_RNDN);
        mpfr_exp(r18588, r18587, MPFR_RNDN);
        mpfr_log(r18589, r18588, MPFR_RNDN);
        mpfr_mul(r18590, r18589, r18589, MPFR_RNDN); mpfr_mul(r18590, r18590, r18589, MPFR_RNDN);
        mpfr_div(r18591, r18590, r18576, MPFR_RNDN);
        ;
        mpfr_div(r18593, r18592, r18576, MPFR_RNDN);
        mpfr_fma(r18594, r18586, r18591, r18593, MPFR_RNDN);
        mpfr_mul(r18595, r18575, r18575, MPFR_RNDN);
        mpfr_exp(r18596, r18595, MPFR_RNDN);
        mpfr_div(r18597, r18594, r18596, MPFR_RNDN);
        mpfr_sub(r18598, r18571, r18597, MPFR_RNDN);
        return mpfr_get_d(r18598, MPFR_RNDN);
}