#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 r26739 = x;
        float r26740 = exp(r26739);
        float r26741 = 1;
        float r26742 = r26740 - r26741;
        float r26743 = r26740 / r26742;
        return r26743;
}

double f_id(double x) {
        double r26744 = x;
        double r26745 = exp(r26744);
        double r26746 = 1;
        double r26747 = r26745 - r26746;
        double r26748 = r26745 / r26747;
        return r26748;
}


double f_of(float x) {
        float r26749 = 1;
        float r26750 = x;
        float r26751 = -r26750;
        float r26752 = exp(r26751);
        float r26753 = r26749 - r26752;
        float r26754 = r26749 / r26753;
        float r26755 = -82.15533673966395;
        bool r26756 = r26754 <= r26755;
        float r26757 = 1/2;
        float r26758 = r26749 / r26750;
        float r26759 = 1/12;
        float r26760 = r26759 * r26750;
        float r26761 = r26758 + r26760;
        float r26762 = r26757 + r26761;
        float r26763 = 1.0;
        bool r26764 = r26754 <= r26763;
        float r26765 = r26764 ? r26754 : r26762;
        float r26766 = r26756 ? r26762 : r26765;
        return r26766;
}

double f_od(double x) {
        double r26767 = 1;
        double r26768 = x;
        double r26769 = -r26768;
        double r26770 = exp(r26769);
        double r26771 = r26767 - r26770;
        double r26772 = r26767 / r26771;
        double r26773 = -82.15533673966395;
        bool r26774 = r26772 <= r26773;
        double r26775 = 1/2;
        double r26776 = r26767 / r26768;
        double r26777 = 1/12;
        double r26778 = r26777 * r26768;
        double r26779 = r26776 + r26778;
        double r26780 = r26775 + r26779;
        double r26781 = 1.0;
        bool r26782 = r26772 <= r26781;
        double r26783 = r26782 ? r26772 : r26780;
        double r26784 = r26774 ? r26780 : r26783;
        return r26784;
}

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 r26785, r26786, r26787, r26788, r26789;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init_set_str(r26787, "1", 10, MPFR_RNDN);
        mpfr_init(r26788);
        mpfr_init(r26789);
}

double f_im(double x) {
        mpfr_set_d(r26785, x, MPFR_RNDN);
        mpfr_exp(r26786, r26785, MPFR_RNDN);
        ;
        mpfr_sub(r26788, r26786, r26787, MPFR_RNDN);
        mpfr_div(r26789, r26786, r26788, MPFR_RNDN);
        return mpfr_get_d(r26789, MPFR_RNDN);
}

static mpfr_t r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26790, "1", 10, MPFR_RNDN);
        mpfr_init(r26791);
        mpfr_init(r26792);
        mpfr_init(r26793);
        mpfr_init(r26794);
        mpfr_init(r26795);
        mpfr_init_set_str(r26796, "-82.15533673966395", 10, MPFR_RNDN);
        mpfr_init(r26797);
        mpfr_init_set_str(r26798, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26799);
        mpfr_init_set_str(r26800, "1/12", 10, MPFR_RNDN);
        mpfr_init(r26801);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init_set_str(r26804, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26805);
        mpfr_init(r26806);
        mpfr_init(r26807);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r26791, x, MPFR_RNDN);
        mpfr_neg(r26792, r26791, MPFR_RNDN);
        mpfr_exp(r26793, r26792, MPFR_RNDN);
        mpfr_sub(r26794, r26790, r26793, MPFR_RNDN);
        mpfr_div(r26795, r26790, r26794, MPFR_RNDN);
        ;
        mpfr_set_si(r26797, mpfr_cmp(r26795, r26796) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26799, r26790, r26791, MPFR_RNDN);
        ;
        mpfr_mul(r26801, r26800, r26791, MPFR_RNDN);
        mpfr_add(r26802, r26799, r26801, MPFR_RNDN);
        mpfr_add(r26803, r26798, r26802, MPFR_RNDN);
        ;
        mpfr_set_si(r26805, mpfr_cmp(r26795, r26804) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r26805, MPFR_RNDN)) { mpfr_set(r26806, r26795, MPFR_RNDN); } else { mpfr_set(r26806, r26803, MPFR_RNDN); };
        if (mpfr_get_si(r26797, MPFR_RNDN)) { mpfr_set(r26807, r26803, MPFR_RNDN); } else { mpfr_set(r26807, r26806, MPFR_RNDN); };
        return mpfr_get_d(r26807, MPFR_RNDN);
}

static mpfr_t r26808, r26809, r26810, r26811, r26812, r26813, r26814, r26815, r26816, r26817, r26818, r26819, r26820, r26821, r26822, r26823, r26824, r26825;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26808, "1", 10, MPFR_RNDN);
        mpfr_init(r26809);
        mpfr_init(r26810);
        mpfr_init(r26811);
        mpfr_init(r26812);
        mpfr_init(r26813);
        mpfr_init_set_str(r26814, "-82.15533673966395", 10, MPFR_RNDN);
        mpfr_init(r26815);
        mpfr_init_set_str(r26816, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26817);
        mpfr_init_set_str(r26818, "1/12", 10, MPFR_RNDN);
        mpfr_init(r26819);
        mpfr_init(r26820);
        mpfr_init(r26821);
        mpfr_init_set_str(r26822, "1.0", 10, MPFR_RNDN);
        mpfr_init(r26823);
        mpfr_init(r26824);
        mpfr_init(r26825);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r26809, x, MPFR_RNDN);
        mpfr_neg(r26810, r26809, MPFR_RNDN);
        mpfr_exp(r26811, r26810, MPFR_RNDN);
        mpfr_sub(r26812, r26808, r26811, MPFR_RNDN);
        mpfr_div(r26813, r26808, r26812, MPFR_RNDN);
        ;
        mpfr_set_si(r26815, mpfr_cmp(r26813, r26814) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26817, r26808, r26809, MPFR_RNDN);
        ;
        mpfr_mul(r26819, r26818, r26809, MPFR_RNDN);
        mpfr_add(r26820, r26817, r26819, MPFR_RNDN);
        mpfr_add(r26821, r26816, r26820, MPFR_RNDN);
        ;
        mpfr_set_si(r26823, mpfr_cmp(r26813, r26822) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r26823, MPFR_RNDN)) { mpfr_set(r26824, r26813, MPFR_RNDN); } else { mpfr_set(r26824, r26821, MPFR_RNDN); };
        if (mpfr_get_si(r26815, MPFR_RNDN)) { mpfr_set(r26825, r26821, MPFR_RNDN); } else { mpfr_set(r26825, r26824, MPFR_RNDN); };
        return mpfr_get_d(r26825, MPFR_RNDN);
}