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

char *name = "expm1 (example 3.7)";

double f_if(float x) {
        float r35186 = x;
        float r35187 = exp(r35186);
        float r35188 = 1.0f;
        float r35189 = r35187 - r35188;
        return r35189;
}

double f_id(double x) {
        double r35190 = x;
        double r35191 = exp(r35190);
        double r35192 = 1.0;
        double r35193 = r35191 - r35192;
        return r35193;
}


double f_of(float x) {
        float r35194 = x;
        float r35195 = r35194 * r35194;
        float r35196 = 0.16666666666666666f;
        float r35197 = r35196 * r35194;
        float r35198 = 0.5f;
        float r35199 = r35197 + r35198;
        float r35200 = r35195 * r35199;
        float r35201 = r35194 + r35200;
        return r35201;
}

double f_od(double x) {
        double r35202 = x;
        double r35203 = r35202 * r35202;
        double r35204 = 0.16666666666666666;
        double r35205 = r35204 * r35202;
        double r35206 = 0.5;
        double r35207 = r35205 + r35206;
        double r35208 = r35203 * r35207;
        double r35209 = r35202 + r35208;
        return r35209;
}

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 r35210, r35211, r35212, r35213;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r35210);
        mpfr_init(r35211);
        mpfr_init_set_str(r35212, "1", 10, MPFR_RNDN);
        mpfr_init(r35213);
}

double f_im(double x) {
        mpfr_set_d(r35210, x, MPFR_RNDN);
        mpfr_exp(r35211, r35210, MPFR_RNDN);
        ;
        mpfr_sub(r35213, r35211, r35212, MPFR_RNDN);
        return mpfr_get_d(r35213, MPFR_RNDN);
}

static mpfr_t r35214, r35215, r35216, r35217, r35218, r35219, r35220, r35221;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r35214);
        mpfr_init(r35215);
        mpfr_init_set_str(r35216, "1/6", 10, MPFR_RNDN);
        mpfr_init(r35217);
        mpfr_init_set_str(r35218, "1/2", 10, MPFR_RNDN);
        mpfr_init(r35219);
        mpfr_init(r35220);
        mpfr_init(r35221);
}

double f_fm(double x) {
        mpfr_set_d(r35214, x, MPFR_RNDN);
        mpfr_sqr(r35215, r35214, MPFR_RNDN);
        ;
        mpfr_mul(r35217, r35216, r35214, MPFR_RNDN);
        ;
        mpfr_add(r35219, r35217, r35218, MPFR_RNDN);
        mpfr_mul(r35220, r35215, r35219, MPFR_RNDN);
        mpfr_add(r35221, r35214, r35220, MPFR_RNDN);
        return mpfr_get_d(r35221, MPFR_RNDN);
}

static mpfr_t r35222, r35223, r35224, r35225, r35226, r35227, r35228, r35229;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r35222);
        mpfr_init(r35223);
        mpfr_init_set_str(r35224, "1/6", 10, MPFR_RNDN);
        mpfr_init(r35225);
        mpfr_init_set_str(r35226, "1/2", 10, MPFR_RNDN);
        mpfr_init(r35227);
        mpfr_init(r35228);
        mpfr_init(r35229);
}

double f_dm(double x) {
        mpfr_set_d(r35222, x, MPFR_RNDN);
        mpfr_sqr(r35223, r35222, MPFR_RNDN);
        ;
        mpfr_mul(r35225, r35224, r35222, MPFR_RNDN);
        ;
        mpfr_add(r35227, r35225, r35226, MPFR_RNDN);
        mpfr_mul(r35228, r35223, r35227, MPFR_RNDN);
        mpfr_add(r35229, r35222, r35228, MPFR_RNDN);
        return mpfr_get_d(r35229, MPFR_RNDN);
}