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

char *name = "expax (section 3.5)";

double f_if(float a, float x) {
        float r10828 = a;
        float r10829 = x;
        float r10830 = r10828 * r10829;
        float r10831 = exp(r10830);
        float r10832 = 1;
        float r10833 = r10831 - r10832;
        return r10833;
}

double f_id(double a, double x) {
        double r10834 = a;
        double r10835 = x;
        double r10836 = r10834 * r10835;
        double r10837 = exp(r10836);
        double r10838 = 1;
        double r10839 = r10837 - r10838;
        return r10839;
}


double f_of(float a, float x) {
        float r10840 = a;
        float r10841 = x;
        float r10842 = r10840 * r10841;
        float r10843 = expm1(r10842);
        return r10843;
}

double f_od(double a, double x) {
        double r10844 = a;
        double r10845 = x;
        double r10846 = r10844 * r10845;
        double r10847 = expm1(r10846);
        return r10847;
}

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 r10848, r10849, r10850, r10851, r10852, r10853;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10848);
        mpfr_init(r10849);
        mpfr_init(r10850);
        mpfr_init(r10851);
        mpfr_init_set_str(r10852, "1", 10, MPFR_RNDN);
        mpfr_init(r10853);
}

double f_im(double a, double x) {
        mpfr_set_d(r10848, a, MPFR_RNDN);
        mpfr_set_d(r10849, x, MPFR_RNDN);
        mpfr_mul(r10850, r10848, r10849, MPFR_RNDN);
        mpfr_exp(r10851, r10850, MPFR_RNDN);
        ;
        mpfr_sub(r10853, r10851, r10852, MPFR_RNDN);
        return mpfr_get_d(r10853, MPFR_RNDN);
}

static mpfr_t r10854, r10855, r10856, r10857;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10854);
        mpfr_init(r10855);
        mpfr_init(r10856);
        mpfr_init(r10857);
}

double f_fm(double a, double x) {
        mpfr_set_d(r10854, a, MPFR_RNDN);
        mpfr_set_d(r10855, x, MPFR_RNDN);
        mpfr_mul(r10856, r10854, r10855, MPFR_RNDN);
        mpfr_expm1(r10857, r10856, MPFR_RNDN);
        return mpfr_get_d(r10857, MPFR_RNDN);
}

static mpfr_t r10858, r10859, r10860, r10861;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10858);
        mpfr_init(r10859);
        mpfr_init(r10860);
        mpfr_init(r10861);
}

double f_dm(double a, double x) {
        mpfr_set_d(r10858, a, MPFR_RNDN);
        mpfr_set_d(r10859, x, MPFR_RNDN);
        mpfr_mul(r10860, r10858, r10859, MPFR_RNDN);
        mpfr_expm1(r10861, r10860, MPFR_RNDN);
        return mpfr_get_d(r10861, MPFR_RNDN);
}