#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 r39109 = 1;
        float r39110 = x;
        float r39111 = r39109 + r39110;
        float r39112 = log(r39111);
        return r39112;
}

double f_id(double x) {
        double r39113 = 1;
        double r39114 = x;
        double r39115 = r39113 + r39114;
        double r39116 = log(r39115);
        return r39116;
}


double f_of(float x) {
        float r39117 = x;
        float r39118 = 0.00010417168329262867;
        bool r39119 = r39117 <= r39118;
        float r39120 = r39117 * r39117;
        float r39121 = 1/3;
        float r39122 = r39117 * r39121;
        float r39123 = 1/2;
        float r39124 = r39122 - r39123;
        float r39125 = r39120 * r39124;
        float r39126 = r39117 + r39125;
        float r39127 = 1;
        float r39128 = r39127 + r39117;
        float r39129 = log(r39128);
        float r39130 = r39119 ? r39126 : r39129;
        return r39130;
}

double f_od(double x) {
        double r39131 = x;
        double r39132 = 0.00010417168329262867;
        bool r39133 = r39131 <= r39132;
        double r39134 = r39131 * r39131;
        double r39135 = 1/3;
        double r39136 = r39131 * r39135;
        double r39137 = 1/2;
        double r39138 = r39136 - r39137;
        double r39139 = r39134 * r39138;
        double r39140 = r39131 + r39139;
        double r39141 = 1;
        double r39142 = r39141 + r39131;
        double r39143 = log(r39142);
        double r39144 = r39133 ? r39140 : r39143;
        return r39144;
}

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 r39145, r39146, r39147, r39148;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r39145, "1", 10, MPFR_RNDN);
        mpfr_init(r39146);
        mpfr_init(r39147);
        mpfr_init(r39148);
}

double f_im(double x) {
        ;
        mpfr_set_d(r39146, x, MPFR_RNDN);
        mpfr_add(r39147, r39145, r39146, MPFR_RNDN);
        mpfr_log(r39148, r39147, MPFR_RNDN);
        return mpfr_get_d(r39148, MPFR_RNDN);
}

static mpfr_t r39149, r39150, r39151, r39152, r39153, r39154, r39155, r39156, r39157, r39158, r39159, r39160, r39161, r39162;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r39149);
        mpfr_init_set_str(r39150, "0.00010417168329262867", 10, MPFR_RNDN);
        mpfr_init(r39151);
        mpfr_init(r39152);
        mpfr_init_set_str(r39153, "1/3", 10, MPFR_RNDN);
        mpfr_init(r39154);
        mpfr_init_set_str(r39155, "1/2", 10, MPFR_RNDN);
        mpfr_init(r39156);
        mpfr_init(r39157);
        mpfr_init(r39158);
        mpfr_init_set_str(r39159, "1", 10, MPFR_RNDN);
        mpfr_init(r39160);
        mpfr_init(r39161);
        mpfr_init(r39162);
}

double f_fm(double x) {
        mpfr_set_d(r39149, x, MPFR_RNDN);
        ;
        mpfr_set_si(r39151, mpfr_cmp(r39149, r39150) <= 0, MPFR_RNDN);
        mpfr_mul(r39152, r39149, r39149, MPFR_RNDN);
        ;
        mpfr_mul(r39154, r39149, r39153, MPFR_RNDN);
        ;
        mpfr_sub(r39156, r39154, r39155, MPFR_RNDN);
        mpfr_mul(r39157, r39152, r39156, MPFR_RNDN);
        mpfr_add(r39158, r39149, r39157, MPFR_RNDN);
        ;
        mpfr_add(r39160, r39159, r39149, MPFR_RNDN);
        mpfr_log(r39161, r39160, MPFR_RNDN);
        if (mpfr_get_si(r39151, MPFR_RNDN)) { mpfr_set(r39162, r39158, MPFR_RNDN); } else { mpfr_set(r39162, r39161, MPFR_RNDN); };
        return mpfr_get_d(r39162, MPFR_RNDN);
}

static mpfr_t r39163, r39164, r39165, r39166, r39167, r39168, r39169, r39170, r39171, r39172, r39173, r39174, r39175, r39176;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r39163);
        mpfr_init_set_str(r39164, "0.00010417168329262867", 10, MPFR_RNDN);
        mpfr_init(r39165);
        mpfr_init(r39166);
        mpfr_init_set_str(r39167, "1/3", 10, MPFR_RNDN);
        mpfr_init(r39168);
        mpfr_init_set_str(r39169, "1/2", 10, MPFR_RNDN);
        mpfr_init(r39170);
        mpfr_init(r39171);
        mpfr_init(r39172);
        mpfr_init_set_str(r39173, "1", 10, MPFR_RNDN);
        mpfr_init(r39174);
        mpfr_init(r39175);
        mpfr_init(r39176);
}

double f_dm(double x) {
        mpfr_set_d(r39163, x, MPFR_RNDN);
        ;
        mpfr_set_si(r39165, mpfr_cmp(r39163, r39164) <= 0, MPFR_RNDN);
        mpfr_mul(r39166, r39163, r39163, MPFR_RNDN);
        ;
        mpfr_mul(r39168, r39163, r39167, MPFR_RNDN);
        ;
        mpfr_sub(r39170, r39168, r39169, MPFR_RNDN);
        mpfr_mul(r39171, r39166, r39170, MPFR_RNDN);
        mpfr_add(r39172, r39163, r39171, MPFR_RNDN);
        ;
        mpfr_add(r39174, r39173, r39163, MPFR_RNDN);
        mpfr_log(r39175, r39174, MPFR_RNDN);
        if (mpfr_get_si(r39165, MPFR_RNDN)) { mpfr_set(r39176, r39172, MPFR_RNDN); } else { mpfr_set(r39176, r39175, MPFR_RNDN); };
        return mpfr_get_d(r39176, MPFR_RNDN);
}