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

char *name = "NMSE problem 3.3.2";

double f_if(float x, float eps) {
        float r16188 = x;
        float r16189 = eps;
        float r16190 = r16188 + r16189;
        float r16191 = tan(r16190);
        float r16192 = tan(r16188);
        float r16193 = r16191 - r16192;
        return r16193;
}

double f_id(double x, double eps) {
        double r16194 = x;
        double r16195 = eps;
        double r16196 = r16194 + r16195;
        double r16197 = tan(r16196);
        double r16198 = tan(r16194);
        double r16199 = r16197 - r16198;
        return r16199;
}


double f_of(float x, float eps) {
        float r16200 = eps;
        float r16201 = -0.00010876831220230088f;
        bool r16202 = r16200 <= r16201;
        float r16203 = x;
        float r16204 = r16203 + r16200;
        float r16205 = tan(r16204);
        float r16206 = tan(r16203);
        float r16207 = cbrt(r16206);
        float r16208 = r16207 * (r16207 * r16207);
        float r16209 = r16205 - r16208;
        float r16210 = 3.691105754377588e-10f;
        bool r16211 = r16200 <= r16210;
        float r16212 = 3.0f;
        float r16213 = pow(r16200, r16212);
        float r16214 = r16203 * r16203;
        float r16215 = r16213 * r16214;
        float r16216 = 4.0f;
        float r16217 = pow(r16200, r16216);
        float r16218 = pow(r16203, r16212);
        float r16219 = r16217 * r16218;
        float r16220 = r16215 + r16219;
        float r16221 = r16200 + r16220;
        float r16222 = expm1(r16206);
        float r16223 = cbrt(r16222);
        float r16224 = r16223 * (r16223 * r16223);
        float r16225 = log1p(r16224);
        float r16226 = r16205 - r16225;
        float r16227 = r16211 ? r16221 : r16226;
        float r16228 = r16202 ? r16209 : r16227;
        return r16228;
}

double f_od(double x, double eps) {
        double r16229 = eps;
        double r16230 = -0.00010876831220230088;
        bool r16231 = r16229 <= r16230;
        double r16232 = x;
        double r16233 = r16232 + r16229;
        double r16234 = tan(r16233);
        double r16235 = tan(r16232);
        double r16236 = cbrt(r16235);
        double r16237 = r16236 * (r16236 * r16236);
        double r16238 = r16234 - r16237;
        double r16239 = 3.691105754377588e-10;
        bool r16240 = r16229 <= r16239;
        double r16241 = 3.0;
        double r16242 = pow(r16229, r16241);
        double r16243 = r16232 * r16232;
        double r16244 = r16242 * r16243;
        double r16245 = 4.0;
        double r16246 = pow(r16229, r16245);
        double r16247 = pow(r16232, r16241);
        double r16248 = r16246 * r16247;
        double r16249 = r16244 + r16248;
        double r16250 = r16229 + r16249;
        double r16251 = expm1(r16235);
        double r16252 = cbrt(r16251);
        double r16253 = r16252 * (r16252 * r16252);
        double r16254 = log1p(r16253);
        double r16255 = r16234 - r16254;
        double r16256 = r16240 ? r16250 : r16255;
        double r16257 = r16231 ? r16238 : r16256;
        return r16257;
}

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 r16258, r16259, r16260, r16261, r16262, r16263;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r16258);
        mpfr_init(r16259);
        mpfr_init(r16260);
        mpfr_init(r16261);
        mpfr_init(r16262);
        mpfr_init(r16263);
}

double f_im(double x, double eps) {
        mpfr_set_d(r16258, x, MPFR_RNDN);
        mpfr_set_d(r16259, eps, MPFR_RNDN);
        mpfr_add(r16260, r16258, r16259, MPFR_RNDN);
        mpfr_tan(r16261, r16260, MPFR_RNDN);
        mpfr_tan(r16262, r16258, MPFR_RNDN);
        mpfr_sub(r16263, r16261, r16262, MPFR_RNDN);
        return mpfr_get_d(r16263, MPFR_RNDN);
}

static mpfr_t r16264, r16265, r16266, r16267, r16268, r16269, r16270, r16271, r16272, r16273, r16274, r16275, r16276, r16277, r16278, r16279, r16280, r16281, r16282, r16283, r16284, r16285, r16286, r16287, r16288, r16289, r16290, r16291, r16292;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16264);
        mpfr_init_set_str(r16265, "-0.00010876831f0", 10, MPFR_RNDN);
        mpfr_init(r16266);
        mpfr_init(r16267);
        mpfr_init(r16268);
        mpfr_init(r16269);
        mpfr_init(r16270);
        mpfr_init(r16271);
        mpfr_init(r16272);
        mpfr_init(r16273);
        mpfr_init_set_str(r16274, "3.6911058f-10", 10, MPFR_RNDN);
        mpfr_init(r16275);
        mpfr_init_set_str(r16276, "3", 10, MPFR_RNDN);
        mpfr_init(r16277);
        mpfr_init(r16278);
        mpfr_init(r16279);
        mpfr_init_set_str(r16280, "4", 10, MPFR_RNDN);
        mpfr_init(r16281);
        mpfr_init(r16282);
        mpfr_init(r16283);
        mpfr_init(r16284);
        mpfr_init(r16285);
        mpfr_init(r16286);
        mpfr_init(r16287);
        mpfr_init(r16288);
        mpfr_init(r16289);
        mpfr_init(r16290);
        mpfr_init(r16291);
        mpfr_init(r16292);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r16264, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r16266, mpfr_cmp(r16264, r16265) <= 0, MPFR_RNDN);
        mpfr_set_d(r16267, x, MPFR_RNDN);
        mpfr_add(r16268, r16267, r16264, MPFR_RNDN);
        mpfr_tan(r16269, r16268, MPFR_RNDN);
        mpfr_tan(r16270, r16267, MPFR_RNDN);
        mpfr_cbrt(r16271, r16270, MPFR_RNDN);
        mpfr_mul(r16272, r16271, r16271, MPFR_RNDN); mpfr_mul(r16272, r16272, r16271, MPFR_RNDN);
        mpfr_sub(r16273, r16269, r16272, MPFR_RNDN);
        ;
        mpfr_set_si(r16275, mpfr_cmp(r16264, r16274) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r16277, r16264, r16276, MPFR_RNDN);
        mpfr_sqr(r16278, r16267, MPFR_RNDN);
        mpfr_mul(r16279, r16277, r16278, MPFR_RNDN);
        ;
        mpfr_pow(r16281, r16264, r16280, MPFR_RNDN);
        mpfr_pow(r16282, r16267, r16276, MPFR_RNDN);
        mpfr_mul(r16283, r16281, r16282, MPFR_RNDN);
        mpfr_add(r16284, r16279, r16283, MPFR_RNDN);
        mpfr_add(r16285, r16264, r16284, MPFR_RNDN);
        mpfr_expm1(r16286, r16270, MPFR_RNDN);
        mpfr_cbrt(r16287, r16286, MPFR_RNDN);
        mpfr_mul(r16288, r16287, r16287, MPFR_RNDN); mpfr_mul(r16288, r16288, r16287, MPFR_RNDN);
        mpfr_log1p(r16289, r16288, MPFR_RNDN);
        mpfr_sub(r16290, r16269, r16289, MPFR_RNDN);
        if (mpfr_get_si(r16275, MPFR_RNDN)) { mpfr_set(r16291, r16285, MPFR_RNDN); } else { mpfr_set(r16291, r16290, MPFR_RNDN); };
        if (mpfr_get_si(r16266, MPFR_RNDN)) { mpfr_set(r16292, r16273, MPFR_RNDN); } else { mpfr_set(r16292, r16291, MPFR_RNDN); };
        return mpfr_get_d(r16292, MPFR_RNDN);
}

static mpfr_t r16293, r16294, r16295, r16296, r16297, r16298, r16299, r16300, r16301, r16302, r16303, r16304, r16305, r16306, r16307, r16308, r16309, r16310, r16311, r16312, r16313, r16314, r16315, r16316, r16317, r16318, r16319, r16320, r16321;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r16293);
        mpfr_init_set_str(r16294, "-0.00010876831f0", 10, MPFR_RNDN);
        mpfr_init(r16295);
        mpfr_init(r16296);
        mpfr_init(r16297);
        mpfr_init(r16298);
        mpfr_init(r16299);
        mpfr_init(r16300);
        mpfr_init(r16301);
        mpfr_init(r16302);
        mpfr_init_set_str(r16303, "3.6911058f-10", 10, MPFR_RNDN);
        mpfr_init(r16304);
        mpfr_init_set_str(r16305, "3", 10, MPFR_RNDN);
        mpfr_init(r16306);
        mpfr_init(r16307);
        mpfr_init(r16308);
        mpfr_init_set_str(r16309, "4", 10, MPFR_RNDN);
        mpfr_init(r16310);
        mpfr_init(r16311);
        mpfr_init(r16312);
        mpfr_init(r16313);
        mpfr_init(r16314);
        mpfr_init(r16315);
        mpfr_init(r16316);
        mpfr_init(r16317);
        mpfr_init(r16318);
        mpfr_init(r16319);
        mpfr_init(r16320);
        mpfr_init(r16321);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r16293, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r16295, mpfr_cmp(r16293, r16294) <= 0, MPFR_RNDN);
        mpfr_set_d(r16296, x, MPFR_RNDN);
        mpfr_add(r16297, r16296, r16293, MPFR_RNDN);
        mpfr_tan(r16298, r16297, MPFR_RNDN);
        mpfr_tan(r16299, r16296, MPFR_RNDN);
        mpfr_cbrt(r16300, r16299, MPFR_RNDN);
        mpfr_mul(r16301, r16300, r16300, MPFR_RNDN); mpfr_mul(r16301, r16301, r16300, MPFR_RNDN);
        mpfr_sub(r16302, r16298, r16301, MPFR_RNDN);
        ;
        mpfr_set_si(r16304, mpfr_cmp(r16293, r16303) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r16306, r16293, r16305, MPFR_RNDN);
        mpfr_sqr(r16307, r16296, MPFR_RNDN);
        mpfr_mul(r16308, r16306, r16307, MPFR_RNDN);
        ;
        mpfr_pow(r16310, r16293, r16309, MPFR_RNDN);
        mpfr_pow(r16311, r16296, r16305, MPFR_RNDN);
        mpfr_mul(r16312, r16310, r16311, MPFR_RNDN);
        mpfr_add(r16313, r16308, r16312, MPFR_RNDN);
        mpfr_add(r16314, r16293, r16313, MPFR_RNDN);
        mpfr_expm1(r16315, r16299, MPFR_RNDN);
        mpfr_cbrt(r16316, r16315, MPFR_RNDN);
        mpfr_mul(r16317, r16316, r16316, MPFR_RNDN); mpfr_mul(r16317, r16317, r16316, MPFR_RNDN);
        mpfr_log1p(r16318, r16317, MPFR_RNDN);
        mpfr_sub(r16319, r16298, r16318, MPFR_RNDN);
        if (mpfr_get_si(r16304, MPFR_RNDN)) { mpfr_set(r16320, r16314, MPFR_RNDN); } else { mpfr_set(r16320, r16319, MPFR_RNDN); };
        if (mpfr_get_si(r16295, MPFR_RNDN)) { mpfr_set(r16321, r16302, MPFR_RNDN); } else { mpfr_set(r16321, r16320, MPFR_RNDN); };
        return mpfr_get_d(r16321, MPFR_RNDN);
}