#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 r18465 = 1.0f;
        float r18466 = 0.3275911f;
        float r18467 = x;
        float r18468 = fabs(r18467);
        float r18469 = r18466 * r18468;
        float r18470 = r18465 + r18469;
        float r18471 = r18465 / r18470;
        float r18472 = 0.254829592f;
        float r18473 = -0.284496736f;
        float r18474 = 1.421413741f;
        float r18475 = -1.453152027f;
        float r18476 = 1.061405429f;
        float r18477 = r18471 * r18476;
        float r18478 = r18475 + r18477;
        float r18479 = r18471 * r18478;
        float r18480 = r18474 + r18479;
        float r18481 = r18471 * r18480;
        float r18482 = r18473 + r18481;
        float r18483 = r18471 * r18482;
        float r18484 = r18472 + r18483;
        float r18485 = r18471 * r18484;
        float r18486 = r18468 * r18468;
        float r18487 = -r18486;
        float r18488 = exp(r18487);
        float r18489 = r18485 * r18488;
        float r18490 = r18465 - r18489;
        return r18490;
}

double f_id(double x) {
        double r18491 = 1.0;
        double r18492 = 0.3275911;
        double r18493 = x;
        double r18494 = fabs(r18493);
        double r18495 = r18492 * r18494;
        double r18496 = r18491 + r18495;
        double r18497 = r18491 / r18496;
        double r18498 = 0.254829592;
        double r18499 = -0.284496736;
        double r18500 = 1.421413741;
        double r18501 = -1.453152027;
        double r18502 = 1.061405429;
        double r18503 = r18497 * r18502;
        double r18504 = r18501 + r18503;
        double r18505 = r18497 * r18504;
        double r18506 = r18500 + r18505;
        double r18507 = r18497 * r18506;
        double r18508 = r18499 + r18507;
        double r18509 = r18497 * r18508;
        double r18510 = r18498 + r18509;
        double r18511 = r18497 * r18510;
        double r18512 = r18494 * r18494;
        double r18513 = -r18512;
        double r18514 = exp(r18513);
        double r18515 = r18511 * r18514;
        double r18516 = r18491 - r18515;
        return r18516;
}


double f_of(float x) {
        float r18517 = 1.0f;
        float r18518 = 1.061405429f;
        float r18519 = 0.3275911f;
        float r18520 = x;
        float r18521 = fabs(r18520);
        float r18522 = fma(r18519, r18521, r18517);
        float r18523 = r18518 / r18522;
        float r18524 = -1.453152027f;
        float r18525 = r18523 + r18524;
        float r18526 = r18517 / r18522;
        float r18527 = r18526 / r18522;
        float r18528 = 1.421413741f;
        float r18529 = r18528 / r18522;
        float r18530 = -0.284496736f;
        float r18531 = r18529 + r18530;
        float r18532 = fma(r18525, r18527, r18531);
        float r18533 = 0.254829592f;
        float r18534 = r18533 / r18522;
        float r18535 = fma(r18532, r18527, r18534);
        float r18536 = r18521 * r18521;
        float r18537 = exp(r18536);
        float r18538 = r18535 / r18537;
        float r18539 = r18517 - r18538;
        float r18540 = exp(r18539);
        float r18541 = log(r18540);
        return r18541;
}

double f_od(double x) {
        double r18542 = 1.0;
        double r18543 = 1.061405429;
        double r18544 = 0.3275911;
        double r18545 = x;
        double r18546 = fabs(r18545);
        double r18547 = fma(r18544, r18546, r18542);
        double r18548 = r18543 / r18547;
        double r18549 = -1.453152027;
        double r18550 = r18548 + r18549;
        double r18551 = r18542 / r18547;
        double r18552 = r18551 / r18547;
        double r18553 = 1.421413741;
        double r18554 = r18553 / r18547;
        double r18555 = -0.284496736;
        double r18556 = r18554 + r18555;
        double r18557 = fma(r18550, r18552, r18556);
        double r18558 = 0.254829592;
        double r18559 = r18558 / r18547;
        double r18560 = fma(r18557, r18552, r18559);
        double r18561 = r18546 * r18546;
        double r18562 = exp(r18561);
        double r18563 = r18560 / r18562;
        double r18564 = r18542 - r18563;
        double r18565 = exp(r18564);
        double r18566 = log(r18565);
        return r18566;
}

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 r18567, r18568, r18569, r18570, r18571, r18572, r18573, r18574, r18575, r18576, r18577, r18578, r18579, r18580, r18581, r18582, r18583, r18584, r18585, r18586, r18587, r18588, r18589, r18590, r18591, r18592;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18567, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18568, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18569);
        mpfr_init(r18570);
        mpfr_init(r18571);
        mpfr_init(r18572);
        mpfr_init(r18573);
        mpfr_init_set_str(r18574, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18575, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18576, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18577, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18578, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18579);
        mpfr_init(r18580);
        mpfr_init(r18581);
        mpfr_init(r18582);
        mpfr_init(r18583);
        mpfr_init(r18584);
        mpfr_init(r18585);
        mpfr_init(r18586);
        mpfr_init(r18587);
        mpfr_init(r18588);
        mpfr_init(r18589);
        mpfr_init(r18590);
        mpfr_init(r18591);
        mpfr_init(r18592);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18569, x, MPFR_RNDN);
        mpfr_abs(r18570, r18569, MPFR_RNDN);
        mpfr_mul(r18571, r18568, r18570, MPFR_RNDN);
        mpfr_add(r18572, r18567, r18571, MPFR_RNDN);
        mpfr_div(r18573, r18567, r18572, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18579, r18573, r18578, MPFR_RNDN);
        mpfr_add(r18580, r18577, r18579, MPFR_RNDN);
        mpfr_mul(r18581, r18573, r18580, MPFR_RNDN);
        mpfr_add(r18582, r18576, r18581, MPFR_RNDN);
        mpfr_mul(r18583, r18573, r18582, MPFR_RNDN);
        mpfr_add(r18584, r18575, r18583, MPFR_RNDN);
        mpfr_mul(r18585, r18573, r18584, MPFR_RNDN);
        mpfr_add(r18586, r18574, r18585, MPFR_RNDN);
        mpfr_mul(r18587, r18573, r18586, MPFR_RNDN);
        mpfr_mul(r18588, r18570, r18570, MPFR_RNDN);
        mpfr_neg(r18589, r18588, MPFR_RNDN);
        mpfr_exp(r18590, r18589, MPFR_RNDN);
        mpfr_mul(r18591, r18587, r18590, MPFR_RNDN);
        mpfr_sub(r18592, r18567, r18591, MPFR_RNDN);
        return mpfr_get_d(r18592, MPFR_RNDN);
}

static mpfr_t r18593, r18594, r18595, r18596, r18597, r18598, r18599, r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609, r18610, r18611, r18612, r18613, r18614, r18615, r18616, r18617;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18593, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18594, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18595, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18596);
        mpfr_init(r18597);
        mpfr_init(r18598);
        mpfr_init(r18599);
        mpfr_init_set_str(r18600, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18601);
        mpfr_init(r18602);
        mpfr_init(r18603);
        mpfr_init_set_str(r18604, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18605);
        mpfr_init_set_str(r18606, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18607);
        mpfr_init(r18608);
        mpfr_init_set_str(r18609, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18610);
        mpfr_init(r18611);
        mpfr_init(r18612);
        mpfr_init(r18613);
        mpfr_init(r18614);
        mpfr_init(r18615);
        mpfr_init(r18616);
        mpfr_init(r18617);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18596, x, MPFR_RNDN);
        mpfr_abs(r18597, r18596, MPFR_RNDN);
        mpfr_fma(r18598, r18595, r18597, r18593, MPFR_RNDN);
        mpfr_div(r18599, r18594, r18598, MPFR_RNDN);
        ;
        mpfr_add(r18601, r18599, r18600, MPFR_RNDN);
        mpfr_div(r18602, r18593, r18598, MPFR_RNDN);
        mpfr_div(r18603, r18602, r18598, MPFR_RNDN);
        ;
        mpfr_div(r18605, r18604, r18598, MPFR_RNDN);
        ;
        mpfr_add(r18607, r18605, r18606, MPFR_RNDN);
        mpfr_fma(r18608, r18601, r18603, r18607, MPFR_RNDN);
        ;
        mpfr_div(r18610, r18609, r18598, MPFR_RNDN);
        mpfr_fma(r18611, r18608, r18603, r18610, MPFR_RNDN);
        mpfr_mul(r18612, r18597, r18597, MPFR_RNDN);
        mpfr_exp(r18613, r18612, MPFR_RNDN);
        mpfr_div(r18614, r18611, r18613, MPFR_RNDN);
        mpfr_sub(r18615, r18593, r18614, MPFR_RNDN);
        mpfr_exp(r18616, r18615, MPFR_RNDN);
        mpfr_log(r18617, r18616, MPFR_RNDN);
        return mpfr_get_d(r18617, MPFR_RNDN);
}

static mpfr_t r18618, r18619, r18620, r18621, r18622, r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634, r18635, r18636, r18637, r18638, r18639, r18640, r18641, r18642;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18618, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18619, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18620, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18621);
        mpfr_init(r18622);
        mpfr_init(r18623);
        mpfr_init(r18624);
        mpfr_init_set_str(r18625, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18626);
        mpfr_init(r18627);
        mpfr_init(r18628);
        mpfr_init_set_str(r18629, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18630);
        mpfr_init_set_str(r18631, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18632);
        mpfr_init(r18633);
        mpfr_init_set_str(r18634, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18635);
        mpfr_init(r18636);
        mpfr_init(r18637);
        mpfr_init(r18638);
        mpfr_init(r18639);
        mpfr_init(r18640);
        mpfr_init(r18641);
        mpfr_init(r18642);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18621, x, MPFR_RNDN);
        mpfr_abs(r18622, r18621, MPFR_RNDN);
        mpfr_fma(r18623, r18620, r18622, r18618, MPFR_RNDN);
        mpfr_div(r18624, r18619, r18623, MPFR_RNDN);
        ;
        mpfr_add(r18626, r18624, r18625, MPFR_RNDN);
        mpfr_div(r18627, r18618, r18623, MPFR_RNDN);
        mpfr_div(r18628, r18627, r18623, MPFR_RNDN);
        ;
        mpfr_div(r18630, r18629, r18623, MPFR_RNDN);
        ;
        mpfr_add(r18632, r18630, r18631, MPFR_RNDN);
        mpfr_fma(r18633, r18626, r18628, r18632, MPFR_RNDN);
        ;
        mpfr_div(r18635, r18634, r18623, MPFR_RNDN);
        mpfr_fma(r18636, r18633, r18628, r18635, MPFR_RNDN);
        mpfr_mul(r18637, r18622, r18622, MPFR_RNDN);
        mpfr_exp(r18638, r18637, MPFR_RNDN);
        mpfr_div(r18639, r18636, r18638, MPFR_RNDN);
        mpfr_sub(r18640, r18618, r18639, MPFR_RNDN);
        mpfr_exp(r18641, r18640, MPFR_RNDN);
        mpfr_log(r18642, r18641, MPFR_RNDN);
        return mpfr_get_d(r18642, MPFR_RNDN);
}