#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 r9089 = N;
        float r9090 = 1.0f;
        float r9091 = r9089 + r9090;
        float r9092 = log(r9091);
        float r9093 = log(r9089);
        float r9094 = r9092 - r9093;
        return r9094;
}

double f_id(double N) {
        double r9095 = N;
        double r9096 = 1.0;
        double r9097 = r9095 + r9096;
        double r9098 = log(r9097);
        double r9099 = log(r9095);
        double r9100 = r9098 - r9099;
        return r9100;
}


double f_of(float N) {
        float r9101 = 1.0f;
        float r9102 = N;
        float r9103 = r9101 + r9102;
        float r9104 = log(r9103);
        float r9105 = log(r9102);
        float r9106 = r9104 - r9105;
        float r9107 = 3.542537891792339e-05f;
        bool r9108 = r9106 <= r9107;
        float r9109 = 0.3333333333333333f;
        float r9110 = r9109 / r9102;
        float r9111 = 0.5f;
        float r9112 = r9110 - r9111;
        float r9113 = r9102 * r9102;
        float r9114 = r9101 / r9113;
        float r9115 = r9101 / r9102;
        float r9116 = fma(r9112, r9114, r9115);
        float r9117 = r9103 / r9102;
        float r9118 = log(r9117);
        float r9119 = r9108 ? r9116 : r9118;
        return r9119;
}

double f_od(double N) {
        double r9120 = 1.0;
        double r9121 = N;
        double r9122 = r9120 + r9121;
        double r9123 = log(r9122);
        double r9124 = log(r9121);
        double r9125 = r9123 - r9124;
        double r9126 = 3.542537891792339e-05;
        bool r9127 = r9125 <= r9126;
        double r9128 = 0.3333333333333333;
        double r9129 = r9128 / r9121;
        double r9130 = 0.5;
        double r9131 = r9129 - r9130;
        double r9132 = r9121 * r9121;
        double r9133 = r9120 / r9132;
        double r9134 = r9120 / r9121;
        double r9135 = fma(r9131, r9133, r9134);
        double r9136 = r9122 / r9121;
        double r9137 = log(r9136);
        double r9138 = r9127 ? r9135 : r9137;
        return r9138;
}

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 r9139, r9140, r9141, r9142, r9143, r9144;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9139);
        mpfr_init_set_str(r9140, "1", 10, MPFR_RNDN);
        mpfr_init(r9141);
        mpfr_init(r9142);
        mpfr_init(r9143);
        mpfr_init(r9144);
}

double f_im(double N) {
        mpfr_set_d(r9139, N, MPFR_RNDN);
        ;
        mpfr_add(r9141, r9139, r9140, MPFR_RNDN);
        mpfr_log(r9142, r9141, MPFR_RNDN);
        mpfr_log(r9143, r9139, MPFR_RNDN);
        mpfr_sub(r9144, r9142, r9143, MPFR_RNDN);
        return mpfr_get_d(r9144, MPFR_RNDN);
}

static mpfr_t r9145, r9146, r9147, r9148, r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159, r9160, r9161, r9162, r9163;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9145, "1", 10, MPFR_RNDN);
        mpfr_init(r9146);
        mpfr_init(r9147);
        mpfr_init(r9148);
        mpfr_init(r9149);
        mpfr_init(r9150);
        mpfr_init_set_str(r9151, "3.542537891792339e-05", 10, MPFR_RNDN);
        mpfr_init(r9152);
        mpfr_init_set_str(r9153, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9154);
        mpfr_init_set_str(r9155, "1/2", 10, MPFR_RNDN);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
        mpfr_init(r9159);
        mpfr_init(r9160);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init(r9163);
}

double f_fm(double N) {
        ;
        mpfr_set_d(r9146, N, MPFR_RNDN);
        mpfr_add(r9147, r9145, r9146, MPFR_RNDN);
        mpfr_log(r9148, r9147, MPFR_RNDN);
        mpfr_log(r9149, r9146, MPFR_RNDN);
        mpfr_sub(r9150, r9148, r9149, MPFR_RNDN);
        ;
        mpfr_set_si(r9152, mpfr_cmp(r9150, r9151) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9154, r9153, r9146, MPFR_RNDN);
        ;
        mpfr_sub(r9156, r9154, r9155, MPFR_RNDN);
        mpfr_mul(r9157, r9146, r9146, MPFR_RNDN);
        mpfr_div(r9158, r9145, r9157, MPFR_RNDN);
        mpfr_div(r9159, r9145, r9146, MPFR_RNDN);
        mpfr_fma(r9160, r9156, r9158, r9159, MPFR_RNDN);
        mpfr_div(r9161, r9147, r9146, MPFR_RNDN);
        mpfr_log(r9162, r9161, MPFR_RNDN);
        if (mpfr_get_si(r9152, MPFR_RNDN)) { mpfr_set(r9163, r9160, MPFR_RNDN); } else { mpfr_set(r9163, r9162, MPFR_RNDN); };
        return mpfr_get_d(r9163, MPFR_RNDN);
}

static mpfr_t r9164, r9165, r9166, r9167, r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1616);
        mpfr_init_set_str(r9164, "1", 10, MPFR_RNDN);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init_set_str(r9170, "3.542537891792339e-05", 10, MPFR_RNDN);
        mpfr_init(r9171);
        mpfr_init_set_str(r9172, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9173);
        mpfr_init_set_str(r9174, "1/2", 10, MPFR_RNDN);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
}

double f_dm(double N) {
        ;
        mpfr_set_d(r9165, N, MPFR_RNDN);
        mpfr_add(r9166, r9164, r9165, MPFR_RNDN);
        mpfr_log(r9167, r9166, MPFR_RNDN);
        mpfr_log(r9168, r9165, MPFR_RNDN);
        mpfr_sub(r9169, r9167, r9168, MPFR_RNDN);
        ;
        mpfr_set_si(r9171, mpfr_cmp(r9169, r9170) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9173, r9172, r9165, MPFR_RNDN);
        ;
        mpfr_sub(r9175, r9173, r9174, MPFR_RNDN);
        mpfr_mul(r9176, r9165, r9165, MPFR_RNDN);
        mpfr_div(r9177, r9164, r9176, MPFR_RNDN);
        mpfr_div(r9178, r9164, r9165, MPFR_RNDN);
        mpfr_fma(r9179, r9175, r9177, r9178, MPFR_RNDN);
        mpfr_div(r9180, r9166, r9165, MPFR_RNDN);
        mpfr_log(r9181, r9180, MPFR_RNDN);
        if (mpfr_get_si(r9171, MPFR_RNDN)) { mpfr_set(r9182, r9179, MPFR_RNDN); } else { mpfr_set(r9182, r9181, MPFR_RNDN); };
        return mpfr_get_d(r9182, MPFR_RNDN);
}