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

char *name = "exp2 (problem 3.3.7)";

double f_if(float x) {
        float r10397 = x;
        float r10398 = exp(r10397);
        float r10399 = 2.0f;
        float r10400 = r10398 - r10399;
        float r10401 = -r10397;
        float r10402 = exp(r10401);
        float r10403 = r10400 + r10402;
        return r10403;
}

double f_id(double x) {
        double r10404 = x;
        double r10405 = exp(r10404);
        double r10406 = 2.0;
        double r10407 = r10405 - r10406;
        double r10408 = -r10404;
        double r10409 = exp(r10408);
        double r10410 = r10407 + r10409;
        return r10410;
}


double f_of(float x) {
        float r10411 = 0.08333333333333333f;
        float r10412 = x;
        float r10413 = 4.0f;
        float r10414 = pow(r10412, r10413);
        float r10415 = 0.002777777777777778f;
        float r10416 = 6.0f;
        float r10417 = pow(r10412, r10416);
        float r10418 = r10412 * r10412;
        float r10419 = fma(r10415, r10417, r10418);
        float r10420 = fma(r10411, r10414, r10419);
        return r10420;
}

double f_od(double x) {
        double r10421 = 0.08333333333333333;
        double r10422 = x;
        double r10423 = 4.0;
        double r10424 = pow(r10422, r10423);
        double r10425 = 0.002777777777777778;
        double r10426 = 6.0;
        double r10427 = pow(r10422, r10426);
        double r10428 = r10422 * r10422;
        double r10429 = fma(r10425, r10427, r10428);
        double r10430 = fma(r10421, r10424, r10429);
        return r10430;
}

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 r10431, r10432, r10433, r10434, r10435, r10436, r10437;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10431);
        mpfr_init(r10432);
        mpfr_init_set_str(r10433, "2", 10, MPFR_RNDN);
        mpfr_init(r10434);
        mpfr_init(r10435);
        mpfr_init(r10436);
        mpfr_init(r10437);
}

double f_im(double x) {
        mpfr_set_d(r10431, x, MPFR_RNDN);
        mpfr_exp(r10432, r10431, MPFR_RNDN);
        ;
        mpfr_sub(r10434, r10432, r10433, MPFR_RNDN);
        mpfr_neg(r10435, r10431, MPFR_RNDN);
        mpfr_exp(r10436, r10435, MPFR_RNDN);
        mpfr_add(r10437, r10434, r10436, MPFR_RNDN);
        return mpfr_get_d(r10437, MPFR_RNDN);
}

static mpfr_t r10438, r10439, r10440, r10441, r10442, r10443, r10444, r10445, r10446, r10447;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10438, "1/12", 10, MPFR_RNDN);
        mpfr_init(r10439);
        mpfr_init_set_str(r10440, "4", 10, MPFR_RNDN);
        mpfr_init(r10441);
        mpfr_init_set_str(r10442, "1/360", 10, MPFR_RNDN);
        mpfr_init_set_str(r10443, "6", 10, MPFR_RNDN);
        mpfr_init(r10444);
        mpfr_init(r10445);
        mpfr_init(r10446);
        mpfr_init(r10447);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10439, x, MPFR_RNDN);
        ;
        mpfr_pow(r10441, r10439, r10440, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10444, r10439, r10443, MPFR_RNDN);
        mpfr_mul(r10445, r10439, r10439, MPFR_RNDN);
        mpfr_fma(r10446, r10442, r10444, r10445, MPFR_RNDN);
        mpfr_fma(r10447, r10438, r10441, r10446, MPFR_RNDN);
        return mpfr_get_d(r10447, MPFR_RNDN);
}

static mpfr_t r10448, r10449, r10450, r10451, r10452, r10453, r10454, r10455, r10456, r10457;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10448, "1/12", 10, MPFR_RNDN);
        mpfr_init(r10449);
        mpfr_init_set_str(r10450, "4", 10, MPFR_RNDN);
        mpfr_init(r10451);
        mpfr_init_set_str(r10452, "1/360", 10, MPFR_RNDN);
        mpfr_init_set_str(r10453, "6", 10, MPFR_RNDN);
        mpfr_init(r10454);
        mpfr_init(r10455);
        mpfr_init(r10456);
        mpfr_init(r10457);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10449, x, MPFR_RNDN);
        ;
        mpfr_pow(r10451, r10449, r10450, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10454, r10449, r10453, MPFR_RNDN);
        mpfr_mul(r10455, r10449, r10449, MPFR_RNDN);
        mpfr_fma(r10456, r10452, r10454, r10455, MPFR_RNDN);
        mpfr_fma(r10457, r10448, r10451, r10456, MPFR_RNDN);
        return mpfr_get_d(r10457, MPFR_RNDN);
}