#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 r26775 = x;
        float r26776 = exp(r26775);
        float r26777 = 1;
        float r26778 = r26776 - r26777;
        float r26779 = r26776 / r26778;
        return r26779;
}

double f_id(double x) {
        double r26780 = x;
        double r26781 = exp(r26780);
        double r26782 = 1;
        double r26783 = r26781 - r26782;
        double r26784 = r26781 / r26783;
        return r26784;
}


double f_of(float x) {
        float r26785 = 1;
        float r26786 = x;
        float r26787 = -r26786;
        float r26788 = exp(r26787);
        float r26789 = r26785 - r26788;
        float r26790 = r26785 / r26789;
        float r26791 = 1.0308725086867097;
        bool r26792 = r26790 <= r26791;
        float r26793 = 1/2;
        float r26794 = r26785 / r26786;
        float r26795 = 1/12;
        float r26796 = r26795 * r26786;
        float r26797 = r26794 + r26796;
        float r26798 = r26793 + r26797;
        float r26799 = r26792 ? r26790 : r26798;
        return r26799;
}

double f_od(double x) {
        double r26800 = 1;
        double r26801 = x;
        double r26802 = -r26801;
        double r26803 = exp(r26802);
        double r26804 = r26800 - r26803;
        double r26805 = r26800 / r26804;
        double r26806 = 1.0308725086867097;
        bool r26807 = r26805 <= r26806;
        double r26808 = 1/2;
        double r26809 = r26800 / r26801;
        double r26810 = 1/12;
        double r26811 = r26810 * r26801;
        double r26812 = r26809 + r26811;
        double r26813 = r26808 + r26812;
        double r26814 = r26807 ? r26805 : r26813;
        return r26814;
}

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 r26815, r26816, r26817, r26818, r26819;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26815);
        mpfr_init(r26816);
        mpfr_init_set_str(r26817, "1", 10, MPFR_RNDN);
        mpfr_init(r26818);
        mpfr_init(r26819);
}

double f_im(double x) {
        mpfr_set_d(r26815, x, MPFR_RNDN);
        mpfr_exp(r26816, r26815, MPFR_RNDN);
        ;
        mpfr_sub(r26818, r26816, r26817, MPFR_RNDN);
        mpfr_div(r26819, r26816, r26818, MPFR_RNDN);
        return mpfr_get_d(r26819, MPFR_RNDN);
}

static mpfr_t r26820, r26821, r26822, r26823, r26824, r26825, r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26820, "1", 10, MPFR_RNDN);
        mpfr_init(r26821);
        mpfr_init(r26822);
        mpfr_init(r26823);
        mpfr_init(r26824);
        mpfr_init(r26825);
        mpfr_init_set_str(r26826, "1.0308725086867097", 10, MPFR_RNDN);
        mpfr_init(r26827);
        mpfr_init_set_str(r26828, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26829);
        mpfr_init_set_str(r26830, "1/12", 10, MPFR_RNDN);
        mpfr_init(r26831);
        mpfr_init(r26832);
        mpfr_init(r26833);
        mpfr_init(r26834);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r26821, x, MPFR_RNDN);
        mpfr_neg(r26822, r26821, MPFR_RNDN);
        mpfr_exp(r26823, r26822, MPFR_RNDN);
        mpfr_sub(r26824, r26820, r26823, MPFR_RNDN);
        mpfr_div(r26825, r26820, r26824, MPFR_RNDN);
        ;
        mpfr_set_si(r26827, mpfr_cmp(r26825, r26826) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26829, r26820, r26821, MPFR_RNDN);
        ;
        mpfr_mul(r26831, r26830, r26821, MPFR_RNDN);
        mpfr_add(r26832, r26829, r26831, MPFR_RNDN);
        mpfr_add(r26833, r26828, r26832, MPFR_RNDN);
        if (mpfr_get_si(r26827, MPFR_RNDN)) { mpfr_set(r26834, r26825, MPFR_RNDN); } else { mpfr_set(r26834, r26833, MPFR_RNDN); };
        return mpfr_get_d(r26834, MPFR_RNDN);
}

static mpfr_t r26835, r26836, r26837, r26838, r26839, r26840, r26841, r26842, r26843, r26844, r26845, r26846, r26847, r26848, r26849;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26835, "1", 10, MPFR_RNDN);
        mpfr_init(r26836);
        mpfr_init(r26837);
        mpfr_init(r26838);
        mpfr_init(r26839);
        mpfr_init(r26840);
        mpfr_init_set_str(r26841, "1.0308725086867097", 10, MPFR_RNDN);
        mpfr_init(r26842);
        mpfr_init_set_str(r26843, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26844);
        mpfr_init_set_str(r26845, "1/12", 10, MPFR_RNDN);
        mpfr_init(r26846);
        mpfr_init(r26847);
        mpfr_init(r26848);
        mpfr_init(r26849);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r26836, x, MPFR_RNDN);
        mpfr_neg(r26837, r26836, MPFR_RNDN);
        mpfr_exp(r26838, r26837, MPFR_RNDN);
        mpfr_sub(r26839, r26835, r26838, MPFR_RNDN);
        mpfr_div(r26840, r26835, r26839, MPFR_RNDN);
        ;
        mpfr_set_si(r26842, mpfr_cmp(r26840, r26841) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r26844, r26835, r26836, MPFR_RNDN);
        ;
        mpfr_mul(r26846, r26845, r26836, MPFR_RNDN);
        mpfr_add(r26847, r26844, r26846, MPFR_RNDN);
        mpfr_add(r26848, r26843, r26847, MPFR_RNDN);
        if (mpfr_get_si(r26842, MPFR_RNDN)) { mpfr_set(r26849, r26840, MPFR_RNDN); } else { mpfr_set(r26849, r26848, MPFR_RNDN); };
        return mpfr_get_d(r26849, MPFR_RNDN);
}