#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 r18844 = x;
        float r18845 = eps;
        float r18846 = r18844 + r18845;
        float r18847 = cos(r18846);
        float r18848 = cos(r18844);
        float r18849 = r18847 - r18848;
        return r18849;
}

double f_id(double x, double eps) {
        double r18850 = x;
        double r18851 = eps;
        double r18852 = r18850 + r18851;
        double r18853 = cos(r18852);
        double r18854 = cos(r18850);
        double r18855 = r18853 - r18854;
        return r18855;
}


double f_of(float x, float eps) {
        float r18856 = eps;
        float r18857 = -0.11215751184115397f;
        bool r18858 = r18856 <= r18857;
        float r18859 = x;
        float r18860 = cos(r18859);
        float r18861 = cos(r18856);
        float r18862 = r18860 * r18861;
        float r18863 = sin(r18859);
        float r18864 = sin(r18856);
        float r18865 = r18863 * r18864;
        float r18866 = r18862 - r18865;
        float r18867 = r18866 - r18860;
        float r18868 = 0.005563461367400353f;
        bool r18869 = r18856 <= r18868;
        float r18870 = -2.0f;
        float r18871 = 2.0f;
        float r18872 = r18856 / r18871;
        float r18873 = sin(r18872);
        float r18874 = r18859 + r18859;
        float r18875 = r18856 + r18874;
        float r18876 = r18875 / r18871;
        float r18877 = sin(r18876);
        float r18878 = r18873 * r18877;
        float r18879 = r18870 * r18878;
        float r18880 = r18869 ? r18879 : r18867;
        float r18881 = r18858 ? r18867 : r18880;
        return r18881;
}

double f_od(double x, double eps) {
        double r18882 = eps;
        double r18883 = -0.11215751184115397;
        bool r18884 = r18882 <= r18883;
        double r18885 = x;
        double r18886 = cos(r18885);
        double r18887 = cos(r18882);
        double r18888 = r18886 * r18887;
        double r18889 = sin(r18885);
        double r18890 = sin(r18882);
        double r18891 = r18889 * r18890;
        double r18892 = r18888 - r18891;
        double r18893 = r18892 - r18886;
        double r18894 = 0.005563461367400353;
        bool r18895 = r18882 <= r18894;
        double r18896 = -2.0;
        double r18897 = 2.0;
        double r18898 = r18882 / r18897;
        double r18899 = sin(r18898);
        double r18900 = r18885 + r18885;
        double r18901 = r18882 + r18900;
        double r18902 = r18901 / r18897;
        double r18903 = sin(r18902);
        double r18904 = r18899 * r18903;
        double r18905 = r18896 * r18904;
        double r18906 = r18895 ? r18905 : r18893;
        double r18907 = r18884 ? r18893 : r18906;
        return r18907;
}

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 r18908, r18909, r18910, r18911, r18912, r18913;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2192);
        mpfr_init(r18908);
        mpfr_init(r18909);
        mpfr_init(r18910);
        mpfr_init(r18911);
        mpfr_init(r18912);
        mpfr_init(r18913);
}

double f_im(double x, double eps) {
        mpfr_set_d(r18908, x, MPFR_RNDN);
        mpfr_set_d(r18909, eps, MPFR_RNDN);
        mpfr_add(r18910, r18908, r18909, MPFR_RNDN);
        mpfr_cos(r18911, r18910, MPFR_RNDN);
        mpfr_cos(r18912, r18908, MPFR_RNDN);
        mpfr_sub(r18913, r18911, r18912, MPFR_RNDN);
        return mpfr_get_d(r18913, MPFR_RNDN);
}

static mpfr_t r18914, r18915, r18916, r18917, r18918, r18919, r18920, r18921, r18922, r18923, r18924, r18925, r18926, r18927, r18928, r18929, r18930, r18931, r18932, r18933, r18934, r18935, r18936, r18937, r18938, r18939;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r18914);
        mpfr_init_set_str(r18915, "-0.11215751184115397", 10, MPFR_RNDN);
        mpfr_init(r18916);
        mpfr_init(r18917);
        mpfr_init(r18918);
        mpfr_init(r18919);
        mpfr_init(r18920);
        mpfr_init(r18921);
        mpfr_init(r18922);
        mpfr_init(r18923);
        mpfr_init(r18924);
        mpfr_init(r18925);
        mpfr_init_set_str(r18926, "0.005563461367400353", 10, MPFR_RNDN);
        mpfr_init(r18927);
        mpfr_init_set_str(r18928, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r18929, "2", 10, MPFR_RNDN);
        mpfr_init(r18930);
        mpfr_init(r18931);
        mpfr_init(r18932);
        mpfr_init(r18933);
        mpfr_init(r18934);
        mpfr_init(r18935);
        mpfr_init(r18936);
        mpfr_init(r18937);
        mpfr_init(r18938);
        mpfr_init(r18939);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r18914, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r18916, mpfr_cmp(r18914, r18915) <= 0, MPFR_RNDN);
        mpfr_set_d(r18917, x, MPFR_RNDN);
        mpfr_cos(r18918, r18917, MPFR_RNDN);
        mpfr_cos(r18919, r18914, MPFR_RNDN);
        mpfr_mul(r18920, r18918, r18919, MPFR_RNDN);
        mpfr_sin(r18921, r18917, MPFR_RNDN);
        mpfr_sin(r18922, r18914, MPFR_RNDN);
        mpfr_mul(r18923, r18921, r18922, MPFR_RNDN);
        mpfr_sub(r18924, r18920, r18923, MPFR_RNDN);
        mpfr_sub(r18925, r18924, r18918, MPFR_RNDN);
        ;
        mpfr_set_si(r18927, mpfr_cmp(r18914, r18926) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r18930, r18914, r18929, MPFR_RNDN);
        mpfr_sin(r18931, r18930, MPFR_RNDN);
        mpfr_add(r18932, r18917, r18917, MPFR_RNDN);
        mpfr_add(r18933, r18914, r18932, MPFR_RNDN);
        mpfr_div(r18934, r18933, r18929, MPFR_RNDN);
        mpfr_sin(r18935, r18934, MPFR_RNDN);
        mpfr_mul(r18936, r18931, r18935, MPFR_RNDN);
        mpfr_mul(r18937, r18928, r18936, MPFR_RNDN);
        if (mpfr_get_si(r18927, MPFR_RNDN)) { mpfr_set(r18938, r18937, MPFR_RNDN); } else { mpfr_set(r18938, r18925, MPFR_RNDN); };
        if (mpfr_get_si(r18916, MPFR_RNDN)) { mpfr_set(r18939, r18925, MPFR_RNDN); } else { mpfr_set(r18939, r18938, MPFR_RNDN); };
        return mpfr_get_d(r18939, MPFR_RNDN);
}

static mpfr_t r18940, r18941, r18942, r18943, r18944, r18945, r18946, r18947, r18948, r18949, r18950, r18951, r18952, r18953, r18954, r18955, r18956, r18957, r18958, r18959, r18960, r18961, r18962, r18963, r18964, r18965;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r18940);
        mpfr_init_set_str(r18941, "-0.11215751184115397", 10, MPFR_RNDN);
        mpfr_init(r18942);
        mpfr_init(r18943);
        mpfr_init(r18944);
        mpfr_init(r18945);
        mpfr_init(r18946);
        mpfr_init(r18947);
        mpfr_init(r18948);
        mpfr_init(r18949);
        mpfr_init(r18950);
        mpfr_init(r18951);
        mpfr_init_set_str(r18952, "0.005563461367400353", 10, MPFR_RNDN);
        mpfr_init(r18953);
        mpfr_init_set_str(r18954, "-2", 10, MPFR_RNDN);
        mpfr_init_set_str(r18955, "2", 10, MPFR_RNDN);
        mpfr_init(r18956);
        mpfr_init(r18957);
        mpfr_init(r18958);
        mpfr_init(r18959);
        mpfr_init(r18960);
        mpfr_init(r18961);
        mpfr_init(r18962);
        mpfr_init(r18963);
        mpfr_init(r18964);
        mpfr_init(r18965);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r18940, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r18942, mpfr_cmp(r18940, r18941) <= 0, MPFR_RNDN);
        mpfr_set_d(r18943, x, MPFR_RNDN);
        mpfr_cos(r18944, r18943, MPFR_RNDN);
        mpfr_cos(r18945, r18940, MPFR_RNDN);
        mpfr_mul(r18946, r18944, r18945, MPFR_RNDN);
        mpfr_sin(r18947, r18943, MPFR_RNDN);
        mpfr_sin(r18948, r18940, MPFR_RNDN);
        mpfr_mul(r18949, r18947, r18948, MPFR_RNDN);
        mpfr_sub(r18950, r18946, r18949, MPFR_RNDN);
        mpfr_sub(r18951, r18950, r18944, MPFR_RNDN);
        ;
        mpfr_set_si(r18953, mpfr_cmp(r18940, r18952) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_div(r18956, r18940, r18955, MPFR_RNDN);
        mpfr_sin(r18957, r18956, MPFR_RNDN);
        mpfr_add(r18958, r18943, r18943, MPFR_RNDN);
        mpfr_add(r18959, r18940, r18958, MPFR_RNDN);
        mpfr_div(r18960, r18959, r18955, MPFR_RNDN);
        mpfr_sin(r18961, r18960, MPFR_RNDN);
        mpfr_mul(r18962, r18957, r18961, MPFR_RNDN);
        mpfr_mul(r18963, r18954, r18962, MPFR_RNDN);
        if (mpfr_get_si(r18953, MPFR_RNDN)) { mpfr_set(r18964, r18963, MPFR_RNDN); } else { mpfr_set(r18964, r18951, MPFR_RNDN); };
        if (mpfr_get_si(r18942, MPFR_RNDN)) { mpfr_set(r18965, r18951, MPFR_RNDN); } else { mpfr_set(r18965, r18964, MPFR_RNDN); };
        return mpfr_get_d(r18965, MPFR_RNDN);
}

