#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 r27403 = x;
        float r27404 = exp(r27403);
        float r27405 = 2;
        float r27406 = r27404 - r27405;
        float r27407 = -r27403;
        float r27408 = exp(r27407);
        float r27409 = r27406 + r27408;
        return r27409;
}

double f_id(double x) {
        double r27410 = x;
        double r27411 = exp(r27410);
        double r27412 = 2;
        double r27413 = r27411 - r27412;
        double r27414 = -r27410;
        double r27415 = exp(r27414);
        double r27416 = r27413 + r27415;
        return r27416;
}


double f_of(float x) {
        float r27417 = x;
        float r27418 = 2;
        float r27419 = pow(r27417, r27418);
        float r27420 = 1/12;
        float r27421 = 4;
        float r27422 = pow(r27417, r27421);
        float r27423 = r27420 * r27422;
        float r27424 = 1/360;
        float r27425 = 6;
        float r27426 = pow(r27417, r27425);
        float r27427 = r27424 * r27426;
        float r27428 = r27423 + r27427;
        float r27429 = r27419 + r27428;
        return r27429;
}

double f_od(double x) {
        double r27430 = x;
        double r27431 = 2;
        double r27432 = pow(r27430, r27431);
        double r27433 = 1/12;
        double r27434 = 4;
        double r27435 = pow(r27430, r27434);
        double r27436 = r27433 * r27435;
        double r27437 = 1/360;
        double r27438 = 6;
        double r27439 = pow(r27430, r27438);
        double r27440 = r27437 * r27439;
        double r27441 = r27436 + r27440;
        double r27442 = r27432 + r27441;
        return r27442;
}

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 r27443, r27444, r27445, r27446, r27447, r27448, r27449;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27443);
        mpfr_init(r27444);
        mpfr_init_set_str(r27445, "2", 10, MPFR_RNDN);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init(r27449);
}

double f_im(double x) {
        mpfr_set_d(r27443, x, MPFR_RNDN);
        mpfr_exp(r27444, r27443, MPFR_RNDN);
        ;
        mpfr_sub(r27446, r27444, r27445, MPFR_RNDN);
        mpfr_neg(r27447, r27443, MPFR_RNDN);
        mpfr_exp(r27448, r27447, MPFR_RNDN);
        mpfr_add(r27449, r27446, r27448, MPFR_RNDN);
        return mpfr_get_d(r27449, MPFR_RNDN);
}

static mpfr_t r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459, r27460, r27461, r27462;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27450);
        mpfr_init_set_str(r27451, "2", 10, MPFR_RNDN);
        mpfr_init(r27452);
        mpfr_init_set_str(r27453, "1/12", 10, MPFR_RNDN);
        mpfr_init_set_str(r27454, "4", 10, MPFR_RNDN);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init_set_str(r27457, "1/360", 10, MPFR_RNDN);
        mpfr_init_set_str(r27458, "6", 10, MPFR_RNDN);
        mpfr_init(r27459);
        mpfr_init(r27460);
        mpfr_init(r27461);
        mpfr_init(r27462);
}

double f_fm(double x) {
        mpfr_set_d(r27450, x, MPFR_RNDN);
        ;
        mpfr_pow(r27452, r27450, r27451, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27455, r27450, r27454, MPFR_RNDN);
        mpfr_mul(r27456, r27453, r27455, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27459, r27450, r27458, MPFR_RNDN);
        mpfr_mul(r27460, r27457, r27459, MPFR_RNDN);
        mpfr_add(r27461, r27456, r27460, MPFR_RNDN);
        mpfr_add(r27462, r27452, r27461, MPFR_RNDN);
        return mpfr_get_d(r27462, MPFR_RNDN);
}

static mpfr_t r27463, r27464, r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472, r27473, r27474, r27475;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27463);
        mpfr_init_set_str(r27464, "2", 10, MPFR_RNDN);
        mpfr_init(r27465);
        mpfr_init_set_str(r27466, "1/12", 10, MPFR_RNDN);
        mpfr_init_set_str(r27467, "4", 10, MPFR_RNDN);
        mpfr_init(r27468);
        mpfr_init(r27469);
        mpfr_init_set_str(r27470, "1/360", 10, MPFR_RNDN);
        mpfr_init_set_str(r27471, "6", 10, MPFR_RNDN);
        mpfr_init(r27472);
        mpfr_init(r27473);
        mpfr_init(r27474);
        mpfr_init(r27475);
}

double f_dm(double x) {
        mpfr_set_d(r27463, x, MPFR_RNDN);
        ;
        mpfr_pow(r27465, r27463, r27464, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27468, r27463, r27467, MPFR_RNDN);
        mpfr_mul(r27469, r27466, r27468, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27472, r27463, r27471, MPFR_RNDN);
        mpfr_mul(r27473, r27470, r27472, MPFR_RNDN);
        mpfr_add(r27474, r27469, r27473, MPFR_RNDN);
        mpfr_add(r27475, r27465, r27474, MPFR_RNDN);
        return mpfr_get_d(r27475, MPFR_RNDN);
}

