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

char *name = "sintan (problem 3.4.5)";

double f_if(float x) {
        float r21293 = x;
        float r21294 = sin(r21293);
        float r21295 = r21293 - r21294;
        float r21296 = tan(r21293);
        float r21297 = r21293 - r21296;
        float r21298 = r21295 / r21297;
        return r21298;
}

double f_id(double x) {
        double r21299 = x;
        double r21300 = sin(r21299);
        double r21301 = r21299 - r21300;
        double r21302 = tan(r21299);
        double r21303 = r21299 - r21302;
        double r21304 = r21301 / r21303;
        return r21304;
}


double f_of(float x) {
        float r21305 = x;
        float r21306 = -0.03977511220299969;
        bool r21307 = r21305 <= r21306;
        float r21308 = tan(r21305);
        float r21309 = r21305 - r21308;
        float r21310 = r21305 / r21309;
        float r21311 = sin(r21305);
        float r21312 = r21311 / r21309;
        float r21313 = r21310 - r21312;
        float r21314 = 0.03999707453083916;
        bool r21315 = r21305 <= r21314;
        float r21316 = r21305 * r21305;
        float r21317 = 9/40;
        float r21318 = 1/2;
        float r21319 = exp(r21318);
        float r21320 = r21317 / r21319;
        float r21321 = r21316 * r21320;
        float r21322 = 351/22400;
        float r21323 = 4;
        float r21324 = pow(r21305, r21323);
        float r21325 = r21324 / r21319;
        float r21326 = r21322 * r21325;
        float r21327 = r21321 + r21326;
        float r21328 = -r21318;
        float r21329 = exp(r21328);
        float r21330 = r21327 + r21329;
        float r21331 = log(r21330);
        float r21332 = r21315 ? r21331 : r21313;
        float r21333 = r21307 ? r21313 : r21332;
        return r21333;
}

double f_od(double x) {
        double r21334 = x;
        double r21335 = -0.03977511220299969;
        bool r21336 = r21334 <= r21335;
        double r21337 = tan(r21334);
        double r21338 = r21334 - r21337;
        double r21339 = r21334 / r21338;
        double r21340 = sin(r21334);
        double r21341 = r21340 / r21338;
        double r21342 = r21339 - r21341;
        double r21343 = 0.03999707453083916;
        bool r21344 = r21334 <= r21343;
        double r21345 = r21334 * r21334;
        double r21346 = 9/40;
        double r21347 = 1/2;
        double r21348 = exp(r21347);
        double r21349 = r21346 / r21348;
        double r21350 = r21345 * r21349;
        double r21351 = 351/22400;
        double r21352 = 4;
        double r21353 = pow(r21334, r21352);
        double r21354 = r21353 / r21348;
        double r21355 = r21351 * r21354;
        double r21356 = r21350 + r21355;
        double r21357 = -r21347;
        double r21358 = exp(r21357);
        double r21359 = r21356 + r21358;
        double r21360 = log(r21359);
        double r21361 = r21344 ? r21360 : r21342;
        double r21362 = r21336 ? r21342 : r21361;
        return r21362;
}

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 r21363, r21364, r21365, r21366, r21367, r21368;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r21363);
        mpfr_init(r21364);
        mpfr_init(r21365);
        mpfr_init(r21366);
        mpfr_init(r21367);
        mpfr_init(r21368);
}

double f_im(double x) {
        mpfr_set_d(r21363, x, MPFR_RNDN);
        mpfr_sin(r21364, r21363, MPFR_RNDN);
        mpfr_sub(r21365, r21363, r21364, MPFR_RNDN);
        mpfr_tan(r21366, r21363, MPFR_RNDN);
        mpfr_sub(r21367, r21363, r21366, MPFR_RNDN);
        mpfr_div(r21368, r21365, r21367, MPFR_RNDN);
        return mpfr_get_d(r21368, MPFR_RNDN);
}

static mpfr_t r21369, r21370, r21371, r21372, r21373, r21374, r21375, r21376, r21377, r21378, r21379, r21380, r21381, r21382, r21383, r21384, r21385, r21386, r21387, r21388, r21389, r21390, r21391, r21392, r21393, r21394, r21395, r21396, r21397;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r21369);
        mpfr_init_set_str(r21370, "-0.03977511220299969", 10, MPFR_RNDN);
        mpfr_init(r21371);
        mpfr_init(r21372);
        mpfr_init(r21373);
        mpfr_init(r21374);
        mpfr_init(r21375);
        mpfr_init(r21376);
        mpfr_init(r21377);
        mpfr_init_set_str(r21378, "0.03999707453083916", 10, MPFR_RNDN);
        mpfr_init(r21379);
        mpfr_init(r21380);
        mpfr_init_set_str(r21381, "9/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r21382, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21383);
        mpfr_init(r21384);
        mpfr_init(r21385);
        mpfr_init_set_str(r21386, "351/22400", 10, MPFR_RNDN);
        mpfr_init_set_str(r21387, "4", 10, MPFR_RNDN);
        mpfr_init(r21388);
        mpfr_init(r21389);
        mpfr_init(r21390);
        mpfr_init(r21391);
        mpfr_init(r21392);
        mpfr_init(r21393);
        mpfr_init(r21394);
        mpfr_init(r21395);
        mpfr_init(r21396);
        mpfr_init(r21397);
}

double f_fm(double x) {
        mpfr_set_d(r21369, x, MPFR_RNDN);
        ;
        mpfr_set_si(r21371, mpfr_cmp(r21369, r21370) <= 0, MPFR_RNDN);
        mpfr_tan(r21372, r21369, MPFR_RNDN);
        mpfr_sub(r21373, r21369, r21372, MPFR_RNDN);
        mpfr_div(r21374, r21369, r21373, MPFR_RNDN);
        mpfr_sin(r21375, r21369, MPFR_RNDN);
        mpfr_div(r21376, r21375, r21373, MPFR_RNDN);
        mpfr_sub(r21377, r21374, r21376, MPFR_RNDN);
        ;
        mpfr_set_si(r21379, mpfr_cmp(r21369, r21378) <= 0, MPFR_RNDN);
        mpfr_mul(r21380, r21369, r21369, MPFR_RNDN);
        ;
        ;
        mpfr_exp(r21383, r21382, MPFR_RNDN);
        mpfr_div(r21384, r21381, r21383, MPFR_RNDN);
        mpfr_mul(r21385, r21380, r21384, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r21388, r21369, r21387, MPFR_RNDN);
        mpfr_div(r21389, r21388, r21383, MPFR_RNDN);
        mpfr_mul(r21390, r21386, r21389, MPFR_RNDN);
        mpfr_add(r21391, r21385, r21390, MPFR_RNDN);
        mpfr_neg(r21392, r21382, MPFR_RNDN);
        mpfr_exp(r21393, r21392, MPFR_RNDN);
        mpfr_add(r21394, r21391, r21393, MPFR_RNDN);
        mpfr_log(r21395, r21394, MPFR_RNDN);
        if (mpfr_get_si(r21379, MPFR_RNDN)) { mpfr_set(r21396, r21395, MPFR_RNDN); } else { mpfr_set(r21396, r21377, MPFR_RNDN); };
        if (mpfr_get_si(r21371, MPFR_RNDN)) { mpfr_set(r21397, r21377, MPFR_RNDN); } else { mpfr_set(r21397, r21396, MPFR_RNDN); };
        return mpfr_get_d(r21397, MPFR_RNDN);
}

static mpfr_t r21398, r21399, r21400, r21401, r21402, r21403, r21404, r21405, r21406, r21407, r21408, r21409, r21410, r21411, r21412, r21413, r21414, r21415, r21416, r21417, r21418, r21419, r21420, r21421, r21422, r21423, r21424, r21425, r21426;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r21398);
        mpfr_init_set_str(r21399, "-0.03977511220299969", 10, MPFR_RNDN);
        mpfr_init(r21400);
        mpfr_init(r21401);
        mpfr_init(r21402);
        mpfr_init(r21403);
        mpfr_init(r21404);
        mpfr_init(r21405);
        mpfr_init(r21406);
        mpfr_init_set_str(r21407, "0.03999707453083916", 10, MPFR_RNDN);
        mpfr_init(r21408);
        mpfr_init(r21409);
        mpfr_init_set_str(r21410, "9/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r21411, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21412);
        mpfr_init(r21413);
        mpfr_init(r21414);
        mpfr_init_set_str(r21415, "351/22400", 10, MPFR_RNDN);
        mpfr_init_set_str(r21416, "4", 10, MPFR_RNDN);
        mpfr_init(r21417);
        mpfr_init(r21418);
        mpfr_init(r21419);
        mpfr_init(r21420);
        mpfr_init(r21421);
        mpfr_init(r21422);
        mpfr_init(r21423);
        mpfr_init(r21424);
        mpfr_init(r21425);
        mpfr_init(r21426);
}

double f_dm(double x) {
        mpfr_set_d(r21398, x, MPFR_RNDN);
        ;
        mpfr_set_si(r21400, mpfr_cmp(r21398, r21399) <= 0, MPFR_RNDN);
        mpfr_tan(r21401, r21398, MPFR_RNDN);
        mpfr_sub(r21402, r21398, r21401, MPFR_RNDN);
        mpfr_div(r21403, r21398, r21402, MPFR_RNDN);
        mpfr_sin(r21404, r21398, MPFR_RNDN);
        mpfr_div(r21405, r21404, r21402, MPFR_RNDN);
        mpfr_sub(r21406, r21403, r21405, MPFR_RNDN);
        ;
        mpfr_set_si(r21408, mpfr_cmp(r21398, r21407) <= 0, MPFR_RNDN);
        mpfr_mul(r21409, r21398, r21398, MPFR_RNDN);
        ;
        ;
        mpfr_exp(r21412, r21411, MPFR_RNDN);
        mpfr_div(r21413, r21410, r21412, MPFR_RNDN);
        mpfr_mul(r21414, r21409, r21413, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r21417, r21398, r21416, MPFR_RNDN);
        mpfr_div(r21418, r21417, r21412, MPFR_RNDN);
        mpfr_mul(r21419, r21415, r21418, MPFR_RNDN);
        mpfr_add(r21420, r21414, r21419, MPFR_RNDN);
        mpfr_neg(r21421, r21411, MPFR_RNDN);
        mpfr_exp(r21422, r21421, MPFR_RNDN);
        mpfr_add(r21423, r21420, r21422, MPFR_RNDN);
        mpfr_log(r21424, r21423, MPFR_RNDN);
        if (mpfr_get_si(r21408, MPFR_RNDN)) { mpfr_set(r21425, r21424, MPFR_RNDN); } else { mpfr_set(r21425, r21406, MPFR_RNDN); };
        if (mpfr_get_si(r21400, MPFR_RNDN)) { mpfr_set(r21426, r21406, MPFR_RNDN); } else { mpfr_set(r21426, r21425, MPFR_RNDN); };
        return mpfr_get_d(r21426, MPFR_RNDN);
}