#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 r18529 = 1.0f;
        float r18530 = 0.3275911f;
        float r18531 = x;
        float r18532 = fabs(r18531);
        float r18533 = r18530 * r18532;
        float r18534 = r18529 + r18533;
        float r18535 = r18529 / r18534;
        float r18536 = 0.254829592f;
        float r18537 = -0.284496736f;
        float r18538 = 1.421413741f;
        float r18539 = -1.453152027f;
        float r18540 = 1.061405429f;
        float r18541 = r18535 * r18540;
        float r18542 = r18539 + r18541;
        float r18543 = r18535 * r18542;
        float r18544 = r18538 + r18543;
        float r18545 = r18535 * r18544;
        float r18546 = r18537 + r18545;
        float r18547 = r18535 * r18546;
        float r18548 = r18536 + r18547;
        float r18549 = r18535 * r18548;
        float r18550 = r18532 * r18532;
        float r18551 = -r18550;
        float r18552 = exp(r18551);
        float r18553 = r18549 * r18552;
        float r18554 = r18529 - r18553;
        return r18554;
}

double f_id(double x) {
        double r18555 = 1.0;
        double r18556 = 0.3275911;
        double r18557 = x;
        double r18558 = fabs(r18557);
        double r18559 = r18556 * r18558;
        double r18560 = r18555 + r18559;
        double r18561 = r18555 / r18560;
        double r18562 = 0.254829592;
        double r18563 = -0.284496736;
        double r18564 = 1.421413741;
        double r18565 = -1.453152027;
        double r18566 = 1.061405429;
        double r18567 = r18561 * r18566;
        double r18568 = r18565 + r18567;
        double r18569 = r18561 * r18568;
        double r18570 = r18564 + r18569;
        double r18571 = r18561 * r18570;
        double r18572 = r18563 + r18571;
        double r18573 = r18561 * r18572;
        double r18574 = r18562 + r18573;
        double r18575 = r18561 * r18574;
        double r18576 = r18558 * r18558;
        double r18577 = -r18576;
        double r18578 = exp(r18577);
        double r18579 = r18575 * r18578;
        double r18580 = r18555 - r18579;
        return r18580;
}


double f_of(float x) {
        float r18581 = 1.0f;
        float r18582 = 1.061405429f;
        float r18583 = 0.3275911f;
        float r18584 = x;
        float r18585 = fabs(r18584);
        float r18586 = fma(r18583, r18585, r18581);
        float r18587 = r18582 / r18586;
        float r18588 = -1.453152027f;
        float r18589 = r18587 + r18588;
        float r18590 = r18581 / r18586;
        float r18591 = r18590 / r18586;
        float r18592 = 1.421413741f;
        float r18593 = r18592 / r18586;
        float r18594 = -0.284496736f;
        float r18595 = r18593 + r18594;
        float r18596 = fma(r18589, r18591, r18595);
        float r18597 = 0.254829592f;
        float r18598 = r18597 / r18586;
        float r18599 = cbrt(r18598);
        float r18600 = r18599 * (r18599 * r18599);
        float r18601 = fma(r18596, r18591, r18600);
        float r18602 = r18585 * r18585;
        float r18603 = exp(r18602);
        float r18604 = r18601 / r18603;
        float r18605 = r18581 - r18604;
        return r18605;
}

double f_od(double x) {
        double r18606 = 1.0;
        double r18607 = 1.061405429;
        double r18608 = 0.3275911;
        double r18609 = x;
        double r18610 = fabs(r18609);
        double r18611 = fma(r18608, r18610, r18606);
        double r18612 = r18607 / r18611;
        double r18613 = -1.453152027;
        double r18614 = r18612 + r18613;
        double r18615 = r18606 / r18611;
        double r18616 = r18615 / r18611;
        double r18617 = 1.421413741;
        double r18618 = r18617 / r18611;
        double r18619 = -0.284496736;
        double r18620 = r18618 + r18619;
        double r18621 = fma(r18614, r18616, r18620);
        double r18622 = 0.254829592;
        double r18623 = r18622 / r18611;
        double r18624 = cbrt(r18623);
        double r18625 = r18624 * (r18624 * r18624);
        double r18626 = fma(r18621, r18616, r18625);
        double r18627 = r18610 * r18610;
        double r18628 = exp(r18627);
        double r18629 = r18626 / r18628;
        double r18630 = r18606 - r18629;
        return r18630;
}

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 r18631, r18632, r18633, r18634, r18635, r18636, r18637, r18638, r18639, r18640, r18641, r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18631, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18632, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18633);
        mpfr_init(r18634);
        mpfr_init(r18635);
        mpfr_init(r18636);
        mpfr_init(r18637);
        mpfr_init_set_str(r18638, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18639, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18640, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18641, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18642, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18643);
        mpfr_init(r18644);
        mpfr_init(r18645);
        mpfr_init(r18646);
        mpfr_init(r18647);
        mpfr_init(r18648);
        mpfr_init(r18649);
        mpfr_init(r18650);
        mpfr_init(r18651);
        mpfr_init(r18652);
        mpfr_init(r18653);
        mpfr_init(r18654);
        mpfr_init(r18655);
        mpfr_init(r18656);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18633, x, MPFR_RNDN);
        mpfr_abs(r18634, r18633, MPFR_RNDN);
        mpfr_mul(r18635, r18632, r18634, MPFR_RNDN);
        mpfr_add(r18636, r18631, r18635, MPFR_RNDN);
        mpfr_div(r18637, r18631, r18636, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18643, r18637, r18642, MPFR_RNDN);
        mpfr_add(r18644, r18641, r18643, MPFR_RNDN);
        mpfr_mul(r18645, r18637, r18644, MPFR_RNDN);
        mpfr_add(r18646, r18640, r18645, MPFR_RNDN);
        mpfr_mul(r18647, r18637, r18646, MPFR_RNDN);
        mpfr_add(r18648, r18639, r18647, MPFR_RNDN);
        mpfr_mul(r18649, r18637, r18648, MPFR_RNDN);
        mpfr_add(r18650, r18638, r18649, MPFR_RNDN);
        mpfr_mul(r18651, r18637, r18650, MPFR_RNDN);
        mpfr_mul(r18652, r18634, r18634, MPFR_RNDN);
        mpfr_neg(r18653, r18652, MPFR_RNDN);
        mpfr_exp(r18654, r18653, MPFR_RNDN);
        mpfr_mul(r18655, r18651, r18654, MPFR_RNDN);
        mpfr_sub(r18656, r18631, r18655, MPFR_RNDN);
        return mpfr_get_d(r18656, MPFR_RNDN);
}

static mpfr_t r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667, r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18657, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18658, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18659, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18660);
        mpfr_init(r18661);
        mpfr_init(r18662);
        mpfr_init(r18663);
        mpfr_init_set_str(r18664, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18665);
        mpfr_init(r18666);
        mpfr_init(r18667);
        mpfr_init_set_str(r18668, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18669);
        mpfr_init_set_str(r18670, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18671);
        mpfr_init(r18672);
        mpfr_init_set_str(r18673, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18674);
        mpfr_init(r18675);
        mpfr_init(r18676);
        mpfr_init(r18677);
        mpfr_init(r18678);
        mpfr_init(r18679);
        mpfr_init(r18680);
        mpfr_init(r18681);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18660, x, MPFR_RNDN);
        mpfr_abs(r18661, r18660, MPFR_RNDN);
        mpfr_fma(r18662, r18659, r18661, r18657, MPFR_RNDN);
        mpfr_div(r18663, r18658, r18662, MPFR_RNDN);
        ;
        mpfr_add(r18665, r18663, r18664, MPFR_RNDN);
        mpfr_div(r18666, r18657, r18662, MPFR_RNDN);
        mpfr_div(r18667, r18666, r18662, MPFR_RNDN);
        ;
        mpfr_div(r18669, r18668, r18662, MPFR_RNDN);
        ;
        mpfr_add(r18671, r18669, r18670, MPFR_RNDN);
        mpfr_fma(r18672, r18665, r18667, r18671, MPFR_RNDN);
        ;
        mpfr_div(r18674, r18673, r18662, MPFR_RNDN);
        mpfr_cbrt(r18675, r18674, MPFR_RNDN);
        mpfr_mul(r18676, r18675, r18675, MPFR_RNDN); mpfr_mul(r18676, r18676, r18675, MPFR_RNDN);
        mpfr_fma(r18677, r18672, r18667, r18676, MPFR_RNDN);
        mpfr_mul(r18678, r18661, r18661, MPFR_RNDN);
        mpfr_exp(r18679, r18678, MPFR_RNDN);
        mpfr_div(r18680, r18677, r18679, MPFR_RNDN);
        mpfr_sub(r18681, r18657, r18680, MPFR_RNDN);
        return mpfr_get_d(r18681, MPFR_RNDN);
}

static mpfr_t r18682, r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703, r18704, r18705, r18706;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18682, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18683, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18684, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18685);
        mpfr_init(r18686);
        mpfr_init(r18687);
        mpfr_init(r18688);
        mpfr_init_set_str(r18689, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18690);
        mpfr_init(r18691);
        mpfr_init(r18692);
        mpfr_init_set_str(r18693, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18694);
        mpfr_init_set_str(r18695, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18696);
        mpfr_init(r18697);
        mpfr_init_set_str(r18698, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18699);
        mpfr_init(r18700);
        mpfr_init(r18701);
        mpfr_init(r18702);
        mpfr_init(r18703);
        mpfr_init(r18704);
        mpfr_init(r18705);
        mpfr_init(r18706);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18685, x, MPFR_RNDN);
        mpfr_abs(r18686, r18685, MPFR_RNDN);
        mpfr_fma(r18687, r18684, r18686, r18682, MPFR_RNDN);
        mpfr_div(r18688, r18683, r18687, MPFR_RNDN);
        ;
        mpfr_add(r18690, r18688, r18689, MPFR_RNDN);
        mpfr_div(r18691, r18682, r18687, MPFR_RNDN);
        mpfr_div(r18692, r18691, r18687, MPFR_RNDN);
        ;
        mpfr_div(r18694, r18693, r18687, MPFR_RNDN);
        ;
        mpfr_add(r18696, r18694, r18695, MPFR_RNDN);
        mpfr_fma(r18697, r18690, r18692, r18696, MPFR_RNDN);
        ;
        mpfr_div(r18699, r18698, r18687, MPFR_RNDN);
        mpfr_cbrt(r18700, r18699, MPFR_RNDN);
        mpfr_mul(r18701, r18700, r18700, MPFR_RNDN); mpfr_mul(r18701, r18701, r18700, MPFR_RNDN);
        mpfr_fma(r18702, r18697, r18692, r18701, MPFR_RNDN);
        mpfr_mul(r18703, r18686, r18686, MPFR_RNDN);
        mpfr_exp(r18704, r18703, MPFR_RNDN);
        mpfr_div(r18705, r18702, r18704, MPFR_RNDN);
        mpfr_sub(r18706, r18682, r18705, MPFR_RNDN);
        return mpfr_get_d(r18706, MPFR_RNDN);
}