#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 r20219 = x;
        float r20220 = eps;
        float r20221 = r20219 + r20220;
        float r20222 = tan(r20221);
        float r20223 = tan(r20219);
        float r20224 = r20222 - r20223;
        return r20224;
}

double f_id(double x, double eps) {
        double r20225 = x;
        double r20226 = eps;
        double r20227 = r20225 + r20226;
        double r20228 = tan(r20227);
        double r20229 = tan(r20225);
        double r20230 = r20228 - r20229;
        return r20230;
}


double f_of(float x, float eps) {
        float r20231 = eps;
        float r20232 = -1.9479733961868206e-39f;
        bool r20233 = r20231 <= r20232;
        float r20234 = x;
        float r20235 = tan(r20234);
        float r20236 = tan(r20231);
        float r20237 = r20235 + r20236;
        float r20238 = 1.0f;
        float r20239 = sin(r20234);
        float r20240 = sin(r20231);
        float r20241 = r20239 * r20240;
        float r20242 = cos(r20234);
        float r20243 = cos(r20231);
        float r20244 = r20242 * r20243;
        float r20245 = r20241 / r20244;
        float r20246 = r20238 - r20245;
        float r20247 = r20237 / r20246;
        float r20248 = r20247 - r20235;
        float r20249 = 5.496498246559763e-32f;
        bool r20250 = r20231 <= r20249;
        float r20251 = r20234 * r20234;
        float r20252 = r20231 * (r20231 * r20231);
        float r20253 = r20251 * r20252;
        float r20254 = 4.0f;
        float r20255 = pow(r20231, r20254);
        float r20256 = r20234 * (r20234 * r20234);
        float r20257 = r20255 * r20256;
        float r20258 = r20253 + r20257;
        float r20259 = r20258 + r20231;
        float r20260 = r20235 * r20235;
        float r20261 = r20236 * r20236;
        float r20262 = r20260 - r20261;
        float r20263 = r20235 * r20236;
        float r20264 = r20238 - r20263;
        float r20265 = r20235 - r20236;
        float r20266 = r20264 * r20265;
        float r20267 = r20262 / r20266;
        float r20268 = r20267 - r20235;
        float r20269 = r20250 ? r20259 : r20268;
        float r20270 = r20233 ? r20248 : r20269;
        return r20270;
}

double f_od(double x, double eps) {
        double r20271 = eps;
        double r20272 = -1.9479733961868206e-39;
        bool r20273 = r20271 <= r20272;
        double r20274 = x;
        double r20275 = tan(r20274);
        double r20276 = tan(r20271);
        double r20277 = r20275 + r20276;
        double r20278 = 1.0;
        double r20279 = sin(r20274);
        double r20280 = sin(r20271);
        double r20281 = r20279 * r20280;
        double r20282 = cos(r20274);
        double r20283 = cos(r20271);
        double r20284 = r20282 * r20283;
        double r20285 = r20281 / r20284;
        double r20286 = r20278 - r20285;
        double r20287 = r20277 / r20286;
        double r20288 = r20287 - r20275;
        double r20289 = 5.496498246559763e-32;
        bool r20290 = r20271 <= r20289;
        double r20291 = r20274 * r20274;
        double r20292 = r20271 * (r20271 * r20271);
        double r20293 = r20291 * r20292;
        double r20294 = 4.0;
        double r20295 = pow(r20271, r20294);
        double r20296 = r20274 * (r20274 * r20274);
        double r20297 = r20295 * r20296;
        double r20298 = r20293 + r20297;
        double r20299 = r20298 + r20271;
        double r20300 = r20275 * r20275;
        double r20301 = r20276 * r20276;
        double r20302 = r20300 - r20301;
        double r20303 = r20275 * r20276;
        double r20304 = r20278 - r20303;
        double r20305 = r20275 - r20276;
        double r20306 = r20304 * r20305;
        double r20307 = r20302 / r20306;
        double r20308 = r20307 - r20275;
        double r20309 = r20290 ? r20299 : r20308;
        double r20310 = r20273 ? r20288 : r20309;
        return r20310;
}

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 r20311, r20312, r20313, r20314, r20315, r20316;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20311);
        mpfr_init(r20312);
        mpfr_init(r20313);
        mpfr_init(r20314);
        mpfr_init(r20315);
        mpfr_init(r20316);
}

double f_im(double x, double eps) {
        mpfr_set_d(r20311, x, MPFR_RNDN);
        mpfr_set_d(r20312, eps, MPFR_RNDN);
        mpfr_add(r20313, r20311, r20312, MPFR_RNDN);
        mpfr_tan(r20314, r20313, MPFR_RNDN);
        mpfr_tan(r20315, r20311, MPFR_RNDN);
        mpfr_sub(r20316, r20314, r20315, MPFR_RNDN);
        return mpfr_get_d(r20316, MPFR_RNDN);
}

static mpfr_t r20317, r20318, r20319, r20320, r20321, r20322, r20323, r20324, r20325, r20326, r20327, r20328, r20329, r20330, r20331, r20332, r20333, r20334, r20335, r20336, r20337, r20338, r20339, r20340, r20341, r20342, r20343, r20344, r20345, r20346, r20347, r20348, r20349, r20350, r20351, r20352, r20353, r20354, r20355, r20356;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20317);
        mpfr_init_set_str(r20318, "-1.9479733961868206e-39", 10, MPFR_RNDN);
        mpfr_init(r20319);
        mpfr_init(r20320);
        mpfr_init(r20321);
        mpfr_init(r20322);
        mpfr_init(r20323);
        mpfr_init_set_str(r20324, "1", 10, MPFR_RNDN);
        mpfr_init(r20325);
        mpfr_init(r20326);
        mpfr_init(r20327);
        mpfr_init(r20328);
        mpfr_init(r20329);
        mpfr_init(r20330);
        mpfr_init(r20331);
        mpfr_init(r20332);
        mpfr_init(r20333);
        mpfr_init(r20334);
        mpfr_init_set_str(r20335, "5.496498246559763e-32", 10, MPFR_RNDN);
        mpfr_init(r20336);
        mpfr_init(r20337);
        mpfr_init(r20338);
        mpfr_init(r20339);
        mpfr_init_set_str(r20340, "4", 10, MPFR_RNDN);
        mpfr_init(r20341);
        mpfr_init(r20342);
        mpfr_init(r20343);
        mpfr_init(r20344);
        mpfr_init(r20345);
        mpfr_init(r20346);
        mpfr_init(r20347);
        mpfr_init(r20348);
        mpfr_init(r20349);
        mpfr_init(r20350);
        mpfr_init(r20351);
        mpfr_init(r20352);
        mpfr_init(r20353);
        mpfr_init(r20354);
        mpfr_init(r20355);
        mpfr_init(r20356);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r20317, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20319, mpfr_cmp(r20317, r20318) <= 0, MPFR_RNDN);
        mpfr_set_d(r20320, x, MPFR_RNDN);
        mpfr_tan(r20321, r20320, MPFR_RNDN);
        mpfr_tan(r20322, r20317, MPFR_RNDN);
        mpfr_add(r20323, r20321, r20322, MPFR_RNDN);
        ;
        mpfr_sin(r20325, r20320, MPFR_RNDN);
        mpfr_sin(r20326, r20317, MPFR_RNDN);
        mpfr_mul(r20327, r20325, r20326, MPFR_RNDN);
        mpfr_cos(r20328, r20320, MPFR_RNDN);
        mpfr_cos(r20329, r20317, MPFR_RNDN);
        mpfr_mul(r20330, r20328, r20329, MPFR_RNDN);
        mpfr_div(r20331, r20327, r20330, MPFR_RNDN);
        mpfr_sub(r20332, r20324, r20331, MPFR_RNDN);
        mpfr_div(r20333, r20323, r20332, MPFR_RNDN);
        mpfr_sub(r20334, r20333, r20321, MPFR_RNDN);
        ;
        mpfr_set_si(r20336, mpfr_cmp(r20317, r20335) <= 0, MPFR_RNDN);
        mpfr_sqr(r20337, r20320, MPFR_RNDN);
        mpfr_mul(r20338, r20317, r20317, MPFR_RNDN); mpfr_mul(r20338, r20338, r20317, MPFR_RNDN);
        mpfr_mul(r20339, r20337, r20338, MPFR_RNDN);
        ;
        mpfr_pow(r20341, r20317, r20340, MPFR_RNDN);
        mpfr_mul(r20342, r20320, r20320, MPFR_RNDN); mpfr_mul(r20342, r20342, r20320, MPFR_RNDN);
        mpfr_mul(r20343, r20341, r20342, MPFR_RNDN);
        mpfr_add(r20344, r20339, r20343, MPFR_RNDN);
        mpfr_add(r20345, r20344, r20317, MPFR_RNDN);
        mpfr_sqr(r20346, r20321, MPFR_RNDN);
        mpfr_sqr(r20347, r20322, MPFR_RNDN);
        mpfr_sub(r20348, r20346, r20347, MPFR_RNDN);
        mpfr_mul(r20349, r20321, r20322, MPFR_RNDN);
        mpfr_sub(r20350, r20324, r20349, MPFR_RNDN);
        mpfr_sub(r20351, r20321, r20322, MPFR_RNDN);
        mpfr_mul(r20352, r20350, r20351, MPFR_RNDN);
        mpfr_div(r20353, r20348, r20352, MPFR_RNDN);
        mpfr_sub(r20354, r20353, r20321, MPFR_RNDN);
        if (mpfr_get_si(r20336, MPFR_RNDN)) { mpfr_set(r20355, r20345, MPFR_RNDN); } else { mpfr_set(r20355, r20354, MPFR_RNDN); };
        if (mpfr_get_si(r20319, MPFR_RNDN)) { mpfr_set(r20356, r20334, MPFR_RNDN); } else { mpfr_set(r20356, r20355, MPFR_RNDN); };
        return mpfr_get_d(r20356, MPFR_RNDN);
}

static mpfr_t r20357, r20358, r20359, r20360, r20361, r20362, r20363, r20364, r20365, r20366, r20367, r20368, r20369, r20370, r20371, r20372, r20373, r20374, r20375, r20376, r20377, r20378, r20379, r20380, r20381, r20382, r20383, r20384, r20385, r20386, r20387, r20388, r20389, r20390, r20391, r20392, r20393, r20394, r20395, r20396;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20357);
        mpfr_init_set_str(r20358, "-1.9479733961868206e-39", 10, MPFR_RNDN);
        mpfr_init(r20359);
        mpfr_init(r20360);
        mpfr_init(r20361);
        mpfr_init(r20362);
        mpfr_init(r20363);
        mpfr_init_set_str(r20364, "1", 10, MPFR_RNDN);
        mpfr_init(r20365);
        mpfr_init(r20366);
        mpfr_init(r20367);
        mpfr_init(r20368);
        mpfr_init(r20369);
        mpfr_init(r20370);
        mpfr_init(r20371);
        mpfr_init(r20372);
        mpfr_init(r20373);
        mpfr_init(r20374);
        mpfr_init_set_str(r20375, "5.496498246559763e-32", 10, MPFR_RNDN);
        mpfr_init(r20376);
        mpfr_init(r20377);
        mpfr_init(r20378);
        mpfr_init(r20379);
        mpfr_init_set_str(r20380, "4", 10, MPFR_RNDN);
        mpfr_init(r20381);
        mpfr_init(r20382);
        mpfr_init(r20383);
        mpfr_init(r20384);
        mpfr_init(r20385);
        mpfr_init(r20386);
        mpfr_init(r20387);
        mpfr_init(r20388);
        mpfr_init(r20389);
        mpfr_init(r20390);
        mpfr_init(r20391);
        mpfr_init(r20392);
        mpfr_init(r20393);
        mpfr_init(r20394);
        mpfr_init(r20395);
        mpfr_init(r20396);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r20357, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20359, mpfr_cmp(r20357, r20358) <= 0, MPFR_RNDN);
        mpfr_set_d(r20360, x, MPFR_RNDN);
        mpfr_tan(r20361, r20360, MPFR_RNDN);
        mpfr_tan(r20362, r20357, MPFR_RNDN);
        mpfr_add(r20363, r20361, r20362, MPFR_RNDN);
        ;
        mpfr_sin(r20365, r20360, MPFR_RNDN);
        mpfr_sin(r20366, r20357, MPFR_RNDN);
        mpfr_mul(r20367, r20365, r20366, MPFR_RNDN);
        mpfr_cos(r20368, r20360, MPFR_RNDN);
        mpfr_cos(r20369, r20357, MPFR_RNDN);
        mpfr_mul(r20370, r20368, r20369, MPFR_RNDN);
        mpfr_div(r20371, r20367, r20370, MPFR_RNDN);
        mpfr_sub(r20372, r20364, r20371, MPFR_RNDN);
        mpfr_div(r20373, r20363, r20372, MPFR_RNDN);
        mpfr_sub(r20374, r20373, r20361, MPFR_RNDN);
        ;
        mpfr_set_si(r20376, mpfr_cmp(r20357, r20375) <= 0, MPFR_RNDN);
        mpfr_sqr(r20377, r20360, MPFR_RNDN);
        mpfr_mul(r20378, r20357, r20357, MPFR_RNDN); mpfr_mul(r20378, r20378, r20357, MPFR_RNDN);
        mpfr_mul(r20379, r20377, r20378, MPFR_RNDN);
        ;
        mpfr_pow(r20381, r20357, r20380, MPFR_RNDN);
        mpfr_mul(r20382, r20360, r20360, MPFR_RNDN); mpfr_mul(r20382, r20382, r20360, MPFR_RNDN);
        mpfr_mul(r20383, r20381, r20382, MPFR_RNDN);
        mpfr_add(r20384, r20379, r20383, MPFR_RNDN);
        mpfr_add(r20385, r20384, r20357, MPFR_RNDN);
        mpfr_sqr(r20386, r20361, MPFR_RNDN);
        mpfr_sqr(r20387, r20362, MPFR_RNDN);
        mpfr_sub(r20388, r20386, r20387, MPFR_RNDN);
        mpfr_mul(r20389, r20361, r20362, MPFR_RNDN);
        mpfr_sub(r20390, r20364, r20389, MPFR_RNDN);
        mpfr_sub(r20391, r20361, r20362, MPFR_RNDN);
        mpfr_mul(r20392, r20390, r20391, MPFR_RNDN);
        mpfr_div(r20393, r20388, r20392, MPFR_RNDN);
        mpfr_sub(r20394, r20393, r20361, MPFR_RNDN);
        if (mpfr_get_si(r20376, MPFR_RNDN)) { mpfr_set(r20395, r20385, MPFR_RNDN); } else { mpfr_set(r20395, r20394, MPFR_RNDN); };
        if (mpfr_get_si(r20359, MPFR_RNDN)) { mpfr_set(r20396, r20374, MPFR_RNDN); } else { mpfr_set(r20396, r20395, MPFR_RNDN); };
        return mpfr_get_d(r20396, MPFR_RNDN);
}

