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

char *name = "cos2 (problem 3.4.1)";

double f_if(float x) {
        float r20921 = 1;
        float r20922 = x;
        float r20923 = cos(r20922);
        float r20924 = r20921 - r20923;
        float r20925 = r20922 * r20922;
        float r20926 = r20924 / r20925;
        return r20926;
}

double f_id(double x) {
        double r20927 = 1;
        double r20928 = x;
        double r20929 = cos(r20928);
        double r20930 = r20927 - r20929;
        double r20931 = r20928 * r20928;
        double r20932 = r20930 / r20931;
        return r20932;
}


double f_of(float x) {
        float r20933 = x;
        float r20934 = -0.013047469261293217;
        bool r20935 = r20933 <= r20934;
        float r20936 = 2;
        float r20937 = r20933 / r20936;
        float r20938 = tan(r20937);
        float r20939 = sin(r20933);
        float r20940 = r20938 * r20939;
        float r20941 = sqrt(r20940);
        float r20942 = r20941 / r20933;
        float r20943 = r20942 * r20942;
        float r20944 = 0.0003791239906132239;
        bool r20945 = r20933 <= r20944;
        float r20946 = 4;
        float r20947 = pow(r20933, r20946);
        float r20948 = 1/720;
        float r20949 = 1/2;
        float r20950 = fma(r20947, r20948, r20949);
        float r20951 = r20933 * r20933;
        float r20952 = 1/24;
        float r20953 = r20951 * r20952;
        float r20954 = r20950 - r20953;
        float r20955 = -r20939;
        float r20956 = r20955 * r20938;
        float r20957 = -r20951;
        float r20958 = r20956 / r20957;
        float r20959 = r20945 ? r20954 : r20958;
        float r20960 = r20935 ? r20943 : r20959;
        return r20960;
}

double f_od(double x) {
        double r20961 = x;
        double r20962 = -0.013047469261293217;
        bool r20963 = r20961 <= r20962;
        double r20964 = 2;
        double r20965 = r20961 / r20964;
        double r20966 = tan(r20965);
        double r20967 = sin(r20961);
        double r20968 = r20966 * r20967;
        double r20969 = sqrt(r20968);
        double r20970 = r20969 / r20961;
        double r20971 = r20970 * r20970;
        double r20972 = 0.0003791239906132239;
        bool r20973 = r20961 <= r20972;
        double r20974 = 4;
        double r20975 = pow(r20961, r20974);
        double r20976 = 1/720;
        double r20977 = 1/2;
        double r20978 = fma(r20975, r20976, r20977);
        double r20979 = r20961 * r20961;
        double r20980 = 1/24;
        double r20981 = r20979 * r20980;
        double r20982 = r20978 - r20981;
        double r20983 = -r20967;
        double r20984 = r20983 * r20966;
        double r20985 = -r20979;
        double r20986 = r20984 / r20985;
        double r20987 = r20973 ? r20982 : r20986;
        double r20988 = r20963 ? r20971 : r20987;
        return r20988;
}

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 r20989, r20990, r20991, r20992, r20993, r20994;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20989, "1", 10, MPFR_RNDN);
        mpfr_init(r20990);
        mpfr_init(r20991);
        mpfr_init(r20992);
        mpfr_init(r20993);
        mpfr_init(r20994);
}

double f_im(double x) {
        ;
        mpfr_set_d(r20990, x, MPFR_RNDN);
        mpfr_cos(r20991, r20990, MPFR_RNDN);
        mpfr_sub(r20992, r20989, r20991, MPFR_RNDN);
        mpfr_mul(r20993, r20990, r20990, MPFR_RNDN);
        mpfr_div(r20994, r20992, r20993, MPFR_RNDN);
        return mpfr_get_d(r20994, MPFR_RNDN);
}

static mpfr_t r20995, r20996, r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020, r21021, r21022;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r20995);
        mpfr_init_set_str(r20996, "-0.013047469261293217", 10, MPFR_RNDN);
        mpfr_init(r20997);
        mpfr_init_set_str(r20998, "2", 10, MPFR_RNDN);
        mpfr_init(r20999);
        mpfr_init(r21000);
        mpfr_init(r21001);
        mpfr_init(r21002);
        mpfr_init(r21003);
        mpfr_init(r21004);
        mpfr_init(r21005);
        mpfr_init_set_str(r21006, "0.0003791239906132239", 10, MPFR_RNDN);
        mpfr_init(r21007);
        mpfr_init_set_str(r21008, "4", 10, MPFR_RNDN);
        mpfr_init(r21009);
        mpfr_init_set_str(r21010, "1/720", 10, MPFR_RNDN);
        mpfr_init_set_str(r21011, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21012);
        mpfr_init(r21013);
        mpfr_init_set_str(r21014, "1/24", 10, MPFR_RNDN);
        mpfr_init(r21015);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
        mpfr_init(r21019);
        mpfr_init(r21020);
        mpfr_init(r21021);
        mpfr_init(r21022);
}

double f_fm(double x) {
        mpfr_set_d(r20995, x, MPFR_RNDN);
        ;
        mpfr_set_si(r20997, mpfr_cmp(r20995, r20996) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r20999, r20995, r20998, MPFR_RNDN);
        mpfr_tan(r21000, r20999, MPFR_RNDN);
        mpfr_sin(r21001, r20995, MPFR_RNDN);
        mpfr_mul(r21002, r21000, r21001, MPFR_RNDN);
        mpfr_sqrt(r21003, r21002, MPFR_RNDN);
        mpfr_div(r21004, r21003, r20995, MPFR_RNDN);
        mpfr_mul(r21005, r21004, r21004, MPFR_RNDN);
        ;
        mpfr_set_si(r21007, mpfr_cmp(r20995, r21006) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r21009, r20995, r21008, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r21012, r21009, r21010, r21011, MPFR_RNDN);
        mpfr_mul(r21013, r20995, r20995, MPFR_RNDN);
        ;
        mpfr_mul(r21015, r21013, r21014, MPFR_RNDN);
        mpfr_sub(r21016, r21012, r21015, MPFR_RNDN);
        mpfr_neg(r21017, r21001, MPFR_RNDN);
        mpfr_mul(r21018, r21017, r21000, MPFR_RNDN);
        mpfr_neg(r21019, r21013, MPFR_RNDN);
        mpfr_div(r21020, r21018, r21019, MPFR_RNDN);
        if (mpfr_get_si(r21007, MPFR_RNDN)) { mpfr_set(r21021, r21016, MPFR_RNDN); } else { mpfr_set(r21021, r21020, MPFR_RNDN); };
        if (mpfr_get_si(r20997, MPFR_RNDN)) { mpfr_set(r21022, r21005, MPFR_RNDN); } else { mpfr_set(r21022, r21021, MPFR_RNDN); };
        return mpfr_get_d(r21022, MPFR_RNDN);
}

static mpfr_t r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042, r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r21023);
        mpfr_init_set_str(r21024, "-0.013047469261293217", 10, MPFR_RNDN);
        mpfr_init(r21025);
        mpfr_init_set_str(r21026, "2", 10, MPFR_RNDN);
        mpfr_init(r21027);
        mpfr_init(r21028);
        mpfr_init(r21029);
        mpfr_init(r21030);
        mpfr_init(r21031);
        mpfr_init(r21032);
        mpfr_init(r21033);
        mpfr_init_set_str(r21034, "0.0003791239906132239", 10, MPFR_RNDN);
        mpfr_init(r21035);
        mpfr_init_set_str(r21036, "4", 10, MPFR_RNDN);
        mpfr_init(r21037);
        mpfr_init_set_str(r21038, "1/720", 10, MPFR_RNDN);
        mpfr_init_set_str(r21039, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21040);
        mpfr_init(r21041);
        mpfr_init_set_str(r21042, "1/24", 10, MPFR_RNDN);
        mpfr_init(r21043);
        mpfr_init(r21044);
        mpfr_init(r21045);
        mpfr_init(r21046);
        mpfr_init(r21047);
        mpfr_init(r21048);
        mpfr_init(r21049);
        mpfr_init(r21050);
}

double f_dm(double x) {
        mpfr_set_d(r21023, x, MPFR_RNDN);
        ;
        mpfr_set_si(r21025, mpfr_cmp(r21023, r21024) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r21027, r21023, r21026, MPFR_RNDN);
        mpfr_tan(r21028, r21027, MPFR_RNDN);
        mpfr_sin(r21029, r21023, MPFR_RNDN);
        mpfr_mul(r21030, r21028, r21029, MPFR_RNDN);
        mpfr_sqrt(r21031, r21030, MPFR_RNDN);
        mpfr_div(r21032, r21031, r21023, MPFR_RNDN);
        mpfr_mul(r21033, r21032, r21032, MPFR_RNDN);
        ;
        mpfr_set_si(r21035, mpfr_cmp(r21023, r21034) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r21037, r21023, r21036, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r21040, r21037, r21038, r21039, MPFR_RNDN);
        mpfr_mul(r21041, r21023, r21023, MPFR_RNDN);
        ;
        mpfr_mul(r21043, r21041, r21042, MPFR_RNDN);
        mpfr_sub(r21044, r21040, r21043, MPFR_RNDN);
        mpfr_neg(r21045, r21029, MPFR_RNDN);
        mpfr_mul(r21046, r21045, r21028, MPFR_RNDN);
        mpfr_neg(r21047, r21041, MPFR_RNDN);
        mpfr_div(r21048, r21046, r21047, MPFR_RNDN);
        if (mpfr_get_si(r21035, MPFR_RNDN)) { mpfr_set(r21049, r21044, MPFR_RNDN); } else { mpfr_set(r21049, r21048, MPFR_RNDN); };
        if (mpfr_get_si(r21025, MPFR_RNDN)) { mpfr_set(r21050, r21033, MPFR_RNDN); } else { mpfr_set(r21050, r21049, MPFR_RNDN); };
        return mpfr_get_d(r21050, MPFR_RNDN);
}