#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 r25728 = x;
        float r25729 = eps;
        float r25730 = r25728 + r25729;
        float r25731 = cos(r25730);
        float r25732 = cos(r25728);
        float r25733 = r25731 - r25732;
        return r25733;
}

double f_id(double x, double eps) {
        double r25734 = x;
        double r25735 = eps;
        double r25736 = r25734 + r25735;
        double r25737 = cos(r25736);
        double r25738 = cos(r25734);
        double r25739 = r25737 - r25738;
        return r25739;
}


double f_of(float x, float eps) {
        float r25740 = -2;
        float r25741 = eps;
        float r25742 = 2;
        float r25743 = r25741 / r25742;
        float r25744 = sin(r25743);
        float r25745 = x;
        float r25746 = r25741 + r25745;
        float r25747 = r25745 + r25746;
        float r25748 = r25747 / r25742;
        float r25749 = sin(r25748);
        float r25750 = r25744 * r25749;
        float r25751 = r25740 * r25750;
        float r25752 = cbrt(r25751);
        float r25753 = r25752 * r25752;
        float r25754 = r25753 * r25752;
        float r25755 = -1.0181736631411995e-07;
        bool r25756 = r25754 <= r25755;
        float r25757 = cos(r25745);
        float r25758 = cos(r25741);
        float r25759 = r25757 * r25758;
        float r25760 = sin(r25741);
        float r25761 = sin(r25745);
        float r25762 = fma(r25760, r25761, r25757);
        float r25763 = r25759 - r25762;
        float r25764 = 2.8999654517490028e-06;
        bool r25765 = r25754 <= r25764;
        float r25766 = expm1(r25749);
        float r25767 = log1p(r25766);
        float r25768 = r25744 * r25767;
        float r25769 = r25740 * r25768;
        float r25770 = r25761 * r25760;
        float r25771 = r25759 - r25770;
        float r25772 = r25771 - r25757;
        float r25773 = r25765 ? r25769 : r25772;
        float r25774 = r25756 ? r25763 : r25773;
        return r25774;
}

double f_od(double x, double eps) {
        double r25775 = -2;
        double r25776 = eps;
        double r25777 = 2;
        double r25778 = r25776 / r25777;
        double r25779 = sin(r25778);
        double r25780 = x;
        double r25781 = r25776 + r25780;
        double r25782 = r25780 + r25781;
        double r25783 = r25782 / r25777;
        double r25784 = sin(r25783);
        double r25785 = r25779 * r25784;
        double r25786 = r25775 * r25785;
        double r25787 = cbrt(r25786);
        double r25788 = r25787 * r25787;
        double r25789 = r25788 * r25787;
        double r25790 = -1.0181736631411995e-07;
        bool r25791 = r25789 <= r25790;
        double r25792 = cos(r25780);
        double r25793 = cos(r25776);
        double r25794 = r25792 * r25793;
        double r25795 = sin(r25776);
        double r25796 = sin(r25780);
        double r25797 = fma(r25795, r25796, r25792);
        double r25798 = r25794 - r25797;
        double r25799 = 2.8999654517490028e-06;
        bool r25800 = r25789 <= r25799;
        double r25801 = expm1(r25784);
        double r25802 = log1p(r25801);
        double r25803 = r25779 * r25802;
        double r25804 = r25775 * r25803;
        double r25805 = r25796 * r25795;
        double r25806 = r25794 - r25805;
        double r25807 = r25806 - r25792;
        double r25808 = r25800 ? r25804 : r25807;
        double r25809 = r25791 ? r25798 : r25808;
        return r25809;
}

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 r25810, r25811, r25812, r25813, r25814, r25815;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r25810);
        mpfr_init(r25811);
        mpfr_init(r25812);
        mpfr_init(r25813);
        mpfr_init(r25814);
        mpfr_init(r25815);
}

double f_im(double x, double eps) {
        mpfr_set_d(r25810, x, MPFR_RNDN);
        mpfr_set_d(r25811, eps, MPFR_RNDN);
        mpfr_add(r25812, r25810, r25811, MPFR_RNDN);
        mpfr_cos(r25813, r25812, MPFR_RNDN);
        mpfr_cos(r25814, r25810, MPFR_RNDN);
        mpfr_sub(r25815, r25813, r25814, MPFR_RNDN);
        return mpfr_get_d(r25815, MPFR_RNDN);
}

static mpfr_t r25816, r25817, r25818, r25819, r25820, r25821, r25822, r25823, r25824, r25825, r25826, r25827, r25828, r25829, r25830, r25831, r25832, r25833, r25834, r25835, r25836, r25837, r25838, r25839, r25840, r25841, r25842, r25843, r25844, r25845, r25846, r25847, r25848, r25849, r25850;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r25816, "-2", 10, MPFR_RNDN);
        mpfr_init(r25817);
        mpfr_init_set_str(r25818, "2", 10, MPFR_RNDN);
        mpfr_init(r25819);
        mpfr_init(r25820);
        mpfr_init(r25821);
        mpfr_init(r25822);
        mpfr_init(r25823);
        mpfr_init(r25824);
        mpfr_init(r25825);
        mpfr_init(r25826);
        mpfr_init(r25827);
        mpfr_init(r25828);
        mpfr_init(r25829);
        mpfr_init(r25830);
        mpfr_init_set_str(r25831, "-1.0181736631411995e-07", 10, MPFR_RNDN);
        mpfr_init(r25832);
        mpfr_init(r25833);
        mpfr_init(r25834);
        mpfr_init(r25835);
        mpfr_init(r25836);
        mpfr_init(r25837);
        mpfr_init(r25838);
        mpfr_init(r25839);
        mpfr_init_set_str(r25840, "2.8999654517490028e-06", 10, MPFR_RNDN);
        mpfr_init(r25841);
        mpfr_init(r25842);
        mpfr_init(r25843);
        mpfr_init(r25844);
        mpfr_init(r25845);
        mpfr_init(r25846);
        mpfr_init(r25847);
        mpfr_init(r25848);
        mpfr_init(r25849);
        mpfr_init(r25850);
}

double f_fm(double x, double eps) {
        ;
        mpfr_set_d(r25817, eps, MPFR_RNDN);
        ;
        mpfr_div(r25819, r25817, r25818, MPFR_RNDN);
        mpfr_sin(r25820, r25819, MPFR_RNDN);
        mpfr_set_d(r25821, x, MPFR_RNDN);
        mpfr_add(r25822, r25817, r25821, MPFR_RNDN);
        mpfr_add(r25823, r25821, r25822, MPFR_RNDN);
        mpfr_div(r25824, r25823, r25818, MPFR_RNDN);
        mpfr_sin(r25825, r25824, MPFR_RNDN);
        mpfr_mul(r25826, r25820, r25825, MPFR_RNDN);
        mpfr_mul(r25827, r25816, r25826, MPFR_RNDN);
        mpfr_cbrt(r25828, r25827, MPFR_RNDN);
        mpfr_mul(r25829, r25828, r25828, MPFR_RNDN);
        mpfr_mul(r25830, r25829, r25828, MPFR_RNDN);
        ;
        mpfr_set_si(r25832, mpfr_cmp(r25830, r25831) <= 0, MPFR_RNDN);
        mpfr_cos(r25833, r25821, MPFR_RNDN);
        mpfr_cos(r25834, r25817, MPFR_RNDN);
        mpfr_mul(r25835, r25833, r25834, MPFR_RNDN);
        mpfr_sin(r25836, r25817, MPFR_RNDN);
        mpfr_sin(r25837, r25821, MPFR_RNDN);
        mpfr_fma(r25838, r25836, r25837, r25833, MPFR_RNDN);
        mpfr_sub(r25839, r25835, r25838, MPFR_RNDN);
        ;
        mpfr_set_si(r25841, mpfr_cmp(r25830, r25840) <= 0, MPFR_RNDN);
        mpfr_expm1(r25842, r25825, MPFR_RNDN);
        mpfr_log1p(r25843, r25842, MPFR_RNDN);
        mpfr_mul(r25844, r25820, r25843, MPFR_RNDN);
        mpfr_mul(r25845, r25816, r25844, MPFR_RNDN);
        mpfr_mul(r25846, r25837, r25836, MPFR_RNDN);
        mpfr_sub(r25847, r25835, r25846, MPFR_RNDN);
        mpfr_sub(r25848, r25847, r25833, MPFR_RNDN);
        if (mpfr_get_si(r25841, MPFR_RNDN)) { mpfr_set(r25849, r25845, MPFR_RNDN); } else { mpfr_set(r25849, r25848, MPFR_RNDN); };
        if (mpfr_get_si(r25832, MPFR_RNDN)) { mpfr_set(r25850, r25839, MPFR_RNDN); } else { mpfr_set(r25850, r25849, MPFR_RNDN); };
        return mpfr_get_d(r25850, MPFR_RNDN);
}

static mpfr_t r25851, r25852, r25853, r25854, r25855, r25856, r25857, r25858, r25859, r25860, r25861, r25862, r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871, r25872, r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r25851, "-2", 10, MPFR_RNDN);
        mpfr_init(r25852);
        mpfr_init_set_str(r25853, "2", 10, MPFR_RNDN);
        mpfr_init(r25854);
        mpfr_init(r25855);
        mpfr_init(r25856);
        mpfr_init(r25857);
        mpfr_init(r25858);
        mpfr_init(r25859);
        mpfr_init(r25860);
        mpfr_init(r25861);
        mpfr_init(r25862);
        mpfr_init(r25863);
        mpfr_init(r25864);
        mpfr_init(r25865);
        mpfr_init_set_str(r25866, "-1.0181736631411995e-07", 10, MPFR_RNDN);
        mpfr_init(r25867);
        mpfr_init(r25868);
        mpfr_init(r25869);
        mpfr_init(r25870);
        mpfr_init(r25871);
        mpfr_init(r25872);
        mpfr_init(r25873);
        mpfr_init(r25874);
        mpfr_init_set_str(r25875, "2.8999654517490028e-06", 10, MPFR_RNDN);
        mpfr_init(r25876);
        mpfr_init(r25877);
        mpfr_init(r25878);
        mpfr_init(r25879);
        mpfr_init(r25880);
        mpfr_init(r25881);
        mpfr_init(r25882);
        mpfr_init(r25883);
        mpfr_init(r25884);
        mpfr_init(r25885);
}

double f_dm(double x, double eps) {
        ;
        mpfr_set_d(r25852, eps, MPFR_RNDN);
        ;
        mpfr_div(r25854, r25852, r25853, MPFR_RNDN);
        mpfr_sin(r25855, r25854, MPFR_RNDN);
        mpfr_set_d(r25856, x, MPFR_RNDN);
        mpfr_add(r25857, r25852, r25856, MPFR_RNDN);
        mpfr_add(r25858, r25856, r25857, MPFR_RNDN);
        mpfr_div(r25859, r25858, r25853, MPFR_RNDN);
        mpfr_sin(r25860, r25859, MPFR_RNDN);
        mpfr_mul(r25861, r25855, r25860, MPFR_RNDN);
        mpfr_mul(r25862, r25851, r25861, MPFR_RNDN);
        mpfr_cbrt(r25863, r25862, MPFR_RNDN);
        mpfr_mul(r25864, r25863, r25863, MPFR_RNDN);
        mpfr_mul(r25865, r25864, r25863, MPFR_RNDN);
        ;
        mpfr_set_si(r25867, mpfr_cmp(r25865, r25866) <= 0, MPFR_RNDN);
        mpfr_cos(r25868, r25856, MPFR_RNDN);
        mpfr_cos(r25869, r25852, MPFR_RNDN);
        mpfr_mul(r25870, r25868, r25869, MPFR_RNDN);
        mpfr_sin(r25871, r25852, MPFR_RNDN);
        mpfr_sin(r25872, r25856, MPFR_RNDN);
        mpfr_fma(r25873, r25871, r25872, r25868, MPFR_RNDN);
        mpfr_sub(r25874, r25870, r25873, MPFR_RNDN);
        ;
        mpfr_set_si(r25876, mpfr_cmp(r25865, r25875) <= 0, MPFR_RNDN);
        mpfr_expm1(r25877, r25860, MPFR_RNDN);
        mpfr_log1p(r25878, r25877, MPFR_RNDN);
        mpfr_mul(r25879, r25855, r25878, MPFR_RNDN);
        mpfr_mul(r25880, r25851, r25879, MPFR_RNDN);
        mpfr_mul(r25881, r25872, r25871, MPFR_RNDN);
        mpfr_sub(r25882, r25870, r25881, MPFR_RNDN);
        mpfr_sub(r25883, r25882, r25868, MPFR_RNDN);
        if (mpfr_get_si(r25876, MPFR_RNDN)) { mpfr_set(r25884, r25880, MPFR_RNDN); } else { mpfr_set(r25884, r25883, MPFR_RNDN); };
        if (mpfr_get_si(r25867, MPFR_RNDN)) { mpfr_set(r25885, r25874, MPFR_RNDN); } else { mpfr_set(r25885, r25884, MPFR_RNDN); };
        return mpfr_get_d(r25885, MPFR_RNDN);
}