#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 r27545 = x;
        float r27546 = eps;
        float r27547 = r27545 + r27546;
        float r27548 = sin(r27547);
        float r27549 = sin(r27545);
        float r27550 = r27548 - r27549;
        return r27550;
}

double f_id(double x, double eps) {
        double r27551 = x;
        double r27552 = eps;
        double r27553 = r27551 + r27552;
        double r27554 = sin(r27553);
        double r27555 = sin(r27551);
        double r27556 = r27554 - r27555;
        return r27556;
}


double f_of(float x, float eps) {
        float r27557 = 2;
        float r27558 = eps;
        float r27559 = r27558 / r27557;
        float r27560 = sin(r27559);
        float r27561 = x;
        float r27562 = r27558 + r27561;
        float r27563 = r27561 + r27562;
        float r27564 = r27563 / r27557;
        float r27565 = cos(r27564);
        float r27566 = r27560 * r27565;
        float r27567 = cbrt(r27566);
        float r27568 = r27567 * r27567;
        float r27569 = r27568 * r27567;
        float r27570 = r27557 * r27569;
        float r27571 = -3.2534648250179595e-17;
        bool r27572 = r27570 <= r27571;
        float r27573 = sin(r27561);
        float r27574 = cos(r27558);
        float r27575 = r27573 * r27574;
        float r27576 = cos(r27561);
        float r27577 = sin(r27558);
        float r27578 = r27576 * r27577;
        float r27579 = r27575 + r27578;
        float r27580 = r27579 - r27573;
        float r27581 = 3.626326386292518e-07;
        bool r27582 = r27570 <= r27581;
        float r27583 = r27557 * r27566;
        float r27584 = r27578 - r27573;
        float r27585 = r27575 + r27584;
        float r27586 = r27582 ? r27583 : r27585;
        float r27587 = r27572 ? r27580 : r27586;
        return r27587;
}

double f_od(double x, double eps) {
        double r27588 = 2;
        double r27589 = eps;
        double r27590 = r27589 / r27588;
        double r27591 = sin(r27590);
        double r27592 = x;
        double r27593 = r27589 + r27592;
        double r27594 = r27592 + r27593;
        double r27595 = r27594 / r27588;
        double r27596 = cos(r27595);
        double r27597 = r27591 * r27596;
        double r27598 = cbrt(r27597);
        double r27599 = r27598 * r27598;
        double r27600 = r27599 * r27598;
        double r27601 = r27588 * r27600;
        double r27602 = -3.2534648250179595e-17;
        bool r27603 = r27601 <= r27602;
        double r27604 = sin(r27592);
        double r27605 = cos(r27589);
        double r27606 = r27604 * r27605;
        double r27607 = cos(r27592);
        double r27608 = sin(r27589);
        double r27609 = r27607 * r27608;
        double r27610 = r27606 + r27609;
        double r27611 = r27610 - r27604;
        double r27612 = 3.626326386292518e-07;
        bool r27613 = r27601 <= r27612;
        double r27614 = r27588 * r27597;
        double r27615 = r27609 - r27604;
        double r27616 = r27606 + r27615;
        double r27617 = r27613 ? r27614 : r27616;
        double r27618 = r27603 ? r27611 : r27617;
        return r27618;
}

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 r27619, r27620, r27621, r27622, r27623, r27624;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27619);
        mpfr_init(r27620);
        mpfr_init(r27621);
        mpfr_init(r27622);
        mpfr_init(r27623);
        mpfr_init(r27624);
}

double f_im(double x, double eps) {
        mpfr_set_d(r27619, x, MPFR_RNDN);
        mpfr_set_d(r27620, eps, MPFR_RNDN);
        mpfr_add(r27621, r27619, r27620, MPFR_RNDN);
        mpfr_sin(r27622, r27621, MPFR_RNDN);
        mpfr_sin(r27623, r27619, MPFR_RNDN);
        mpfr_sub(r27624, r27622, r27623, MPFR_RNDN);
        return mpfr_get_d(r27624, MPFR_RNDN);
}

static mpfr_t r27625, r27626, r27627, r27628, r27629, r27630, r27631, r27632, r27633, r27634, r27635, r27636, r27637, r27638, r27639, r27640, r27641, r27642, r27643, r27644, r27645, r27646, r27647, r27648, r27649, r27650, r27651, r27652, r27653, r27654, r27655;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r27625, "2", 10, MPFR_RNDN);
        mpfr_init(r27626);
        mpfr_init(r27627);
        mpfr_init(r27628);
        mpfr_init(r27629);
        mpfr_init(r27630);
        mpfr_init(r27631);
        mpfr_init(r27632);
        mpfr_init(r27633);
        mpfr_init(r27634);
        mpfr_init(r27635);
        mpfr_init(r27636);
        mpfr_init(r27637);
        mpfr_init(r27638);
        mpfr_init_set_str(r27639, "-3.2534648250179595e-17", 10, MPFR_RNDN);
        mpfr_init(r27640);
        mpfr_init(r27641);
        mpfr_init(r27642);
        mpfr_init(r27643);
        mpfr_init(r27644);
        mpfr_init(r27645);
        mpfr_init(r27646);
        mpfr_init(r27647);
        mpfr_init(r27648);
        mpfr_init_set_str(r27649, "3.626326386292518e-07", 10, MPFR_RNDN);
        mpfr_init(r27650);
        mpfr_init(r27651);
        mpfr_init(r27652);
        mpfr_init(r27653);
        mpfr_init(r27654);
        mpfr_init(r27655);
}

double f_fm(double x, double eps) {
        ;
        mpfr_set_d(r27626, eps, MPFR_RNDN);
        mpfr_div(r27627, r27626, r27625, MPFR_RNDN);
        mpfr_sin(r27628, r27627, MPFR_RNDN);
        mpfr_set_d(r27629, x, MPFR_RNDN);
        mpfr_add(r27630, r27626, r27629, MPFR_RNDN);
        mpfr_add(r27631, r27629, r27630, MPFR_RNDN);
        mpfr_div(r27632, r27631, r27625, MPFR_RNDN);
        mpfr_cos(r27633, r27632, MPFR_RNDN);
        mpfr_mul(r27634, r27628, r27633, MPFR_RNDN);
        mpfr_cbrt(r27635, r27634, MPFR_RNDN);
        mpfr_mul(r27636, r27635, r27635, MPFR_RNDN);
        mpfr_mul(r27637, r27636, r27635, MPFR_RNDN);
        mpfr_mul(r27638, r27625, r27637, MPFR_RNDN);
        ;
        mpfr_set_si(r27640, mpfr_cmp(r27638, r27639) <= 0, MPFR_RNDN);
        mpfr_sin(r27641, r27629, MPFR_RNDN);
        mpfr_cos(r27642, r27626, MPFR_RNDN);
        mpfr_mul(r27643, r27641, r27642, MPFR_RNDN);
        mpfr_cos(r27644, r27629, MPFR_RNDN);
        mpfr_sin(r27645, r27626, MPFR_RNDN);
        mpfr_mul(r27646, r27644, r27645, MPFR_RNDN);
        mpfr_add(r27647, r27643, r27646, MPFR_RNDN);
        mpfr_sub(r27648, r27647, r27641, MPFR_RNDN);
        ;
        mpfr_set_si(r27650, mpfr_cmp(r27638, r27649) <= 0, MPFR_RNDN);
        mpfr_mul(r27651, r27625, r27634, MPFR_RNDN);
        mpfr_sub(r27652, r27646, r27641, MPFR_RNDN);
        mpfr_add(r27653, r27643, r27652, MPFR_RNDN);
        if (mpfr_get_si(r27650, MPFR_RNDN)) { mpfr_set(r27654, r27651, MPFR_RNDN); } else { mpfr_set(r27654, r27653, MPFR_RNDN); };
        if (mpfr_get_si(r27640, MPFR_RNDN)) { mpfr_set(r27655, r27648, MPFR_RNDN); } else { mpfr_set(r27655, r27654, MPFR_RNDN); };
        return mpfr_get_d(r27655, MPFR_RNDN);
}

static mpfr_t r27656, r27657, r27658, r27659, r27660, r27661, r27662, r27663, r27664, r27665, r27666, r27667, r27668, r27669, r27670, r27671, r27672, r27673, r27674, r27675, r27676, r27677, r27678, r27679, r27680, r27681, r27682, r27683, r27684, r27685, r27686;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r27656, "2", 10, MPFR_RNDN);
        mpfr_init(r27657);
        mpfr_init(r27658);
        mpfr_init(r27659);
        mpfr_init(r27660);
        mpfr_init(r27661);
        mpfr_init(r27662);
        mpfr_init(r27663);
        mpfr_init(r27664);
        mpfr_init(r27665);
        mpfr_init(r27666);
        mpfr_init(r27667);
        mpfr_init(r27668);
        mpfr_init(r27669);
        mpfr_init_set_str(r27670, "-3.2534648250179595e-17", 10, MPFR_RNDN);
        mpfr_init(r27671);
        mpfr_init(r27672);
        mpfr_init(r27673);
        mpfr_init(r27674);
        mpfr_init(r27675);
        mpfr_init(r27676);
        mpfr_init(r27677);
        mpfr_init(r27678);
        mpfr_init(r27679);
        mpfr_init_set_str(r27680, "3.626326386292518e-07", 10, MPFR_RNDN);
        mpfr_init(r27681);
        mpfr_init(r27682);
        mpfr_init(r27683);
        mpfr_init(r27684);
        mpfr_init(r27685);
        mpfr_init(r27686);
}

double f_dm(double x, double eps) {
        ;
        mpfr_set_d(r27657, eps, MPFR_RNDN);
        mpfr_div(r27658, r27657, r27656, MPFR_RNDN);
        mpfr_sin(r27659, r27658, MPFR_RNDN);
        mpfr_set_d(r27660, x, MPFR_RNDN);
        mpfr_add(r27661, r27657, r27660, MPFR_RNDN);
        mpfr_add(r27662, r27660, r27661, MPFR_RNDN);
        mpfr_div(r27663, r27662, r27656, MPFR_RNDN);
        mpfr_cos(r27664, r27663, MPFR_RNDN);
        mpfr_mul(r27665, r27659, r27664, MPFR_RNDN);
        mpfr_cbrt(r27666, r27665, MPFR_RNDN);
        mpfr_mul(r27667, r27666, r27666, MPFR_RNDN);
        mpfr_mul(r27668, r27667, r27666, MPFR_RNDN);
        mpfr_mul(r27669, r27656, r27668, MPFR_RNDN);
        ;
        mpfr_set_si(r27671, mpfr_cmp(r27669, r27670) <= 0, MPFR_RNDN);
        mpfr_sin(r27672, r27660, MPFR_RNDN);
        mpfr_cos(r27673, r27657, MPFR_RNDN);
        mpfr_mul(r27674, r27672, r27673, MPFR_RNDN);
        mpfr_cos(r27675, r27660, MPFR_RNDN);
        mpfr_sin(r27676, r27657, MPFR_RNDN);
        mpfr_mul(r27677, r27675, r27676, MPFR_RNDN);
        mpfr_add(r27678, r27674, r27677, MPFR_RNDN);
        mpfr_sub(r27679, r27678, r27672, MPFR_RNDN);
        ;
        mpfr_set_si(r27681, mpfr_cmp(r27669, r27680) <= 0, MPFR_RNDN);
        mpfr_mul(r27682, r27656, r27665, MPFR_RNDN);
        mpfr_sub(r27683, r27677, r27672, MPFR_RNDN);
        mpfr_add(r27684, r27674, r27683, MPFR_RNDN);
        if (mpfr_get_si(r27681, MPFR_RNDN)) { mpfr_set(r27685, r27682, MPFR_RNDN); } else { mpfr_set(r27685, r27684, MPFR_RNDN); };
        if (mpfr_get_si(r27671, MPFR_RNDN)) { mpfr_set(r27686, r27679, MPFR_RNDN); } else { mpfr_set(r27686, r27685, MPFR_RNDN); };
        return mpfr_get_d(r27686, MPFR_RNDN);
}

