#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 r25990 = x;
        float r25991 = exp(r25990);
        float r25992 = 1;
        float r25993 = r25991 - r25992;
        float r25994 = r25993 / r25990;
        return r25994;
}

double f_id(double x) {
        double r25995 = x;
        double r25996 = exp(r25995);
        double r25997 = 1;
        double r25998 = r25996 - r25997;
        double r25999 = r25998 / r25995;
        return r25999;
}


double f_of(float x) {
        float r26000 = 1/6;
        float r26001 = x;
        float r26002 = 2;
        float r26003 = pow(r26001, r26002);
        float r26004 = r26000 * r26003;
        float r26005 = 1;
        float r26006 = 1/2;
        float r26007 = r26006 * r26001;
        float r26008 = r26005 + r26007;
        float r26009 = r26004 + r26008;
        float r26010 = 1.005934458562277;
        bool r26011 = r26009 <= r26010;
        float r26012 = exp(r26001);
        float r26013 = r26012 / r26001;
        float r26014 = r26005 / r26001;
        float r26015 = r26013 - r26014;
        float r26016 = r26011 ? r26009 : r26015;
        return r26016;
}

double f_od(double x) {
        double r26017 = 1/6;
        double r26018 = x;
        double r26019 = 2;
        double r26020 = pow(r26018, r26019);
        double r26021 = r26017 * r26020;
        double r26022 = 1;
        double r26023 = 1/2;
        double r26024 = r26023 * r26018;
        double r26025 = r26022 + r26024;
        double r26026 = r26021 + r26025;
        double r26027 = 1.005934458562277;
        bool r26028 = r26026 <= r26027;
        double r26029 = exp(r26018);
        double r26030 = r26029 / r26018;
        double r26031 = r26022 / r26018;
        double r26032 = r26030 - r26031;
        double r26033 = r26028 ? r26026 : r26032;
        return r26033;
}

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 r26034, r26035, r26036, r26037, r26038;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26034);
        mpfr_init(r26035);
        mpfr_init_set_str(r26036, "1", 10, MPFR_RNDN);
        mpfr_init(r26037);
        mpfr_init(r26038);
}

double f_im(double x) {
        mpfr_set_d(r26034, x, MPFR_RNDN);
        mpfr_exp(r26035, r26034, MPFR_RNDN);
        ;
        mpfr_sub(r26037, r26035, r26036, MPFR_RNDN);
        mpfr_div(r26038, r26037, r26034, MPFR_RNDN);
        return mpfr_get_d(r26038, MPFR_RNDN);
}

static mpfr_t r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26039, "1/6", 10, MPFR_RNDN);
        mpfr_init(r26040);
        mpfr_init_set_str(r26041, "2", 10, MPFR_RNDN);
        mpfr_init(r26042);
        mpfr_init(r26043);
        mpfr_init_set_str(r26044, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26045, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26046);
        mpfr_init(r26047);
        mpfr_init(r26048);
        mpfr_init_set_str(r26049, "1.005934458562277", 10, MPFR_RNDN);
        mpfr_init(r26050);
        mpfr_init(r26051);
        mpfr_init(r26052);
        mpfr_init(r26053);
        mpfr_init(r26054);
        mpfr_init(r26055);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r26040, x, MPFR_RNDN);
        ;
        mpfr_pow(r26042, r26040, r26041, MPFR_RNDN);
        mpfr_mul(r26043, r26039, r26042, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26046, r26045, r26040, MPFR_RNDN);
        mpfr_add(r26047, r26044, r26046, MPFR_RNDN);
        mpfr_add(r26048, r26043, r26047, MPFR_RNDN);
        ;
        mpfr_set_si(r26050, mpfr_cmp(r26048, r26049) <= 0, MPFR_RNDN);
        mpfr_exp(r26051, r26040, MPFR_RNDN);
        mpfr_div(r26052, r26051, r26040, MPFR_RNDN);
        mpfr_div(r26053, r26044, r26040, MPFR_RNDN);
        mpfr_sub(r26054, r26052, r26053, MPFR_RNDN);
        if (mpfr_get_si(r26050, MPFR_RNDN)) { mpfr_set(r26055, r26048, MPFR_RNDN); } else { mpfr_set(r26055, r26054, MPFR_RNDN); };
        return mpfr_get_d(r26055, MPFR_RNDN);
}

static mpfr_t r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26056, "1/6", 10, MPFR_RNDN);
        mpfr_init(r26057);
        mpfr_init_set_str(r26058, "2", 10, MPFR_RNDN);
        mpfr_init(r26059);
        mpfr_init(r26060);
        mpfr_init_set_str(r26061, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26062, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26063);
        mpfr_init(r26064);
        mpfr_init(r26065);
        mpfr_init_set_str(r26066, "1.005934458562277", 10, MPFR_RNDN);
        mpfr_init(r26067);
        mpfr_init(r26068);
        mpfr_init(r26069);
        mpfr_init(r26070);
        mpfr_init(r26071);
        mpfr_init(r26072);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r26057, x, MPFR_RNDN);
        ;
        mpfr_pow(r26059, r26057, r26058, MPFR_RNDN);
        mpfr_mul(r26060, r26056, r26059, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r26063, r26062, r26057, MPFR_RNDN);
        mpfr_add(r26064, r26061, r26063, MPFR_RNDN);
        mpfr_add(r26065, r26060, r26064, MPFR_RNDN);
        ;
        mpfr_set_si(r26067, mpfr_cmp(r26065, r26066) <= 0, MPFR_RNDN);
        mpfr_exp(r26068, r26057, MPFR_RNDN);
        mpfr_div(r26069, r26068, r26057, MPFR_RNDN);
        mpfr_div(r26070, r26061, r26057, MPFR_RNDN);
        mpfr_sub(r26071, r26069, r26070, MPFR_RNDN);
        if (mpfr_get_si(r26067, MPFR_RNDN)) { mpfr_set(r26072, r26065, MPFR_RNDN); } else { mpfr_set(r26072, r26071, MPFR_RNDN); };
        return mpfr_get_d(r26072, MPFR_RNDN);
}