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

char *name = "NMSE problem 3.3.6";

double f_if(float N) {
        float r26244 = N;
        float r26245 = 1.0f;
        float r26246 = r26244 + r26245;
        float r26247 = log(r26246);
        float r26248 = log(r26244);
        float r26249 = r26247 - r26248;
        return r26249;
}

double f_id(double N) {
        double r26250 = N;
        double r26251 = 1.0;
        double r26252 = r26250 + r26251;
        double r26253 = log(r26252);
        double r26254 = log(r26250);
        double r26255 = r26253 - r26254;
        return r26255;
}


double f_of(float N) {
        float r26256 = N;
        float r26257 = 591618366256079.6f;
        bool r26258 = r26256 <= r26257;
        float r26259 = 1.0f;
        float r26260 = r26256 + r26259;
        float r26261 = r26260 / r26256;
        float r26262 = log(r26261);
        float r26263 = 0.3333333333333333f;
        float r26264 = r26263 / r26256;
        float r26265 = 0.5f;
        float r26266 = r26264 - r26265;
        float r26267 = r26256 * r26256;
        float r26268 = r26266 / r26267;
        float r26269 = r26259 / r26256;
        float r26270 = r26268 + r26269;
        float r26271 = r26258 ? r26262 : r26270;
        return r26271;
}

double f_od(double N) {
        double r26272 = N;
        double r26273 = 591618366256079.6;
        bool r26274 = r26272 <= r26273;
        double r26275 = 1.0;
        double r26276 = r26272 + r26275;
        double r26277 = r26276 / r26272;
        double r26278 = log(r26277);
        double r26279 = 0.3333333333333333;
        double r26280 = r26279 / r26272;
        double r26281 = 0.5;
        double r26282 = r26280 - r26281;
        double r26283 = r26272 * r26272;
        double r26284 = r26282 / r26283;
        double r26285 = r26275 / r26272;
        double r26286 = r26284 + r26285;
        double r26287 = r26274 ? r26278 : r26286;
        return r26287;
}

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 r26288, r26289, r26290, r26291, r26292, r26293;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r26288);
        mpfr_init_set_str(r26289, "1", 10, MPFR_RNDN);
        mpfr_init(r26290);
        mpfr_init(r26291);
        mpfr_init(r26292);
        mpfr_init(r26293);
}

double f_im(double N) {
        mpfr_set_d(r26288, N, MPFR_RNDN);
        ;
        mpfr_add(r26290, r26288, r26289, MPFR_RNDN);
        mpfr_log(r26291, r26290, MPFR_RNDN);
        mpfr_log(r26292, r26288, MPFR_RNDN);
        mpfr_sub(r26293, r26291, r26292, MPFR_RNDN);
        return mpfr_get_d(r26293, MPFR_RNDN);
}

static mpfr_t r26294, r26295, r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r26294);
        mpfr_init_set_str(r26295, "591618366256079.6", 10, MPFR_RNDN);
        mpfr_init(r26296);
        mpfr_init_set_str(r26297, "1", 10, MPFR_RNDN);
        mpfr_init(r26298);
        mpfr_init(r26299);
        mpfr_init(r26300);
        mpfr_init_set_str(r26301, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26302);
        mpfr_init_set_str(r26303, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
}

double f_fm(double N) {
        mpfr_set_d(r26294, N, MPFR_RNDN);
        ;
        mpfr_set_si(r26296, mpfr_cmp(r26294, r26295) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r26298, r26294, r26297, MPFR_RNDN);
        mpfr_div(r26299, r26298, r26294, MPFR_RNDN);
        mpfr_log(r26300, r26299, MPFR_RNDN);
        ;
        mpfr_div(r26302, r26301, r26294, MPFR_RNDN);
        ;
        mpfr_sub(r26304, r26302, r26303, MPFR_RNDN);
        mpfr_sqr(r26305, r26294, MPFR_RNDN);
        mpfr_div(r26306, r26304, r26305, MPFR_RNDN);
        mpfr_div(r26307, r26297, r26294, MPFR_RNDN);
        mpfr_add(r26308, r26306, r26307, MPFR_RNDN);
        if (mpfr_get_si(r26296, MPFR_RNDN)) { mpfr_set(r26309, r26300, MPFR_RNDN); } else { mpfr_set(r26309, r26308, MPFR_RNDN); };
        return mpfr_get_d(r26309, MPFR_RNDN);
}

static mpfr_t r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r26310);
        mpfr_init_set_str(r26311, "591618366256079.6", 10, MPFR_RNDN);
        mpfr_init(r26312);
        mpfr_init_set_str(r26313, "1", 10, MPFR_RNDN);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init_set_str(r26317, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26318);
        mpfr_init_set_str(r26319, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26320);
        mpfr_init(r26321);
        mpfr_init(r26322);
        mpfr_init(r26323);
        mpfr_init(r26324);
        mpfr_init(r26325);
}

double f_dm(double N) {
        mpfr_set_d(r26310, N, MPFR_RNDN);
        ;
        mpfr_set_si(r26312, mpfr_cmp(r26310, r26311) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r26314, r26310, r26313, MPFR_RNDN);
        mpfr_div(r26315, r26314, r26310, MPFR_RNDN);
        mpfr_log(r26316, r26315, MPFR_RNDN);
        ;
        mpfr_div(r26318, r26317, r26310, MPFR_RNDN);
        ;
        mpfr_sub(r26320, r26318, r26319, MPFR_RNDN);
        mpfr_sqr(r26321, r26310, MPFR_RNDN);
        mpfr_div(r26322, r26320, r26321, MPFR_RNDN);
        mpfr_div(r26323, r26313, r26310, MPFR_RNDN);
        mpfr_add(r26324, r26322, r26323, MPFR_RNDN);
        if (mpfr_get_si(r26312, MPFR_RNDN)) { mpfr_set(r26325, r26316, MPFR_RNDN); } else { mpfr_set(r26325, r26324, MPFR_RNDN); };
        return mpfr_get_d(r26325, MPFR_RNDN);
}