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

char *name = "2sin (example 3.3)";

double f_if(float x, float eps) {
        float r20155 = x;
        float r20156 = eps;
        float r20157 = r20155 + r20156;
        float r20158 = sin(r20157);
        float r20159 = sin(r20155);
        float r20160 = r20158 - r20159;
        return r20160;
}

double f_id(double x, double eps) {
        double r20161 = x;
        double r20162 = eps;
        double r20163 = r20161 + r20162;
        double r20164 = sin(r20163);
        double r20165 = sin(r20161);
        double r20166 = r20164 - r20165;
        return r20166;
}


double f_of(float x, float eps) {
        float r20167 = eps;
        float r20168 = -3.645937152382937e+19f;
        bool r20169 = r20167 <= r20168;
        float r20170 = x;
        float r20171 = sin(r20170);
        float r20172 = cos(r20167);
        float r20173 = r20171 * r20172;
        float r20174 = cos(r20170);
        float r20175 = sin(r20167);
        float r20176 = r20174 * r20175;
        float r20177 = r20173 + r20176;
        float r20178 = r20177 - r20171;
        float r20179 = 6.326235572596747e-15f;
        bool r20180 = r20167 <= r20179;
        float r20181 = 2.0f;
        float r20182 = r20167 / r20181;
        float r20183 = sin(r20182);
        float r20184 = r20170 + r20167;
        float r20185 = r20184 + r20170;
        float r20186 = r20185 / r20181;
        float r20187 = cos(r20186);
        float r20188 = r20183 * r20187;
        float r20189 = r20181 * r20188;
        float r20190 = r20180 ? r20189 : r20178;
        float r20191 = r20169 ? r20178 : r20190;
        return r20191;
}

double f_od(double x, double eps) {
        double r20192 = eps;
        double r20193 = -3.645937152382937e+19;
        bool r20194 = r20192 <= r20193;
        double r20195 = x;
        double r20196 = sin(r20195);
        double r20197 = cos(r20192);
        double r20198 = r20196 * r20197;
        double r20199 = cos(r20195);
        double r20200 = sin(r20192);
        double r20201 = r20199 * r20200;
        double r20202 = r20198 + r20201;
        double r20203 = r20202 - r20196;
        double r20204 = 6.326235572596747e-15;
        bool r20205 = r20192 <= r20204;
        double r20206 = 2.0;
        double r20207 = r20192 / r20206;
        double r20208 = sin(r20207);
        double r20209 = r20195 + r20192;
        double r20210 = r20209 + r20195;
        double r20211 = r20210 / r20206;
        double r20212 = cos(r20211);
        double r20213 = r20208 * r20212;
        double r20214 = r20206 * r20213;
        double r20215 = r20205 ? r20214 : r20203;
        double r20216 = r20194 ? r20203 : r20215;
        return r20216;
}

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 r20217, r20218, r20219, r20220, r20221, r20222;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20217);
        mpfr_init(r20218);
        mpfr_init(r20219);
        mpfr_init(r20220);
        mpfr_init(r20221);
        mpfr_init(r20222);
}

double f_im(double x, double eps) {
        mpfr_set_d(r20217, x, MPFR_RNDN);
        mpfr_set_d(r20218, eps, MPFR_RNDN);
        mpfr_add(r20219, r20217, r20218, MPFR_RNDN);
        mpfr_sin(r20220, r20219, MPFR_RNDN);
        mpfr_sin(r20221, r20217, MPFR_RNDN);
        mpfr_sub(r20222, r20220, r20221, MPFR_RNDN);
        return mpfr_get_d(r20222, MPFR_RNDN);
}

static mpfr_t r20223, r20224, r20225, r20226, r20227, r20228, r20229, r20230, r20231, r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242, r20243, r20244, r20245, r20246, r20247;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20223);
        mpfr_init_set_str(r20224, "-3.645937152382937e+19", 10, MPFR_RNDN);
        mpfr_init(r20225);
        mpfr_init(r20226);
        mpfr_init(r20227);
        mpfr_init(r20228);
        mpfr_init(r20229);
        mpfr_init(r20230);
        mpfr_init(r20231);
        mpfr_init(r20232);
        mpfr_init(r20233);
        mpfr_init(r20234);
        mpfr_init_set_str(r20235, "6.326235572596747e-15", 10, MPFR_RNDN);
        mpfr_init(r20236);
        mpfr_init_set_str(r20237, "2", 10, MPFR_RNDN);
        mpfr_init(r20238);
        mpfr_init(r20239);
        mpfr_init(r20240);
        mpfr_init(r20241);
        mpfr_init(r20242);
        mpfr_init(r20243);
        mpfr_init(r20244);
        mpfr_init(r20245);
        mpfr_init(r20246);
        mpfr_init(r20247);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r20223, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20225, mpfr_cmp(r20223, r20224) <= 0, MPFR_RNDN);
        mpfr_set_d(r20226, x, MPFR_RNDN);
        mpfr_sin(r20227, r20226, MPFR_RNDN);
        mpfr_cos(r20228, r20223, MPFR_RNDN);
        mpfr_mul(r20229, r20227, r20228, MPFR_RNDN);
        mpfr_cos(r20230, r20226, MPFR_RNDN);
        mpfr_sin(r20231, r20223, MPFR_RNDN);
        mpfr_mul(r20232, r20230, r20231, MPFR_RNDN);
        mpfr_add(r20233, r20229, r20232, MPFR_RNDN);
        mpfr_sub(r20234, r20233, r20227, MPFR_RNDN);
        ;
        mpfr_set_si(r20236, mpfr_cmp(r20223, r20235) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r20238, r20223, r20237, MPFR_RNDN);
        mpfr_sin(r20239, r20238, MPFR_RNDN);
        mpfr_add(r20240, r20226, r20223, MPFR_RNDN);
        mpfr_add(r20241, r20240, r20226, MPFR_RNDN);
        mpfr_div(r20242, r20241, r20237, MPFR_RNDN);
        mpfr_cos(r20243, r20242, MPFR_RNDN);
        mpfr_mul(r20244, r20239, r20243, MPFR_RNDN);
        mpfr_mul(r20245, r20237, r20244, MPFR_RNDN);
        if (mpfr_get_si(r20236, MPFR_RNDN)) { mpfr_set(r20246, r20245, MPFR_RNDN); } else { mpfr_set(r20246, r20234, MPFR_RNDN); };
        if (mpfr_get_si(r20225, MPFR_RNDN)) { mpfr_set(r20247, r20234, MPFR_RNDN); } else { mpfr_set(r20247, r20246, MPFR_RNDN); };
        return mpfr_get_d(r20247, MPFR_RNDN);
}

static mpfr_t r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267, r20268, r20269, r20270, r20271, r20272;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20248);
        mpfr_init_set_str(r20249, "-3.645937152382937e+19", 10, MPFR_RNDN);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
        mpfr_init(r20258);
        mpfr_init(r20259);
        mpfr_init_set_str(r20260, "6.326235572596747e-15", 10, MPFR_RNDN);
        mpfr_init(r20261);
        mpfr_init_set_str(r20262, "2", 10, MPFR_RNDN);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init(r20265);
        mpfr_init(r20266);
        mpfr_init(r20267);
        mpfr_init(r20268);
        mpfr_init(r20269);
        mpfr_init(r20270);
        mpfr_init(r20271);
        mpfr_init(r20272);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r20248, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20250, mpfr_cmp(r20248, r20249) <= 0, MPFR_RNDN);
        mpfr_set_d(r20251, x, MPFR_RNDN);
        mpfr_sin(r20252, r20251, MPFR_RNDN);
        mpfr_cos(r20253, r20248, MPFR_RNDN);
        mpfr_mul(r20254, r20252, r20253, MPFR_RNDN);
        mpfr_cos(r20255, r20251, MPFR_RNDN);
        mpfr_sin(r20256, r20248, MPFR_RNDN);
        mpfr_mul(r20257, r20255, r20256, MPFR_RNDN);
        mpfr_add(r20258, r20254, r20257, MPFR_RNDN);
        mpfr_sub(r20259, r20258, r20252, MPFR_RNDN);
        ;
        mpfr_set_si(r20261, mpfr_cmp(r20248, r20260) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r20263, r20248, r20262, MPFR_RNDN);
        mpfr_sin(r20264, r20263, MPFR_RNDN);
        mpfr_add(r20265, r20251, r20248, MPFR_RNDN);
        mpfr_add(r20266, r20265, r20251, MPFR_RNDN);
        mpfr_div(r20267, r20266, r20262, MPFR_RNDN);
        mpfr_cos(r20268, r20267, MPFR_RNDN);
        mpfr_mul(r20269, r20264, r20268, MPFR_RNDN);
        mpfr_mul(r20270, r20262, r20269, MPFR_RNDN);
        if (mpfr_get_si(r20261, MPFR_RNDN)) { mpfr_set(r20271, r20270, MPFR_RNDN); } else { mpfr_set(r20271, r20259, MPFR_RNDN); };
        if (mpfr_get_si(r20250, MPFR_RNDN)) { mpfr_set(r20272, r20259, MPFR_RNDN); } else { mpfr_set(r20272, r20271, MPFR_RNDN); };
        return mpfr_get_d(r20272, MPFR_RNDN);
}

