#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 r9113 = N;
        float r9114 = 1.0f;
        float r9115 = r9113 + r9114;
        float r9116 = log(r9115);
        float r9117 = log(r9113);
        float r9118 = r9116 - r9117;
        return r9118;
}

double f_id(double N) {
        double r9119 = N;
        double r9120 = 1.0;
        double r9121 = r9119 + r9120;
        double r9122 = log(r9121);
        double r9123 = log(r9119);
        double r9124 = r9122 - r9123;
        return r9124;
}


double f_of(float N) {
        float r9125 = 1.0f;
        float r9126 = N;
        float r9127 = r9125 + r9126;
        float r9128 = log(r9127);
        float r9129 = log(r9126);
        float r9130 = r9128 - r9129;
        float r9131 = 3.542537891792339e-05f;
        bool r9132 = r9130 <= r9131;
        float r9133 = r9125 / r9126;
        float r9134 = 0.5f;
        float r9135 = 0.3333333333333333f;
        float r9136 = r9135 / r9126;
        float r9137 = r9134 - r9136;
        float r9138 = r9126 * r9126;
        float r9139 = r9137 / r9138;
        float r9140 = r9133 - r9139;
        float r9141 = r9127 / r9126;
        float r9142 = log(r9141);
        float r9143 = r9132 ? r9140 : r9142;
        return r9143;
}

double f_od(double N) {
        double r9144 = 1.0;
        double r9145 = N;
        double r9146 = r9144 + r9145;
        double r9147 = log(r9146);
        double r9148 = log(r9145);
        double r9149 = r9147 - r9148;
        double r9150 = 3.542537891792339e-05;
        bool r9151 = r9149 <= r9150;
        double r9152 = r9144 / r9145;
        double r9153 = 0.5;
        double r9154 = 0.3333333333333333;
        double r9155 = r9154 / r9145;
        double r9156 = r9153 - r9155;
        double r9157 = r9145 * r9145;
        double r9158 = r9156 / r9157;
        double r9159 = r9152 - r9158;
        double r9160 = r9146 / r9145;
        double r9161 = log(r9160);
        double r9162 = r9151 ? r9159 : r9161;
        return r9162;
}

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 r9163, r9164, r9165, r9166, r9167, r9168;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9163);
        mpfr_init_set_str(r9164, "1", 10, MPFR_RNDN);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
        mpfr_init(r9168);
}

double f_im(double N) {
        mpfr_set_d(r9163, N, MPFR_RNDN);
        ;
        mpfr_add(r9165, r9163, r9164, MPFR_RNDN);
        mpfr_log(r9166, r9165, MPFR_RNDN);
        mpfr_log(r9167, r9163, MPFR_RNDN);
        mpfr_sub(r9168, r9166, r9167, MPFR_RNDN);
        return mpfr_get_d(r9168, MPFR_RNDN);
}

static mpfr_t r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185, r9186, r9187;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9169, "1", 10, MPFR_RNDN);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init_set_str(r9175, "3.542537891792339e-05", 10, MPFR_RNDN);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init_set_str(r9178, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r9179, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r9170, N, MPFR_RNDN);
        mpfr_add(r9171, r9169, r9170, MPFR_RNDN);
        mpfr_log(r9172, r9171, MPFR_RNDN);
        mpfr_log(r9173, r9170, MPFR_RNDN);
        mpfr_sub(r9174, r9172, r9173, MPFR_RNDN);
        ;
        mpfr_set_si(r9176, mpfr_cmp(r9174, r9175) <= 0, MPFR_RNDN);
        mpfr_div(r9177, r9169, r9170, MPFR_RNDN);
        ;
        ;
        mpfr_div(r9180, r9179, r9170, MPFR_RNDN);
        mpfr_sub(r9181, r9178, r9180, MPFR_RNDN);
        mpfr_mul(r9182, r9170, r9170, MPFR_RNDN);
        mpfr_div(r9183, r9181, r9182, MPFR_RNDN);
        mpfr_sub(r9184, r9177, r9183, MPFR_RNDN);
        mpfr_div(r9185, r9171, r9170, MPFR_RNDN);
        mpfr_log(r9186, r9185, MPFR_RNDN);
        if (mpfr_get_si(r9176, MPFR_RNDN)) { mpfr_set(r9187, r9184, MPFR_RNDN); } else { mpfr_set(r9187, r9186, MPFR_RNDN); };
        return mpfr_get_d(r9187, MPFR_RNDN);
}

static mpfr_t r9188, r9189, r9190, r9191, r9192, r9193, r9194, r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204, r9205, r9206;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9188, "1", 10, MPFR_RNDN);
        mpfr_init(r9189);
        mpfr_init(r9190);
        mpfr_init(r9191);
        mpfr_init(r9192);
        mpfr_init(r9193);
        mpfr_init_set_str(r9194, "3.542537891792339e-05", 10, MPFR_RNDN);
        mpfr_init(r9195);
        mpfr_init(r9196);
        mpfr_init_set_str(r9197, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r9198, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9199);
        mpfr_init(r9200);
        mpfr_init(r9201);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
        mpfr_init(r9205);
        mpfr_init(r9206);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r9189, N, MPFR_RNDN);
        mpfr_add(r9190, r9188, r9189, MPFR_RNDN);
        mpfr_log(r9191, r9190, MPFR_RNDN);
        mpfr_log(r9192, r9189, MPFR_RNDN);
        mpfr_sub(r9193, r9191, r9192, MPFR_RNDN);
        ;
        mpfr_set_si(r9195, mpfr_cmp(r9193, r9194) <= 0, MPFR_RNDN);
        mpfr_div(r9196, r9188, r9189, MPFR_RNDN);
        ;
        ;
        mpfr_div(r9199, r9198, r9189, MPFR_RNDN);
        mpfr_sub(r9200, r9197, r9199, MPFR_RNDN);
        mpfr_mul(r9201, r9189, r9189, MPFR_RNDN);
        mpfr_div(r9202, r9200, r9201, MPFR_RNDN);
        mpfr_sub(r9203, r9196, r9202, MPFR_RNDN);
        mpfr_div(r9204, r9190, r9189, MPFR_RNDN);
        mpfr_log(r9205, r9204, MPFR_RNDN);
        if (mpfr_get_si(r9195, MPFR_RNDN)) { mpfr_set(r9206, r9203, MPFR_RNDN); } else { mpfr_set(r9206, r9205, MPFR_RNDN); };
        return mpfr_get_d(r9206, MPFR_RNDN);
}