#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 r10519 = 1.0f;
        float r10520 = 0.3275911f;
        float r10521 = x;
        float r10522 = fabs(r10521);
        float r10523 = r10520 * r10522;
        float r10524 = r10519 + r10523;
        float r10525 = r10519 / r10524;
        float r10526 = 0.254829592f;
        float r10527 = -0.284496736f;
        float r10528 = 1.421413741f;
        float r10529 = -1.453152027f;
        float r10530 = 1.061405429f;
        float r10531 = r10525 * r10530;
        float r10532 = r10529 + r10531;
        float r10533 = r10525 * r10532;
        float r10534 = r10528 + r10533;
        float r10535 = r10525 * r10534;
        float r10536 = r10527 + r10535;
        float r10537 = r10525 * r10536;
        float r10538 = r10526 + r10537;
        float r10539 = r10525 * r10538;
        float r10540 = r10522 * r10522;
        float r10541 = -r10540;
        float r10542 = exp(r10541);
        float r10543 = r10539 * r10542;
        float r10544 = r10519 - r10543;
        return r10544;
}

double f_id(double x) {
        double r10545 = 1.0;
        double r10546 = 0.3275911;
        double r10547 = x;
        double r10548 = fabs(r10547);
        double r10549 = r10546 * r10548;
        double r10550 = r10545 + r10549;
        double r10551 = r10545 / r10550;
        double r10552 = 0.254829592;
        double r10553 = -0.284496736;
        double r10554 = 1.421413741;
        double r10555 = -1.453152027;
        double r10556 = 1.061405429;
        double r10557 = r10551 * r10556;
        double r10558 = r10555 + r10557;
        double r10559 = r10551 * r10558;
        double r10560 = r10554 + r10559;
        double r10561 = r10551 * r10560;
        double r10562 = r10553 + r10561;
        double r10563 = r10551 * r10562;
        double r10564 = r10552 + r10563;
        double r10565 = r10551 * r10564;
        double r10566 = r10548 * r10548;
        double r10567 = -r10566;
        double r10568 = exp(r10567);
        double r10569 = r10565 * r10568;
        double r10570 = r10545 - r10569;
        return r10570;
}


double f_of(float x) {
        float r10571 = 1.0f;
        float r10572 = x;
        float r10573 = fabs(r10572);
        float r10574 = 0.3275911f;
        float r10575 = r10573 * r10574;
        float r10576 = r10571 + r10575;
        float r10577 = r10571 / r10576;
        float r10578 = 0.254829592f;
        float r10579 = -0.284496736f;
        float r10580 = r10571 - r10575;
        float r10581 = 1.061405429f;
        float r10582 = r10577 * r10581;
        float r10583 = -1.453152027f;
        float r10584 = r10582 + r10583;
        float r10585 = r10577 * r10584;
        float r10586 = 1.421413741f;
        float r10587 = r10585 + r10586;
        float r10588 = r10580 * r10587;
        float r10589 = -r10574;
        float r10590 = r10589 * r10573;
        float r10591 = fma(r10575, r10590, r10571);
        float r10592 = r10571 / r10591;
        float r10593 = cbrt(r10592);
        float r10594 = r10593 * r10593;
        float r10595 = r10594 * r10593;
        float r10596 = r10588 * r10595;
        float r10597 = r10579 + r10596;
        float r10598 = r10597 * r10577;
        float r10599 = r10578 + r10598;
        float r10600 = r10577 * r10599;
        float r10601 = -r10573;
        float r10602 = r10601 * r10573;
        float r10603 = exp(r10602);
        float r10604 = r10600 * r10603;
        float r10605 = r10571 - r10604;
        return r10605;
}

double f_od(double x) {
        double r10606 = 1.0;
        double r10607 = x;
        double r10608 = fabs(r10607);
        double r10609 = 0.3275911;
        double r10610 = r10608 * r10609;
        double r10611 = r10606 + r10610;
        double r10612 = r10606 / r10611;
        double r10613 = 0.254829592;
        double r10614 = -0.284496736;
        double r10615 = r10606 - r10610;
        double r10616 = 1.061405429;
        double r10617 = r10612 * r10616;
        double r10618 = -1.453152027;
        double r10619 = r10617 + r10618;
        double r10620 = r10612 * r10619;
        double r10621 = 1.421413741;
        double r10622 = r10620 + r10621;
        double r10623 = r10615 * r10622;
        double r10624 = -r10609;
        double r10625 = r10624 * r10608;
        double r10626 = fma(r10610, r10625, r10606);
        double r10627 = r10606 / r10626;
        double r10628 = cbrt(r10627);
        double r10629 = r10628 * r10628;
        double r10630 = r10629 * r10628;
        double r10631 = r10623 * r10630;
        double r10632 = r10614 + r10631;
        double r10633 = r10632 * r10612;
        double r10634 = r10613 + r10633;
        double r10635 = r10612 * r10634;
        double r10636 = -r10608;
        double r10637 = r10636 * r10608;
        double r10638 = exp(r10637);
        double r10639 = r10635 * r10638;
        double r10640 = r10606 - r10639;
        return r10640;
}

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 r10641, r10642, r10643, r10644, r10645, r10646, r10647, r10648, r10649, r10650, r10651, r10652, r10653, r10654, r10655, r10656, r10657, r10658, r10659, r10660, r10661, r10662, r10663, r10664, r10665, r10666;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10641, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r10642, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r10643);
        mpfr_init(r10644);
        mpfr_init(r10645);
        mpfr_init(r10646);
        mpfr_init(r10647);
        mpfr_init_set_str(r10648, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r10649, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r10650, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r10651, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r10652, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r10653);
        mpfr_init(r10654);
        mpfr_init(r10655);
        mpfr_init(r10656);
        mpfr_init(r10657);
        mpfr_init(r10658);
        mpfr_init(r10659);
        mpfr_init(r10660);
        mpfr_init(r10661);
        mpfr_init(r10662);
        mpfr_init(r10663);
        mpfr_init(r10664);
        mpfr_init(r10665);
        mpfr_init(r10666);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r10643, x, MPFR_RNDN);
        mpfr_abs(r10644, r10643, MPFR_RNDN);
        mpfr_mul(r10645, r10642, r10644, MPFR_RNDN);
        mpfr_add(r10646, r10641, r10645, MPFR_RNDN);
        mpfr_div(r10647, r10641, r10646, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r10653, r10647, r10652, MPFR_RNDN);
        mpfr_add(r10654, r10651, r10653, MPFR_RNDN);
        mpfr_mul(r10655, r10647, r10654, MPFR_RNDN);
        mpfr_add(r10656, r10650, r10655, MPFR_RNDN);
        mpfr_mul(r10657, r10647, r10656, MPFR_RNDN);
        mpfr_add(r10658, r10649, r10657, MPFR_RNDN);
        mpfr_mul(r10659, r10647, r10658, MPFR_RNDN);
        mpfr_add(r10660, r10648, r10659, MPFR_RNDN);
        mpfr_mul(r10661, r10647, r10660, MPFR_RNDN);
        mpfr_mul(r10662, r10644, r10644, MPFR_RNDN);
        mpfr_neg(r10663, r10662, MPFR_RNDN);
        mpfr_exp(r10664, r10663, MPFR_RNDN);
        mpfr_mul(r10665, r10661, r10664, MPFR_RNDN);
        mpfr_sub(r10666, r10641, r10665, MPFR_RNDN);
        return mpfr_get_d(r10666, MPFR_RNDN);
}

static mpfr_t r10667, r10668, r10669, r10670, r10671, r10672, r10673, r10674, r10675, r10676, r10677, r10678, r10679, r10680, r10681, r10682, r10683, r10684, r10685, r10686, r10687, r10688, r10689, r10690, r10691, r10692, r10693, r10694, r10695, r10696, r10697, r10698, r10699, r10700, r10701;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10667, "1", 10, MPFR_RNDN);
        mpfr_init(r10668);
        mpfr_init(r10669);
        mpfr_init_set_str(r10670, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r10671);
        mpfr_init(r10672);
        mpfr_init(r10673);
        mpfr_init_set_str(r10674, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r10675, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r10676);
        mpfr_init_set_str(r10677, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r10678);
        mpfr_init_set_str(r10679, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r10680);
        mpfr_init(r10681);
        mpfr_init_set_str(r10682, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r10683);
        mpfr_init(r10684);
        mpfr_init(r10685);
        mpfr_init(r10686);
        mpfr_init(r10687);
        mpfr_init(r10688);
        mpfr_init(r10689);
        mpfr_init(r10690);
        mpfr_init(r10691);
        mpfr_init(r10692);
        mpfr_init(r10693);
        mpfr_init(r10694);
        mpfr_init(r10695);
        mpfr_init(r10696);
        mpfr_init(r10697);
        mpfr_init(r10698);
        mpfr_init(r10699);
        mpfr_init(r10700);
        mpfr_init(r10701);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10668, x, MPFR_RNDN);
        mpfr_abs(r10669, r10668, MPFR_RNDN);
        ;
        mpfr_mul(r10671, r10669, r10670, MPFR_RNDN);
        mpfr_add(r10672, r10667, r10671, MPFR_RNDN);
        mpfr_div(r10673, r10667, r10672, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10676, r10667, r10671, MPFR_RNDN);
        ;
        mpfr_mul(r10678, r10673, r10677, MPFR_RNDN);
        ;
        mpfr_add(r10680, r10678, r10679, MPFR_RNDN);
        mpfr_mul(r10681, r10673, r10680, MPFR_RNDN);
        ;
        mpfr_add(r10683, r10681, r10682, MPFR_RNDN);
        mpfr_mul(r10684, r10676, r10683, MPFR_RNDN);
        mpfr_neg(r10685, r10670, MPFR_RNDN);
        mpfr_mul(r10686, r10685, r10669, MPFR_RNDN);
        mpfr_fma(r10687, r10671, r10686, r10667, MPFR_RNDN);
        mpfr_div(r10688, r10667, r10687, MPFR_RNDN);
        mpfr_cbrt(r10689, r10688, MPFR_RNDN);
        mpfr_mul(r10690, r10689, r10689, MPFR_RNDN);
        mpfr_mul(r10691, r10690, r10689, MPFR_RNDN);
        mpfr_mul(r10692, r10684, r10691, MPFR_RNDN);
        mpfr_add(r10693, r10675, r10692, MPFR_RNDN);
        mpfr_mul(r10694, r10693, r10673, MPFR_RNDN);
        mpfr_add(r10695, r10674, r10694, MPFR_RNDN);
        mpfr_mul(r10696, r10673, r10695, MPFR_RNDN);
        mpfr_neg(r10697, r10669, MPFR_RNDN);
        mpfr_mul(r10698, r10697, r10669, MPFR_RNDN);
        mpfr_exp(r10699, r10698, MPFR_RNDN);
        mpfr_mul(r10700, r10696, r10699, MPFR_RNDN);
        mpfr_sub(r10701, r10667, r10700, MPFR_RNDN);
        return mpfr_get_d(r10701, MPFR_RNDN);
}

static mpfr_t r10702, r10703, r10704, r10705, r10706, r10707, r10708, r10709, r10710, r10711, r10712, r10713, r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721, r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729, r10730, r10731, r10732, r10733, r10734, r10735, r10736;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10702, "1", 10, MPFR_RNDN);
        mpfr_init(r10703);
        mpfr_init(r10704);
        mpfr_init_set_str(r10705, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r10706);
        mpfr_init(r10707);
        mpfr_init(r10708);
        mpfr_init_set_str(r10709, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r10710, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r10711);
        mpfr_init_set_str(r10712, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r10713);
        mpfr_init_set_str(r10714, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r10715);
        mpfr_init(r10716);
        mpfr_init_set_str(r10717, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r10718);
        mpfr_init(r10719);
        mpfr_init(r10720);
        mpfr_init(r10721);
        mpfr_init(r10722);
        mpfr_init(r10723);
        mpfr_init(r10724);
        mpfr_init(r10725);
        mpfr_init(r10726);
        mpfr_init(r10727);
        mpfr_init(r10728);
        mpfr_init(r10729);
        mpfr_init(r10730);
        mpfr_init(r10731);
        mpfr_init(r10732);
        mpfr_init(r10733);
        mpfr_init(r10734);
        mpfr_init(r10735);
        mpfr_init(r10736);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10703, x, MPFR_RNDN);
        mpfr_abs(r10704, r10703, MPFR_RNDN);
        ;
        mpfr_mul(r10706, r10704, r10705, MPFR_RNDN);
        mpfr_add(r10707, r10702, r10706, MPFR_RNDN);
        mpfr_div(r10708, r10702, r10707, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10711, r10702, r10706, MPFR_RNDN);
        ;
        mpfr_mul(r10713, r10708, r10712, MPFR_RNDN);
        ;
        mpfr_add(r10715, r10713, r10714, MPFR_RNDN);
        mpfr_mul(r10716, r10708, r10715, MPFR_RNDN);
        ;
        mpfr_add(r10718, r10716, r10717, MPFR_RNDN);
        mpfr_mul(r10719, r10711, r10718, MPFR_RNDN);
        mpfr_neg(r10720, r10705, MPFR_RNDN);
        mpfr_mul(r10721, r10720, r10704, MPFR_RNDN);
        mpfr_fma(r10722, r10706, r10721, r10702, MPFR_RNDN);
        mpfr_div(r10723, r10702, r10722, MPFR_RNDN);
        mpfr_cbrt(r10724, r10723, MPFR_RNDN);
        mpfr_mul(r10725, r10724, r10724, MPFR_RNDN);
        mpfr_mul(r10726, r10725, r10724, MPFR_RNDN);
        mpfr_mul(r10727, r10719, r10726, MPFR_RNDN);
        mpfr_add(r10728, r10710, r10727, MPFR_RNDN);
        mpfr_mul(r10729, r10728, r10708, MPFR_RNDN);
        mpfr_add(r10730, r10709, r10729, MPFR_RNDN);
        mpfr_mul(r10731, r10708, r10730, MPFR_RNDN);
        mpfr_neg(r10732, r10704, MPFR_RNDN);
        mpfr_mul(r10733, r10732, r10704, MPFR_RNDN);
        mpfr_exp(r10734, r10733, MPFR_RNDN);
        mpfr_mul(r10735, r10731, r10734, MPFR_RNDN);
        mpfr_sub(r10736, r10702, r10735, MPFR_RNDN);
        return mpfr_get_d(r10736, MPFR_RNDN);
}