#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 r20329 = x;
        float r20330 = eps;
        float r20331 = r20329 + r20330;
        float r20332 = cos(r20331);
        float r20333 = cos(r20329);
        float r20334 = r20332 - r20333;
        return r20334;
}

double f_id(double x, double eps) {
        double r20335 = x;
        double r20336 = eps;
        double r20337 = r20335 + r20336;
        double r20338 = cos(r20337);
        double r20339 = cos(r20335);
        double r20340 = r20338 - r20339;
        return r20340;
}


double f_of(float x, float eps) {
        float r20341 = eps;
        float r20342 = -4.937900603750566e-17;
        bool r20343 = r20341 <= r20342;
        float r20344 = x;
        float r20345 = cos(r20344);
        float r20346 = cos(r20341);
        float r20347 = r20345 * r20346;
        float r20348 = sin(r20344);
        float r20349 = sin(r20341);
        float r20350 = r20348 * r20349;
        float r20351 = r20347 - r20350;
        float r20352 = r20351 - r20345;
        float r20353 = 0.802860764679358;
        bool r20354 = r20341 <= r20353;
        float r20355 = -2;
        float r20356 = 2;
        float r20357 = r20341 / r20356;
        float r20358 = sin(r20357);
        float r20359 = r20344 + r20344;
        float r20360 = r20341 + r20359;
        float r20361 = r20360 / r20356;
        float r20362 = sin(r20361);
        float r20363 = r20358 * r20362;
        float r20364 = r20355 * r20363;
        float r20365 = r20354 ? r20364 : r20352;
        float r20366 = r20343 ? r20352 : r20365;
        return r20366;
}

double f_od(double x, double eps) {
        double r20367 = eps;
        double r20368 = -4.937900603750566e-17;
        bool r20369 = r20367 <= r20368;
        double r20370 = x;
        double r20371 = cos(r20370);
        double r20372 = cos(r20367);
        double r20373 = r20371 * r20372;
        double r20374 = sin(r20370);
        double r20375 = sin(r20367);
        double r20376 = r20374 * r20375;
        double r20377 = r20373 - r20376;
        double r20378 = r20377 - r20371;
        double r20379 = 0.802860764679358;
        bool r20380 = r20367 <= r20379;
        double r20381 = -2;
        double r20382 = 2;
        double r20383 = r20367 / r20382;
        double r20384 = sin(r20383);
        double r20385 = r20370 + r20370;
        double r20386 = r20367 + r20385;
        double r20387 = r20386 / r20382;
        double r20388 = sin(r20387);
        double r20389 = r20384 * r20388;
        double r20390 = r20381 * r20389;
        double r20391 = r20380 ? r20390 : r20378;
        double r20392 = r20369 ? r20378 : r20391;
        return r20392;
}

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 r20393, r20394, r20395, r20396, r20397, r20398;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20393);
        mpfr_init(r20394);
        mpfr_init(r20395);
        mpfr_init(r20396);
        mpfr_init(r20397);
        mpfr_init(r20398);
}

double f_im(double x, double eps) {
        mpfr_set_d(r20393, x, MPFR_RNDN);
        mpfr_set_d(r20394, eps, MPFR_RNDN);
        mpfr_add(r20395, r20393, r20394, MPFR_RNDN);
        mpfr_cos(r20396, r20395, MPFR_RNDN);
        mpfr_cos(r20397, r20393, MPFR_RNDN);
        mpfr_sub(r20398, r20396, r20397, MPFR_RNDN);
        return mpfr_get_d(r20398, MPFR_RNDN);
}

static mpfr_t r20399, r20400, r20401, r20402, r20403, r20404, r20405, r20406, r20407, r20408, r20409, r20410, r20411, r20412, r20413, r20414, r20415, r20416, r20417, r20418, r20419, r20420, r20421, r20422, r20423, r20424;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20399);
        mpfr_init_set_str(r20400, "-4.937900603750566e-17", 10, MPFR_RNDN);
        mpfr_init(r20401);
        mpfr_init(r20402);
        mpfr_init(r20403);
        mpfr_init(r20404);
        mpfr_init(r20405);
        mpfr_init(r20406);
        mpfr_init(r20407);
        mpfr_init(r20408);
        mpfr_init(r20409);
        mpfr_init(r20410);
        mpfr_init_set_str(r20411, "0.802860764679358", 10, MPFR_RNDN);
        mpfr_init(r20412);
        mpfr_init_set_str(r20413, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r20414, "2", 10, MPFR_RNDN);
        mpfr_init(r20415);
        mpfr_init(r20416);
        mpfr_init(r20417);
        mpfr_init(r20418);
        mpfr_init(r20419);
        mpfr_init(r20420);
        mpfr_init(r20421);
        mpfr_init(r20422);
        mpfr_init(r20423);
        mpfr_init(r20424);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r20399, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20401, mpfr_cmp(r20399, r20400) <= 0, MPFR_RNDN);
        mpfr_set_d(r20402, x, MPFR_RNDN);
        mpfr_cos(r20403, r20402, MPFR_RNDN);
        mpfr_cos(r20404, r20399, MPFR_RNDN);
        mpfr_mul(r20405, r20403, r20404, MPFR_RNDN);
        mpfr_sin(r20406, r20402, MPFR_RNDN);
        mpfr_sin(r20407, r20399, MPFR_RNDN);
        mpfr_mul(r20408, r20406, r20407, MPFR_RNDN);
        mpfr_sub(r20409, r20405, r20408, MPFR_RNDN);
        mpfr_sub(r20410, r20409, r20403, MPFR_RNDN);
        ;
        mpfr_set_si(r20412, mpfr_cmp(r20399, r20411) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r20415, r20399, r20414, MPFR_RNDN);
        mpfr_sin(r20416, r20415, MPFR_RNDN);
        mpfr_add(r20417, r20402, r20402, MPFR_RNDN);
        mpfr_add(r20418, r20399, r20417, MPFR_RNDN);
        mpfr_div(r20419, r20418, r20414, MPFR_RNDN);
        mpfr_sin(r20420, r20419, MPFR_RNDN);
        mpfr_mul(r20421, r20416, r20420, MPFR_RNDN);
        mpfr_mul(r20422, r20413, r20421, MPFR_RNDN);
        if (mpfr_get_si(r20412, MPFR_RNDN)) { mpfr_set(r20423, r20422, MPFR_RNDN); } else { mpfr_set(r20423, r20410, MPFR_RNDN); };
        if (mpfr_get_si(r20401, MPFR_RNDN)) { mpfr_set(r20424, r20410, MPFR_RNDN); } else { mpfr_set(r20424, r20423, MPFR_RNDN); };
        return mpfr_get_d(r20424, MPFR_RNDN);
}

static mpfr_t r20425, r20426, r20427, r20428, r20429, r20430, r20431, r20432, r20433, r20434, r20435, r20436, r20437, r20438, r20439, r20440, r20441, r20442, r20443, r20444, r20445, r20446, r20447, r20448, r20449, r20450;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20425);
        mpfr_init_set_str(r20426, "-4.937900603750566e-17", 10, MPFR_RNDN);
        mpfr_init(r20427);
        mpfr_init(r20428);
        mpfr_init(r20429);
        mpfr_init(r20430);
        mpfr_init(r20431);
        mpfr_init(r20432);
        mpfr_init(r20433);
        mpfr_init(r20434);
        mpfr_init(r20435);
        mpfr_init(r20436);
        mpfr_init_set_str(r20437, "0.802860764679358", 10, MPFR_RNDN);
        mpfr_init(r20438);
        mpfr_init_set_str(r20439, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r20440, "2", 10, MPFR_RNDN);
        mpfr_init(r20441);
        mpfr_init(r20442);
        mpfr_init(r20443);
        mpfr_init(r20444);
        mpfr_init(r20445);
        mpfr_init(r20446);
        mpfr_init(r20447);
        mpfr_init(r20448);
        mpfr_init(r20449);
        mpfr_init(r20450);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r20425, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20427, mpfr_cmp(r20425, r20426) <= 0, MPFR_RNDN);
        mpfr_set_d(r20428, x, MPFR_RNDN);
        mpfr_cos(r20429, r20428, MPFR_RNDN);
        mpfr_cos(r20430, r20425, MPFR_RNDN);
        mpfr_mul(r20431, r20429, r20430, MPFR_RNDN);
        mpfr_sin(r20432, r20428, MPFR_RNDN);
        mpfr_sin(r20433, r20425, MPFR_RNDN);
        mpfr_mul(r20434, r20432, r20433, MPFR_RNDN);
        mpfr_sub(r20435, r20431, r20434, MPFR_RNDN);
        mpfr_sub(r20436, r20435, r20429, MPFR_RNDN);
        ;
        mpfr_set_si(r20438, mpfr_cmp(r20425, r20437) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r20441, r20425, r20440, MPFR_RNDN);
        mpfr_sin(r20442, r20441, MPFR_RNDN);
        mpfr_add(r20443, r20428, r20428, MPFR_RNDN);
        mpfr_add(r20444, r20425, r20443, MPFR_RNDN);
        mpfr_div(r20445, r20444, r20440, MPFR_RNDN);
        mpfr_sin(r20446, r20445, MPFR_RNDN);
        mpfr_mul(r20447, r20442, r20446, MPFR_RNDN);
        mpfr_mul(r20448, r20439, r20447, MPFR_RNDN);
        if (mpfr_get_si(r20438, MPFR_RNDN)) { mpfr_set(r20449, r20448, MPFR_RNDN); } else { mpfr_set(r20449, r20436, MPFR_RNDN); };
        if (mpfr_get_si(r20427, MPFR_RNDN)) { mpfr_set(r20450, r20436, MPFR_RNDN); } else { mpfr_set(r20450, r20449, MPFR_RNDN); };
        return mpfr_get_d(r20450, MPFR_RNDN);
}