#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5";

double f_if(float x) {
        float r17560 = 1.0f;
        float r17561 = atan2(1.0, 0.0);
        float r17562 = sqrt(r17561);
        float r17563 = r17560 / r17562;
        float r17564 = 2.0f;
        float r17565 = x;
        float r17566 = fabs(r17565);
        float r17567 = r17564 * r17566;
        float r17568 = 3.0f;
        float r17569 = r17564 / r17568;
        float r17570 = r17566 * r17566;
        float r17571 = r17570 * r17566;
        float r17572 = r17569 * r17571;
        float r17573 = r17567 + r17572;
        float r17574 = 5.0f;
        float r17575 = r17560 / r17574;
        float r17576 = r17571 * r17566;
        float r17577 = r17576 * r17566;
        float r17578 = r17575 * r17577;
        float r17579 = r17573 + r17578;
        float r17580 = 21.0f;
        float r17581 = r17560 / r17580;
        float r17582 = r17577 * r17566;
        float r17583 = r17582 * r17566;
        float r17584 = r17581 * r17583;
        float r17585 = r17579 + r17584;
        float r17586 = r17563 * r17585;
        float r17587 = fabs(r17586);
        return r17587;
}

double f_id(double x) {
        double r17588 = 1.0;
        double r17589 = atan2(1.0, 0.0);
        double r17590 = sqrt(r17589);
        double r17591 = r17588 / r17590;
        double r17592 = 2.0;
        double r17593 = x;
        double r17594 = fabs(r17593);
        double r17595 = r17592 * r17594;
        double r17596 = 3.0;
        double r17597 = r17592 / r17596;
        double r17598 = r17594 * r17594;
        double r17599 = r17598 * r17594;
        double r17600 = r17597 * r17599;
        double r17601 = r17595 + r17600;
        double r17602 = 5.0;
        double r17603 = r17588 / r17602;
        double r17604 = r17599 * r17594;
        double r17605 = r17604 * r17594;
        double r17606 = r17603 * r17605;
        double r17607 = r17601 + r17606;
        double r17608 = 21.0;
        double r17609 = r17588 / r17608;
        double r17610 = r17605 * r17594;
        double r17611 = r17610 * r17594;
        double r17612 = r17609 * r17611;
        double r17613 = r17607 + r17612;
        double r17614 = r17591 * r17613;
        double r17615 = fabs(r17614);
        return r17615;
}


double f_of(float x) {
        float r17616 = 2.0f;
        float r17617 = x;
        float r17618 = fabs(r17617);
        float r17619 = r17616 * r17618;
        float r17620 = 0.6666666865348816f;
        float r17621 = 3.0f;
        float r17622 = pow(r17618, r17621);
        float r17623 = r17620 * r17622;
        float r17624 = r17619 + r17623;
        float r17625 = r17618 * (r17618 * r17618);
        float r17626 = r17625 * r17625;
        float r17627 = 21.0f;
        float r17628 = r17627 / r17618;
        float r17629 = r17626 / r17628;
        float r17630 = r17618 * r17618;
        float r17631 = r17625 * r17630;
        float r17632 = 5.0f;
        float r17633 = r17631 / r17632;
        float r17634 = r17629 + r17633;
        float r17635 = r17624 + r17634;
        float r17636 = atan2(1.0, 0.0);
        float r17637 = sqrt(r17636);
        float r17638 = r17635 / r17637;
        float r17639 = fabs(r17638);
        return r17639;
}

double f_od(double x) {
        double r17640 = 2.0;
        double r17641 = x;
        double r17642 = fabs(r17641);
        double r17643 = r17640 * r17642;
        double r17644 = 0.6666666865348816;
        double r17645 = 3.0;
        double r17646 = pow(r17642, r17645);
        double r17647 = r17644 * r17646;
        double r17648 = r17643 + r17647;
        double r17649 = r17642 * (r17642 * r17642);
        double r17650 = r17649 * r17649;
        double r17651 = 21.0;
        double r17652 = r17651 / r17642;
        double r17653 = r17650 / r17652;
        double r17654 = r17642 * r17642;
        double r17655 = r17649 * r17654;
        double r17656 = 5.0;
        double r17657 = r17655 / r17656;
        double r17658 = r17653 + r17657;
        double r17659 = r17648 + r17658;
        double r17660 = atan2(1.0, 0.0);
        double r17661 = sqrt(r17660);
        double r17662 = r17659 / r17661;
        double r17663 = fabs(r17662);
        return r17663;
}

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 r17664, r17665, r17666, r17667, r17668, r17669, r17670, r17671, r17672, r17673, r17674, r17675, r17676, r17677, r17678, r17679, r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17664, "1", 10, MPFR_RNDN);
        mpfr_init(r17665);
        mpfr_init(r17666);
        mpfr_init(r17667);
        mpfr_init_set_str(r17668, "2", 10, MPFR_RNDN);
        mpfr_init(r17669);
        mpfr_init(r17670);
        mpfr_init(r17671);
        mpfr_init_set_str(r17672, "3", 10, MPFR_RNDN);
        mpfr_init(r17673);
        mpfr_init(r17674);
        mpfr_init(r17675);
        mpfr_init(r17676);
        mpfr_init(r17677);
        mpfr_init_set_str(r17678, "5", 10, MPFR_RNDN);
        mpfr_init(r17679);
        mpfr_init(r17680);
        mpfr_init(r17681);
        mpfr_init(r17682);
        mpfr_init(r17683);
        mpfr_init_set_str(r17684, "21", 10, MPFR_RNDN);
        mpfr_init(r17685);
        mpfr_init(r17686);
        mpfr_init(r17687);
        mpfr_init(r17688);
        mpfr_init(r17689);
        mpfr_init(r17690);
        mpfr_init(r17691);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r17665, MPFR_RNDN);
        mpfr_sqrt(r17666, r17665, MPFR_RNDN);
        mpfr_div(r17667, r17664, r17666, MPFR_RNDN);
        ;
        mpfr_set_d(r17669, x, MPFR_RNDN);
        mpfr_abs(r17670, r17669, MPFR_RNDN);
        mpfr_mul(r17671, r17668, r17670, MPFR_RNDN);
        ;
        mpfr_div(r17673, r17668, r17672, MPFR_RNDN);
        mpfr_mul(r17674, r17670, r17670, MPFR_RNDN);
        mpfr_mul(r17675, r17674, r17670, MPFR_RNDN);
        mpfr_mul(r17676, r17673, r17675, MPFR_RNDN);
        mpfr_add(r17677, r17671, r17676, MPFR_RNDN);
        ;
        mpfr_div(r17679, r17664, r17678, MPFR_RNDN);
        mpfr_mul(r17680, r17675, r17670, MPFR_RNDN);
        mpfr_mul(r17681, r17680, r17670, MPFR_RNDN);
        mpfr_mul(r17682, r17679, r17681, MPFR_RNDN);
        mpfr_add(r17683, r17677, r17682, MPFR_RNDN);
        ;
        mpfr_div(r17685, r17664, r17684, MPFR_RNDN);
        mpfr_mul(r17686, r17681, r17670, MPFR_RNDN);
        mpfr_mul(r17687, r17686, r17670, MPFR_RNDN);
        mpfr_mul(r17688, r17685, r17687, MPFR_RNDN);
        mpfr_add(r17689, r17683, r17688, MPFR_RNDN);
        mpfr_mul(r17690, r17667, r17689, MPFR_RNDN);
        mpfr_abs(r17691, r17690, MPFR_RNDN);
        return mpfr_get_d(r17691, MPFR_RNDN);
}

static mpfr_t r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703, r17704, r17705, r17706, r17707, r17708, r17709, r17710, r17711, r17712, r17713, r17714, r17715;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17692, "2", 10, MPFR_RNDN);
        mpfr_init(r17693);
        mpfr_init(r17694);
        mpfr_init(r17695);
        mpfr_init_set_str(r17696, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r17697, "3", 10, MPFR_RNDN);
        mpfr_init(r17698);
        mpfr_init(r17699);
        mpfr_init(r17700);
        mpfr_init(r17701);
        mpfr_init(r17702);
        mpfr_init_set_str(r17703, "21", 10, MPFR_RNDN);
        mpfr_init(r17704);
        mpfr_init(r17705);
        mpfr_init(r17706);
        mpfr_init(r17707);
        mpfr_init_set_str(r17708, "5", 10, MPFR_RNDN);
        mpfr_init(r17709);
        mpfr_init(r17710);
        mpfr_init(r17711);
        mpfr_init(r17712);
        mpfr_init(r17713);
        mpfr_init(r17714);
        mpfr_init(r17715);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r17693, x, MPFR_RNDN);
        mpfr_abs(r17694, r17693, MPFR_RNDN);
        mpfr_mul(r17695, r17692, r17694, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r17698, r17694, r17697, MPFR_RNDN);
        mpfr_mul(r17699, r17696, r17698, MPFR_RNDN);
        mpfr_add(r17700, r17695, r17699, MPFR_RNDN);
        mpfr_mul(r17701, r17694, r17694, MPFR_RNDN); mpfr_mul(r17701, r17701, r17694, MPFR_RNDN);
        mpfr_sqr(r17702, r17701, MPFR_RNDN);
        ;
        mpfr_div(r17704, r17703, r17694, MPFR_RNDN);
        mpfr_div(r17705, r17702, r17704, MPFR_RNDN);
        mpfr_sqr(r17706, r17694, MPFR_RNDN);
        mpfr_mul(r17707, r17701, r17706, MPFR_RNDN);
        ;
        mpfr_div(r17709, r17707, r17708, MPFR_RNDN);
        mpfr_add(r17710, r17705, r17709, MPFR_RNDN);
        mpfr_add(r17711, r17700, r17710, MPFR_RNDN);
        mpfr_const_pi(r17712, MPFR_RNDN);
        mpfr_sqrt(r17713, r17712, MPFR_RNDN);
        mpfr_div(r17714, r17711, r17713, MPFR_RNDN);
        mpfr_abs(r17715, r17714, MPFR_RNDN);
        return mpfr_get_d(r17715, MPFR_RNDN);
}

static mpfr_t r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732, r17733, r17734, r17735, r17736, r17737, r17738, r17739;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17716, "2", 10, MPFR_RNDN);
        mpfr_init(r17717);
        mpfr_init(r17718);
        mpfr_init(r17719);
        mpfr_init_set_str(r17720, "2/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r17721, "3", 10, MPFR_RNDN);
        mpfr_init(r17722);
        mpfr_init(r17723);
        mpfr_init(r17724);
        mpfr_init(r17725);
        mpfr_init(r17726);
        mpfr_init_set_str(r17727, "21", 10, MPFR_RNDN);
        mpfr_init(r17728);
        mpfr_init(r17729);
        mpfr_init(r17730);
        mpfr_init(r17731);
        mpfr_init_set_str(r17732, "5", 10, MPFR_RNDN);
        mpfr_init(r17733);
        mpfr_init(r17734);
        mpfr_init(r17735);
        mpfr_init(r17736);
        mpfr_init(r17737);
        mpfr_init(r17738);
        mpfr_init(r17739);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r17717, x, MPFR_RNDN);
        mpfr_abs(r17718, r17717, MPFR_RNDN);
        mpfr_mul(r17719, r17716, r17718, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r17722, r17718, r17721, MPFR_RNDN);
        mpfr_mul(r17723, r17720, r17722, MPFR_RNDN);
        mpfr_add(r17724, r17719, r17723, MPFR_RNDN);
        mpfr_mul(r17725, r17718, r17718, MPFR_RNDN); mpfr_mul(r17725, r17725, r17718, MPFR_RNDN);
        mpfr_sqr(r17726, r17725, MPFR_RNDN);
        ;
        mpfr_div(r17728, r17727, r17718, MPFR_RNDN);
        mpfr_div(r17729, r17726, r17728, MPFR_RNDN);
        mpfr_sqr(r17730, r17718, MPFR_RNDN);
        mpfr_mul(r17731, r17725, r17730, MPFR_RNDN);
        ;
        mpfr_div(r17733, r17731, r17732, MPFR_RNDN);
        mpfr_add(r17734, r17729, r17733, MPFR_RNDN);
        mpfr_add(r17735, r17724, r17734, MPFR_RNDN);
        mpfr_const_pi(r17736, MPFR_RNDN);
        mpfr_sqrt(r17737, r17736, MPFR_RNDN);
        mpfr_div(r17738, r17735, r17737, MPFR_RNDN);
        mpfr_abs(r17739, r17738, MPFR_RNDN);
        return mpfr_get_d(r17739, MPFR_RNDN);
}