#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 r27409 = x;
        float r27410 = eps;
        float r27411 = r27409 + r27410;
        float r27412 = sin(r27411);
        float r27413 = sin(r27409);
        float r27414 = r27412 - r27413;
        return r27414;
}

double f_id(double x, double eps) {
        double r27415 = x;
        double r27416 = eps;
        double r27417 = r27415 + r27416;
        double r27418 = sin(r27417);
        double r27419 = sin(r27415);
        double r27420 = r27418 - r27419;
        return r27420;
}


double f_of(float x, float eps) {
        float r27421 = x;
        float r27422 = sin(r27421);
        float r27423 = eps;
        float r27424 = cos(r27423);
        float r27425 = r27422 * r27424;
        float r27426 = cos(r27421);
        float r27427 = sin(r27423);
        float r27428 = r27426 * r27427;
        float r27429 = r27428 - r27422;
        float r27430 = r27425 + r27429;
        float r27431 = -0.09882177259665836;
        bool r27432 = r27430 <= r27431;
        float r27433 = 2.856382425847474e-10;
        bool r27434 = r27430 <= r27433;
        float r27435 = 2;
        float r27436 = r27423 / r27435;
        float r27437 = sin(r27436);
        float r27438 = r27423 + r27421;
        float r27439 = r27421 + r27438;
        float r27440 = r27439 / r27435;
        float r27441 = cos(r27440);
        float r27442 = r27437 * r27441;
        float r27443 = r27435 * r27442;
        float r27444 = r27425 + r27428;
        float r27445 = r27444 - r27422;
        float r27446 = r27434 ? r27443 : r27445;
        float r27447 = r27432 ? r27430 : r27446;
        return r27447;
}

double f_od(double x, double eps) {
        double r27448 = x;
        double r27449 = sin(r27448);
        double r27450 = eps;
        double r27451 = cos(r27450);
        double r27452 = r27449 * r27451;
        double r27453 = cos(r27448);
        double r27454 = sin(r27450);
        double r27455 = r27453 * r27454;
        double r27456 = r27455 - r27449;
        double r27457 = r27452 + r27456;
        double r27458 = -0.09882177259665836;
        bool r27459 = r27457 <= r27458;
        double r27460 = 2.856382425847474e-10;
        bool r27461 = r27457 <= r27460;
        double r27462 = 2;
        double r27463 = r27450 / r27462;
        double r27464 = sin(r27463);
        double r27465 = r27450 + r27448;
        double r27466 = r27448 + r27465;
        double r27467 = r27466 / r27462;
        double r27468 = cos(r27467);
        double r27469 = r27464 * r27468;
        double r27470 = r27462 * r27469;
        double r27471 = r27452 + r27455;
        double r27472 = r27471 - r27449;
        double r27473 = r27461 ? r27470 : r27472;
        double r27474 = r27459 ? r27457 : r27473;
        return r27474;
}

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 r27475, r27476, r27477, r27478, r27479, r27480;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27475);
        mpfr_init(r27476);
        mpfr_init(r27477);
        mpfr_init(r27478);
        mpfr_init(r27479);
        mpfr_init(r27480);
}

double f_im(double x, double eps) {
        mpfr_set_d(r27475, x, MPFR_RNDN);
        mpfr_set_d(r27476, eps, MPFR_RNDN);
        mpfr_add(r27477, r27475, r27476, MPFR_RNDN);
        mpfr_sin(r27478, r27477, MPFR_RNDN);
        mpfr_sin(r27479, r27475, MPFR_RNDN);
        mpfr_sub(r27480, r27478, r27479, MPFR_RNDN);
        return mpfr_get_d(r27480, MPFR_RNDN);
}

static mpfr_t r27481, r27482, r27483, r27484, r27485, r27486, r27487, r27488, r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500, r27501, r27502, r27503, r27504, r27505, r27506, r27507;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27481);
        mpfr_init(r27482);
        mpfr_init(r27483);
        mpfr_init(r27484);
        mpfr_init(r27485);
        mpfr_init(r27486);
        mpfr_init(r27487);
        mpfr_init(r27488);
        mpfr_init(r27489);
        mpfr_init(r27490);
        mpfr_init_set_str(r27491, "-0.09882177259665836", 10, MPFR_RNDN);
        mpfr_init(r27492);
        mpfr_init_set_str(r27493, "2.856382425847474e-10", 10, MPFR_RNDN);
        mpfr_init(r27494);
        mpfr_init_set_str(r27495, "2", 10, MPFR_RNDN);
        mpfr_init(r27496);
        mpfr_init(r27497);
        mpfr_init(r27498);
        mpfr_init(r27499);
        mpfr_init(r27500);
        mpfr_init(r27501);
        mpfr_init(r27502);
        mpfr_init(r27503);
        mpfr_init(r27504);
        mpfr_init(r27505);
        mpfr_init(r27506);
        mpfr_init(r27507);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r27481, x, MPFR_RNDN);
        mpfr_sin(r27482, r27481, MPFR_RNDN);
        mpfr_set_d(r27483, eps, MPFR_RNDN);
        mpfr_cos(r27484, r27483, MPFR_RNDN);
        mpfr_mul(r27485, r27482, r27484, MPFR_RNDN);
        mpfr_cos(r27486, r27481, MPFR_RNDN);
        mpfr_sin(r27487, r27483, MPFR_RNDN);
        mpfr_mul(r27488, r27486, r27487, MPFR_RNDN);
        mpfr_sub(r27489, r27488, r27482, MPFR_RNDN);
        mpfr_add(r27490, r27485, r27489, MPFR_RNDN);
        ;
        mpfr_set_si(r27492, mpfr_cmp(r27490, r27491) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27494, mpfr_cmp(r27490, r27493) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27496, r27483, r27495, MPFR_RNDN);
        mpfr_sin(r27497, r27496, MPFR_RNDN);
        mpfr_add(r27498, r27483, r27481, MPFR_RNDN);
        mpfr_add(r27499, r27481, r27498, MPFR_RNDN);
        mpfr_div(r27500, r27499, r27495, MPFR_RNDN);
        mpfr_cos(r27501, r27500, MPFR_RNDN);
        mpfr_mul(r27502, r27497, r27501, MPFR_RNDN);
        mpfr_mul(r27503, r27495, r27502, MPFR_RNDN);
        mpfr_add(r27504, r27485, r27488, MPFR_RNDN);
        mpfr_sub(r27505, r27504, r27482, MPFR_RNDN);
        if (mpfr_get_si(r27494, MPFR_RNDN)) { mpfr_set(r27506, r27503, MPFR_RNDN); } else { mpfr_set(r27506, r27505, MPFR_RNDN); };
        if (mpfr_get_si(r27492, MPFR_RNDN)) { mpfr_set(r27507, r27490, MPFR_RNDN); } else { mpfr_set(r27507, r27506, MPFR_RNDN); };
        return mpfr_get_d(r27507, MPFR_RNDN);
}

static mpfr_t r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517, r27518, r27519, r27520, r27521, r27522, r27523, r27524, r27525, r27526, r27527, r27528, r27529, r27530, r27531, r27532, r27533, r27534;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27508);
        mpfr_init(r27509);
        mpfr_init(r27510);
        mpfr_init(r27511);
        mpfr_init(r27512);
        mpfr_init(r27513);
        mpfr_init(r27514);
        mpfr_init(r27515);
        mpfr_init(r27516);
        mpfr_init(r27517);
        mpfr_init_set_str(r27518, "-0.09882177259665836", 10, MPFR_RNDN);
        mpfr_init(r27519);
        mpfr_init_set_str(r27520, "2.856382425847474e-10", 10, MPFR_RNDN);
        mpfr_init(r27521);
        mpfr_init_set_str(r27522, "2", 10, MPFR_RNDN);
        mpfr_init(r27523);
        mpfr_init(r27524);
        mpfr_init(r27525);
        mpfr_init(r27526);
        mpfr_init(r27527);
        mpfr_init(r27528);
        mpfr_init(r27529);
        mpfr_init(r27530);
        mpfr_init(r27531);
        mpfr_init(r27532);
        mpfr_init(r27533);
        mpfr_init(r27534);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r27508, x, MPFR_RNDN);
        mpfr_sin(r27509, r27508, MPFR_RNDN);
        mpfr_set_d(r27510, eps, MPFR_RNDN);
        mpfr_cos(r27511, r27510, MPFR_RNDN);
        mpfr_mul(r27512, r27509, r27511, MPFR_RNDN);
        mpfr_cos(r27513, r27508, MPFR_RNDN);
        mpfr_sin(r27514, r27510, MPFR_RNDN);
        mpfr_mul(r27515, r27513, r27514, MPFR_RNDN);
        mpfr_sub(r27516, r27515, r27509, MPFR_RNDN);
        mpfr_add(r27517, r27512, r27516, MPFR_RNDN);
        ;
        mpfr_set_si(r27519, mpfr_cmp(r27517, r27518) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27521, mpfr_cmp(r27517, r27520) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27523, r27510, r27522, MPFR_RNDN);
        mpfr_sin(r27524, r27523, MPFR_RNDN);
        mpfr_add(r27525, r27510, r27508, MPFR_RNDN);
        mpfr_add(r27526, r27508, r27525, MPFR_RNDN);
        mpfr_div(r27527, r27526, r27522, MPFR_RNDN);
        mpfr_cos(r27528, r27527, MPFR_RNDN);
        mpfr_mul(r27529, r27524, r27528, MPFR_RNDN);
        mpfr_mul(r27530, r27522, r27529, MPFR_RNDN);
        mpfr_add(r27531, r27512, r27515, MPFR_RNDN);
        mpfr_sub(r27532, r27531, r27509, MPFR_RNDN);
        if (mpfr_get_si(r27521, MPFR_RNDN)) { mpfr_set(r27533, r27530, MPFR_RNDN); } else { mpfr_set(r27533, r27532, MPFR_RNDN); };
        if (mpfr_get_si(r27519, MPFR_RNDN)) { mpfr_set(r27534, r27517, MPFR_RNDN); } else { mpfr_set(r27534, r27533, MPFR_RNDN); };
        return mpfr_get_d(r27534, MPFR_RNDN);
}

