#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 r27453 = 1;
        float r27454 = 0.3275911;
        float r27455 = x;
        float r27456 = fabs(r27455);
        float r27457 = r27454 * r27456;
        float r27458 = r27453 + r27457;
        float r27459 = r27453 / r27458;
        float r27460 = 0.254829592;
        float r27461 = -0.284496736;
        float r27462 = 1.421413741;
        float r27463 = -1.453152027;
        float r27464 = 1.061405429;
        float r27465 = r27459 * r27464;
        float r27466 = r27463 + r27465;
        float r27467 = r27459 * r27466;
        float r27468 = r27462 + r27467;
        float r27469 = r27459 * r27468;
        float r27470 = r27461 + r27469;
        float r27471 = r27459 * r27470;
        float r27472 = r27460 + r27471;
        float r27473 = r27459 * r27472;
        float r27474 = r27456 * r27456;
        float r27475 = -r27474;
        float r27476 = exp(r27475);
        float r27477 = r27473 * r27476;
        float r27478 = r27453 - r27477;
        return r27478;
}

double f_id(double x) {
        double r27479 = 1;
        double r27480 = 0.3275911;
        double r27481 = x;
        double r27482 = fabs(r27481);
        double r27483 = r27480 * r27482;
        double r27484 = r27479 + r27483;
        double r27485 = r27479 / r27484;
        double r27486 = 0.254829592;
        double r27487 = -0.284496736;
        double r27488 = 1.421413741;
        double r27489 = -1.453152027;
        double r27490 = 1.061405429;
        double r27491 = r27485 * r27490;
        double r27492 = r27489 + r27491;
        double r27493 = r27485 * r27492;
        double r27494 = r27488 + r27493;
        double r27495 = r27485 * r27494;
        double r27496 = r27487 + r27495;
        double r27497 = r27485 * r27496;
        double r27498 = r27486 + r27497;
        double r27499 = r27485 * r27498;
        double r27500 = r27482 * r27482;
        double r27501 = -r27500;
        double r27502 = exp(r27501);
        double r27503 = r27499 * r27502;
        double r27504 = r27479 - r27503;
        return r27504;
}


double f_of(float x) {
        float r27505 = 1;
        float r27506 = -r27505;
        float r27507 = x;
        float r27508 = fabs(r27507);
        float r27509 = r27508 * r27508;
        float r27510 = exp(r27509);
        float r27511 = r27506 / r27510;
        float r27512 = 1.061405429;
        float r27513 = 0.3275911;
        float r27514 = fma(r27508, r27513, r27505);
        float r27515 = r27512 / r27514;
        float r27516 = -1.453152027;
        float r27517 = r27515 + r27516;
        float r27518 = r27505 / r27514;
        float r27519 = 1.421413741;
        float r27520 = fma(r27517, r27518, r27519);
        float r27521 = -0.284496736;
        float r27522 = fma(r27520, r27518, r27521);
        float r27523 = fma(r27513, r27508, r27505);
        float r27524 = r27523 * r27523;
        float r27525 = r27505 / r27524;
        float r27526 = cbrt(r27525);
        float r27527 = r27526 * r27526;
        float r27528 = r27505 / r27523;
        float r27529 = r27528 / r27523;
        float r27530 = cbrt(r27529);
        float r27531 = r27527 * r27530;
        float r27532 = 0.254829592;
        float r27533 = r27532 / r27514;
        float r27534 = fma(r27522, r27531, r27533);
        float r27535 = fma(r27511, r27534, r27505);
        return r27535;
}

double f_od(double x) {
        double r27536 = 1;
        double r27537 = -r27536;
        double r27538 = x;
        double r27539 = fabs(r27538);
        double r27540 = r27539 * r27539;
        double r27541 = exp(r27540);
        double r27542 = r27537 / r27541;
        double r27543 = 1.061405429;
        double r27544 = 0.3275911;
        double r27545 = fma(r27539, r27544, r27536);
        double r27546 = r27543 / r27545;
        double r27547 = -1.453152027;
        double r27548 = r27546 + r27547;
        double r27549 = r27536 / r27545;
        double r27550 = 1.421413741;
        double r27551 = fma(r27548, r27549, r27550);
        double r27552 = -0.284496736;
        double r27553 = fma(r27551, r27549, r27552);
        double r27554 = fma(r27544, r27539, r27536);
        double r27555 = r27554 * r27554;
        double r27556 = r27536 / r27555;
        double r27557 = cbrt(r27556);
        double r27558 = r27557 * r27557;
        double r27559 = r27536 / r27554;
        double r27560 = r27559 / r27554;
        double r27561 = cbrt(r27560);
        double r27562 = r27558 * r27561;
        double r27563 = 0.254829592;
        double r27564 = r27563 / r27545;
        double r27565 = fma(r27553, r27562, r27564);
        double r27566 = fma(r27542, r27565, r27536);
        return r27566;
}

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 r27567, r27568, r27569, r27570, r27571, r27572, r27573, r27574, r27575, r27576, r27577, r27578, r27579, r27580, r27581, r27582, r27583, r27584, r27585, r27586, r27587, r27588, r27589, r27590, r27591, r27592;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27567, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27568, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r27569);
        mpfr_init(r27570);
        mpfr_init(r27571);
        mpfr_init(r27572);
        mpfr_init(r27573);
        mpfr_init_set_str(r27574, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r27575, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r27576, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r27577, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r27578, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r27579);
        mpfr_init(r27580);
        mpfr_init(r27581);
        mpfr_init(r27582);
        mpfr_init(r27583);
        mpfr_init(r27584);
        mpfr_init(r27585);
        mpfr_init(r27586);
        mpfr_init(r27587);
        mpfr_init(r27588);
        mpfr_init(r27589);
        mpfr_init(r27590);
        mpfr_init(r27591);
        mpfr_init(r27592);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r27569, x, MPFR_RNDN);
        mpfr_abs(r27570, r27569, MPFR_RNDN);
        mpfr_mul(r27571, r27568, r27570, MPFR_RNDN);
        mpfr_add(r27572, r27567, r27571, MPFR_RNDN);
        mpfr_div(r27573, r27567, r27572, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r27579, r27573, r27578, MPFR_RNDN);
        mpfr_add(r27580, r27577, r27579, MPFR_RNDN);
        mpfr_mul(r27581, r27573, r27580, MPFR_RNDN);
        mpfr_add(r27582, r27576, r27581, MPFR_RNDN);
        mpfr_mul(r27583, r27573, r27582, MPFR_RNDN);
        mpfr_add(r27584, r27575, r27583, MPFR_RNDN);
        mpfr_mul(r27585, r27573, r27584, MPFR_RNDN);
        mpfr_add(r27586, r27574, r27585, MPFR_RNDN);
        mpfr_mul(r27587, r27573, r27586, MPFR_RNDN);
        mpfr_mul(r27588, r27570, r27570, MPFR_RNDN);
        mpfr_neg(r27589, r27588, MPFR_RNDN);
        mpfr_exp(r27590, r27589, MPFR_RNDN);
        mpfr_mul(r27591, r27587, r27590, MPFR_RNDN);
        mpfr_sub(r27592, r27567, r27591, MPFR_RNDN);
        return mpfr_get_d(r27592, MPFR_RNDN);
}

static mpfr_t r27593, r27594, r27595, r27596, r27597, r27598, r27599, r27600, r27601, r27602, r27603, r27604, r27605, r27606, r27607, r27608, r27609, r27610, r27611, r27612, r27613, r27614, r27615, r27616, r27617, r27618, r27619, r27620, r27621, r27622, r27623;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27593, "1", 10, MPFR_RNDN);
        mpfr_init(r27594);
        mpfr_init(r27595);
        mpfr_init(r27596);
        mpfr_init(r27597);
        mpfr_init(r27598);
        mpfr_init(r27599);
        mpfr_init_set_str(r27600, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r27601, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r27602);
        mpfr_init(r27603);
        mpfr_init_set_str(r27604, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r27605);
        mpfr_init(r27606);
        mpfr_init_set_str(r27607, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r27608);
        mpfr_init_set_str(r27609, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r27610);
        mpfr_init(r27611);
        mpfr_init(r27612);
        mpfr_init(r27613);
        mpfr_init(r27614);
        mpfr_init(r27615);
        mpfr_init(r27616);
        mpfr_init(r27617);
        mpfr_init(r27618);
        mpfr_init(r27619);
        mpfr_init_set_str(r27620, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r27621);
        mpfr_init(r27622);
        mpfr_init(r27623);
}

double f_fm(double x) {
        ;
        mpfr_neg(r27594, r27593, MPFR_RNDN);
        mpfr_set_d(r27595, x, MPFR_RNDN);
        mpfr_abs(r27596, r27595, MPFR_RNDN);
        mpfr_mul(r27597, r27596, r27596, MPFR_RNDN);
        mpfr_exp(r27598, r27597, MPFR_RNDN);
        mpfr_div(r27599, r27594, r27598, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r27602, r27596, r27601, r27593, MPFR_RNDN);
        mpfr_div(r27603, r27600, r27602, MPFR_RNDN);
        ;
        mpfr_add(r27605, r27603, r27604, MPFR_RNDN);
        mpfr_div(r27606, r27593, r27602, MPFR_RNDN);
        ;
        mpfr_fma(r27608, r27605, r27606, r27607, MPFR_RNDN);
        ;
        mpfr_fma(r27610, r27608, r27606, r27609, MPFR_RNDN);
        mpfr_fma(r27611, r27601, r27596, r27593, MPFR_RNDN);
        mpfr_mul(r27612, r27611, r27611, MPFR_RNDN);
        mpfr_div(r27613, r27593, r27612, MPFR_RNDN);
        mpfr_cbrt(r27614, r27613, MPFR_RNDN);
        mpfr_mul(r27615, r27614, r27614, MPFR_RNDN);
        mpfr_div(r27616, r27593, r27611, MPFR_RNDN);
        mpfr_div(r27617, r27616, r27611, MPFR_RNDN);
        mpfr_cbrt(r27618, r27617, MPFR_RNDN);
        mpfr_mul(r27619, r27615, r27618, MPFR_RNDN);
        ;
        mpfr_div(r27621, r27620, r27602, MPFR_RNDN);
        mpfr_fma(r27622, r27610, r27619, r27621, MPFR_RNDN);
        mpfr_fma(r27623, r27599, r27622, r27593, MPFR_RNDN);
        return mpfr_get_d(r27623, MPFR_RNDN);
}

static mpfr_t r27624, r27625, r27626, r27627, r27628, r27629, r27630, r27631, r27632, r27633, r27634, r27635, r27636, r27637, r27638, r27639, r27640, r27641, r27642, r27643, r27644, r27645, r27646, r27647, r27648, r27649, r27650, r27651, r27652, r27653, r27654;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27624, "1", 10, MPFR_RNDN);
        mpfr_init(r27625);
        mpfr_init(r27626);
        mpfr_init(r27627);
        mpfr_init(r27628);
        mpfr_init(r27629);
        mpfr_init(r27630);
        mpfr_init_set_str(r27631, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r27632, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r27633);
        mpfr_init(r27634);
        mpfr_init_set_str(r27635, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r27636);
        mpfr_init(r27637);
        mpfr_init_set_str(r27638, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r27639);
        mpfr_init_set_str(r27640, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r27641);
        mpfr_init(r27642);
        mpfr_init(r27643);
        mpfr_init(r27644);
        mpfr_init(r27645);
        mpfr_init(r27646);
        mpfr_init(r27647);
        mpfr_init(r27648);
        mpfr_init(r27649);
        mpfr_init(r27650);
        mpfr_init_set_str(r27651, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r27652);
        mpfr_init(r27653);
        mpfr_init(r27654);
}

double f_dm(double x) {
        ;
        mpfr_neg(r27625, r27624, MPFR_RNDN);
        mpfr_set_d(r27626, x, MPFR_RNDN);
        mpfr_abs(r27627, r27626, MPFR_RNDN);
        mpfr_mul(r27628, r27627, r27627, MPFR_RNDN);
        mpfr_exp(r27629, r27628, MPFR_RNDN);
        mpfr_div(r27630, r27625, r27629, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r27633, r27627, r27632, r27624, MPFR_RNDN);
        mpfr_div(r27634, r27631, r27633, MPFR_RNDN);
        ;
        mpfr_add(r27636, r27634, r27635, MPFR_RNDN);
        mpfr_div(r27637, r27624, r27633, MPFR_RNDN);
        ;
        mpfr_fma(r27639, r27636, r27637, r27638, MPFR_RNDN);
        ;
        mpfr_fma(r27641, r27639, r27637, r27640, MPFR_RNDN);
        mpfr_fma(r27642, r27632, r27627, r27624, MPFR_RNDN);
        mpfr_mul(r27643, r27642, r27642, MPFR_RNDN);
        mpfr_div(r27644, r27624, r27643, MPFR_RNDN);
        mpfr_cbrt(r27645, r27644, MPFR_RNDN);
        mpfr_mul(r27646, r27645, r27645, MPFR_RNDN);
        mpfr_div(r27647, r27624, r27642, MPFR_RNDN);
        mpfr_div(r27648, r27647, r27642, MPFR_RNDN);
        mpfr_cbrt(r27649, r27648, MPFR_RNDN);
        mpfr_mul(r27650, r27646, r27649, MPFR_RNDN);
        ;
        mpfr_div(r27652, r27651, r27633, MPFR_RNDN);
        mpfr_fma(r27653, r27641, r27650, r27652, MPFR_RNDN);
        mpfr_fma(r27654, r27630, r27653, r27624, MPFR_RNDN);
        return mpfr_get_d(r27654, MPFR_RNDN);
}