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

char *name = "NMSE Section 6.1 mentioned, A";

double f_if(float x, float eps) {
        float r26195 = 1.0f;
        float r26196 = eps;
        float r26197 = r26195 / r26196;
        float r26198 = r26195 + r26197;
        float r26199 = r26195 - r26196;
        float r26200 = x;
        float r26201 = r26199 * r26200;
        float r26202 = -r26201;
        float r26203 = exp(r26202);
        float r26204 = r26198 * r26203;
        float r26205 = r26197 - r26195;
        float r26206 = r26195 + r26196;
        float r26207 = r26206 * r26200;
        float r26208 = -r26207;
        float r26209 = exp(r26208);
        float r26210 = r26205 * r26209;
        float r26211 = r26204 - r26210;
        float r26212 = 2.0f;
        float r26213 = r26211 / r26212;
        return r26213;
}

double f_id(double x, double eps) {
        double r26214 = 1.0;
        double r26215 = eps;
        double r26216 = r26214 / r26215;
        double r26217 = r26214 + r26216;
        double r26218 = r26214 - r26215;
        double r26219 = x;
        double r26220 = r26218 * r26219;
        double r26221 = -r26220;
        double r26222 = exp(r26221);
        double r26223 = r26217 * r26222;
        double r26224 = r26216 - r26214;
        double r26225 = r26214 + r26215;
        double r26226 = r26225 * r26219;
        double r26227 = -r26226;
        double r26228 = exp(r26227);
        double r26229 = r26224 * r26228;
        double r26230 = r26223 - r26229;
        double r26231 = 2.0;
        double r26232 = r26230 / r26231;
        return r26232;
}


double f_of(float x, float eps) {
        float r26233 = x;
        float r26234 = 921224087.6603982f;
        bool r26235 = r26233 <= r26234;
        float r26236 = 1.0f;
        float r26237 = eps;
        float r26238 = r26236 / r26237;
        float r26239 = r26236 + r26238;
        float r26240 = exp(r26233);
        float r26241 = r26236 - r26237;
        float r26242 = pow(r26240, r26241);
        float r26243 = r26239 / r26242;
        float r26244 = cbrt(r26243);
        float r26245 = r26244 * (r26244 * r26244);
        float r26246 = r26238 - r26236;
        float r26247 = r26236 + r26237;
        float r26248 = r26247 * r26233;
        float r26249 = -r26248;
        float r26250 = exp(r26249);
        float r26251 = r26246 * r26250;
        float r26252 = r26245 - r26251;
        float r26253 = 2.0f;
        float r26254 = r26252 / r26253;
        float r26255 = r26235 ? r26236 : r26254;
        return r26255;
}

double f_od(double x, double eps) {
        double r26256 = x;
        double r26257 = 921224087.6603982;
        bool r26258 = r26256 <= r26257;
        double r26259 = 1.0;
        double r26260 = eps;
        double r26261 = r26259 / r26260;
        double r26262 = r26259 + r26261;
        double r26263 = exp(r26256);
        double r26264 = r26259 - r26260;
        double r26265 = pow(r26263, r26264);
        double r26266 = r26262 / r26265;
        double r26267 = cbrt(r26266);
        double r26268 = r26267 * (r26267 * r26267);
        double r26269 = r26261 - r26259;
        double r26270 = r26259 + r26260;
        double r26271 = r26270 * r26256;
        double r26272 = -r26271;
        double r26273 = exp(r26272);
        double r26274 = r26269 * r26273;
        double r26275 = r26268 - r26274;
        double r26276 = 2.0;
        double r26277 = r26275 / r26276;
        double r26278 = r26258 ? r26259 : r26277;
        return r26278;
}

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 r26279, r26280, r26281, r26282, r26283, r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r26279, "1", 10, MPFR_RNDN);
        mpfr_init(r26280);
        mpfr_init(r26281);
        mpfr_init(r26282);
        mpfr_init(r26283);
        mpfr_init(r26284);
        mpfr_init(r26285);
        mpfr_init(r26286);
        mpfr_init(r26287);
        mpfr_init(r26288);
        mpfr_init(r26289);
        mpfr_init(r26290);
        mpfr_init(r26291);
        mpfr_init(r26292);
        mpfr_init(r26293);
        mpfr_init(r26294);
        mpfr_init(r26295);
        mpfr_init_set_str(r26296, "2", 10, MPFR_RNDN);
        mpfr_init(r26297);
}

double f_im(double x, double eps) {
        ;
        mpfr_set_d(r26280, eps, MPFR_RNDN);
        mpfr_div(r26281, r26279, r26280, MPFR_RNDN);
        mpfr_add(r26282, r26279, r26281, MPFR_RNDN);
        mpfr_sub(r26283, r26279, r26280, MPFR_RNDN);
        mpfr_set_d(r26284, x, MPFR_RNDN);
        mpfr_mul(r26285, r26283, r26284, MPFR_RNDN);
        mpfr_neg(r26286, r26285, MPFR_RNDN);
        mpfr_exp(r26287, r26286, MPFR_RNDN);
        mpfr_mul(r26288, r26282, r26287, MPFR_RNDN);
        mpfr_sub(r26289, r26281, r26279, MPFR_RNDN);
        mpfr_add(r26290, r26279, r26280, MPFR_RNDN);
        mpfr_mul(r26291, r26290, r26284, MPFR_RNDN);
        mpfr_neg(r26292, r26291, MPFR_RNDN);
        mpfr_exp(r26293, r26292, MPFR_RNDN);
        mpfr_mul(r26294, r26289, r26293, MPFR_RNDN);
        mpfr_sub(r26295, r26288, r26294, MPFR_RNDN);
        ;
        mpfr_div(r26297, r26295, r26296, MPFR_RNDN);
        return mpfr_get_d(r26297, MPFR_RNDN);
}

static mpfr_t r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r26298);
        mpfr_init_set_str(r26299, "921224087.6603982", 10, MPFR_RNDN);
        mpfr_init(r26300);
        mpfr_init_set_str(r26301, "1", 10, MPFR_RNDN);
        mpfr_init(r26302);
        mpfr_init(r26303);
        mpfr_init(r26304);
        mpfr_init(r26305);
        mpfr_init(r26306);
        mpfr_init(r26307);
        mpfr_init(r26308);
        mpfr_init(r26309);
        mpfr_init(r26310);
        mpfr_init(r26311);
        mpfr_init(r26312);
        mpfr_init(r26313);
        mpfr_init(r26314);
        mpfr_init(r26315);
        mpfr_init(r26316);
        mpfr_init(r26317);
        mpfr_init_set_str(r26318, "2", 10, MPFR_RNDN);
        mpfr_init(r26319);
        mpfr_init(r26320);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r26298, x, MPFR_RNDN);
        ;
        mpfr_set_si(r26300, mpfr_cmp(r26298, r26299) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r26302, eps, MPFR_RNDN);
        mpfr_div(r26303, r26301, r26302, MPFR_RNDN);
        mpfr_add(r26304, r26301, r26303, MPFR_RNDN);
        mpfr_exp(r26305, r26298, MPFR_RNDN);
        mpfr_sub(r26306, r26301, r26302, MPFR_RNDN);
        mpfr_pow(r26307, r26305, r26306, MPFR_RNDN);
        mpfr_div(r26308, r26304, r26307, MPFR_RNDN);
        mpfr_cbrt(r26309, r26308, MPFR_RNDN);
        mpfr_mul(r26310, r26309, r26309, MPFR_RNDN); mpfr_mul(r26310, r26310, r26309, MPFR_RNDN);
        mpfr_sub(r26311, r26303, r26301, MPFR_RNDN);
        mpfr_add(r26312, r26301, r26302, MPFR_RNDN);
        mpfr_mul(r26313, r26312, r26298, MPFR_RNDN);
        mpfr_neg(r26314, r26313, MPFR_RNDN);
        mpfr_exp(r26315, r26314, MPFR_RNDN);
        mpfr_mul(r26316, r26311, r26315, MPFR_RNDN);
        mpfr_sub(r26317, r26310, r26316, MPFR_RNDN);
        ;
        mpfr_div(r26319, r26317, r26318, MPFR_RNDN);
        if (mpfr_get_si(r26300, MPFR_RNDN)) { mpfr_set(r26320, r26301, MPFR_RNDN); } else { mpfr_set(r26320, r26319, MPFR_RNDN); };
        return mpfr_get_d(r26320, MPFR_RNDN);
}

static mpfr_t r26321, r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r26321);
        mpfr_init_set_str(r26322, "921224087.6603982", 10, MPFR_RNDN);
        mpfr_init(r26323);
        mpfr_init_set_str(r26324, "1", 10, MPFR_RNDN);
        mpfr_init(r26325);
        mpfr_init(r26326);
        mpfr_init(r26327);
        mpfr_init(r26328);
        mpfr_init(r26329);
        mpfr_init(r26330);
        mpfr_init(r26331);
        mpfr_init(r26332);
        mpfr_init(r26333);
        mpfr_init(r26334);
        mpfr_init(r26335);
        mpfr_init(r26336);
        mpfr_init(r26337);
        mpfr_init(r26338);
        mpfr_init(r26339);
        mpfr_init(r26340);
        mpfr_init_set_str(r26341, "2", 10, MPFR_RNDN);
        mpfr_init(r26342);
        mpfr_init(r26343);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r26321, x, MPFR_RNDN);
        ;
        mpfr_set_si(r26323, mpfr_cmp(r26321, r26322) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r26325, eps, MPFR_RNDN);
        mpfr_div(r26326, r26324, r26325, MPFR_RNDN);
        mpfr_add(r26327, r26324, r26326, MPFR_RNDN);
        mpfr_exp(r26328, r26321, MPFR_RNDN);
        mpfr_sub(r26329, r26324, r26325, MPFR_RNDN);
        mpfr_pow(r26330, r26328, r26329, MPFR_RNDN);
        mpfr_div(r26331, r26327, r26330, MPFR_RNDN);
        mpfr_cbrt(r26332, r26331, MPFR_RNDN);
        mpfr_mul(r26333, r26332, r26332, MPFR_RNDN); mpfr_mul(r26333, r26333, r26332, MPFR_RNDN);
        mpfr_sub(r26334, r26326, r26324, MPFR_RNDN);
        mpfr_add(r26335, r26324, r26325, MPFR_RNDN);
        mpfr_mul(r26336, r26335, r26321, MPFR_RNDN);
        mpfr_neg(r26337, r26336, MPFR_RNDN);
        mpfr_exp(r26338, r26337, MPFR_RNDN);
        mpfr_mul(r26339, r26334, r26338, MPFR_RNDN);
        mpfr_sub(r26340, r26333, r26339, MPFR_RNDN);
        ;
        mpfr_div(r26342, r26340, r26341, MPFR_RNDN);
        if (mpfr_get_si(r26323, MPFR_RNDN)) { mpfr_set(r26343, r26324, MPFR_RNDN); } else { mpfr_set(r26343, r26342, MPFR_RNDN); };
        return mpfr_get_d(r26343, MPFR_RNDN);
}