#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 r20344 = N;
        float r20345 = 1;
        float r20346 = r20344 + r20345;
        float r20347 = log(r20346);
        float r20348 = log(r20344);
        float r20349 = r20347 - r20348;
        return r20349;
}

double f_id(double N) {
        double r20350 = N;
        double r20351 = 1;
        double r20352 = r20350 + r20351;
        double r20353 = log(r20352);
        double r20354 = log(r20350);
        double r20355 = r20353 - r20354;
        return r20355;
}


double f_of(float N) {
        float r20356 = N;
        float r20357 = 7752.741142320907;
        bool r20358 = r20356 <= r20357;
        float r20359 = 1;
        float r20360 = r20356 + r20359;
        float r20361 = r20360 / r20356;
        float r20362 = log(r20361);
        float r20363 = r20359 / r20356;
        float r20364 = 1/3;
        float r20365 = r20364 / r20356;
        float r20366 = 1/2;
        float r20367 = r20365 - r20366;
        float r20368 = r20356 * r20356;
        float r20369 = r20367 / r20368;
        float r20370 = r20363 + r20369;
        float r20371 = r20358 ? r20362 : r20370;
        return r20371;
}

double f_od(double N) {
        double r20372 = N;
        double r20373 = 7752.741142320907;
        bool r20374 = r20372 <= r20373;
        double r20375 = 1;
        double r20376 = r20372 + r20375;
        double r20377 = r20376 / r20372;
        double r20378 = log(r20377);
        double r20379 = r20375 / r20372;
        double r20380 = 1/3;
        double r20381 = r20380 / r20372;
        double r20382 = 1/2;
        double r20383 = r20381 - r20382;
        double r20384 = r20372 * r20372;
        double r20385 = r20383 / r20384;
        double r20386 = r20379 + r20385;
        double r20387 = r20374 ? r20378 : r20386;
        return r20387;
}

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 r20388, r20389, r20390, r20391, r20392, r20393;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20388);
        mpfr_init_set_str(r20389, "1", 10, MPFR_RNDN);
        mpfr_init(r20390);
        mpfr_init(r20391);
        mpfr_init(r20392);
        mpfr_init(r20393);
}

double f_im(double N) {
        mpfr_set_d(r20388, N, MPFR_RNDN);
        ;
        mpfr_add(r20390, r20388, r20389, MPFR_RNDN);
        mpfr_log(r20391, r20390, MPFR_RNDN);
        mpfr_log(r20392, r20388, MPFR_RNDN);
        mpfr_sub(r20393, r20391, r20392, MPFR_RNDN);
        return mpfr_get_d(r20393, MPFR_RNDN);
}

static mpfr_t r20394, r20395, r20396, r20397, r20398, r20399, r20400, r20401, r20402, r20403, r20404, r20405, r20406, r20407, r20408, r20409;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20394);
        mpfr_init_set_str(r20395, "7752.741142320907", 10, MPFR_RNDN);
        mpfr_init(r20396);
        mpfr_init_set_str(r20397, "1", 10, MPFR_RNDN);
        mpfr_init(r20398);
        mpfr_init(r20399);
        mpfr_init(r20400);
        mpfr_init(r20401);
        mpfr_init_set_str(r20402, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20403);
        mpfr_init_set_str(r20404, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20405);
        mpfr_init(r20406);
        mpfr_init(r20407);
        mpfr_init(r20408);
        mpfr_init(r20409);
}

double f_fm(double N) {
        mpfr_set_d(r20394, N, MPFR_RNDN);
        ;
        mpfr_set_si(r20396, mpfr_cmp(r20394, r20395) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r20398, r20394, r20397, MPFR_RNDN);
        mpfr_div(r20399, r20398, r20394, MPFR_RNDN);
        mpfr_log(r20400, r20399, MPFR_RNDN);
        mpfr_div(r20401, r20397, r20394, MPFR_RNDN);
        ;
        mpfr_div(r20403, r20402, r20394, MPFR_RNDN);
        ;
        mpfr_sub(r20405, r20403, r20404, MPFR_RNDN);
        mpfr_mul(r20406, r20394, r20394, MPFR_RNDN);
        mpfr_div(r20407, r20405, r20406, MPFR_RNDN);
        mpfr_add(r20408, r20401, r20407, MPFR_RNDN);
        if (mpfr_get_si(r20396, MPFR_RNDN)) { mpfr_set(r20409, r20400, MPFR_RNDN); } else { mpfr_set(r20409, r20408, MPFR_RNDN); };
        return mpfr_get_d(r20409, MPFR_RNDN);
}

static mpfr_t r20410, r20411, r20412, r20413, r20414, r20415, r20416, r20417, r20418, r20419, r20420, r20421, r20422, r20423, r20424, r20425;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20410);
        mpfr_init_set_str(r20411, "7752.741142320907", 10, MPFR_RNDN);
        mpfr_init(r20412);
        mpfr_init_set_str(r20413, "1", 10, MPFR_RNDN);
        mpfr_init(r20414);
        mpfr_init(r20415);
        mpfr_init(r20416);
        mpfr_init(r20417);
        mpfr_init_set_str(r20418, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20419);
        mpfr_init_set_str(r20420, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20421);
        mpfr_init(r20422);
        mpfr_init(r20423);
        mpfr_init(r20424);
        mpfr_init(r20425);
}

double f_dm(double N) {
        mpfr_set_d(r20410, N, MPFR_RNDN);
        ;
        mpfr_set_si(r20412, mpfr_cmp(r20410, r20411) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r20414, r20410, r20413, MPFR_RNDN);
        mpfr_div(r20415, r20414, r20410, MPFR_RNDN);
        mpfr_log(r20416, r20415, MPFR_RNDN);
        mpfr_div(r20417, r20413, r20410, MPFR_RNDN);
        ;
        mpfr_div(r20419, r20418, r20410, MPFR_RNDN);
        ;
        mpfr_sub(r20421, r20419, r20420, MPFR_RNDN);
        mpfr_mul(r20422, r20410, r20410, MPFR_RNDN);
        mpfr_div(r20423, r20421, r20422, MPFR_RNDN);
        mpfr_add(r20424, r20417, r20423, MPFR_RNDN);
        if (mpfr_get_si(r20412, MPFR_RNDN)) { mpfr_set(r20425, r20416, MPFR_RNDN); } else { mpfr_set(r20425, r20424, MPFR_RNDN); };
        return mpfr_get_d(r20425, MPFR_RNDN);
}