#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 r20221 = N;
        float r20222 = 1;
        float r20223 = r20221 + r20222;
        float r20224 = log(r20223);
        float r20225 = log(r20221);
        float r20226 = r20224 - r20225;
        return r20226;
}

double f_id(double N) {
        double r20227 = N;
        double r20228 = 1;
        double r20229 = r20227 + r20228;
        double r20230 = log(r20229);
        double r20231 = log(r20227);
        double r20232 = r20230 - r20231;
        return r20232;
}


double f_of(float N) {
        float r20233 = N;
        float r20234 = 7531.5475611680995;
        bool r20235 = r20233 <= r20234;
        float r20236 = 1;
        float r20237 = r20233 + r20236;
        float r20238 = r20237 / r20233;
        float r20239 = log(r20238);
        float r20240 = r20236 / r20233;
        float r20241 = 1/3;
        float r20242 = r20241 / r20233;
        float r20243 = 1/2;
        float r20244 = r20242 - r20243;
        float r20245 = r20233 * r20233;
        float r20246 = r20244 / r20245;
        float r20247 = r20240 + r20246;
        float r20248 = r20235 ? r20239 : r20247;
        return r20248;
}

double f_od(double N) {
        double r20249 = N;
        double r20250 = 7531.5475611680995;
        bool r20251 = r20249 <= r20250;
        double r20252 = 1;
        double r20253 = r20249 + r20252;
        double r20254 = r20253 / r20249;
        double r20255 = log(r20254);
        double r20256 = r20252 / r20249;
        double r20257 = 1/3;
        double r20258 = r20257 / r20249;
        double r20259 = 1/2;
        double r20260 = r20258 - r20259;
        double r20261 = r20249 * r20249;
        double r20262 = r20260 / r20261;
        double r20263 = r20256 + r20262;
        double r20264 = r20251 ? r20255 : r20263;
        return r20264;
}

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 r20265, r20266, r20267, r20268, r20269, r20270;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20265);
        mpfr_init_set_str(r20266, "1", 10, MPFR_RNDN);
        mpfr_init(r20267);
        mpfr_init(r20268);
        mpfr_init(r20269);
        mpfr_init(r20270);
}

double f_im(double N) {
        mpfr_set_d(r20265, N, MPFR_RNDN);
        ;
        mpfr_add(r20267, r20265, r20266, MPFR_RNDN);
        mpfr_log(r20268, r20267, MPFR_RNDN);
        mpfr_log(r20269, r20265, MPFR_RNDN);
        mpfr_sub(r20270, r20268, r20269, MPFR_RNDN);
        return mpfr_get_d(r20270, MPFR_RNDN);
}

static mpfr_t r20271, r20272, r20273, r20274, r20275, r20276, r20277, r20278, r20279, r20280, r20281, r20282, r20283, r20284, r20285, r20286;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20271);
        mpfr_init_set_str(r20272, "7531.5475611680995", 10, MPFR_RNDN);
        mpfr_init(r20273);
        mpfr_init_set_str(r20274, "1", 10, MPFR_RNDN);
        mpfr_init(r20275);
        mpfr_init(r20276);
        mpfr_init(r20277);
        mpfr_init(r20278);
        mpfr_init_set_str(r20279, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20280);
        mpfr_init_set_str(r20281, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20282);
        mpfr_init(r20283);
        mpfr_init(r20284);
        mpfr_init(r20285);
        mpfr_init(r20286);
}

double f_fm(double N) {
        mpfr_set_d(r20271, N, MPFR_RNDN);
        ;
        mpfr_set_si(r20273, mpfr_cmp(r20271, r20272) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r20275, r20271, r20274, MPFR_RNDN);
        mpfr_div(r20276, r20275, r20271, MPFR_RNDN);
        mpfr_log(r20277, r20276, MPFR_RNDN);
        mpfr_div(r20278, r20274, r20271, MPFR_RNDN);
        ;
        mpfr_div(r20280, r20279, r20271, MPFR_RNDN);
        ;
        mpfr_sub(r20282, r20280, r20281, MPFR_RNDN);
        mpfr_mul(r20283, r20271, r20271, MPFR_RNDN);
        mpfr_div(r20284, r20282, r20283, MPFR_RNDN);
        mpfr_add(r20285, r20278, r20284, MPFR_RNDN);
        if (mpfr_get_si(r20273, MPFR_RNDN)) { mpfr_set(r20286, r20277, MPFR_RNDN); } else { mpfr_set(r20286, r20285, MPFR_RNDN); };
        return mpfr_get_d(r20286, MPFR_RNDN);
}

static mpfr_t r20287, r20288, r20289, r20290, r20291, r20292, r20293, r20294, r20295, r20296, r20297, r20298, r20299, r20300, r20301, r20302;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20287);
        mpfr_init_set_str(r20288, "7531.5475611680995", 10, MPFR_RNDN);
        mpfr_init(r20289);
        mpfr_init_set_str(r20290, "1", 10, MPFR_RNDN);
        mpfr_init(r20291);
        mpfr_init(r20292);
        mpfr_init(r20293);
        mpfr_init(r20294);
        mpfr_init_set_str(r20295, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20296);
        mpfr_init_set_str(r20297, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20298);
        mpfr_init(r20299);
        mpfr_init(r20300);
        mpfr_init(r20301);
        mpfr_init(r20302);
}

double f_dm(double N) {
        mpfr_set_d(r20287, N, MPFR_RNDN);
        ;
        mpfr_set_si(r20289, mpfr_cmp(r20287, r20288) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r20291, r20287, r20290, MPFR_RNDN);
        mpfr_div(r20292, r20291, r20287, MPFR_RNDN);
        mpfr_log(r20293, r20292, MPFR_RNDN);
        mpfr_div(r20294, r20290, r20287, MPFR_RNDN);
        ;
        mpfr_div(r20296, r20295, r20287, MPFR_RNDN);
        ;
        mpfr_sub(r20298, r20296, r20297, MPFR_RNDN);
        mpfr_mul(r20299, r20287, r20287, MPFR_RNDN);
        mpfr_div(r20300, r20298, r20299, MPFR_RNDN);
        mpfr_add(r20301, r20294, r20300, MPFR_RNDN);
        if (mpfr_get_si(r20289, MPFR_RNDN)) { mpfr_set(r20302, r20293, MPFR_RNDN); } else { mpfr_set(r20302, r20301, MPFR_RNDN); };
        return mpfr_get_d(r20302, MPFR_RNDN);
}