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

char *name = "2cos (problem 3.3.5)";

double f_if(float x, float eps) {
        float r21331 = x;
        float r21332 = eps;
        float r21333 = r21331 + r21332;
        float r21334 = cos(r21333);
        float r21335 = cos(r21331);
        float r21336 = r21334 - r21335;
        return r21336;
}

double f_id(double x, double eps) {
        double r21337 = x;
        double r21338 = eps;
        double r21339 = r21337 + r21338;
        double r21340 = cos(r21339);
        double r21341 = cos(r21337);
        double r21342 = r21340 - r21341;
        return r21342;
}


double f_of(float x, float eps) {
        float r21343 = -2;
        float r21344 = eps;
        float r21345 = 2;
        float r21346 = r21344 / r21345;
        float r21347 = sin(r21346);
        float r21348 = x;
        float r21349 = fma(r21345, r21348, r21344);
        float r21350 = r21349 / r21345;
        float r21351 = sin(r21350);
        float r21352 = cbrt(r21351);
        float r21353 = r21352 * r21352;
        float r21354 = r21347 * r21353;
        float r21355 = r21354 * r21352;
        float r21356 = r21343 * r21355;
        float r21357 = -4.960559738300167e-06;
        bool r21358 = r21356 <= r21357;
        float r21359 = cos(r21348);
        float r21360 = cos(r21344);
        float r21361 = r21359 * r21360;
        float r21362 = sin(r21348);
        float r21363 = sin(r21344);
        float r21364 = r21362 * r21363;
        float r21365 = r21361 - r21364;
        float r21366 = r21365 - r21359;
        float r21367 = 0.014495001146888918;
        bool r21368 = r21356 <= r21367;
        float r21369 = r21347 * r21351;
        float r21370 = r21343 * r21369;
        float r21371 = fma(r21363, r21362, r21359);
        float r21372 = r21361 - r21371;
        float r21373 = r21368 ? r21370 : r21372;
        float r21374 = r21358 ? r21366 : r21373;
        return r21374;
}

double f_od(double x, double eps) {
        double r21375 = -2;
        double r21376 = eps;
        double r21377 = 2;
        double r21378 = r21376 / r21377;
        double r21379 = sin(r21378);
        double r21380 = x;
        double r21381 = fma(r21377, r21380, r21376);
        double r21382 = r21381 / r21377;
        double r21383 = sin(r21382);
        double r21384 = cbrt(r21383);
        double r21385 = r21384 * r21384;
        double r21386 = r21379 * r21385;
        double r21387 = r21386 * r21384;
        double r21388 = r21375 * r21387;
        double r21389 = -4.960559738300167e-06;
        bool r21390 = r21388 <= r21389;
        double r21391 = cos(r21380);
        double r21392 = cos(r21376);
        double r21393 = r21391 * r21392;
        double r21394 = sin(r21380);
        double r21395 = sin(r21376);
        double r21396 = r21394 * r21395;
        double r21397 = r21393 - r21396;
        double r21398 = r21397 - r21391;
        double r21399 = 0.014495001146888918;
        bool r21400 = r21388 <= r21399;
        double r21401 = r21379 * r21383;
        double r21402 = r21375 * r21401;
        double r21403 = fma(r21395, r21394, r21391);
        double r21404 = r21393 - r21403;
        double r21405 = r21400 ? r21402 : r21404;
        double r21406 = r21390 ? r21398 : r21405;
        return r21406;
}

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 r21407, r21408, r21409, r21410, r21411, r21412;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r21407);
        mpfr_init(r21408);
        mpfr_init(r21409);
        mpfr_init(r21410);
        mpfr_init(r21411);
        mpfr_init(r21412);
}

double f_im(double x, double eps) {
        mpfr_set_d(r21407, x, MPFR_RNDN);
        mpfr_set_d(r21408, eps, MPFR_RNDN);
        mpfr_add(r21409, r21407, r21408, MPFR_RNDN);
        mpfr_cos(r21410, r21409, MPFR_RNDN);
        mpfr_cos(r21411, r21407, MPFR_RNDN);
        mpfr_sub(r21412, r21410, r21411, MPFR_RNDN);
        return mpfr_get_d(r21412, MPFR_RNDN);
}

static mpfr_t r21413, r21414, r21415, r21416, r21417, r21418, r21419, r21420, r21421, r21422, r21423, r21424, r21425, r21426, r21427, r21428, r21429, r21430, r21431, r21432, r21433, r21434, r21435, r21436, r21437, r21438, r21439, r21440, r21441, r21442, r21443, r21444;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r21413, "-2", 10, MPFR_RNDN);
        mpfr_init(r21414);
        mpfr_init_set_str(r21415, "2", 10, MPFR_RNDN);
        mpfr_init(r21416);
        mpfr_init(r21417);
        mpfr_init(r21418);
        mpfr_init(r21419);
        mpfr_init(r21420);
        mpfr_init(r21421);
        mpfr_init(r21422);
        mpfr_init(r21423);
        mpfr_init(r21424);
        mpfr_init(r21425);
        mpfr_init(r21426);
        mpfr_init_set_str(r21427, "-4.960559738300167e-06", 10, MPFR_RNDN);
        mpfr_init(r21428);
        mpfr_init(r21429);
        mpfr_init(r21430);
        mpfr_init(r21431);
        mpfr_init(r21432);
        mpfr_init(r21433);
        mpfr_init(r21434);
        mpfr_init(r21435);
        mpfr_init(r21436);
        mpfr_init_set_str(r21437, "0.014495001146888918", 10, MPFR_RNDN);
        mpfr_init(r21438);
        mpfr_init(r21439);
        mpfr_init(r21440);
        mpfr_init(r21441);
        mpfr_init(r21442);
        mpfr_init(r21443);
        mpfr_init(r21444);
}

double f_fm(double x, double eps) {
        ;
        mpfr_set_d(r21414, eps, MPFR_RNDN);
        ;
        mpfr_div(r21416, r21414, r21415, MPFR_RNDN);
        mpfr_sin(r21417, r21416, MPFR_RNDN);
        mpfr_set_d(r21418, x, MPFR_RNDN);
        mpfr_fma(r21419, r21415, r21418, r21414, MPFR_RNDN);
        mpfr_div(r21420, r21419, r21415, MPFR_RNDN);
        mpfr_sin(r21421, r21420, MPFR_RNDN);
        mpfr_cbrt(r21422, r21421, MPFR_RNDN);
        mpfr_mul(r21423, r21422, r21422, MPFR_RNDN);
        mpfr_mul(r21424, r21417, r21423, MPFR_RNDN);
        mpfr_mul(r21425, r21424, r21422, MPFR_RNDN);
        mpfr_mul(r21426, r21413, r21425, MPFR_RNDN);
        ;
        mpfr_set_si(r21428, mpfr_cmp(r21426, r21427) <= 0, MPFR_RNDN);
        mpfr_cos(r21429, r21418, MPFR_RNDN);
        mpfr_cos(r21430, r21414, MPFR_RNDN);
        mpfr_mul(r21431, r21429, r21430, MPFR_RNDN);
        mpfr_sin(r21432, r21418, MPFR_RNDN);
        mpfr_sin(r21433, r21414, MPFR_RNDN);
        mpfr_mul(r21434, r21432, r21433, MPFR_RNDN);
        mpfr_sub(r21435, r21431, r21434, MPFR_RNDN);
        mpfr_sub(r21436, r21435, r21429, MPFR_RNDN);
        ;
        mpfr_set_si(r21438, mpfr_cmp(r21426, r21437) <= 0, MPFR_RNDN);
        mpfr_mul(r21439, r21417, r21421, MPFR_RNDN);
        mpfr_mul(r21440, r21413, r21439, MPFR_RNDN);
        mpfr_fma(r21441, r21433, r21432, r21429, MPFR_RNDN);
        mpfr_sub(r21442, r21431, r21441, MPFR_RNDN);
        if (mpfr_get_si(r21438, MPFR_RNDN)) { mpfr_set(r21443, r21440, MPFR_RNDN); } else { mpfr_set(r21443, r21442, MPFR_RNDN); };
        if (mpfr_get_si(r21428, MPFR_RNDN)) { mpfr_set(r21444, r21436, MPFR_RNDN); } else { mpfr_set(r21444, r21443, MPFR_RNDN); };
        return mpfr_get_d(r21444, MPFR_RNDN);
}

static mpfr_t r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459, r21460, r21461, r21462, r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r21445, "-2", 10, MPFR_RNDN);
        mpfr_init(r21446);
        mpfr_init_set_str(r21447, "2", 10, MPFR_RNDN);
        mpfr_init(r21448);
        mpfr_init(r21449);
        mpfr_init(r21450);
        mpfr_init(r21451);
        mpfr_init(r21452);
        mpfr_init(r21453);
        mpfr_init(r21454);
        mpfr_init(r21455);
        mpfr_init(r21456);
        mpfr_init(r21457);
        mpfr_init(r21458);
        mpfr_init_set_str(r21459, "-4.960559738300167e-06", 10, MPFR_RNDN);
        mpfr_init(r21460);
        mpfr_init(r21461);
        mpfr_init(r21462);
        mpfr_init(r21463);
        mpfr_init(r21464);
        mpfr_init(r21465);
        mpfr_init(r21466);
        mpfr_init(r21467);
        mpfr_init(r21468);
        mpfr_init_set_str(r21469, "0.014495001146888918", 10, MPFR_RNDN);
        mpfr_init(r21470);
        mpfr_init(r21471);
        mpfr_init(r21472);
        mpfr_init(r21473);
        mpfr_init(r21474);
        mpfr_init(r21475);
        mpfr_init(r21476);
}

double f_dm(double x, double eps) {
        ;
        mpfr_set_d(r21446, eps, MPFR_RNDN);
        ;
        mpfr_div(r21448, r21446, r21447, MPFR_RNDN);
        mpfr_sin(r21449, r21448, MPFR_RNDN);
        mpfr_set_d(r21450, x, MPFR_RNDN);
        mpfr_fma(r21451, r21447, r21450, r21446, MPFR_RNDN);
        mpfr_div(r21452, r21451, r21447, MPFR_RNDN);
        mpfr_sin(r21453, r21452, MPFR_RNDN);
        mpfr_cbrt(r21454, r21453, MPFR_RNDN);
        mpfr_mul(r21455, r21454, r21454, MPFR_RNDN);
        mpfr_mul(r21456, r21449, r21455, MPFR_RNDN);
        mpfr_mul(r21457, r21456, r21454, MPFR_RNDN);
        mpfr_mul(r21458, r21445, r21457, MPFR_RNDN);
        ;
        mpfr_set_si(r21460, mpfr_cmp(r21458, r21459) <= 0, MPFR_RNDN);
        mpfr_cos(r21461, r21450, MPFR_RNDN);
        mpfr_cos(r21462, r21446, MPFR_RNDN);
        mpfr_mul(r21463, r21461, r21462, MPFR_RNDN);
        mpfr_sin(r21464, r21450, MPFR_RNDN);
        mpfr_sin(r21465, r21446, MPFR_RNDN);
        mpfr_mul(r21466, r21464, r21465, MPFR_RNDN);
        mpfr_sub(r21467, r21463, r21466, MPFR_RNDN);
        mpfr_sub(r21468, r21467, r21461, MPFR_RNDN);
        ;
        mpfr_set_si(r21470, mpfr_cmp(r21458, r21469) <= 0, MPFR_RNDN);
        mpfr_mul(r21471, r21449, r21453, MPFR_RNDN);
        mpfr_mul(r21472, r21445, r21471, MPFR_RNDN);
        mpfr_fma(r21473, r21465, r21464, r21461, MPFR_RNDN);
        mpfr_sub(r21474, r21463, r21473, MPFR_RNDN);
        if (mpfr_get_si(r21470, MPFR_RNDN)) { mpfr_set(r21475, r21472, MPFR_RNDN); } else { mpfr_set(r21475, r21474, MPFR_RNDN); };
        if (mpfr_get_si(r21460, MPFR_RNDN)) { mpfr_set(r21476, r21468, MPFR_RNDN); } else { mpfr_set(r21476, r21475, MPFR_RNDN); };
        return mpfr_get_d(r21476, MPFR_RNDN);
}