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

char *name = "2log (problem 3.3.6)";

double f_if(float N) {
        float r20899 = N;
        float r20900 = 1;
        float r20901 = r20899 + r20900;
        float r20902 = log(r20901);
        float r20903 = log(r20899);
        float r20904 = r20902 - r20903;
        return r20904;
}

double f_id(double N) {
        double r20905 = N;
        double r20906 = 1;
        double r20907 = r20905 + r20906;
        double r20908 = log(r20907);
        double r20909 = log(r20905);
        double r20910 = r20908 - r20909;
        return r20910;
}


double f_of(float N) {
        float r20911 = N;
        float r20912 = 1;
        float r20913 = r20911 + r20912;
        float r20914 = log(r20913);
        float r20915 = log(r20911);
        float r20916 = r20914 - r20915;
        float r20917 = 1.5458009983149168e-08;
        bool r20918 = r20916 <= r20917;
        float r20919 = r20912 / r20911;
        float r20920 = log(r20912);
        float r20921 = r20919 - r20920;
        float r20922 = 1/2;
        float r20923 = r20922 / r20911;
        float r20924 = r20923 / r20911;
        float r20925 = r20921 - r20924;
        float r20926 = 3;
        float r20927 = pow(r20911, r20926);
        float r20928 = r20912 + r20927;
        float r20929 = log(r20928);
        float r20930 = r20911 * r20911;
        float r20931 = r20912 - r20911;
        float r20932 = r20930 + r20931;
        float r20933 = log(r20932);
        float r20934 = r20929 - r20933;
        float r20935 = r20934 - r20915;
        float r20936 = r20918 ? r20925 : r20935;
        return r20936;
}

double f_od(double N) {
        double r20937 = N;
        double r20938 = 1;
        double r20939 = r20937 + r20938;
        double r20940 = log(r20939);
        double r20941 = log(r20937);
        double r20942 = r20940 - r20941;
        double r20943 = 1.5458009983149168e-08;
        bool r20944 = r20942 <= r20943;
        double r20945 = r20938 / r20937;
        double r20946 = log(r20938);
        double r20947 = r20945 - r20946;
        double r20948 = 1/2;
        double r20949 = r20948 / r20937;
        double r20950 = r20949 / r20937;
        double r20951 = r20947 - r20950;
        double r20952 = 3;
        double r20953 = pow(r20937, r20952);
        double r20954 = r20938 + r20953;
        double r20955 = log(r20954);
        double r20956 = r20937 * r20937;
        double r20957 = r20938 - r20937;
        double r20958 = r20956 + r20957;
        double r20959 = log(r20958);
        double r20960 = r20955 - r20959;
        double r20961 = r20960 - r20941;
        double r20962 = r20944 ? r20951 : r20961;
        return r20962;
}

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 r20963, r20964, r20965, r20966, r20967, r20968;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20963);
        mpfr_init_set_str(r20964, "1", 10, MPFR_RNDN);
        mpfr_init(r20965);
        mpfr_init(r20966);
        mpfr_init(r20967);
        mpfr_init(r20968);
}

double f_im(double N) {
        mpfr_set_d(r20963, N, MPFR_RNDN);
        ;
        mpfr_add(r20965, r20963, r20964, MPFR_RNDN);
        mpfr_log(r20966, r20965, MPFR_RNDN);
        mpfr_log(r20967, r20963, MPFR_RNDN);
        mpfr_sub(r20968, r20966, r20967, MPFR_RNDN);
        return mpfr_get_d(r20968, MPFR_RNDN);
}

static mpfr_t r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976, r20977, r20978, r20979, r20980, r20981, r20982, r20983, r20984, r20985, r20986, r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20969);
        mpfr_init_set_str(r20970, "1", 10, MPFR_RNDN);
        mpfr_init(r20971);
        mpfr_init(r20972);
        mpfr_init(r20973);
        mpfr_init(r20974);
        mpfr_init_set_str(r20975, "1.5458009983149168e-08", 10, MPFR_RNDN);
        mpfr_init(r20976);
        mpfr_init(r20977);
        mpfr_init(r20978);
        mpfr_init(r20979);
        mpfr_init_set_str(r20980, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20981);
        mpfr_init(r20982);
        mpfr_init(r20983);
        mpfr_init_set_str(r20984, "3", 10, MPFR_RNDN);
        mpfr_init(r20985);
        mpfr_init(r20986);
        mpfr_init(r20987);
        mpfr_init(r20988);
        mpfr_init(r20989);
        mpfr_init(r20990);
        mpfr_init(r20991);
        mpfr_init(r20992);
        mpfr_init(r20993);
        mpfr_init(r20994);
}

double f_fm(double N) {
        mpfr_set_d(r20969, N, MPFR_RNDN);
        ;
        mpfr_add(r20971, r20969, r20970, MPFR_RNDN);
        mpfr_log(r20972, r20971, MPFR_RNDN);
        mpfr_log(r20973, r20969, MPFR_RNDN);
        mpfr_sub(r20974, r20972, r20973, MPFR_RNDN);
        ;
        mpfr_set_si(r20976, mpfr_cmp(r20974, r20975) <= 0, MPFR_RNDN);
        mpfr_div(r20977, r20970, r20969, MPFR_RNDN);
        mpfr_log(r20978, r20970, MPFR_RNDN);
        mpfr_sub(r20979, r20977, r20978, MPFR_RNDN);
        ;
        mpfr_div(r20981, r20980, r20969, MPFR_RNDN);
        mpfr_div(r20982, r20981, r20969, MPFR_RNDN);
        mpfr_sub(r20983, r20979, r20982, MPFR_RNDN);
        ;
        mpfr_pow(r20985, r20969, r20984, MPFR_RNDN);
        mpfr_add(r20986, r20970, r20985, MPFR_RNDN);
        mpfr_log(r20987, r20986, MPFR_RNDN);
        mpfr_mul(r20988, r20969, r20969, MPFR_RNDN);
        mpfr_sub(r20989, r20970, r20969, MPFR_RNDN);
        mpfr_add(r20990, r20988, r20989, MPFR_RNDN);
        mpfr_log(r20991, r20990, MPFR_RNDN);
        mpfr_sub(r20992, r20987, r20991, MPFR_RNDN);
        mpfr_sub(r20993, r20992, r20973, MPFR_RNDN);
        if (mpfr_get_si(r20976, MPFR_RNDN)) { mpfr_set(r20994, r20983, MPFR_RNDN); } else { mpfr_set(r20994, 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;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20995);
        mpfr_init_set_str(r20996, "1", 10, MPFR_RNDN);
        mpfr_init(r20997);
        mpfr_init(r20998);
        mpfr_init(r20999);
        mpfr_init(r21000);
        mpfr_init_set_str(r21001, "1.5458009983149168e-08", 10, MPFR_RNDN);
        mpfr_init(r21002);
        mpfr_init(r21003);
        mpfr_init(r21004);
        mpfr_init(r21005);
        mpfr_init_set_str(r21006, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21007);
        mpfr_init(r21008);
        mpfr_init(r21009);
        mpfr_init_set_str(r21010, "3", 10, MPFR_RNDN);
        mpfr_init(r21011);
        mpfr_init(r21012);
        mpfr_init(r21013);
        mpfr_init(r21014);
        mpfr_init(r21015);
        mpfr_init(r21016);
        mpfr_init(r21017);
        mpfr_init(r21018);
        mpfr_init(r21019);
        mpfr_init(r21020);
}

double f_dm(double N) {
        mpfr_set_d(r20995, N, MPFR_RNDN);
        ;
        mpfr_add(r20997, r20995, r20996, MPFR_RNDN);
        mpfr_log(r20998, r20997, MPFR_RNDN);
        mpfr_log(r20999, r20995, MPFR_RNDN);
        mpfr_sub(r21000, r20998, r20999, MPFR_RNDN);
        ;
        mpfr_set_si(r21002, mpfr_cmp(r21000, r21001) <= 0, MPFR_RNDN);
        mpfr_div(r21003, r20996, r20995, MPFR_RNDN);
        mpfr_log(r21004, r20996, MPFR_RNDN);
        mpfr_sub(r21005, r21003, r21004, MPFR_RNDN);
        ;
        mpfr_div(r21007, r21006, r20995, MPFR_RNDN);
        mpfr_div(r21008, r21007, r20995, MPFR_RNDN);
        mpfr_sub(r21009, r21005, r21008, MPFR_RNDN);
        ;
        mpfr_pow(r21011, r20995, r21010, MPFR_RNDN);
        mpfr_add(r21012, r20996, r21011, MPFR_RNDN);
        mpfr_log(r21013, r21012, MPFR_RNDN);
        mpfr_mul(r21014, r20995, r20995, MPFR_RNDN);
        mpfr_sub(r21015, r20996, r20995, MPFR_RNDN);
        mpfr_add(r21016, r21014, r21015, MPFR_RNDN);
        mpfr_log(r21017, r21016, MPFR_RNDN);
        mpfr_sub(r21018, r21013, r21017, MPFR_RNDN);
        mpfr_sub(r21019, r21018, r20999, MPFR_RNDN);
        if (mpfr_get_si(r21002, MPFR_RNDN)) { mpfr_set(r21020, r21009, MPFR_RNDN); } else { mpfr_set(r21020, r21019, MPFR_RNDN); };
        return mpfr_get_d(r21020, MPFR_RNDN);
}