#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 r17444 = 1.0f;
        float r17445 = atan2(1.0, 0.0);
        float r17446 = sqrt(r17445);
        float r17447 = r17444 / r17446;
        float r17448 = 2.0f;
        float r17449 = x;
        float r17450 = fabs(r17449);
        float r17451 = r17448 * r17450;
        float r17452 = 3.0f;
        float r17453 = r17448 / r17452;
        float r17454 = r17450 * r17450;
        float r17455 = r17454 * r17450;
        float r17456 = r17453 * r17455;
        float r17457 = r17451 + r17456;
        float r17458 = 5.0f;
        float r17459 = r17444 / r17458;
        float r17460 = r17455 * r17450;
        float r17461 = r17460 * r17450;
        float r17462 = r17459 * r17461;
        float r17463 = r17457 + r17462;
        float r17464 = 21.0f;
        float r17465 = r17444 / r17464;
        float r17466 = r17461 * r17450;
        float r17467 = r17466 * r17450;
        float r17468 = r17465 * r17467;
        float r17469 = r17463 + r17468;
        float r17470 = r17447 * r17469;
        float r17471 = fabs(r17470);
        return r17471;
}

double f_id(double x) {
        double r17472 = 1.0;
        double r17473 = atan2(1.0, 0.0);
        double r17474 = sqrt(r17473);
        double r17475 = r17472 / r17474;
        double r17476 = 2.0;
        double r17477 = x;
        double r17478 = fabs(r17477);
        double r17479 = r17476 * r17478;
        double r17480 = 3.0;
        double r17481 = r17476 / r17480;
        double r17482 = r17478 * r17478;
        double r17483 = r17482 * r17478;
        double r17484 = r17481 * r17483;
        double r17485 = r17479 + r17484;
        double r17486 = 5.0;
        double r17487 = r17472 / r17486;
        double r17488 = r17483 * r17478;
        double r17489 = r17488 * r17478;
        double r17490 = r17487 * r17489;
        double r17491 = r17485 + r17490;
        double r17492 = 21.0;
        double r17493 = r17472 / r17492;
        double r17494 = r17489 * r17478;
        double r17495 = r17494 * r17478;
        double r17496 = r17493 * r17495;
        double r17497 = r17491 + r17496;
        double r17498 = r17475 * r17497;
        double r17499 = fabs(r17498);
        return r17499;
}


double f_of(float x) {
        float r17500 = x;
        float r17501 = fabs(r17500);
        float r17502 = 5.0f;
        float r17503 = r17501 / r17502;
        float r17504 = r17501 * (r17501 * r17501);
        float r17505 = r17503 * r17504;
        float r17506 = 2.0f;
        float r17507 = 3.0f;
        float r17508 = r17506 / r17507;
        float r17509 = r17506 * r17501;
        float r17510 = fma(r17508, r17504, r17509);
        float r17511 = fma(r17505, r17501, r17510);
        float r17512 = 0.0476190485060215f;
        float r17513 = r17501 * r17501;
        float r17514 = r17513 * (r17513 * r17513);
        float r17515 = r17514 * r17501;
        float r17516 = r17512 * r17515;
        float r17517 = r17511 + r17516;
        float r17518 = atan2(1.0, 0.0);
        float r17519 = sqrt(r17518);
        float r17520 = r17517 / r17519;
        float r17521 = fabs(r17520);
        return r17521;
}

double f_od(double x) {
        double r17522 = x;
        double r17523 = fabs(r17522);
        double r17524 = 5.0;
        double r17525 = r17523 / r17524;
        double r17526 = r17523 * (r17523 * r17523);
        double r17527 = r17525 * r17526;
        double r17528 = 2.0;
        double r17529 = 3.0;
        double r17530 = r17528 / r17529;
        double r17531 = r17528 * r17523;
        double r17532 = fma(r17530, r17526, r17531);
        double r17533 = fma(r17527, r17523, r17532);
        double r17534 = 0.0476190485060215;
        double r17535 = r17523 * r17523;
        double r17536 = r17535 * (r17535 * r17535);
        double r17537 = r17536 * r17523;
        double r17538 = r17534 * r17537;
        double r17539 = r17533 + r17538;
        double r17540 = atan2(1.0, 0.0);
        double r17541 = sqrt(r17540);
        double r17542 = r17539 / r17541;
        double r17543 = fabs(r17542);
        return r17543;
}

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 r17544, r17545, r17546, r17547, r17548, r17549, r17550, r17551, r17552, r17553, r17554, r17555, r17556, r17557, r17558, r17559, r17560, r17561, r17562, r17563, r17564, r17565, r17566, r17567, r17568, r17569, r17570, r17571;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17544, "1", 10, MPFR_RNDN);
        mpfr_init(r17545);
        mpfr_init(r17546);
        mpfr_init(r17547);
        mpfr_init_set_str(r17548, "2", 10, MPFR_RNDN);
        mpfr_init(r17549);
        mpfr_init(r17550);
        mpfr_init(r17551);
        mpfr_init_set_str(r17552, "3", 10, MPFR_RNDN);
        mpfr_init(r17553);
        mpfr_init(r17554);
        mpfr_init(r17555);
        mpfr_init(r17556);
        mpfr_init(r17557);
        mpfr_init_set_str(r17558, "5", 10, MPFR_RNDN);
        mpfr_init(r17559);
        mpfr_init(r17560);
        mpfr_init(r17561);
        mpfr_init(r17562);
        mpfr_init(r17563);
        mpfr_init_set_str(r17564, "21", 10, MPFR_RNDN);
        mpfr_init(r17565);
        mpfr_init(r17566);
        mpfr_init(r17567);
        mpfr_init(r17568);
        mpfr_init(r17569);
        mpfr_init(r17570);
        mpfr_init(r17571);
}

double f_im(double x) {
        ;
        mpfr_const_pi(r17545, MPFR_RNDN);
        mpfr_sqrt(r17546, r17545, MPFR_RNDN);
        mpfr_div(r17547, r17544, r17546, MPFR_RNDN);
        ;
        mpfr_set_d(r17549, x, MPFR_RNDN);
        mpfr_abs(r17550, r17549, MPFR_RNDN);
        mpfr_mul(r17551, r17548, r17550, MPFR_RNDN);
        ;
        mpfr_div(r17553, r17548, r17552, MPFR_RNDN);
        mpfr_mul(r17554, r17550, r17550, MPFR_RNDN);
        mpfr_mul(r17555, r17554, r17550, MPFR_RNDN);
        mpfr_mul(r17556, r17553, r17555, MPFR_RNDN);
        mpfr_add(r17557, r17551, r17556, MPFR_RNDN);
        ;
        mpfr_div(r17559, r17544, r17558, MPFR_RNDN);
        mpfr_mul(r17560, r17555, r17550, MPFR_RNDN);
        mpfr_mul(r17561, r17560, r17550, MPFR_RNDN);
        mpfr_mul(r17562, r17559, r17561, MPFR_RNDN);
        mpfr_add(r17563, r17557, r17562, MPFR_RNDN);
        ;
        mpfr_div(r17565, r17544, r17564, MPFR_RNDN);
        mpfr_mul(r17566, r17561, r17550, MPFR_RNDN);
        mpfr_mul(r17567, r17566, r17550, MPFR_RNDN);
        mpfr_mul(r17568, r17565, r17567, MPFR_RNDN);
        mpfr_add(r17569, r17563, r17568, MPFR_RNDN);
        mpfr_mul(r17570, r17547, r17569, MPFR_RNDN);
        mpfr_abs(r17571, r17570, MPFR_RNDN);
        return mpfr_get_d(r17571, MPFR_RNDN);
}

static mpfr_t r17572, r17573, r17574, r17575, r17576, r17577, r17578, r17579, r17580, r17581, r17582, r17583, r17584, r17585, r17586, r17587, r17588, r17589, r17590, r17591, r17592, r17593;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17572);
        mpfr_init(r17573);
        mpfr_init_set_str(r17574, "5", 10, MPFR_RNDN);
        mpfr_init(r17575);
        mpfr_init(r17576);
        mpfr_init(r17577);
        mpfr_init_set_str(r17578, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r17579, "3", 10, MPFR_RNDN);
        mpfr_init(r17580);
        mpfr_init(r17581);
        mpfr_init(r17582);
        mpfr_init(r17583);
        mpfr_init_set_str(r17584, "1/21", 10, MPFR_RNDN);
        mpfr_init(r17585);
        mpfr_init(r17586);
        mpfr_init(r17587);
        mpfr_init(r17588);
        mpfr_init(r17589);
        mpfr_init(r17590);
        mpfr_init(r17591);
        mpfr_init(r17592);
        mpfr_init(r17593);
}

double f_fm(double x) {
        mpfr_set_d(r17572, x, MPFR_RNDN);
        mpfr_abs(r17573, r17572, MPFR_RNDN);
        ;
        mpfr_div(r17575, r17573, r17574, MPFR_RNDN);
        mpfr_mul(r17576, r17573, r17573, MPFR_RNDN); mpfr_mul(r17576, r17576, r17573, MPFR_RNDN);
        mpfr_mul(r17577, r17575, r17576, MPFR_RNDN);
        ;
        ;
        mpfr_div(r17580, r17578, r17579, MPFR_RNDN);
        mpfr_mul(r17581, r17578, r17573, MPFR_RNDN);
        mpfr_fma(r17582, r17580, r17576, r17581, MPFR_RNDN);
        mpfr_fma(r17583, r17577, r17573, r17582, MPFR_RNDN);
        ;
        mpfr_sqr(r17585, r17573, MPFR_RNDN);
        mpfr_mul(r17586, r17585, r17585, MPFR_RNDN); mpfr_mul(r17586, r17586, r17585, MPFR_RNDN);
        mpfr_mul(r17587, r17586, r17573, MPFR_RNDN);
        mpfr_mul(r17588, r17584, r17587, MPFR_RNDN);
        mpfr_add(r17589, r17583, r17588, MPFR_RNDN);
        mpfr_const_pi(r17590, MPFR_RNDN);
        mpfr_sqrt(r17591, r17590, MPFR_RNDN);
        mpfr_div(r17592, r17589, r17591, MPFR_RNDN);
        mpfr_abs(r17593, r17592, MPFR_RNDN);
        return mpfr_get_d(r17593, MPFR_RNDN);
}

static mpfr_t r17594, r17595, r17596, r17597, r17598, r17599, r17600, r17601, r17602, r17603, r17604, r17605, r17606, r17607, r17608, r17609, r17610, r17611, r17612, r17613, r17614, r17615;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17594);
        mpfr_init(r17595);
        mpfr_init_set_str(r17596, "5", 10, MPFR_RNDN);
        mpfr_init(r17597);
        mpfr_init(r17598);
        mpfr_init(r17599);
        mpfr_init_set_str(r17600, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r17601, "3", 10, MPFR_RNDN);
        mpfr_init(r17602);
        mpfr_init(r17603);
        mpfr_init(r17604);
        mpfr_init(r17605);
        mpfr_init_set_str(r17606, "1/21", 10, MPFR_RNDN);
        mpfr_init(r17607);
        mpfr_init(r17608);
        mpfr_init(r17609);
        mpfr_init(r17610);
        mpfr_init(r17611);
        mpfr_init(r17612);
        mpfr_init(r17613);
        mpfr_init(r17614);
        mpfr_init(r17615);
}

double f_dm(double x) {
        mpfr_set_d(r17594, x, MPFR_RNDN);
        mpfr_abs(r17595, r17594, MPFR_RNDN);
        ;
        mpfr_div(r17597, r17595, r17596, MPFR_RNDN);
        mpfr_mul(r17598, r17595, r17595, MPFR_RNDN); mpfr_mul(r17598, r17598, r17595, MPFR_RNDN);
        mpfr_mul(r17599, r17597, r17598, MPFR_RNDN);
        ;
        ;
        mpfr_div(r17602, r17600, r17601, MPFR_RNDN);
        mpfr_mul(r17603, r17600, r17595, MPFR_RNDN);
        mpfr_fma(r17604, r17602, r17598, r17603, MPFR_RNDN);
        mpfr_fma(r17605, r17599, r17595, r17604, MPFR_RNDN);
        ;
        mpfr_sqr(r17607, r17595, MPFR_RNDN);
        mpfr_mul(r17608, r17607, r17607, MPFR_RNDN); mpfr_mul(r17608, r17608, r17607, MPFR_RNDN);
        mpfr_mul(r17609, r17608, r17595, MPFR_RNDN);
        mpfr_mul(r17610, r17606, r17609, MPFR_RNDN);
        mpfr_add(r17611, r17605, r17610, MPFR_RNDN);
        mpfr_const_pi(r17612, MPFR_RNDN);
        mpfr_sqrt(r17613, r17612, MPFR_RNDN);
        mpfr_div(r17614, r17611, r17613, MPFR_RNDN);
        mpfr_abs(r17615, r17614, MPFR_RNDN);
        return mpfr_get_d(r17615, MPFR_RNDN);
}