#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 r23425 = x;
        float r23426 = eps;
        float r23427 = r23425 + r23426;
        float r23428 = cos(r23427);
        float r23429 = cos(r23425);
        float r23430 = r23428 - r23429;
        return r23430;
}

double f_id(double x, double eps) {
        double r23431 = x;
        double r23432 = eps;
        double r23433 = r23431 + r23432;
        double r23434 = cos(r23433);
        double r23435 = cos(r23431);
        double r23436 = r23434 - r23435;
        return r23436;
}


double f_of(float x, float eps) {
        float r23437 = x;
        float r23438 = cos(r23437);
        float r23439 = eps;
        float r23440 = cos(r23439);
        float r23441 = r23438 * r23440;
        float r23442 = sin(r23437);
        float r23443 = sin(r23439);
        float r23444 = r23442 * r23443;
        float r23445 = r23441 - r23444;
        float r23446 = r23445 - r23438;
        float r23447 = -0.01443769579639198;
        bool r23448 = r23446 <= r23447;
        float r23449 = 0.004408463031335852;
        bool r23450 = r23446 <= r23449;
        float r23451 = -2;
        float r23452 = 2;
        float r23453 = r23439 / r23452;
        float r23454 = sin(r23453);
        float r23455 = r23439 + r23437;
        float r23456 = r23437 + r23455;
        float r23457 = r23456 / r23452;
        float r23458 = sin(r23457);
        float r23459 = r23454 * r23458;
        float r23460 = r23451 * r23459;
        float r23461 = r23450 ? r23460 : r23446;
        float r23462 = r23448 ? r23446 : r23461;
        return r23462;
}

double f_od(double x, double eps) {
        double r23463 = x;
        double r23464 = cos(r23463);
        double r23465 = eps;
        double r23466 = cos(r23465);
        double r23467 = r23464 * r23466;
        double r23468 = sin(r23463);
        double r23469 = sin(r23465);
        double r23470 = r23468 * r23469;
        double r23471 = r23467 - r23470;
        double r23472 = r23471 - r23464;
        double r23473 = -0.01443769579639198;
        bool r23474 = r23472 <= r23473;
        double r23475 = 0.004408463031335852;
        bool r23476 = r23472 <= r23475;
        double r23477 = -2;
        double r23478 = 2;
        double r23479 = r23465 / r23478;
        double r23480 = sin(r23479);
        double r23481 = r23465 + r23463;
        double r23482 = r23463 + r23481;
        double r23483 = r23482 / r23478;
        double r23484 = sin(r23483);
        double r23485 = r23480 * r23484;
        double r23486 = r23477 * r23485;
        double r23487 = r23476 ? r23486 : r23472;
        double r23488 = r23474 ? r23472 : r23487;
        return r23488;
}

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 r23489, r23490, r23491, r23492, r23493, r23494;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23489);
        mpfr_init(r23490);
        mpfr_init(r23491);
        mpfr_init(r23492);
        mpfr_init(r23493);
        mpfr_init(r23494);
}

double f_im(double x, double eps) {
        mpfr_set_d(r23489, x, MPFR_RNDN);
        mpfr_set_d(r23490, eps, MPFR_RNDN);
        mpfr_add(r23491, r23489, r23490, MPFR_RNDN);
        mpfr_cos(r23492, r23491, MPFR_RNDN);
        mpfr_cos(r23493, r23489, MPFR_RNDN);
        mpfr_sub(r23494, r23492, r23493, MPFR_RNDN);
        return mpfr_get_d(r23494, MPFR_RNDN);
}

static mpfr_t r23495, r23496, r23497, r23498, r23499, r23500, r23501, r23502, r23503, r23504, r23505, r23506, r23507, r23508, r23509, r23510, r23511, r23512, r23513, r23514, r23515, r23516, r23517, r23518, r23519, r23520;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23495);
        mpfr_init(r23496);
        mpfr_init(r23497);
        mpfr_init(r23498);
        mpfr_init(r23499);
        mpfr_init(r23500);
        mpfr_init(r23501);
        mpfr_init(r23502);
        mpfr_init(r23503);
        mpfr_init(r23504);
        mpfr_init_set_str(r23505, "-0.01443769579639198", 10, MPFR_RNDN);
        mpfr_init(r23506);
        mpfr_init_set_str(r23507, "0.004408463031335852", 10, MPFR_RNDN);
        mpfr_init(r23508);
        mpfr_init_set_str(r23509, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r23510, "2", 10, MPFR_RNDN);
        mpfr_init(r23511);
        mpfr_init(r23512);
        mpfr_init(r23513);
        mpfr_init(r23514);
        mpfr_init(r23515);
        mpfr_init(r23516);
        mpfr_init(r23517);
        mpfr_init(r23518);
        mpfr_init(r23519);
        mpfr_init(r23520);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r23495, x, MPFR_RNDN);
        mpfr_cos(r23496, r23495, MPFR_RNDN);
        mpfr_set_d(r23497, eps, MPFR_RNDN);
        mpfr_cos(r23498, r23497, MPFR_RNDN);
        mpfr_mul(r23499, r23496, r23498, MPFR_RNDN);
        mpfr_sin(r23500, r23495, MPFR_RNDN);
        mpfr_sin(r23501, r23497, MPFR_RNDN);
        mpfr_mul(r23502, r23500, r23501, MPFR_RNDN);
        mpfr_sub(r23503, r23499, r23502, MPFR_RNDN);
        mpfr_sub(r23504, r23503, r23496, MPFR_RNDN);
        ;
        mpfr_set_si(r23506, mpfr_cmp(r23504, r23505) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r23508, mpfr_cmp(r23504, r23507) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r23511, r23497, r23510, MPFR_RNDN);
        mpfr_sin(r23512, r23511, MPFR_RNDN);
        mpfr_add(r23513, r23497, r23495, MPFR_RNDN);
        mpfr_add(r23514, r23495, r23513, MPFR_RNDN);
        mpfr_div(r23515, r23514, r23510, MPFR_RNDN);
        mpfr_sin(r23516, r23515, MPFR_RNDN);
        mpfr_mul(r23517, r23512, r23516, MPFR_RNDN);
        mpfr_mul(r23518, r23509, r23517, MPFR_RNDN);
        if (mpfr_get_si(r23508, MPFR_RNDN)) { mpfr_set(r23519, r23518, MPFR_RNDN); } else { mpfr_set(r23519, r23504, MPFR_RNDN); };
        if (mpfr_get_si(r23506, MPFR_RNDN)) { mpfr_set(r23520, r23504, MPFR_RNDN); } else { mpfr_set(r23520, r23519, MPFR_RNDN); };
        return mpfr_get_d(r23520, MPFR_RNDN);
}

static mpfr_t r23521, r23522, r23523, r23524, r23525, r23526, r23527, r23528, r23529, r23530, r23531, r23532, r23533, r23534, r23535, r23536, r23537, r23538, r23539, r23540, r23541, r23542, r23543, r23544, r23545, r23546;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r23521);
        mpfr_init(r23522);
        mpfr_init(r23523);
        mpfr_init(r23524);
        mpfr_init(r23525);
        mpfr_init(r23526);
        mpfr_init(r23527);
        mpfr_init(r23528);
        mpfr_init(r23529);
        mpfr_init(r23530);
        mpfr_init_set_str(r23531, "-0.01443769579639198", 10, MPFR_RNDN);
        mpfr_init(r23532);
        mpfr_init_set_str(r23533, "0.004408463031335852", 10, MPFR_RNDN);
        mpfr_init(r23534);
        mpfr_init_set_str(r23535, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r23536, "2", 10, MPFR_RNDN);
        mpfr_init(r23537);
        mpfr_init(r23538);
        mpfr_init(r23539);
        mpfr_init(r23540);
        mpfr_init(r23541);
        mpfr_init(r23542);
        mpfr_init(r23543);
        mpfr_init(r23544);
        mpfr_init(r23545);
        mpfr_init(r23546);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r23521, x, MPFR_RNDN);
        mpfr_cos(r23522, r23521, MPFR_RNDN);
        mpfr_set_d(r23523, eps, MPFR_RNDN);
        mpfr_cos(r23524, r23523, MPFR_RNDN);
        mpfr_mul(r23525, r23522, r23524, MPFR_RNDN);
        mpfr_sin(r23526, r23521, MPFR_RNDN);
        mpfr_sin(r23527, r23523, MPFR_RNDN);
        mpfr_mul(r23528, r23526, r23527, MPFR_RNDN);
        mpfr_sub(r23529, r23525, r23528, MPFR_RNDN);
        mpfr_sub(r23530, r23529, r23522, MPFR_RNDN);
        ;
        mpfr_set_si(r23532, mpfr_cmp(r23530, r23531) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r23534, mpfr_cmp(r23530, r23533) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r23537, r23523, r23536, MPFR_RNDN);
        mpfr_sin(r23538, r23537, MPFR_RNDN);
        mpfr_add(r23539, r23523, r23521, MPFR_RNDN);
        mpfr_add(r23540, r23521, r23539, MPFR_RNDN);
        mpfr_div(r23541, r23540, r23536, MPFR_RNDN);
        mpfr_sin(r23542, r23541, MPFR_RNDN);
        mpfr_mul(r23543, r23538, r23542, MPFR_RNDN);
        mpfr_mul(r23544, r23535, r23543, MPFR_RNDN);
        if (mpfr_get_si(r23534, MPFR_RNDN)) { mpfr_set(r23545, r23544, MPFR_RNDN); } else { mpfr_set(r23545, r23530, MPFR_RNDN); };
        if (mpfr_get_si(r23532, MPFR_RNDN)) { mpfr_set(r23546, r23530, MPFR_RNDN); } else { mpfr_set(r23546, r23545, MPFR_RNDN); };
        return mpfr_get_d(r23546, MPFR_RNDN);
}