#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 r18541 = 1.0f;
        float r18542 = 0.3275911f;
        float r18543 = x;
        float r18544 = fabs(r18543);
        float r18545 = r18542 * r18544;
        float r18546 = r18541 + r18545;
        float r18547 = r18541 / r18546;
        float r18548 = 0.254829592f;
        float r18549 = -0.284496736f;
        float r18550 = 1.421413741f;
        float r18551 = -1.453152027f;
        float r18552 = 1.061405429f;
        float r18553 = r18547 * r18552;
        float r18554 = r18551 + r18553;
        float r18555 = r18547 * r18554;
        float r18556 = r18550 + r18555;
        float r18557 = r18547 * r18556;
        float r18558 = r18549 + r18557;
        float r18559 = r18547 * r18558;
        float r18560 = r18548 + r18559;
        float r18561 = r18547 * r18560;
        float r18562 = r18544 * r18544;
        float r18563 = -r18562;
        float r18564 = exp(r18563);
        float r18565 = r18561 * r18564;
        float r18566 = r18541 - r18565;
        return r18566;
}

double f_id(double x) {
        double r18567 = 1.0;
        double r18568 = 0.3275911;
        double r18569 = x;
        double r18570 = fabs(r18569);
        double r18571 = r18568 * r18570;
        double r18572 = r18567 + r18571;
        double r18573 = r18567 / r18572;
        double r18574 = 0.254829592;
        double r18575 = -0.284496736;
        double r18576 = 1.421413741;
        double r18577 = -1.453152027;
        double r18578 = 1.061405429;
        double r18579 = r18573 * r18578;
        double r18580 = r18577 + r18579;
        double r18581 = r18573 * r18580;
        double r18582 = r18576 + r18581;
        double r18583 = r18573 * r18582;
        double r18584 = r18575 + r18583;
        double r18585 = r18573 * r18584;
        double r18586 = r18574 + r18585;
        double r18587 = r18573 * r18586;
        double r18588 = r18570 * r18570;
        double r18589 = -r18588;
        double r18590 = exp(r18589);
        double r18591 = r18587 * r18590;
        double r18592 = r18567 - r18591;
        return r18592;
}


double f_of(float x) {
        float r18593 = 1.0f;
        float r18594 = 1.061405429f;
        float r18595 = 0.3275911f;
        float r18596 = x;
        float r18597 = fabs(r18596);
        float r18598 = fma(r18595, r18597, r18593);
        float r18599 = r18594 / r18598;
        float r18600 = -1.453152027f;
        float r18601 = r18599 + r18600;
        float r18602 = r18593 / r18598;
        float r18603 = r18602 / r18598;
        float r18604 = 1.421413741f;
        float r18605 = r18604 / r18598;
        float r18606 = -0.284496736f;
        float r18607 = r18605 + r18606;
        float r18608 = fma(r18601, r18603, r18607);
        float r18609 = 0.254829592f;
        float r18610 = r18609 / r18598;
        float r18611 = cbrt(r18610);
        float r18612 = r18611 * (r18611 * r18611);
        float r18613 = fma(r18608, r18603, r18612);
        float r18614 = r18597 * r18597;
        float r18615 = exp(r18614);
        float r18616 = r18613 / r18615;
        float r18617 = r18593 - r18616;
        return r18617;
}

double f_od(double x) {
        double r18618 = 1.0;
        double r18619 = 1.061405429;
        double r18620 = 0.3275911;
        double r18621 = x;
        double r18622 = fabs(r18621);
        double r18623 = fma(r18620, r18622, r18618);
        double r18624 = r18619 / r18623;
        double r18625 = -1.453152027;
        double r18626 = r18624 + r18625;
        double r18627 = r18618 / r18623;
        double r18628 = r18627 / r18623;
        double r18629 = 1.421413741;
        double r18630 = r18629 / r18623;
        double r18631 = -0.284496736;
        double r18632 = r18630 + r18631;
        double r18633 = fma(r18626, r18628, r18632);
        double r18634 = 0.254829592;
        double r18635 = r18634 / r18623;
        double r18636 = cbrt(r18635);
        double r18637 = r18636 * (r18636 * r18636);
        double r18638 = fma(r18633, r18628, r18637);
        double r18639 = r18622 * r18622;
        double r18640 = exp(r18639);
        double r18641 = r18638 / r18640;
        double r18642 = r18618 - r18641;
        return r18642;
}

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 r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667, r18668;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18643, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18644, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18645);
        mpfr_init(r18646);
        mpfr_init(r18647);
        mpfr_init(r18648);
        mpfr_init(r18649);
        mpfr_init_set_str(r18650, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18651, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18652, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18653, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18654, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18655);
        mpfr_init(r18656);
        mpfr_init(r18657);
        mpfr_init(r18658);
        mpfr_init(r18659);
        mpfr_init(r18660);
        mpfr_init(r18661);
        mpfr_init(r18662);
        mpfr_init(r18663);
        mpfr_init(r18664);
        mpfr_init(r18665);
        mpfr_init(r18666);
        mpfr_init(r18667);
        mpfr_init(r18668);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18645, x, MPFR_RNDN);
        mpfr_abs(r18646, r18645, MPFR_RNDN);
        mpfr_mul(r18647, r18644, r18646, MPFR_RNDN);
        mpfr_add(r18648, r18643, r18647, MPFR_RNDN);
        mpfr_div(r18649, r18643, r18648, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18655, r18649, r18654, MPFR_RNDN);
        mpfr_add(r18656, r18653, r18655, MPFR_RNDN);
        mpfr_mul(r18657, r18649, r18656, MPFR_RNDN);
        mpfr_add(r18658, r18652, r18657, MPFR_RNDN);
        mpfr_mul(r18659, r18649, r18658, MPFR_RNDN);
        mpfr_add(r18660, r18651, r18659, MPFR_RNDN);
        mpfr_mul(r18661, r18649, r18660, MPFR_RNDN);
        mpfr_add(r18662, r18650, r18661, MPFR_RNDN);
        mpfr_mul(r18663, r18649, r18662, MPFR_RNDN);
        mpfr_mul(r18664, r18646, r18646, MPFR_RNDN);
        mpfr_neg(r18665, r18664, MPFR_RNDN);
        mpfr_exp(r18666, r18665, MPFR_RNDN);
        mpfr_mul(r18667, r18663, r18666, MPFR_RNDN);
        mpfr_sub(r18668, r18643, r18667, MPFR_RNDN);
        return mpfr_get_d(r18668, MPFR_RNDN);
}

static mpfr_t r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681, r18682, r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18669, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18670, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18671, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18672);
        mpfr_init(r18673);
        mpfr_init(r18674);
        mpfr_init(r18675);
        mpfr_init_set_str(r18676, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18677);
        mpfr_init(r18678);
        mpfr_init(r18679);
        mpfr_init_set_str(r18680, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18681);
        mpfr_init_set_str(r18682, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18683);
        mpfr_init(r18684);
        mpfr_init_set_str(r18685, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18686);
        mpfr_init(r18687);
        mpfr_init(r18688);
        mpfr_init(r18689);
        mpfr_init(r18690);
        mpfr_init(r18691);
        mpfr_init(r18692);
        mpfr_init(r18693);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18672, x, MPFR_RNDN);
        mpfr_abs(r18673, r18672, MPFR_RNDN);
        mpfr_fma(r18674, r18671, r18673, r18669, MPFR_RNDN);
        mpfr_div(r18675, r18670, r18674, MPFR_RNDN);
        ;
        mpfr_add(r18677, r18675, r18676, MPFR_RNDN);
        mpfr_div(r18678, r18669, r18674, MPFR_RNDN);
        mpfr_div(r18679, r18678, r18674, MPFR_RNDN);
        ;
        mpfr_div(r18681, r18680, r18674, MPFR_RNDN);
        ;
        mpfr_add(r18683, r18681, r18682, MPFR_RNDN);
        mpfr_fma(r18684, r18677, r18679, r18683, MPFR_RNDN);
        ;
        mpfr_div(r18686, r18685, r18674, MPFR_RNDN);
        mpfr_cbrt(r18687, r18686, MPFR_RNDN);
        mpfr_mul(r18688, r18687, r18687, MPFR_RNDN); mpfr_mul(r18688, r18688, r18687, MPFR_RNDN);
        mpfr_fma(r18689, r18684, r18679, r18688, MPFR_RNDN);
        mpfr_mul(r18690, r18673, r18673, MPFR_RNDN);
        mpfr_exp(r18691, r18690, MPFR_RNDN);
        mpfr_div(r18692, r18689, r18691, MPFR_RNDN);
        mpfr_sub(r18693, r18669, r18692, MPFR_RNDN);
        return mpfr_get_d(r18693, MPFR_RNDN);
}

static mpfr_t r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703, r18704, r18705, r18706, r18707, r18708, r18709, r18710, r18711, r18712, r18713, r18714, r18715, r18716, r18717, r18718;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18694, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18695, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18696, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18697);
        mpfr_init(r18698);
        mpfr_init(r18699);
        mpfr_init(r18700);
        mpfr_init_set_str(r18701, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18702);
        mpfr_init(r18703);
        mpfr_init(r18704);
        mpfr_init_set_str(r18705, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18706);
        mpfr_init_set_str(r18707, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18708);
        mpfr_init(r18709);
        mpfr_init_set_str(r18710, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18711);
        mpfr_init(r18712);
        mpfr_init(r18713);
        mpfr_init(r18714);
        mpfr_init(r18715);
        mpfr_init(r18716);
        mpfr_init(r18717);
        mpfr_init(r18718);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18697, x, MPFR_RNDN);
        mpfr_abs(r18698, r18697, MPFR_RNDN);
        mpfr_fma(r18699, r18696, r18698, r18694, MPFR_RNDN);
        mpfr_div(r18700, r18695, r18699, MPFR_RNDN);
        ;
        mpfr_add(r18702, r18700, r18701, MPFR_RNDN);
        mpfr_div(r18703, r18694, r18699, MPFR_RNDN);
        mpfr_div(r18704, r18703, r18699, MPFR_RNDN);
        ;
        mpfr_div(r18706, r18705, r18699, MPFR_RNDN);
        ;
        mpfr_add(r18708, r18706, r18707, MPFR_RNDN);
        mpfr_fma(r18709, r18702, r18704, r18708, MPFR_RNDN);
        ;
        mpfr_div(r18711, r18710, r18699, MPFR_RNDN);
        mpfr_cbrt(r18712, r18711, MPFR_RNDN);
        mpfr_mul(r18713, r18712, r18712, MPFR_RNDN); mpfr_mul(r18713, r18713, r18712, MPFR_RNDN);
        mpfr_fma(r18714, r18709, r18704, r18713, MPFR_RNDN);
        mpfr_mul(r18715, r18698, r18698, MPFR_RNDN);
        mpfr_exp(r18716, r18715, MPFR_RNDN);
        mpfr_div(r18717, r18714, r18716, MPFR_RNDN);
        mpfr_sub(r18718, r18694, r18717, MPFR_RNDN);
        return mpfr_get_d(r18718, MPFR_RNDN);
}