#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 r18561 = 1.0f;
        float r18562 = 0.32759109139442444f;
        float r18563 = x;
        float r18564 = fabs(r18563);
        float r18565 = r18562 * r18564;
        float r18566 = r18561 + r18565;
        float r18567 = r18561 / r18566;
        float r18568 = 0.2548295855522156f;
        float r18569 = -0.2844967246055603f;
        float r18570 = 1.421413779258728f;
        float r18571 = -1.453152060508728f;
        float r18572 = 1.0614054203033447f;
        float r18573 = r18567 * r18572;
        float r18574 = r18571 + r18573;
        float r18575 = r18567 * r18574;
        float r18576 = r18570 + r18575;
        float r18577 = r18567 * r18576;
        float r18578 = r18569 + r18577;
        float r18579 = r18567 * r18578;
        float r18580 = r18568 + r18579;
        float r18581 = r18567 * r18580;
        float r18582 = r18564 * r18564;
        float r18583 = -r18582;
        float r18584 = exp(r18583);
        float r18585 = r18581 * r18584;
        float r18586 = r18561 - r18585;
        return r18586;
}

double f_id(double x) {
        double r18587 = 1.0;
        double r18588 = 0.32759109139442444;
        double r18589 = x;
        double r18590 = fabs(r18589);
        double r18591 = r18588 * r18590;
        double r18592 = r18587 + r18591;
        double r18593 = r18587 / r18592;
        double r18594 = 0.2548295855522156;
        double r18595 = -0.2844967246055603;
        double r18596 = 1.421413779258728;
        double r18597 = -1.453152060508728;
        double r18598 = 1.0614054203033447;
        double r18599 = r18593 * r18598;
        double r18600 = r18597 + r18599;
        double r18601 = r18593 * r18600;
        double r18602 = r18596 + r18601;
        double r18603 = r18593 * r18602;
        double r18604 = r18595 + r18603;
        double r18605 = r18593 * r18604;
        double r18606 = r18594 + r18605;
        double r18607 = r18593 * r18606;
        double r18608 = r18590 * r18590;
        double r18609 = -r18608;
        double r18610 = exp(r18609);
        double r18611 = r18607 * r18610;
        double r18612 = r18587 - r18611;
        return r18612;
}


double f_of(float x) {
        float r18613 = 1.0f;
        float r18614 = 1.0614054203033447f;
        float r18615 = 0.32759109139442444f;
        float r18616 = x;
        float r18617 = fabs(r18616);
        float r18618 = fma(r18615, r18617, r18613);
        float r18619 = r18614 / r18618;
        float r18620 = -1.453152060508728f;
        float r18621 = r18619 + r18620;
        float r18622 = r18613 / r18618;
        float r18623 = r18622 / r18618;
        float r18624 = 1.421413779258728f;
        float r18625 = r18624 / r18618;
        float r18626 = -0.2844967246055603f;
        float r18627 = fma(r18613, r18625, r18626);
        float r18628 = fma(r18621, r18623, r18627);
        float r18629 = 0.2548295855522156f;
        float r18630 = r18629 / r18618;
        float r18631 = fma(r18628, r18623, r18630);
        float r18632 = r18617 * r18617;
        float r18633 = exp(r18632);
        float r18634 = r18631 / r18633;
        float r18635 = r18613 - r18634;
        float r18636 = log(r18635);
        float r18637 = exp(r18636);
        return r18637;
}

double f_od(double x) {
        double r18638 = 1.0;
        double r18639 = 1.0614054203033447;
        double r18640 = 0.32759109139442444;
        double r18641 = x;
        double r18642 = fabs(r18641);
        double r18643 = fma(r18640, r18642, r18638);
        double r18644 = r18639 / r18643;
        double r18645 = -1.453152060508728;
        double r18646 = r18644 + r18645;
        double r18647 = r18638 / r18643;
        double r18648 = r18647 / r18643;
        double r18649 = 1.421413779258728;
        double r18650 = r18649 / r18643;
        double r18651 = -0.2844967246055603;
        double r18652 = fma(r18638, r18650, r18651);
        double r18653 = fma(r18646, r18648, r18652);
        double r18654 = 0.2548295855522156;
        double r18655 = r18654 / r18643;
        double r18656 = fma(r18653, r18648, r18655);
        double r18657 = r18642 * r18642;
        double r18658 = exp(r18657);
        double r18659 = r18656 / r18658;
        double r18660 = r18638 - r18659;
        double r18661 = log(r18660);
        double r18662 = exp(r18661);
        return r18662;
}

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 r18663, r18664, r18665, r18666, r18667, r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681, r18682, r18683, r18684, r18685, r18686, r18687, r18688;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18663, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18664, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18665);
        mpfr_init(r18666);
        mpfr_init(r18667);
        mpfr_init(r18668);
        mpfr_init(r18669);
        mpfr_init_set_str(r18670, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r18671, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r18672, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r18673, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r18674, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r18675);
        mpfr_init(r18676);
        mpfr_init(r18677);
        mpfr_init(r18678);
        mpfr_init(r18679);
        mpfr_init(r18680);
        mpfr_init(r18681);
        mpfr_init(r18682);
        mpfr_init(r18683);
        mpfr_init(r18684);
        mpfr_init(r18685);
        mpfr_init(r18686);
        mpfr_init(r18687);
        mpfr_init(r18688);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r18665, x, MPFR_RNDN);
        mpfr_abs(r18666, r18665, MPFR_RNDN);
        mpfr_mul(r18667, r18664, r18666, MPFR_RNDN);
        mpfr_add(r18668, r18663, r18667, MPFR_RNDN);
        mpfr_div(r18669, r18663, r18668, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r18675, r18669, r18674, MPFR_RNDN);
        mpfr_add(r18676, r18673, r18675, MPFR_RNDN);
        mpfr_mul(r18677, r18669, r18676, MPFR_RNDN);
        mpfr_add(r18678, r18672, r18677, MPFR_RNDN);
        mpfr_mul(r18679, r18669, r18678, MPFR_RNDN);
        mpfr_add(r18680, r18671, r18679, MPFR_RNDN);
        mpfr_mul(r18681, r18669, r18680, MPFR_RNDN);
        mpfr_add(r18682, r18670, r18681, MPFR_RNDN);
        mpfr_mul(r18683, r18669, r18682, MPFR_RNDN);
        mpfr_mul(r18684, r18666, r18666, MPFR_RNDN);
        mpfr_neg(r18685, r18684, MPFR_RNDN);
        mpfr_exp(r18686, r18685, MPFR_RNDN);
        mpfr_mul(r18687, r18683, r18686, MPFR_RNDN);
        mpfr_sub(r18688, r18663, r18687, MPFR_RNDN);
        return mpfr_get_d(r18688, MPFR_RNDN);
}

static mpfr_t r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703, r18704, r18705, r18706, r18707, r18708, r18709, r18710, r18711, r18712, r18713;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18689, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18690, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18691, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18692);
        mpfr_init(r18693);
        mpfr_init(r18694);
        mpfr_init(r18695);
        mpfr_init_set_str(r18696, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18697);
        mpfr_init(r18698);
        mpfr_init(r18699);
        mpfr_init_set_str(r18700, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18701);
        mpfr_init_set_str(r18702, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18703);
        mpfr_init(r18704);
        mpfr_init_set_str(r18705, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18706);
        mpfr_init(r18707);
        mpfr_init(r18708);
        mpfr_init(r18709);
        mpfr_init(r18710);
        mpfr_init(r18711);
        mpfr_init(r18712);
        mpfr_init(r18713);
}

double f_fm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18692, x, MPFR_RNDN);
        mpfr_abs(r18693, r18692, MPFR_RNDN);
        mpfr_fma(r18694, r18691, r18693, r18689, MPFR_RNDN);
        mpfr_div(r18695, r18690, r18694, MPFR_RNDN);
        ;
        mpfr_add(r18697, r18695, r18696, MPFR_RNDN);
        mpfr_div(r18698, r18689, r18694, MPFR_RNDN);
        mpfr_div(r18699, r18698, r18694, MPFR_RNDN);
        ;
        mpfr_div(r18701, r18700, r18694, MPFR_RNDN);
        ;
        mpfr_fma(r18703, r18689, r18701, r18702, MPFR_RNDN);
        mpfr_fma(r18704, r18697, r18699, r18703, MPFR_RNDN);
        ;
        mpfr_div(r18706, r18705, r18694, MPFR_RNDN);
        mpfr_fma(r18707, r18704, r18699, r18706, MPFR_RNDN);
        mpfr_mul(r18708, r18693, r18693, MPFR_RNDN);
        mpfr_exp(r18709, r18708, MPFR_RNDN);
        mpfr_div(r18710, r18707, r18709, MPFR_RNDN);
        mpfr_sub(r18711, r18689, r18710, MPFR_RNDN);
        mpfr_log(r18712, r18711, MPFR_RNDN);
        mpfr_exp(r18713, r18712, MPFR_RNDN);
        return mpfr_get_d(r18713, MPFR_RNDN);
}

static mpfr_t r18714, r18715, r18716, r18717, r18718, r18719, r18720, r18721, r18722, r18723, r18724, r18725, r18726, r18727, r18728, r18729, r18730, r18731, r18732, r18733, r18734, r18735, r18736, r18737, r18738;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r18714, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r18715, "1.061405429", 10, MPFR_RNDN);
        mpfr_init_set_str(r18716, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r18717);
        mpfr_init(r18718);
        mpfr_init(r18719);
        mpfr_init(r18720);
        mpfr_init_set_str(r18721, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r18722);
        mpfr_init(r18723);
        mpfr_init(r18724);
        mpfr_init_set_str(r18725, "1.421413741", 10, MPFR_RNDN);
        mpfr_init(r18726);
        mpfr_init_set_str(r18727, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init(r18728);
        mpfr_init(r18729);
        mpfr_init_set_str(r18730, "0.254829592", 10, MPFR_RNDN);
        mpfr_init(r18731);
        mpfr_init(r18732);
        mpfr_init(r18733);
        mpfr_init(r18734);
        mpfr_init(r18735);
        mpfr_init(r18736);
        mpfr_init(r18737);
        mpfr_init(r18738);
}

double f_dm(double x) {
        ;
        ;
        ;
        mpfr_set_d(r18717, x, MPFR_RNDN);
        mpfr_abs(r18718, r18717, MPFR_RNDN);
        mpfr_fma(r18719, r18716, r18718, r18714, MPFR_RNDN);
        mpfr_div(r18720, r18715, r18719, MPFR_RNDN);
        ;
        mpfr_add(r18722, r18720, r18721, MPFR_RNDN);
        mpfr_div(r18723, r18714, r18719, MPFR_RNDN);
        mpfr_div(r18724, r18723, r18719, MPFR_RNDN);
        ;
        mpfr_div(r18726, r18725, r18719, MPFR_RNDN);
        ;
        mpfr_fma(r18728, r18714, r18726, r18727, MPFR_RNDN);
        mpfr_fma(r18729, r18722, r18724, r18728, MPFR_RNDN);
        ;
        mpfr_div(r18731, r18730, r18719, MPFR_RNDN);
        mpfr_fma(r18732, r18729, r18724, r18731, MPFR_RNDN);
        mpfr_mul(r18733, r18718, r18718, MPFR_RNDN);
        mpfr_exp(r18734, r18733, MPFR_RNDN);
        mpfr_div(r18735, r18732, r18734, MPFR_RNDN);
        mpfr_sub(r18736, r18714, r18735, MPFR_RNDN);
        mpfr_log(r18737, r18736, MPFR_RNDN);
        mpfr_exp(r18738, r18737, MPFR_RNDN);
        return mpfr_get_d(r18738, MPFR_RNDN);
}