#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 r22048 = 1;
        float r22049 = eps;
        float r22050 = r22048 / r22049;
        float r22051 = r22048 + r22050;
        float r22052 = r22048 - r22049;
        float r22053 = x;
        float r22054 = r22052 * r22053;
        float r22055 = -r22054;
        float r22056 = exp(r22055);
        float r22057 = r22051 * r22056;
        float r22058 = r22050 - r22048;
        float r22059 = r22048 + r22049;
        float r22060 = r22059 * r22053;
        float r22061 = -r22060;
        float r22062 = exp(r22061);
        float r22063 = r22058 * r22062;
        float r22064 = r22057 - r22063;
        float r22065 = 2;
        float r22066 = r22064 / r22065;
        return r22066;
}

double f_id(double x, double eps) {
        double r22067 = 1;
        double r22068 = eps;
        double r22069 = r22067 / r22068;
        double r22070 = r22067 + r22069;
        double r22071 = r22067 - r22068;
        double r22072 = x;
        double r22073 = r22071 * r22072;
        double r22074 = -r22073;
        double r22075 = exp(r22074);
        double r22076 = r22070 * r22075;
        double r22077 = r22069 - r22067;
        double r22078 = r22067 + r22068;
        double r22079 = r22078 * r22072;
        double r22080 = -r22079;
        double r22081 = exp(r22080);
        double r22082 = r22077 * r22081;
        double r22083 = r22076 - r22082;
        double r22084 = 2;
        double r22085 = r22083 / r22084;
        return r22085;
}


double f_of(float x, float eps) {
        float r22086 = 2/3;
        float r22087 = x;
        float r22088 = r22086 * r22087;
        float r22089 = r22087 * r22087;
        float r22090 = 2;
        float r22091 = r22090 - r22089;
        float r22092 = fma(r22088, r22089, r22091);
        float r22093 = r22092 / r22090;
        float r22094 = 1405.389636089941;
        bool r22095 = r22093 <= r22094;
        float r22096 = 1;
        float r22097 = eps;
        float r22098 = r22096 / r22097;
        float r22099 = r22096 + r22098;
        float r22100 = r22096 - r22097;
        float r22101 = -r22087;
        float r22102 = r22100 * r22101;
        float r22103 = exp(r22102);
        float r22104 = r22096 - r22098;
        float r22105 = fma(r22097, r22087, r22087);
        float r22106 = exp(r22105);
        float r22107 = r22104 / r22106;
        float r22108 = fma(r22099, r22103, r22107);
        float r22109 = r22108 / r22090;
        float r22110 = r22095 ? r22093 : r22109;
        return r22110;
}

double f_od(double x, double eps) {
        double r22111 = 2/3;
        double r22112 = x;
        double r22113 = r22111 * r22112;
        double r22114 = r22112 * r22112;
        double r22115 = 2;
        double r22116 = r22115 - r22114;
        double r22117 = fma(r22113, r22114, r22116);
        double r22118 = r22117 / r22115;
        double r22119 = 1405.389636089941;
        bool r22120 = r22118 <= r22119;
        double r22121 = 1;
        double r22122 = eps;
        double r22123 = r22121 / r22122;
        double r22124 = r22121 + r22123;
        double r22125 = r22121 - r22122;
        double r22126 = -r22112;
        double r22127 = r22125 * r22126;
        double r22128 = exp(r22127);
        double r22129 = r22121 - r22123;
        double r22130 = fma(r22122, r22112, r22112);
        double r22131 = exp(r22130);
        double r22132 = r22129 / r22131;
        double r22133 = fma(r22124, r22128, r22132);
        double r22134 = r22133 / r22115;
        double r22135 = r22120 ? r22118 : r22134;
        return r22135;
}

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 r22136, r22137, r22138, r22139, r22140, r22141, r22142, r22143, r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22136, "1", 10, MPFR_RNDN);
        mpfr_init(r22137);
        mpfr_init(r22138);
        mpfr_init(r22139);
        mpfr_init(r22140);
        mpfr_init(r22141);
        mpfr_init(r22142);
        mpfr_init(r22143);
        mpfr_init(r22144);
        mpfr_init(r22145);
        mpfr_init(r22146);
        mpfr_init(r22147);
        mpfr_init(r22148);
        mpfr_init(r22149);
        mpfr_init(r22150);
        mpfr_init(r22151);
        mpfr_init(r22152);
        mpfr_init_set_str(r22153, "2", 10, MPFR_RNDN);
        mpfr_init(r22154);
}

double f_im(double x, double eps) {
        ;
        mpfr_set_d(r22137, eps, MPFR_RNDN);
        mpfr_div(r22138, r22136, r22137, MPFR_RNDN);
        mpfr_add(r22139, r22136, r22138, MPFR_RNDN);
        mpfr_sub(r22140, r22136, r22137, MPFR_RNDN);
        mpfr_set_d(r22141, x, MPFR_RNDN);
        mpfr_mul(r22142, r22140, r22141, MPFR_RNDN);
        mpfr_neg(r22143, r22142, MPFR_RNDN);
        mpfr_exp(r22144, r22143, MPFR_RNDN);
        mpfr_mul(r22145, r22139, r22144, MPFR_RNDN);
        mpfr_sub(r22146, r22138, r22136, MPFR_RNDN);
        mpfr_add(r22147, r22136, r22137, MPFR_RNDN);
        mpfr_mul(r22148, r22147, r22141, MPFR_RNDN);
        mpfr_neg(r22149, r22148, MPFR_RNDN);
        mpfr_exp(r22150, r22149, MPFR_RNDN);
        mpfr_mul(r22151, r22146, r22150, MPFR_RNDN);
        mpfr_sub(r22152, r22145, r22151, MPFR_RNDN);
        ;
        mpfr_div(r22154, r22152, r22153, MPFR_RNDN);
        return mpfr_get_d(r22154, MPFR_RNDN);
}

static mpfr_t r22155, r22156, r22157, r22158, r22159, r22160, r22161, r22162, r22163, r22164, r22165, r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173, r22174, r22175, r22176, r22177, r22178, r22179;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22155, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22156);
        mpfr_init(r22157);
        mpfr_init(r22158);
        mpfr_init_set_str(r22159, "2", 10, MPFR_RNDN);
        mpfr_init(r22160);
        mpfr_init(r22161);
        mpfr_init(r22162);
        mpfr_init_set_str(r22163, "1405.389636089941", 10, MPFR_RNDN);
        mpfr_init(r22164);
        mpfr_init_set_str(r22165, "1", 10, MPFR_RNDN);
        mpfr_init(r22166);
        mpfr_init(r22167);
        mpfr_init(r22168);
        mpfr_init(r22169);
        mpfr_init(r22170);
        mpfr_init(r22171);
        mpfr_init(r22172);
        mpfr_init(r22173);
        mpfr_init(r22174);
        mpfr_init(r22175);
        mpfr_init(r22176);
        mpfr_init(r22177);
        mpfr_init(r22178);
        mpfr_init(r22179);
}

double f_fm(double x, double eps) {
        ;
        mpfr_set_d(r22156, x, MPFR_RNDN);
        mpfr_mul(r22157, r22155, r22156, MPFR_RNDN);
        mpfr_mul(r22158, r22156, r22156, MPFR_RNDN);
        ;
        mpfr_sub(r22160, r22159, r22158, MPFR_RNDN);
        mpfr_fma(r22161, r22157, r22158, r22160, MPFR_RNDN);
        mpfr_div(r22162, r22161, r22159, MPFR_RNDN);
        ;
        mpfr_set_si(r22164, mpfr_cmp(r22162, r22163) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r22166, eps, MPFR_RNDN);
        mpfr_div(r22167, r22165, r22166, MPFR_RNDN);
        mpfr_add(r22168, r22165, r22167, MPFR_RNDN);
        mpfr_sub(r22169, r22165, r22166, MPFR_RNDN);
        mpfr_neg(r22170, r22156, MPFR_RNDN);
        mpfr_mul(r22171, r22169, r22170, MPFR_RNDN);
        mpfr_exp(r22172, r22171, MPFR_RNDN);
        mpfr_sub(r22173, r22165, r22167, MPFR_RNDN);
        mpfr_fma(r22174, r22166, r22156, r22156, MPFR_RNDN);
        mpfr_exp(r22175, r22174, MPFR_RNDN);
        mpfr_div(r22176, r22173, r22175, MPFR_RNDN);
        mpfr_fma(r22177, r22168, r22172, r22176, MPFR_RNDN);
        mpfr_div(r22178, r22177, r22159, MPFR_RNDN);
        if (mpfr_get_si(r22164, MPFR_RNDN)) { mpfr_set(r22179, r22162, MPFR_RNDN); } else { mpfr_set(r22179, r22178, MPFR_RNDN); };
        return mpfr_get_d(r22179, MPFR_RNDN);
}

static mpfr_t r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187, r22188, r22189, r22190, r22191, r22192, r22193, r22194, r22195, r22196, r22197, r22198, r22199, r22200, r22201, r22202, r22203, r22204;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22180, "2/3", 10, MPFR_RNDN);
        mpfr_init(r22181);
        mpfr_init(r22182);
        mpfr_init(r22183);
        mpfr_init_set_str(r22184, "2", 10, MPFR_RNDN);
        mpfr_init(r22185);
        mpfr_init(r22186);
        mpfr_init(r22187);
        mpfr_init_set_str(r22188, "1405.389636089941", 10, MPFR_RNDN);
        mpfr_init(r22189);
        mpfr_init_set_str(r22190, "1", 10, MPFR_RNDN);
        mpfr_init(r22191);
        mpfr_init(r22192);
        mpfr_init(r22193);
        mpfr_init(r22194);
        mpfr_init(r22195);
        mpfr_init(r22196);
        mpfr_init(r22197);
        mpfr_init(r22198);
        mpfr_init(r22199);
        mpfr_init(r22200);
        mpfr_init(r22201);
        mpfr_init(r22202);
        mpfr_init(r22203);
        mpfr_init(r22204);
}

double f_dm(double x, double eps) {
        ;
        mpfr_set_d(r22181, x, MPFR_RNDN);
        mpfr_mul(r22182, r22180, r22181, MPFR_RNDN);
        mpfr_mul(r22183, r22181, r22181, MPFR_RNDN);
        ;
        mpfr_sub(r22185, r22184, r22183, MPFR_RNDN);
        mpfr_fma(r22186, r22182, r22183, r22185, MPFR_RNDN);
        mpfr_div(r22187, r22186, r22184, MPFR_RNDN);
        ;
        mpfr_set_si(r22189, mpfr_cmp(r22187, r22188) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r22191, eps, MPFR_RNDN);
        mpfr_div(r22192, r22190, r22191, MPFR_RNDN);
        mpfr_add(r22193, r22190, r22192, MPFR_RNDN);
        mpfr_sub(r22194, r22190, r22191, MPFR_RNDN);
        mpfr_neg(r22195, r22181, MPFR_RNDN);
        mpfr_mul(r22196, r22194, r22195, MPFR_RNDN);
        mpfr_exp(r22197, r22196, MPFR_RNDN);
        mpfr_sub(r22198, r22190, r22192, MPFR_RNDN);
        mpfr_fma(r22199, r22191, r22181, r22181, MPFR_RNDN);
        mpfr_exp(r22200, r22199, MPFR_RNDN);
        mpfr_div(r22201, r22198, r22200, MPFR_RNDN);
        mpfr_fma(r22202, r22193, r22197, r22201, MPFR_RNDN);
        mpfr_div(r22203, r22202, r22184, MPFR_RNDN);
        if (mpfr_get_si(r22189, MPFR_RNDN)) { mpfr_set(r22204, r22187, MPFR_RNDN); } else { mpfr_set(r22204, r22203, MPFR_RNDN); };
        return mpfr_get_d(r22204, MPFR_RNDN);
}