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

char *name = "expq2 (section 3.11)";

double f_if(float x) {
        float r35099 = x;
        float r35100 = exp(r35099);
        float r35101 = 1.0f;
        float r35102 = r35100 - r35101;
        float r35103 = r35100 / r35102;
        return r35103;
}

double f_id(double x) {
        double r35104 = x;
        double r35105 = exp(r35104);
        double r35106 = 1.0;
        double r35107 = r35105 - r35106;
        double r35108 = r35105 / r35107;
        return r35108;
}


double f_of(float x) {
        float r35109 = x;
        float r35110 = -5.458497096453193e-06f;
        bool r35111 = r35109 <= r35110;
        float r35112 = 1.0f;
        float r35113 = -r35109;
        float r35114 = exp(r35113);
        float r35115 = r35112 - r35114;
        float r35116 = r35112 / r35115;
        float r35117 = 15.465634390463869f;
        bool r35118 = r35109 <= r35117;
        float r35119 = r35112 / r35109;
        float r35120 = 0.08333333333333333f;
        float r35121 = r35120 * r35109;
        float r35122 = 0.5f;
        float r35123 = r35121 + r35122;
        float r35124 = r35119 + r35123;
        float r35125 = r35118 ? r35124 : r35116;
        float r35126 = r35111 ? r35116 : r35125;
        return r35126;
}

double f_od(double x) {
        double r35127 = x;
        double r35128 = -5.458497096453193e-06;
        bool r35129 = r35127 <= r35128;
        double r35130 = 1.0;
        double r35131 = -r35127;
        double r35132 = exp(r35131);
        double r35133 = r35130 - r35132;
        double r35134 = r35130 / r35133;
        double r35135 = 15.465634390463869;
        bool r35136 = r35127 <= r35135;
        double r35137 = r35130 / r35127;
        double r35138 = 0.08333333333333333;
        double r35139 = r35138 * r35127;
        double r35140 = 0.5;
        double r35141 = r35139 + r35140;
        double r35142 = r35137 + r35141;
        double r35143 = r35136 ? r35142 : r35134;
        double r35144 = r35129 ? r35134 : r35143;
        return r35144;
}

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 r35145, r35146, r35147, r35148, r35149;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r35145);
        mpfr_init(r35146);
        mpfr_init_set_str(r35147, "1", 10, MPFR_RNDN);
        mpfr_init(r35148);
        mpfr_init(r35149);
}

double f_im(double x) {
        mpfr_set_d(r35145, x, MPFR_RNDN);
        mpfr_exp(r35146, r35145, MPFR_RNDN);
        ;
        mpfr_sub(r35148, r35146, r35147, MPFR_RNDN);
        mpfr_div(r35149, r35146, r35148, MPFR_RNDN);
        return mpfr_get_d(r35149, MPFR_RNDN);
}

static mpfr_t r35150, r35151, r35152, r35153, r35154, r35155, r35156, r35157, r35158, r35159, r35160, r35161, r35162, r35163, r35164, r35165, r35166, r35167;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r35150);
        mpfr_init_set_str(r35151, "-5.458497096453193e-06", 10, MPFR_RNDN);
        mpfr_init(r35152);
        mpfr_init_set_str(r35153, "1", 10, MPFR_RNDN);
        mpfr_init(r35154);
        mpfr_init(r35155);
        mpfr_init(r35156);
        mpfr_init(r35157);
        mpfr_init_set_str(r35158, "15.465634390463869", 10, MPFR_RNDN);
        mpfr_init(r35159);
        mpfr_init(r35160);
        mpfr_init_set_str(r35161, "1/12", 10, MPFR_RNDN);
        mpfr_init(r35162);
        mpfr_init_set_str(r35163, "1/2", 10, MPFR_RNDN);
        mpfr_init(r35164);
        mpfr_init(r35165);
        mpfr_init(r35166);
        mpfr_init(r35167);
}

double f_fm(double x) {
        mpfr_set_d(r35150, x, MPFR_RNDN);
        ;
        mpfr_set_si(r35152, mpfr_cmp(r35150, r35151) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r35154, r35150, MPFR_RNDN);
        mpfr_exp(r35155, r35154, MPFR_RNDN);
        mpfr_sub(r35156, r35153, r35155, MPFR_RNDN);
        mpfr_div(r35157, r35153, r35156, MPFR_RNDN);
        ;
        mpfr_set_si(r35159, mpfr_cmp(r35150, r35158) <= 0, MPFR_RNDN);
        mpfr_div(r35160, r35153, r35150, MPFR_RNDN);
        ;
        mpfr_mul(r35162, r35161, r35150, MPFR_RNDN);
        ;
        mpfr_add(r35164, r35162, r35163, MPFR_RNDN);
        mpfr_add(r35165, r35160, r35164, MPFR_RNDN);
        if (mpfr_get_si(r35159, MPFR_RNDN)) { mpfr_set(r35166, r35165, MPFR_RNDN); } else { mpfr_set(r35166, r35157, MPFR_RNDN); };
        if (mpfr_get_si(r35152, MPFR_RNDN)) { mpfr_set(r35167, r35157, MPFR_RNDN); } else { mpfr_set(r35167, r35166, MPFR_RNDN); };
        return mpfr_get_d(r35167, MPFR_RNDN);
}

static mpfr_t r35168, r35169, r35170, r35171, r35172, r35173, r35174, r35175, r35176, r35177, r35178, r35179, r35180, r35181, r35182, r35183, r35184, r35185;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r35168);
        mpfr_init_set_str(r35169, "-5.458497096453193e-06", 10, MPFR_RNDN);
        mpfr_init(r35170);
        mpfr_init_set_str(r35171, "1", 10, MPFR_RNDN);
        mpfr_init(r35172);
        mpfr_init(r35173);
        mpfr_init(r35174);
        mpfr_init(r35175);
        mpfr_init_set_str(r35176, "15.465634390463869", 10, MPFR_RNDN);
        mpfr_init(r35177);
        mpfr_init(r35178);
        mpfr_init_set_str(r35179, "1/12", 10, MPFR_RNDN);
        mpfr_init(r35180);
        mpfr_init_set_str(r35181, "1/2", 10, MPFR_RNDN);
        mpfr_init(r35182);
        mpfr_init(r35183);
        mpfr_init(r35184);
        mpfr_init(r35185);
}

double f_dm(double x) {
        mpfr_set_d(r35168, x, MPFR_RNDN);
        ;
        mpfr_set_si(r35170, mpfr_cmp(r35168, r35169) <= 0, MPFR_RNDN);
        ;
        mpfr_neg(r35172, r35168, MPFR_RNDN);
        mpfr_exp(r35173, r35172, MPFR_RNDN);
        mpfr_sub(r35174, r35171, r35173, MPFR_RNDN);
        mpfr_div(r35175, r35171, r35174, MPFR_RNDN);
        ;
        mpfr_set_si(r35177, mpfr_cmp(r35168, r35176) <= 0, MPFR_RNDN);
        mpfr_div(r35178, r35171, r35168, MPFR_RNDN);
        ;
        mpfr_mul(r35180, r35179, r35168, MPFR_RNDN);
        ;
        mpfr_add(r35182, r35180, r35181, MPFR_RNDN);
        mpfr_add(r35183, r35178, r35182, MPFR_RNDN);
        if (mpfr_get_si(r35177, MPFR_RNDN)) { mpfr_set(r35184, r35183, MPFR_RNDN); } else { mpfr_set(r35184, r35175, MPFR_RNDN); };
        if (mpfr_get_si(r35170, MPFR_RNDN)) { mpfr_set(r35185, r35175, MPFR_RNDN); } else { mpfr_set(r35185, r35184, MPFR_RNDN); };
        return mpfr_get_d(r35185, MPFR_RNDN);
}