#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 r18445 = 1.0f;
        float r18446 = 0.32759109139442444f;
        float r18447 = x;
        float r18448 = fabs(r18447);
        float r18449 = r18446 * r18448;
        float r18450 = r18445 + r18449;
        float r18451 = r18445 / r18450;
        float r18452 = 0.2548295855522156f;
        float r18453 = -0.2844967246055603f;
        float r18454 = 1.421413779258728f;
        float r18455 = -1.453152060508728f;
        float r18456 = 1.0614054203033447f;
        float r18457 = r18451 * r18456;
        float r18458 = r18455 + r18457;
        float r18459 = r18451 * r18458;
        float r18460 = r18454 + r18459;
        float r18461 = r18451 * r18460;
        float r18462 = r18453 + r18461;
        float r18463 = r18451 * r18462;
        float r18464 = r18452 + r18463;
        float r18465 = r18451 * r18464;
        float r18466 = r18448 * r18448;
        float r18467 = -r18466;
        float r18468 = exp(r18467);
        float r18469 = r18465 * r18468;
        float r18470 = r18445 - r18469;
        return r18470;
}

double f_id(double x) {
        double r18471 = 1.0;
        double r18472 = 0.32759109139442444;
        double r18473 = x;
        double r18474 = fabs(r18473);
        double r18475 = r18472 * r18474;
        double r18476 = r18471 + r18475;
        double r18477 = r18471 / r18476;
        double r18478 = 0.2548295855522156;
        double r18479 = -0.2844967246055603;
        double r18480 = 1.421413779258728;
        double r18481 = -1.453152060508728;
        double r18482 = 1.0614054203033447;
        double r18483 = r18477 * r18482;
        double r18484 = r18481 + r18483;
        double r18485 = r18477 * r18484;
        double r18486 = r18480 + r18485;
        double r18487 = r18477 * r18486;
        double r18488 = r18479 + r18487;
        double r18489 = r18477 * r18488;
        double r18490 = r18478 + r18489;
        double r18491 = r18477 * r18490;
        double r18492 = r18474 * r18474;
        double r18493 = -r18492;
        double r18494 = exp(r18493);
        double r18495 = r18491 * r18494;
        double r18496 = r18471 - r18495;
        return r18496;
}


double f_of(float x) {
        float r18497 = 1.0f;
        float r18498 = 1.0614054203033447f;
        float r18499 = 0.32759109139442444f;
        float r18500 = x;
        float r18501 = fabs(r18500);
        float r18502 = fma(r18499, r18501, r18497);
        float r18503 = r18498 / r18502;
        float r18504 = -1.453152060508728f;
        float r18505 = r18503 + r18504;
        float r18506 = r18497 / r18502;
        float r18507 = r18506 / r18502;
        float r18508 = 1.421413779258728f;
        float r18509 = r18508 / r18502;
        float r18510 = -0.2844967246055603f;
        float r18511 = r18509 + r18510;
        float r18512 = fma(r18505, r18507, r18511);
        float r18513 = 0.2548295855522156f;
        float r18514 = r18513 / r18502;
        float r18515 = fma(r18512, r18507, r18514);
        float r18516 = r18501 * r18501;
        float r18517 = exp(r18516);
        float r18518 = r18515 / r18517;
        float r18519 = r18497 - r18518;
        float r18520 = cbrt(r18519);
        float r18521 = r18520 * (r18520 * r18520);
        return r18521;
}

double f_od(double x) {
        double r18522 = 1.0;
        double r18523 = 1.0614054203033447;
        double r18524 = 0.32759109139442444;
        double r18525 = x;
        double r18526 = fabs(r18525);
        double r18527 = fma(r18524, r18526, r18522);
        double r18528 = r18523 / r18527;
        double r18529 = -1.453152060508728;
        double r18530 = r18528 + r18529;
        double r18531 = r18522 / r18527;
        double r18532 = r18531 / r18527;
        double r18533 = 1.421413779258728;
        double r18534 = r18533 / r18527;
        double r18535 = -0.2844967246055603;
        double r18536 = r18534 + r18535;
        double r18537 = fma(r18530, r18532, r18536);
        double r18538 = 0.2548295855522156;
        double r18539 = r18538 / r18527;
        double r18540 = fma(r18537, r18532, r18539);
        double r18541 = r18526 * r18526;
        double r18542 = exp(r18541);
        double r18543 = r18540 / r18542;
        double r18544 = r18522 - r18543;
        double r18545 = cbrt(r18544);
        double r18546 = r18545 * (r18545 * r18545);
        return r18546;
}

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 r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567, r18568, r18569, r18570, r18571, r18572;

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

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18549, x, MPFR_RNDN);
        mpfr_abs(r18550, r18549, MPFR_RNDN);
        mpfr_mul(r18551, r18548, r18550, MPFR_RNDN);
        mpfr_add(r18552, r18547, r18551, MPFR_RNDN);
        mpfr_div(r18553, r18547, r18552, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18559, r18553, r18558, MPFR_RNDN);
        mpfr_add(r18560, r18557, r18559, MPFR_RNDN);
        mpfr_mul(r18561, r18553, r18560, MPFR_RNDN);
        mpfr_add(r18562, r18556, r18561, MPFR_RNDN);
        mpfr_mul(r18563, r18553, r18562, MPFR_RNDN);
        mpfr_add(r18564, r18555, r18563, MPFR_RNDN);
        mpfr_mul(r18565, r18553, r18564, MPFR_RNDN);
        mpfr_add(r18566, r18554, r18565, MPFR_RNDN);
        mpfr_mul(r18567, r18553, r18566, MPFR_RNDN);
        mpfr_mul(r18568, r18550, r18550, MPFR_RNDN);
        mpfr_neg(r18569, r18568, MPFR_RNDN);
        mpfr_exp(r18570, r18569, MPFR_RNDN);
        mpfr_mul(r18571, r18567, r18570, MPFR_RNDN);
        mpfr_sub(r18572, r18547, r18571, MPFR_RNDN);
        return mpfr_get_d(r18572, MPFR_RNDN);
}

static mpfr_t 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;

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

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

static mpfr_t r18598, r18599, r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609, r18610, r18611, r18612, r18613, r18614, r18615, r18616, r18617, r18618, r18619, r18620, r18621, r18622;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18598, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18599, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18600, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18601);
        mpfr_init(r18602);
        mpfr_init(r18603);
        mpfr_init(r18604);
        mpfr_init_set_str(r18605, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18606);
        mpfr_init(r18607);
        mpfr_init(r18608);
        mpfr_init_set_str(r18609, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18610);
        mpfr_init_set_str(r18611, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18612);
        mpfr_init(r18613);
        mpfr_init_set_str(r18614, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18615);
        mpfr_init(r18616);
        mpfr_init(r18617);
        mpfr_init(r18618);
        mpfr_init(r18619);
        mpfr_init(r18620);
        mpfr_init(r18621);
        mpfr_init(r18622);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18601, x, MPFR_RNDN);
        mpfr_abs(r18602, r18601, MPFR_RNDN);
        mpfr_fma(r18603, r18600, r18602, r18598, MPFR_RNDN);
        mpfr_div(r18604, r18599, r18603, MPFR_RNDN);
        ;
        mpfr_add(r18606, r18604, r18605, MPFR_RNDN);
        mpfr_div(r18607, r18598, r18603, MPFR_RNDN);
        mpfr_div(r18608, r18607, r18603, MPFR_RNDN);
        ;
        mpfr_div(r18610, r18609, r18603, MPFR_RNDN);
        ;
        mpfr_add(r18612, r18610, r18611, MPFR_RNDN);
        mpfr_fma(r18613, r18606, r18608, r18612, MPFR_RNDN);
        ;
        mpfr_div(r18615, r18614, r18603, MPFR_RNDN);
        mpfr_fma(r18616, r18613, r18608, r18615, MPFR_RNDN);
        mpfr_mul(r18617, r18602, r18602, MPFR_RNDN);
        mpfr_exp(r18618, r18617, MPFR_RNDN);
        mpfr_div(r18619, r18616, r18618, MPFR_RNDN);
        mpfr_sub(r18620, r18598, r18619, MPFR_RNDN);
        mpfr_cbrt(r18621, r18620, MPFR_RNDN);
        mpfr_mul(r18622, r18621, r18621, MPFR_RNDN); mpfr_mul(r18622, r18622, r18621, MPFR_RNDN);
        return mpfr_get_d(r18622, MPFR_RNDN);
}