#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 r27405 = x;
        float r27406 = eps;
        float r27407 = r27405 + r27406;
        float r27408 = sin(r27407);
        float r27409 = sin(r27405);
        float r27410 = r27408 - r27409;
        return r27410;
}

double f_id(double x, double eps) {
        double r27411 = x;
        double r27412 = eps;
        double r27413 = r27411 + r27412;
        double r27414 = sin(r27413);
        double r27415 = sin(r27411);
        double r27416 = r27414 - r27415;
        return r27416;
}


double f_of(float x, float eps) {
        float r27417 = eps;
        float r27418 = -1.223641831715565e-08;
        bool r27419 = r27417 <= r27418;
        float r27420 = 6.378020267029492e-12;
        bool r27421 = r27417 <= r27420;
        float r27422 = !r27421;
        bool r27423 = r27419 || r27422;
        float r27424 = sin(r27417);
        float r27425 = x;
        float r27426 = cos(r27425);
        float r27427 = r27424 * r27426;
        float r27428 = cos(r27417);
        float r27429 = sin(r27425);
        float r27430 = r27428 * r27429;
        float r27431 = r27427 + r27430;
        float r27432 = r27431 - r27429;
        float r27433 = 2;
        float r27434 = r27417 + r27425;
        float r27435 = r27425 + r27434;
        float r27436 = r27435 / r27433;
        float r27437 = cos(r27436);
        float r27438 = r27417 / r27433;
        float r27439 = sin(r27438);
        float r27440 = r27437 * r27439;
        float r27441 = r27433 * r27440;
        float r27442 = r27423 ? r27432 : r27441;
        return r27442;
}

double f_od(double x, double eps) {
        double r27443 = eps;
        double r27444 = -1.223641831715565e-08;
        bool r27445 = r27443 <= r27444;
        double r27446 = 6.378020267029492e-12;
        bool r27447 = r27443 <= r27446;
        double r27448 = !r27447;
        bool r27449 = r27445 || r27448;
        double r27450 = sin(r27443);
        double r27451 = x;
        double r27452 = cos(r27451);
        double r27453 = r27450 * r27452;
        double r27454 = cos(r27443);
        double r27455 = sin(r27451);
        double r27456 = r27454 * r27455;
        double r27457 = r27453 + r27456;
        double r27458 = r27457 - r27455;
        double r27459 = 2;
        double r27460 = r27443 + r27451;
        double r27461 = r27451 + r27460;
        double r27462 = r27461 / r27459;
        double r27463 = cos(r27462);
        double r27464 = r27443 / r27459;
        double r27465 = sin(r27464);
        double r27466 = r27463 * r27465;
        double r27467 = r27459 * r27466;
        double r27468 = r27449 ? r27458 : r27467;
        return r27468;
}

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 r27469, r27470, r27471, r27472, r27473, r27474;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27469);
        mpfr_init(r27470);
        mpfr_init(r27471);
        mpfr_init(r27472);
        mpfr_init(r27473);
        mpfr_init(r27474);
}

double f_im(double x, double eps) {
        mpfr_set_d(r27469, x, MPFR_RNDN);
        mpfr_set_d(r27470, eps, MPFR_RNDN);
        mpfr_add(r27471, r27469, r27470, MPFR_RNDN);
        mpfr_sin(r27472, r27471, MPFR_RNDN);
        mpfr_sin(r27473, r27469, MPFR_RNDN);
        mpfr_sub(r27474, r27472, r27473, MPFR_RNDN);
        return mpfr_get_d(r27474, MPFR_RNDN);
}

static mpfr_t r27475, r27476, r27477, r27478, r27479, r27480, r27481, r27482, r27483, r27484, r27485, r27486, r27487, r27488, r27489, r27490, r27491, r27492, r27493, r27494, r27495, r27496, r27497, r27498, r27499, r27500;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27475);
        mpfr_init_set_str(r27476, "-1.223641831715565e-08", 10, MPFR_RNDN);
        mpfr_init(r27477);
        mpfr_init_set_str(r27478, "6.378020267029492e-12", 10, MPFR_RNDN);
        mpfr_init(r27479);
        mpfr_init(r27480);
        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, "2", 10, MPFR_RNDN);
        mpfr_init(r27492);
        mpfr_init(r27493);
        mpfr_init(r27494);
        mpfr_init(r27495);
        mpfr_init(r27496);
        mpfr_init(r27497);
        mpfr_init(r27498);
        mpfr_init(r27499);
        mpfr_init(r27500);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r27475, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r27477, mpfr_cmp(r27475, r27476) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27479, mpfr_cmp(r27475, r27478) <= 0, MPFR_RNDN);
        mpfr_set_si(r27480, !mpfr_get_si(r27479, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27481, mpfr_get_si(r27477, MPFR_RNDN) || mpfr_get_si(r27480, MPFR_RNDN), MPFR_RNDN);
        mpfr_sin(r27482, r27475, MPFR_RNDN);
        mpfr_set_d(r27483, x, MPFR_RNDN);
        mpfr_cos(r27484, r27483, MPFR_RNDN);
        mpfr_mul(r27485, r27482, r27484, MPFR_RNDN);
        mpfr_cos(r27486, r27475, MPFR_RNDN);
        mpfr_sin(r27487, r27483, MPFR_RNDN);
        mpfr_mul(r27488, r27486, r27487, MPFR_RNDN);
        mpfr_add(r27489, r27485, r27488, MPFR_RNDN);
        mpfr_sub(r27490, r27489, r27487, MPFR_RNDN);
        ;
        mpfr_add(r27492, r27475, r27483, MPFR_RNDN);
        mpfr_add(r27493, r27483, r27492, MPFR_RNDN);
        mpfr_div(r27494, r27493, r27491, MPFR_RNDN);
        mpfr_cos(r27495, r27494, MPFR_RNDN);
        mpfr_div(r27496, r27475, r27491, MPFR_RNDN);
        mpfr_sin(r27497, r27496, MPFR_RNDN);
        mpfr_mul(r27498, r27495, r27497, MPFR_RNDN);
        mpfr_mul(r27499, r27491, r27498, MPFR_RNDN);
        if (mpfr_get_si(r27481, MPFR_RNDN)) { mpfr_set(r27500, r27490, MPFR_RNDN); } else { mpfr_set(r27500, r27499, MPFR_RNDN); };
        return mpfr_get_d(r27500, MPFR_RNDN);
}

static mpfr_t r27501, r27502, r27503, r27504, r27505, r27506, r27507, r27508, r27509, r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517, r27518, r27519, r27520, r27521, r27522, r27523, r27524, r27525, r27526;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27501);
        mpfr_init_set_str(r27502, "-1.223641831715565e-08", 10, MPFR_RNDN);
        mpfr_init(r27503);
        mpfr_init_set_str(r27504, "6.378020267029492e-12", 10, MPFR_RNDN);
        mpfr_init(r27505);
        mpfr_init(r27506);
        mpfr_init(r27507);
        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_set_str(r27517, "2", 10, MPFR_RNDN);
        mpfr_init(r27518);
        mpfr_init(r27519);
        mpfr_init(r27520);
        mpfr_init(r27521);
        mpfr_init(r27522);
        mpfr_init(r27523);
        mpfr_init(r27524);
        mpfr_init(r27525);
        mpfr_init(r27526);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r27501, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r27503, mpfr_cmp(r27501, r27502) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r27505, mpfr_cmp(r27501, r27504) <= 0, MPFR_RNDN);
        mpfr_set_si(r27506, !mpfr_get_si(r27505, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r27507, mpfr_get_si(r27503, MPFR_RNDN) || mpfr_get_si(r27506, MPFR_RNDN), MPFR_RNDN);
        mpfr_sin(r27508, r27501, MPFR_RNDN);
        mpfr_set_d(r27509, x, MPFR_RNDN);
        mpfr_cos(r27510, r27509, MPFR_RNDN);
        mpfr_mul(r27511, r27508, r27510, MPFR_RNDN);
        mpfr_cos(r27512, r27501, MPFR_RNDN);
        mpfr_sin(r27513, r27509, MPFR_RNDN);
        mpfr_mul(r27514, r27512, r27513, MPFR_RNDN);
        mpfr_add(r27515, r27511, r27514, MPFR_RNDN);
        mpfr_sub(r27516, r27515, r27513, MPFR_RNDN);
        ;
        mpfr_add(r27518, r27501, r27509, MPFR_RNDN);
        mpfr_add(r27519, r27509, r27518, MPFR_RNDN);
        mpfr_div(r27520, r27519, r27517, MPFR_RNDN);
        mpfr_cos(r27521, r27520, MPFR_RNDN);
        mpfr_div(r27522, r27501, r27517, MPFR_RNDN);
        mpfr_sin(r27523, r27522, MPFR_RNDN);
        mpfr_mul(r27524, r27521, r27523, MPFR_RNDN);
        mpfr_mul(r27525, r27517, r27524, MPFR_RNDN);
        if (mpfr_get_si(r27507, MPFR_RNDN)) { mpfr_set(r27526, r27516, MPFR_RNDN); } else { mpfr_set(r27526, r27525, MPFR_RNDN); };
        return mpfr_get_d(r27526, MPFR_RNDN);
}

