#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 r24270 = a;
        float r24271 = x;
        float r24272 = r24270 * r24271;
        float r24273 = exp(r24272);
        float r24274 = 1;
        float r24275 = r24273 - r24274;
        return r24275;
}

double f_id(double a, double x) {
        double r24276 = a;
        double r24277 = x;
        double r24278 = r24276 * r24277;
        double r24279 = exp(r24278);
        double r24280 = 1;
        double r24281 = r24279 - r24280;
        return r24281;
}


double f_of(float a, float x) {
        float r24282 = a;
        float r24283 = x;
        float r24284 = r24282 * r24283;
        float r24285 = exp(r24284);
        float r24286 = 1;
        float r24287 = r24285 - r24286;
        float r24288 = -1.6691627015567744e-07;
        bool r24289 = r24287 <= r24288;
        float r24290 = r24284 * r24284;
        float r24291 = 1/6;
        float r24292 = r24284 * r24291;
        float r24293 = 1/2;
        float r24294 = r24292 + r24293;
        float r24295 = r24290 * r24294;
        float r24296 = r24295 + r24284;
        float r24297 = r24289 ? r24287 : r24296;
        return r24297;
}

double f_od(double a, double x) {
        double r24298 = a;
        double r24299 = x;
        double r24300 = r24298 * r24299;
        double r24301 = exp(r24300);
        double r24302 = 1;
        double r24303 = r24301 - r24302;
        double r24304 = -1.6691627015567744e-07;
        bool r24305 = r24303 <= r24304;
        double r24306 = r24300 * r24300;
        double r24307 = 1/6;
        double r24308 = r24300 * r24307;
        double r24309 = 1/2;
        double r24310 = r24308 + r24309;
        double r24311 = r24306 * r24310;
        double r24312 = r24311 + r24300;
        double r24313 = r24305 ? r24303 : r24312;
        return r24313;
}

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 r24314, r24315, r24316, r24317, r24318, r24319;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r24314);
        mpfr_init(r24315);
        mpfr_init(r24316);
        mpfr_init(r24317);
        mpfr_init_set_str(r24318, "1", 10, MPFR_RNDN);
        mpfr_init(r24319);
}

double f_im(double a, double x) {
        mpfr_set_d(r24314, a, MPFR_RNDN);
        mpfr_set_d(r24315, x, MPFR_RNDN);
        mpfr_mul(r24316, r24314, r24315, MPFR_RNDN);
        mpfr_exp(r24317, r24316, MPFR_RNDN);
        ;
        mpfr_sub(r24319, r24317, r24318, MPFR_RNDN);
        return mpfr_get_d(r24319, MPFR_RNDN);
}

static mpfr_t r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332, r24333, r24334, r24335;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r24320);
        mpfr_init(r24321);
        mpfr_init(r24322);
        mpfr_init(r24323);
        mpfr_init_set_str(r24324, "1", 10, MPFR_RNDN);
        mpfr_init(r24325);
        mpfr_init_set_str(r24326, "-1.6691627015567744e-07", 10, MPFR_RNDN);
        mpfr_init(r24327);
        mpfr_init(r24328);
        mpfr_init_set_str(r24329, "1/6", 10, MPFR_RNDN);
        mpfr_init(r24330);
        mpfr_init_set_str(r24331, "1/2", 10, MPFR_RNDN);
        mpfr_init(r24332);
        mpfr_init(r24333);
        mpfr_init(r24334);
        mpfr_init(r24335);
}

double f_fm(double a, double x) {
        mpfr_set_d(r24320, a, MPFR_RNDN);
        mpfr_set_d(r24321, x, MPFR_RNDN);
        mpfr_mul(r24322, r24320, r24321, MPFR_RNDN);
        mpfr_exp(r24323, r24322, MPFR_RNDN);
        ;
        mpfr_sub(r24325, r24323, r24324, MPFR_RNDN);
        ;
        mpfr_set_si(r24327, mpfr_cmp(r24325, r24326) <= 0, MPFR_RNDN);
        mpfr_mul(r24328, r24322, r24322, MPFR_RNDN);
        ;
        mpfr_mul(r24330, r24322, r24329, MPFR_RNDN);
        ;
        mpfr_add(r24332, r24330, r24331, MPFR_RNDN);
        mpfr_mul(r24333, r24328, r24332, MPFR_RNDN);
        mpfr_add(r24334, r24333, r24322, MPFR_RNDN);
        if (mpfr_get_si(r24327, MPFR_RNDN)) { mpfr_set(r24335, r24325, MPFR_RNDN); } else { mpfr_set(r24335, r24334, MPFR_RNDN); };
        return mpfr_get_d(r24335, MPFR_RNDN);
}

static mpfr_t r24336, r24337, r24338, r24339, r24340, r24341, r24342, r24343, r24344, r24345, r24346, r24347, r24348, r24349, r24350, r24351;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r24336);
        mpfr_init(r24337);
        mpfr_init(r24338);
        mpfr_init(r24339);
        mpfr_init_set_str(r24340, "1", 10, MPFR_RNDN);
        mpfr_init(r24341);
        mpfr_init_set_str(r24342, "-1.6691627015567744e-07", 10, MPFR_RNDN);
        mpfr_init(r24343);
        mpfr_init(r24344);
        mpfr_init_set_str(r24345, "1/6", 10, MPFR_RNDN);
        mpfr_init(r24346);
        mpfr_init_set_str(r24347, "1/2", 10, MPFR_RNDN);
        mpfr_init(r24348);
        mpfr_init(r24349);
        mpfr_init(r24350);
        mpfr_init(r24351);
}

double f_dm(double a, double x) {
        mpfr_set_d(r24336, a, MPFR_RNDN);
        mpfr_set_d(r24337, x, MPFR_RNDN);
        mpfr_mul(r24338, r24336, r24337, MPFR_RNDN);
        mpfr_exp(r24339, r24338, MPFR_RNDN);
        ;
        mpfr_sub(r24341, r24339, r24340, MPFR_RNDN);
        ;
        mpfr_set_si(r24343, mpfr_cmp(r24341, r24342) <= 0, MPFR_RNDN);
        mpfr_mul(r24344, r24338, r24338, MPFR_RNDN);
        ;
        mpfr_mul(r24346, r24338, r24345, MPFR_RNDN);
        ;
        mpfr_add(r24348, r24346, r24347, MPFR_RNDN);
        mpfr_mul(r24349, r24344, r24348, MPFR_RNDN);
        mpfr_add(r24350, r24349, r24338, MPFR_RNDN);
        if (mpfr_get_si(r24343, MPFR_RNDN)) { mpfr_set(r24351, r24341, MPFR_RNDN); } else { mpfr_set(r24351, r24350, MPFR_RNDN); };
        return mpfr_get_d(r24351, MPFR_RNDN);
}