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

char *name = "NMSE example 3.3";

double f_if(float x, float eps) {
        float r17516 = x;
        float r17517 = eps;
        float r17518 = r17516 + r17517;
        float r17519 = sin(r17518);
        float r17520 = sin(r17516);
        float r17521 = r17519 - r17520;
        return r17521;
}

double f_id(double x, double eps) {
        double r17522 = x;
        double r17523 = eps;
        double r17524 = r17522 + r17523;
        double r17525 = sin(r17524);
        double r17526 = sin(r17522);
        double r17527 = r17525 - r17526;
        return r17527;
}


double f_of(float x, float eps) {
        float r17528 = eps;
        float r17529 = -1.2324208637780791e-05f;
        bool r17530 = r17528 <= r17529;
        float r17531 = x;
        float r17532 = sin(r17531);
        float r17533 = cos(r17528);
        float r17534 = r17532 * r17533;
        float r17535 = cos(r17531);
        float r17536 = sin(r17528);
        float r17537 = r17535 * r17536;
        float r17538 = r17537 - r17532;
        float r17539 = r17534 + r17538;
        float r17540 = 0.15577836889185903f;
        bool r17541 = r17528 <= r17540;
        float r17542 = 0.16666666666666666f;
        float r17543 = r17528 * (r17528 * r17528);
        float r17544 = r17542 * r17543;
        float r17545 = r17528 - r17544;
        float r17546 = r17545 * r17535;
        float r17547 = 0.5f;
        float r17548 = r17532 * r17547;
        float r17549 = r17528 * r17528;
        float r17550 = r17548 * r17549;
        float r17551 = r17546 - r17550;
        float r17552 = r17541 ? r17551 : r17539;
        float r17553 = r17530 ? r17539 : r17552;
        return r17553;
}

double f_od(double x, double eps) {
        double r17554 = eps;
        double r17555 = -1.2324208637780791e-05;
        bool r17556 = r17554 <= r17555;
        double r17557 = x;
        double r17558 = sin(r17557);
        double r17559 = cos(r17554);
        double r17560 = r17558 * r17559;
        double r17561 = cos(r17557);
        double r17562 = sin(r17554);
        double r17563 = r17561 * r17562;
        double r17564 = r17563 - r17558;
        double r17565 = r17560 + r17564;
        double r17566 = 0.15577836889185903;
        bool r17567 = r17554 <= r17566;
        double r17568 = 0.16666666666666666;
        double r17569 = r17554 * (r17554 * r17554);
        double r17570 = r17568 * r17569;
        double r17571 = r17554 - r17570;
        double r17572 = r17571 * r17561;
        double r17573 = 0.5;
        double r17574 = r17558 * r17573;
        double r17575 = r17554 * r17554;
        double r17576 = r17574 * r17575;
        double r17577 = r17572 - r17576;
        double r17578 = r17567 ? r17577 : r17565;
        double r17579 = r17556 ? r17565 : r17578;
        return r17579;
}

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 r17580, r17581, r17582, r17583, r17584, r17585;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17580);
        mpfr_init(r17581);
        mpfr_init(r17582);
        mpfr_init(r17583);
        mpfr_init(r17584);
        mpfr_init(r17585);
}

double f_im(double x, double eps) {
        mpfr_set_d(r17580, x, MPFR_RNDN);
        mpfr_set_d(r17581, eps, MPFR_RNDN);
        mpfr_add(r17582, r17580, r17581, MPFR_RNDN);
        mpfr_sin(r17583, r17582, MPFR_RNDN);
        mpfr_sin(r17584, r17580, MPFR_RNDN);
        mpfr_sub(r17585, r17583, r17584, MPFR_RNDN);
        return mpfr_get_d(r17585, MPFR_RNDN);
}

static mpfr_t r17586, r17587, r17588, r17589, r17590, r17591, r17592, r17593, r17594, r17595, r17596, r17597, r17598, r17599, r17600, r17601, r17602, r17603, r17604, r17605, r17606, r17607, r17608, r17609, r17610, r17611;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17586);
        mpfr_init_set_str(r17587, "-1.2324208637780791e-05", 10, MPFR_RNDN);
        mpfr_init(r17588);
        mpfr_init(r17589);
        mpfr_init(r17590);
        mpfr_init(r17591);
        mpfr_init(r17592);
        mpfr_init(r17593);
        mpfr_init(r17594);
        mpfr_init(r17595);
        mpfr_init(r17596);
        mpfr_init(r17597);
        mpfr_init_set_str(r17598, "0.15577836889185903", 10, MPFR_RNDN);
        mpfr_init(r17599);
        mpfr_init_set_str(r17600, "1/6", 10, MPFR_RNDN);
        mpfr_init(r17601);
        mpfr_init(r17602);
        mpfr_init(r17603);
        mpfr_init(r17604);
        mpfr_init_set_str(r17605, "1/2", 10, MPFR_RNDN);
        mpfr_init(r17606);
        mpfr_init(r17607);
        mpfr_init(r17608);
        mpfr_init(r17609);
        mpfr_init(r17610);
        mpfr_init(r17611);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r17586, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r17588, mpfr_cmp(r17586, r17587) <= 0, MPFR_RNDN);
        mpfr_set_d(r17589, x, MPFR_RNDN);
        mpfr_sin(r17590, r17589, MPFR_RNDN);
        mpfr_cos(r17591, r17586, MPFR_RNDN);
        mpfr_mul(r17592, r17590, r17591, MPFR_RNDN);
        mpfr_cos(r17593, r17589, MPFR_RNDN);
        mpfr_sin(r17594, r17586, MPFR_RNDN);
        mpfr_mul(r17595, r17593, r17594, MPFR_RNDN);
        mpfr_sub(r17596, r17595, r17590, MPFR_RNDN);
        mpfr_add(r17597, r17592, r17596, MPFR_RNDN);
        ;
        mpfr_set_si(r17599, mpfr_cmp(r17586, r17598) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r17601, r17586, r17586, MPFR_RNDN); mpfr_mul(r17601, r17601, r17586, MPFR_RNDN);
        mpfr_mul(r17602, r17600, r17601, MPFR_RNDN);
        mpfr_sub(r17603, r17586, r17602, MPFR_RNDN);
        mpfr_mul(r17604, r17603, r17593, MPFR_RNDN);
        ;
        mpfr_mul(r17606, r17590, r17605, MPFR_RNDN);
        mpfr_sqr(r17607, r17586, MPFR_RNDN);
        mpfr_mul(r17608, r17606, r17607, MPFR_RNDN);
        mpfr_sub(r17609, r17604, r17608, MPFR_RNDN);
        if (mpfr_get_si(r17599, MPFR_RNDN)) { mpfr_set(r17610, r17609, MPFR_RNDN); } else { mpfr_set(r17610, r17597, MPFR_RNDN); };
        if (mpfr_get_si(r17588, MPFR_RNDN)) { mpfr_set(r17611, r17597, MPFR_RNDN); } else { mpfr_set(r17611, r17610, MPFR_RNDN); };
        return mpfr_get_d(r17611, MPFR_RNDN);
}

static mpfr_t r17612, r17613, r17614, r17615, r17616, r17617, r17618, r17619, r17620, r17621, r17622, r17623, r17624, r17625, r17626, r17627, r17628, r17629, r17630, r17631, r17632, r17633, r17634, r17635, r17636, r17637;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17612);
        mpfr_init_set_str(r17613, "-1.2324208637780791e-05", 10, MPFR_RNDN);
        mpfr_init(r17614);
        mpfr_init(r17615);
        mpfr_init(r17616);
        mpfr_init(r17617);
        mpfr_init(r17618);
        mpfr_init(r17619);
        mpfr_init(r17620);
        mpfr_init(r17621);
        mpfr_init(r17622);
        mpfr_init(r17623);
        mpfr_init_set_str(r17624, "0.15577836889185903", 10, MPFR_RNDN);
        mpfr_init(r17625);
        mpfr_init_set_str(r17626, "1/6", 10, MPFR_RNDN);
        mpfr_init(r17627);
        mpfr_init(r17628);
        mpfr_init(r17629);
        mpfr_init(r17630);
        mpfr_init_set_str(r17631, "1/2", 10, MPFR_RNDN);
        mpfr_init(r17632);
        mpfr_init(r17633);
        mpfr_init(r17634);
        mpfr_init(r17635);
        mpfr_init(r17636);
        mpfr_init(r17637);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r17612, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r17614, mpfr_cmp(r17612, r17613) <= 0, MPFR_RNDN);
        mpfr_set_d(r17615, x, MPFR_RNDN);
        mpfr_sin(r17616, r17615, MPFR_RNDN);
        mpfr_cos(r17617, r17612, MPFR_RNDN);
        mpfr_mul(r17618, r17616, r17617, MPFR_RNDN);
        mpfr_cos(r17619, r17615, MPFR_RNDN);
        mpfr_sin(r17620, r17612, MPFR_RNDN);
        mpfr_mul(r17621, r17619, r17620, MPFR_RNDN);
        mpfr_sub(r17622, r17621, r17616, MPFR_RNDN);
        mpfr_add(r17623, r17618, r17622, MPFR_RNDN);
        ;
        mpfr_set_si(r17625, mpfr_cmp(r17612, r17624) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r17627, r17612, r17612, MPFR_RNDN); mpfr_mul(r17627, r17627, r17612, MPFR_RNDN);
        mpfr_mul(r17628, r17626, r17627, MPFR_RNDN);
        mpfr_sub(r17629, r17612, r17628, MPFR_RNDN);
        mpfr_mul(r17630, r17629, r17619, MPFR_RNDN);
        ;
        mpfr_mul(r17632, r17616, r17631, MPFR_RNDN);
        mpfr_sqr(r17633, r17612, MPFR_RNDN);
        mpfr_mul(r17634, r17632, r17633, MPFR_RNDN);
        mpfr_sub(r17635, r17630, r17634, MPFR_RNDN);
        if (mpfr_get_si(r17625, MPFR_RNDN)) { mpfr_set(r17636, r17635, MPFR_RNDN); } else { mpfr_set(r17636, r17623, MPFR_RNDN); };
        if (mpfr_get_si(r17614, MPFR_RNDN)) { mpfr_set(r17637, r17623, MPFR_RNDN); } else { mpfr_set(r17637, r17636, MPFR_RNDN); };
        return mpfr_get_d(r17637, MPFR_RNDN);
}