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

char *name = "Kahan's exp quotient";

double f_if(float x) {
        float r9101 = x;
        float r9102 = exp(r9101);
        float r9103 = 1.0f;
        float r9104 = r9102 - r9103;
        float r9105 = r9104 / r9101;
        return r9105;
}

double f_id(double x) {
        double r9106 = x;
        double r9107 = exp(r9106);
        double r9108 = 1.0;
        double r9109 = r9107 - r9108;
        double r9110 = r9109 / r9106;
        return r9110;
}


double f_of(float x) {
        float r9111 = x;
        float r9112 = 0.16666666666666666f;
        float r9113 = r9111 * r9112;
        float r9114 = 0.5f;
        float r9115 = r9113 + r9114;
        float r9116 = r9111 * r9111;
        float r9117 = r9115 * r9116;
        float r9118 = r9111 + r9117;
        float r9119 = r9118 / r9111;
        float r9120 = 1.0187894117188008f;
        bool r9121 = r9119 <= r9120;
        float r9122 = exp(r9111);
        float r9123 = r9122 * r9122;
        float r9124 = r9123 * r9123;
        float r9125 = 1.0f;
        float r9126 = r9124 - r9125;
        float r9127 = r9125 + r9123;
        float r9128 = cbrt(r9127);
        float r9129 = r9128 * r9128;
        float r9130 = r9126 / r9129;
        float r9131 = r9130 / r9128;
        float r9132 = r9122 + r9125;
        float r9133 = r9111 * r9132;
        float r9134 = r9131 / r9133;
        float r9135 = r9121 ? r9119 : r9134;
        return r9135;
}

double f_od(double x) {
        double r9136 = x;
        double r9137 = 0.16666666666666666;
        double r9138 = r9136 * r9137;
        double r9139 = 0.5;
        double r9140 = r9138 + r9139;
        double r9141 = r9136 * r9136;
        double r9142 = r9140 * r9141;
        double r9143 = r9136 + r9142;
        double r9144 = r9143 / r9136;
        double r9145 = 1.0187894117188008;
        bool r9146 = r9144 <= r9145;
        double r9147 = exp(r9136);
        double r9148 = r9147 * r9147;
        double r9149 = r9148 * r9148;
        double r9150 = 1.0;
        double r9151 = r9149 - r9150;
        double r9152 = r9150 + r9148;
        double r9153 = cbrt(r9152);
        double r9154 = r9153 * r9153;
        double r9155 = r9151 / r9154;
        double r9156 = r9155 / r9153;
        double r9157 = r9147 + r9150;
        double r9158 = r9136 * r9157;
        double r9159 = r9156 / r9158;
        double r9160 = r9146 ? r9144 : r9159;
        return r9160;
}

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 r9161, r9162, r9163, r9164, r9165;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init_set_str(r9163, "1", 10, MPFR_RNDN);
        mpfr_init(r9164);
        mpfr_init(r9165);
}

double f_im(double x) {
        mpfr_set_d(r9161, x, MPFR_RNDN);
        mpfr_exp(r9162, r9161, MPFR_RNDN);
        ;
        mpfr_sub(r9164, r9162, r9163, MPFR_RNDN);
        mpfr_div(r9165, r9164, r9161, MPFR_RNDN);
        return mpfr_get_d(r9165, MPFR_RNDN);
}

static mpfr_t r9166, r9167, r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185, r9186, r9187, r9188, r9189, r9190;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9166);
        mpfr_init_set_str(r9167, "1/6", 10, MPFR_RNDN);
        mpfr_init(r9168);
        mpfr_init_set_str(r9169, "1/2", 10, MPFR_RNDN);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init_set_str(r9175, "1.0187894117188008", 10, MPFR_RNDN);
        mpfr_init(r9176);
        mpfr_init(r9177);
        mpfr_init(r9178);
        mpfr_init(r9179);
        mpfr_init_set_str(r9180, "1", 10, MPFR_RNDN);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
        mpfr_init(r9189);
        mpfr_init(r9190);
}

double f_fm(double x) {
        mpfr_set_d(r9166, x, MPFR_RNDN);
        ;
        mpfr_mul(r9168, r9166, r9167, MPFR_RNDN);
        ;
        mpfr_add(r9170, r9168, r9169, MPFR_RNDN);
        mpfr_mul(r9171, r9166, r9166, MPFR_RNDN);
        mpfr_mul(r9172, r9170, r9171, MPFR_RNDN);
        mpfr_add(r9173, r9166, r9172, MPFR_RNDN);
        mpfr_div(r9174, r9173, r9166, MPFR_RNDN);
        ;
        mpfr_set_si(r9176, mpfr_cmp(r9174, r9175) <= 0, MPFR_RNDN);
        mpfr_exp(r9177, r9166, MPFR_RNDN);
        mpfr_mul(r9178, r9177, r9177, MPFR_RNDN);
        mpfr_mul(r9179, r9178, r9178, MPFR_RNDN);
        ;
        mpfr_sub(r9181, r9179, r9180, MPFR_RNDN);
        mpfr_add(r9182, r9180, r9178, MPFR_RNDN);
        mpfr_cbrt(r9183, r9182, MPFR_RNDN);
        mpfr_mul(r9184, r9183, r9183, MPFR_RNDN);
        mpfr_div(r9185, r9181, r9184, MPFR_RNDN);
        mpfr_div(r9186, r9185, r9183, MPFR_RNDN);
        mpfr_add(r9187, r9177, r9180, MPFR_RNDN);
        mpfr_mul(r9188, r9166, r9187, MPFR_RNDN);
        mpfr_div(r9189, r9186, r9188, MPFR_RNDN);
        if (mpfr_get_si(r9176, MPFR_RNDN)) { mpfr_set(r9190, r9174, MPFR_RNDN); } else { mpfr_set(r9190, r9189, MPFR_RNDN); };
        return mpfr_get_d(r9190, MPFR_RNDN);
}

static mpfr_t r9191, r9192, r9193, r9194, r9195, r9196, r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204, r9205, r9206, r9207, r9208, r9209, r9210, r9211, r9212, r9213, r9214, r9215;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9191);
        mpfr_init_set_str(r9192, "1/6", 10, MPFR_RNDN);
        mpfr_init(r9193);
        mpfr_init_set_str(r9194, "1/2", 10, MPFR_RNDN);
        mpfr_init(r9195);
        mpfr_init(r9196);
        mpfr_init(r9197);
        mpfr_init(r9198);
        mpfr_init(r9199);
        mpfr_init_set_str(r9200, "1.0187894117188008", 10, MPFR_RNDN);
        mpfr_init(r9201);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
        mpfr_init_set_str(r9205, "1", 10, MPFR_RNDN);
        mpfr_init(r9206);
        mpfr_init(r9207);
        mpfr_init(r9208);
        mpfr_init(r9209);
        mpfr_init(r9210);
        mpfr_init(r9211);
        mpfr_init(r9212);
        mpfr_init(r9213);
        mpfr_init(r9214);
        mpfr_init(r9215);
}

double f_dm(double x) {
        mpfr_set_d(r9191, x, MPFR_RNDN);
        ;
        mpfr_mul(r9193, r9191, r9192, MPFR_RNDN);
        ;
        mpfr_add(r9195, r9193, r9194, MPFR_RNDN);
        mpfr_mul(r9196, r9191, r9191, MPFR_RNDN);
        mpfr_mul(r9197, r9195, r9196, MPFR_RNDN);
        mpfr_add(r9198, r9191, r9197, MPFR_RNDN);
        mpfr_div(r9199, r9198, r9191, MPFR_RNDN);
        ;
        mpfr_set_si(r9201, mpfr_cmp(r9199, r9200) <= 0, MPFR_RNDN);
        mpfr_exp(r9202, r9191, MPFR_RNDN);
        mpfr_mul(r9203, r9202, r9202, MPFR_RNDN);
        mpfr_mul(r9204, r9203, r9203, MPFR_RNDN);
        ;
        mpfr_sub(r9206, r9204, r9205, MPFR_RNDN);
        mpfr_add(r9207, r9205, r9203, MPFR_RNDN);
        mpfr_cbrt(r9208, r9207, MPFR_RNDN);
        mpfr_mul(r9209, r9208, r9208, MPFR_RNDN);
        mpfr_div(r9210, r9206, r9209, MPFR_RNDN);
        mpfr_div(r9211, r9210, r9208, MPFR_RNDN);
        mpfr_add(r9212, r9202, r9205, MPFR_RNDN);
        mpfr_mul(r9213, r9191, r9212, MPFR_RNDN);
        mpfr_div(r9214, r9211, r9213, MPFR_RNDN);
        if (mpfr_get_si(r9201, MPFR_RNDN)) { mpfr_set(r9215, r9199, MPFR_RNDN); } else { mpfr_set(r9215, r9214, MPFR_RNDN); };
        return mpfr_get_d(r9215, MPFR_RNDN);
}