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

char *name = "2tan (problem 3.3.2)";

double f_if(float x, float eps) {
        float r24226 = x;
        float r24227 = eps;
        float r24228 = r24226 + r24227;
        float r24229 = tan(r24228);
        float r24230 = tan(r24226);
        float r24231 = r24229 - r24230;
        return r24231;
}

double f_id(double x, double eps) {
        double r24232 = x;
        double r24233 = eps;
        double r24234 = r24232 + r24233;
        double r24235 = tan(r24234);
        double r24236 = tan(r24232);
        double r24237 = r24235 - r24236;
        return r24237;
}


double f_of(float x, float eps) {
        float r24238 = x;
        float r24239 = eps;
        float r24240 = r24238 * r24239;
        float r24241 = fma(r24240, r24239, r24239);
        float r24242 = fma(r24240, r24241, r24239);
        float r24243 = -3.1068469719484584e-28;
        bool r24244 = r24242 <= r24243;
        float r24245 = cos(r24238);
        float r24246 = tan(r24238);
        float r24247 = tan(r24239);
        float r24248 = r24246 + r24247;
        float r24249 = sin(r24238);
        float r24250 = r24247 * r24246;
        float r24251 = -r24249;
        float r24252 = fma(r24249, r24250, r24251);
        float r24253 = fma(r24245, r24248, r24252);
        float r24254 = r24247 * r24245;
        float r24255 = -r24246;
        float r24256 = fma(r24254, r24255, r24245);
        float r24257 = r24253 / r24256;
        float r24258 = 5.858083040047354e-17;
        bool r24259 = r24242 <= r24258;
        float r24260 = r24259 ? r24242 : r24257;
        float r24261 = r24244 ? r24257 : r24260;
        return r24261;
}

double f_od(double x, double eps) {
        double r24262 = x;
        double r24263 = eps;
        double r24264 = r24262 * r24263;
        double r24265 = fma(r24264, r24263, r24263);
        double r24266 = fma(r24264, r24265, r24263);
        double r24267 = -3.1068469719484584e-28;
        bool r24268 = r24266 <= r24267;
        double r24269 = cos(r24262);
        double r24270 = tan(r24262);
        double r24271 = tan(r24263);
        double r24272 = r24270 + r24271;
        double r24273 = sin(r24262);
        double r24274 = r24271 * r24270;
        double r24275 = -r24273;
        double r24276 = fma(r24273, r24274, r24275);
        double r24277 = fma(r24269, r24272, r24276);
        double r24278 = r24271 * r24269;
        double r24279 = -r24270;
        double r24280 = fma(r24278, r24279, r24269);
        double r24281 = r24277 / r24280;
        double r24282 = 5.858083040047354e-17;
        bool r24283 = r24266 <= r24282;
        double r24284 = r24283 ? r24266 : r24281;
        double r24285 = r24268 ? r24281 : r24284;
        return r24285;
}

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 r24286, r24287, r24288, r24289, r24290, r24291;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r24286);
        mpfr_init(r24287);
        mpfr_init(r24288);
        mpfr_init(r24289);
        mpfr_init(r24290);
        mpfr_init(r24291);
}

double f_im(double x, double eps) {
        mpfr_set_d(r24286, x, MPFR_RNDN);
        mpfr_set_d(r24287, eps, MPFR_RNDN);
        mpfr_add(r24288, r24286, r24287, MPFR_RNDN);
        mpfr_tan(r24289, r24288, MPFR_RNDN);
        mpfr_tan(r24290, r24286, MPFR_RNDN);
        mpfr_sub(r24291, r24289, r24290, MPFR_RNDN);
        return mpfr_get_d(r24291, MPFR_RNDN);
}

static mpfr_t r24292, r24293, r24294, r24295, r24296, r24297, r24298, r24299, r24300, r24301, r24302, r24303, r24304, r24305, r24306, r24307, r24308, r24309, r24310, r24311, r24312, r24313, r24314, r24315;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r24292);
        mpfr_init(r24293);
        mpfr_init(r24294);
        mpfr_init(r24295);
        mpfr_init(r24296);
        mpfr_init_set_str(r24297, "-3.1068469719484584e-28", 10, MPFR_RNDN);
        mpfr_init(r24298);
        mpfr_init(r24299);
        mpfr_init(r24300);
        mpfr_init(r24301);
        mpfr_init(r24302);
        mpfr_init(r24303);
        mpfr_init(r24304);
        mpfr_init(r24305);
        mpfr_init(r24306);
        mpfr_init(r24307);
        mpfr_init(r24308);
        mpfr_init(r24309);
        mpfr_init(r24310);
        mpfr_init(r24311);
        mpfr_init_set_str(r24312, "5.858083040047354e-17", 10, MPFR_RNDN);
        mpfr_init(r24313);
        mpfr_init(r24314);
        mpfr_init(r24315);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r24292, x, MPFR_RNDN);
        mpfr_set_d(r24293, eps, MPFR_RNDN);
        mpfr_mul(r24294, r24292, r24293, MPFR_RNDN);
        mpfr_fma(r24295, r24294, r24293, r24293, MPFR_RNDN);
        mpfr_fma(r24296, r24294, r24295, r24293, MPFR_RNDN);
        ;
        mpfr_set_si(r24298, mpfr_cmp(r24296, r24297) <= 0, MPFR_RNDN);
        mpfr_cos(r24299, r24292, MPFR_RNDN);
        mpfr_tan(r24300, r24292, MPFR_RNDN);
        mpfr_tan(r24301, r24293, MPFR_RNDN);
        mpfr_add(r24302, r24300, r24301, MPFR_RNDN);
        mpfr_sin(r24303, r24292, MPFR_RNDN);
        mpfr_mul(r24304, r24301, r24300, MPFR_RNDN);
        mpfr_neg(r24305, r24303, MPFR_RNDN);
        mpfr_fma(r24306, r24303, r24304, r24305, MPFR_RNDN);
        mpfr_fma(r24307, r24299, r24302, r24306, MPFR_RNDN);
        mpfr_mul(r24308, r24301, r24299, MPFR_RNDN);
        mpfr_neg(r24309, r24300, MPFR_RNDN);
        mpfr_fma(r24310, r24308, r24309, r24299, MPFR_RNDN);
        mpfr_div(r24311, r24307, r24310, MPFR_RNDN);
        ;
        mpfr_set_si(r24313, mpfr_cmp(r24296, r24312) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24313, MPFR_RNDN)) { mpfr_set(r24314, r24296, MPFR_RNDN); } else { mpfr_set(r24314, r24311, MPFR_RNDN); };
        if (mpfr_get_si(r24298, MPFR_RNDN)) { mpfr_set(r24315, r24311, MPFR_RNDN); } else { mpfr_set(r24315, r24314, MPFR_RNDN); };
        return mpfr_get_d(r24315, MPFR_RNDN);
}

static mpfr_t r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332, r24333, r24334, r24335, r24336, r24337, r24338, r24339;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r24316);
        mpfr_init(r24317);
        mpfr_init(r24318);
        mpfr_init(r24319);
        mpfr_init(r24320);
        mpfr_init_set_str(r24321, "-3.1068469719484584e-28", 10, MPFR_RNDN);
        mpfr_init(r24322);
        mpfr_init(r24323);
        mpfr_init(r24324);
        mpfr_init(r24325);
        mpfr_init(r24326);
        mpfr_init(r24327);
        mpfr_init(r24328);
        mpfr_init(r24329);
        mpfr_init(r24330);
        mpfr_init(r24331);
        mpfr_init(r24332);
        mpfr_init(r24333);
        mpfr_init(r24334);
        mpfr_init(r24335);
        mpfr_init_set_str(r24336, "5.858083040047354e-17", 10, MPFR_RNDN);
        mpfr_init(r24337);
        mpfr_init(r24338);
        mpfr_init(r24339);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r24316, x, MPFR_RNDN);
        mpfr_set_d(r24317, eps, MPFR_RNDN);
        mpfr_mul(r24318, r24316, r24317, MPFR_RNDN);
        mpfr_fma(r24319, r24318, r24317, r24317, MPFR_RNDN);
        mpfr_fma(r24320, r24318, r24319, r24317, MPFR_RNDN);
        ;
        mpfr_set_si(r24322, mpfr_cmp(r24320, r24321) <= 0, MPFR_RNDN);
        mpfr_cos(r24323, r24316, MPFR_RNDN);
        mpfr_tan(r24324, r24316, MPFR_RNDN);
        mpfr_tan(r24325, r24317, MPFR_RNDN);
        mpfr_add(r24326, r24324, r24325, MPFR_RNDN);
        mpfr_sin(r24327, r24316, MPFR_RNDN);
        mpfr_mul(r24328, r24325, r24324, MPFR_RNDN);
        mpfr_neg(r24329, r24327, MPFR_RNDN);
        mpfr_fma(r24330, r24327, r24328, r24329, MPFR_RNDN);
        mpfr_fma(r24331, r24323, r24326, r24330, MPFR_RNDN);
        mpfr_mul(r24332, r24325, r24323, MPFR_RNDN);
        mpfr_neg(r24333, r24324, MPFR_RNDN);
        mpfr_fma(r24334, r24332, r24333, r24323, MPFR_RNDN);
        mpfr_div(r24335, r24331, r24334, MPFR_RNDN);
        ;
        mpfr_set_si(r24337, mpfr_cmp(r24320, r24336) <= 0, MPFR_RNDN);
        if (mpfr_get_si(r24337, MPFR_RNDN)) { mpfr_set(r24338, r24320, MPFR_RNDN); } else { mpfr_set(r24338, r24335, MPFR_RNDN); };
        if (mpfr_get_si(r24322, MPFR_RNDN)) { mpfr_set(r24339, r24335, MPFR_RNDN); } else { mpfr_set(r24339, r24338, MPFR_RNDN); };
        return mpfr_get_d(r24339, MPFR_RNDN);
}