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

char *name = "ln(1 + x)";

double f_if(float x) {
        float r25407 = 1;
        float r25408 = x;
        float r25409 = r25407 + r25408;
        float r25410 = log(r25409);
        return r25410;
}

double f_id(double x) {
        double r25411 = 1;
        double r25412 = x;
        double r25413 = r25411 + r25412;
        double r25414 = log(r25413);
        return r25414;
}


double f_of(float x) {
        float r25415 = 1;
        float r25416 = x;
        float r25417 = r25415 + r25416;
        float r25418 = log(r25417);
        float r25419 = 5.063949259907442e-12;
        bool r25420 = r25418 <= r25419;
        float r25421 = 1/2;
        float r25422 = 1/3;
        float r25423 = r25416 * r25422;
        float r25424 = r25421 - r25423;
        float r25425 = r25416 * r25416;
        float r25426 = r25424 * r25425;
        float r25427 = r25416 - r25426;
        float r25428 = r25420 ? r25427 : r25418;
        return r25428;
}

double f_od(double x) {
        double r25429 = 1;
        double r25430 = x;
        double r25431 = r25429 + r25430;
        double r25432 = log(r25431);
        double r25433 = 5.063949259907442e-12;
        bool r25434 = r25432 <= r25433;
        double r25435 = 1/2;
        double r25436 = 1/3;
        double r25437 = r25430 * r25436;
        double r25438 = r25435 - r25437;
        double r25439 = r25430 * r25430;
        double r25440 = r25438 * r25439;
        double r25441 = r25430 - r25440;
        double r25442 = r25434 ? r25441 : r25432;
        return r25442;
}

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 r25443, r25444, r25445, r25446;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r25443, "1", 10, MPFR_RNDN);
        mpfr_init(r25444);
        mpfr_init(r25445);
        mpfr_init(r25446);
}

double f_im(double x) {
        ;
        mpfr_set_d(r25444, x, MPFR_RNDN);
        mpfr_add(r25445, r25443, r25444, MPFR_RNDN);
        mpfr_log(r25446, r25445, MPFR_RNDN);
        return mpfr_get_d(r25446, MPFR_RNDN);
}

static mpfr_t r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r25447, "1", 10, MPFR_RNDN);
        mpfr_init(r25448);
        mpfr_init(r25449);
        mpfr_init(r25450);
        mpfr_init_set_str(r25451, "5.063949259907442e-12", 10, MPFR_RNDN);
        mpfr_init(r25452);
        mpfr_init_set_str(r25453, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r25454, "1/3", 10, MPFR_RNDN);
        mpfr_init(r25455);
        mpfr_init(r25456);
        mpfr_init(r25457);
        mpfr_init(r25458);
        mpfr_init(r25459);
        mpfr_init(r25460);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r25448, x, MPFR_RNDN);
        mpfr_add(r25449, r25447, r25448, MPFR_RNDN);
        mpfr_log(r25450, r25449, MPFR_RNDN);
        ;
        mpfr_set_si(r25452, mpfr_cmp(r25450, r25451) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r25455, r25448, r25454, MPFR_RNDN);
        mpfr_sub(r25456, r25453, r25455, MPFR_RNDN);
        mpfr_mul(r25457, r25448, r25448, MPFR_RNDN);
        mpfr_mul(r25458, r25456, r25457, MPFR_RNDN);
        mpfr_sub(r25459, r25448, r25458, MPFR_RNDN);
        if (mpfr_get_si(r25452, MPFR_RNDN)) { mpfr_set(r25460, r25459, MPFR_RNDN); } else { mpfr_set(r25460, r25450, MPFR_RNDN); };
        return mpfr_get_d(r25460, MPFR_RNDN);
}

static mpfr_t r25461, r25462, r25463, r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r25461, "1", 10, MPFR_RNDN);
        mpfr_init(r25462);
        mpfr_init(r25463);
        mpfr_init(r25464);
        mpfr_init_set_str(r25465, "5.063949259907442e-12", 10, MPFR_RNDN);
        mpfr_init(r25466);
        mpfr_init_set_str(r25467, "1/2", 10, MPFR_RNDN);
        mpfr_init_set_str(r25468, "1/3", 10, MPFR_RNDN);
        mpfr_init(r25469);
        mpfr_init(r25470);
        mpfr_init(r25471);
        mpfr_init(r25472);
        mpfr_init(r25473);
        mpfr_init(r25474);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r25462, x, MPFR_RNDN);
        mpfr_add(r25463, r25461, r25462, MPFR_RNDN);
        mpfr_log(r25464, r25463, MPFR_RNDN);
        ;
        mpfr_set_si(r25466, mpfr_cmp(r25464, r25465) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r25469, r25462, r25468, MPFR_RNDN);
        mpfr_sub(r25470, r25467, r25469, MPFR_RNDN);
        mpfr_mul(r25471, r25462, r25462, MPFR_RNDN);
        mpfr_mul(r25472, r25470, r25471, MPFR_RNDN);
        mpfr_sub(r25473, r25462, r25472, MPFR_RNDN);
        if (mpfr_get_si(r25466, MPFR_RNDN)) { mpfr_set(r25474, r25473, MPFR_RNDN); } else { mpfr_set(r25474, r25464, MPFR_RNDN); };
        return mpfr_get_d(r25474, MPFR_RNDN);
}