#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "sintan (problem 3.4.5)";

double f_if(float x) {
        float r19677 = x;
        float r19678 = sin(r19677);
        float r19679 = r19677 - r19678;
        float r19680 = tan(r19677);
        float r19681 = r19677 - r19680;
        float r19682 = r19679 / r19681;
        return r19682;
}

double f_id(double x) {
        double r19683 = x;
        double r19684 = sin(r19683);
        double r19685 = r19683 - r19684;
        double r19686 = tan(r19683);
        double r19687 = r19683 - r19686;
        double r19688 = r19685 / r19687;
        return r19688;
}


double f_of(float x) {
        float r19689 = x;
        float r19690 = -0.02874884769890613;
        bool r19691 = r19689 <= r19690;
        float r19692 = sin(r19689);
        float r19693 = r19689 - r19692;
        float r19694 = tan(r19689);
        float r19695 = r19689 - r19694;
        float r19696 = r19693 / r19695;
        float r19697 = 0.03136174862029996;
        bool r19698 = r19689 <= r19697;
        float r19699 = 9/40;
        float r19700 = 2;
        float r19701 = pow(r19689, r19700);
        float r19702 = r19699 * r19701;
        float r19703 = 1/2;
        float r19704 = 27/2800;
        float r19705 = 4;
        float r19706 = pow(r19689, r19705);
        float r19707 = r19704 * r19706;
        float r19708 = r19703 + r19707;
        float r19709 = r19702 - r19708;
        float r19710 = r19698 ? r19709 : r19696;
        float r19711 = r19691 ? r19696 : r19710;
        return r19711;
}

double f_od(double x) {
        double r19712 = x;
        double r19713 = -0.02874884769890613;
        bool r19714 = r19712 <= r19713;
        double r19715 = sin(r19712);
        double r19716 = r19712 - r19715;
        double r19717 = tan(r19712);
        double r19718 = r19712 - r19717;
        double r19719 = r19716 / r19718;
        double r19720 = 0.03136174862029996;
        bool r19721 = r19712 <= r19720;
        double r19722 = 9/40;
        double r19723 = 2;
        double r19724 = pow(r19712, r19723);
        double r19725 = r19722 * r19724;
        double r19726 = 1/2;
        double r19727 = 27/2800;
        double r19728 = 4;
        double r19729 = pow(r19712, r19728);
        double r19730 = r19727 * r19729;
        double r19731 = r19726 + r19730;
        double r19732 = r19725 - r19731;
        double r19733 = r19721 ? r19732 : r19719;
        double r19734 = r19714 ? r19719 : r19733;
        return r19734;
}

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 r19735, r19736, r19737, r19738, r19739, r19740;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19735);
        mpfr_init(r19736);
        mpfr_init(r19737);
        mpfr_init(r19738);
        mpfr_init(r19739);
        mpfr_init(r19740);
}

double f_im(double x) {
        mpfr_set_d(r19735, x, MPFR_RNDN);
        mpfr_sin(r19736, r19735, MPFR_RNDN);
        mpfr_sub(r19737, r19735, r19736, MPFR_RNDN);
        mpfr_tan(r19738, r19735, MPFR_RNDN);
        mpfr_sub(r19739, r19735, r19738, MPFR_RNDN);
        mpfr_div(r19740, r19737, r19739, MPFR_RNDN);
        return mpfr_get_d(r19740, MPFR_RNDN);
}

static mpfr_t r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19741);
        mpfr_init_set_str(r19742, "-0.02874884769890613", 10, MPFR_RNDN);
        mpfr_init(r19743);
        mpfr_init(r19744);
        mpfr_init(r19745);
        mpfr_init(r19746);
        mpfr_init(r19747);
        mpfr_init(r19748);
        mpfr_init_set_str(r19749, "0.03136174862029996", 10, MPFR_RNDN);
        mpfr_init(r19750);
        mpfr_init_set_str(r19751, "9/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r19752, "2", 10, MPFR_RNDN);
        mpfr_init(r19753);
        mpfr_init(r19754);
        mpfr_init_set_str(r19755, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r19756, "27/2800", 10, MPFR_RNDN);
        mpfr_init_set_str(r19757, "4", 10, MPFR_RNDN);
        mpfr_init(r19758);
        mpfr_init(r19759);
        mpfr_init(r19760);
        mpfr_init(r19761);
        mpfr_init(r19762);
        mpfr_init(r19763);
}

double f_fm(double x) {
        mpfr_set_d(r19741, x, MPFR_RNDN);
        ;
        mpfr_set_si(r19743, mpfr_cmp(r19741, r19742) <= 0, MPFR_RNDN);
        mpfr_sin(r19744, r19741, MPFR_RNDN);
        mpfr_sub(r19745, r19741, r19744, MPFR_RNDN);
        mpfr_tan(r19746, r19741, MPFR_RNDN);
        mpfr_sub(r19747, r19741, r19746, MPFR_RNDN);
        mpfr_div(r19748, r19745, r19747, MPFR_RNDN);
        ;
        mpfr_set_si(r19750, mpfr_cmp(r19741, r19749) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19753, r19741, r19752, MPFR_RNDN);
        mpfr_mul(r19754, r19751, r19753, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_pow(r19758, r19741, r19757, MPFR_RNDN);
        mpfr_mul(r19759, r19756, r19758, MPFR_RNDN);
        mpfr_add(r19760, r19755, r19759, MPFR_RNDN);
        mpfr_sub(r19761, r19754, r19760, MPFR_RNDN);
        if (mpfr_get_si(r19750, MPFR_RNDN)) { mpfr_set(r19762, r19761, MPFR_RNDN); } else { mpfr_set(r19762, r19748, MPFR_RNDN); };
        if (mpfr_get_si(r19743, MPFR_RNDN)) { mpfr_set(r19763, r19748, MPFR_RNDN); } else { mpfr_set(r19763, r19762, MPFR_RNDN); };
        return mpfr_get_d(r19763, MPFR_RNDN);
}

static mpfr_t r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r19764);
        mpfr_init_set_str(r19765, "-0.02874884769890613", 10, MPFR_RNDN);
        mpfr_init(r19766);
        mpfr_init(r19767);
        mpfr_init(r19768);
        mpfr_init(r19769);
        mpfr_init(r19770);
        mpfr_init(r19771);
        mpfr_init_set_str(r19772, "0.03136174862029996", 10, MPFR_RNDN);
        mpfr_init(r19773);
        mpfr_init_set_str(r19774, "9/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r19775, "2", 10, MPFR_RNDN);
        mpfr_init(r19776);
        mpfr_init(r19777);
        mpfr_init_set_str(r19778, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r19779, "27/2800", 10, MPFR_RNDN);
        mpfr_init_set_str(r19780, "4", 10, MPFR_RNDN);
        mpfr_init(r19781);
        mpfr_init(r19782);
        mpfr_init(r19783);
        mpfr_init(r19784);
        mpfr_init(r19785);
        mpfr_init(r19786);
}

double f_dm(double x) {
        mpfr_set_d(r19764, x, MPFR_RNDN);
        ;
        mpfr_set_si(r19766, mpfr_cmp(r19764, r19765) <= 0, MPFR_RNDN);
        mpfr_sin(r19767, r19764, MPFR_RNDN);
        mpfr_sub(r19768, r19764, r19767, MPFR_RNDN);
        mpfr_tan(r19769, r19764, MPFR_RNDN);
        mpfr_sub(r19770, r19764, r19769, MPFR_RNDN);
        mpfr_div(r19771, r19768, r19770, MPFR_RNDN);
        ;
        mpfr_set_si(r19773, mpfr_cmp(r19764, r19772) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r19776, r19764, r19775, MPFR_RNDN);
        mpfr_mul(r19777, r19774, r19776, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_pow(r19781, r19764, r19780, MPFR_RNDN);
        mpfr_mul(r19782, r19779, r19781, MPFR_RNDN);
        mpfr_add(r19783, r19778, r19782, MPFR_RNDN);
        mpfr_sub(r19784, r19777, r19783, MPFR_RNDN);
        if (mpfr_get_si(r19773, MPFR_RNDN)) { mpfr_set(r19785, r19784, MPFR_RNDN); } else { mpfr_set(r19785, r19771, MPFR_RNDN); };
        if (mpfr_get_si(r19766, MPFR_RNDN)) { mpfr_set(r19786, r19771, MPFR_RNDN); } else { mpfr_set(r19786, r19785, MPFR_RNDN); };
        return mpfr_get_d(r19786, MPFR_RNDN);
}