#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 r26933 = x;
        float r26934 = eps;
        float r26935 = r26933 + r26934;
        float r26936 = tan(r26935);
        float r26937 = tan(r26933);
        float r26938 = r26936 - r26937;
        return r26938;
}

double f_id(double x, double eps) {
        double r26939 = x;
        double r26940 = eps;
        double r26941 = r26939 + r26940;
        double r26942 = tan(r26941);
        double r26943 = tan(r26939);
        double r26944 = r26942 - r26943;
        return r26944;
}


double f_of(float x, float eps) {
        float r26945 = x;
        float r26946 = sin(r26945);
        float r26947 = cos(r26945);
        float r26948 = r26946 / r26947;
        float r26949 = r26948 * r26948;
        float r26950 = 1;
        float r26951 = eps;
        float r26952 = sin(r26951);
        float r26953 = r26952 / r26947;
        float r26954 = r26953 * r26953;
        float r26955 = cos(r26951);
        float r26956 = r26946 / r26955;
        float r26957 = r26956 * r26956;
        float r26958 = r26954 * r26957;
        float r26959 = r26950 - r26958;
        float r26960 = r26949 / r26959;
        float r26961 = r26952 / r26955;
        float r26962 = r26961 / r26959;
        float r26963 = fma(r26960, r26961, r26962);
        float r26964 = fma(r26961, r26961, r26950);
        float r26965 = r26964 * r26948;
        float r26966 = r26965 / r26959;
        float r26967 = r26966 - r26948;
        float r26968 = r26963 + r26967;
        return r26968;
}

double f_od(double x, double eps) {
        double r26969 = x;
        double r26970 = sin(r26969);
        double r26971 = cos(r26969);
        double r26972 = r26970 / r26971;
        double r26973 = r26972 * r26972;
        double r26974 = 1;
        double r26975 = eps;
        double r26976 = sin(r26975);
        double r26977 = r26976 / r26971;
        double r26978 = r26977 * r26977;
        double r26979 = cos(r26975);
        double r26980 = r26970 / r26979;
        double r26981 = r26980 * r26980;
        double r26982 = r26978 * r26981;
        double r26983 = r26974 - r26982;
        double r26984 = r26973 / r26983;
        double r26985 = r26976 / r26979;
        double r26986 = r26985 / r26983;
        double r26987 = fma(r26984, r26985, r26986);
        double r26988 = fma(r26985, r26985, r26974);
        double r26989 = r26988 * r26972;
        double r26990 = r26989 / r26983;
        double r26991 = r26990 - r26972;
        double r26992 = r26987 + r26991;
        return r26992;
}

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 r26993, r26994, r26995, r26996, r26997, r26998;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r26993);
        mpfr_init(r26994);
        mpfr_init(r26995);
        mpfr_init(r26996);
        mpfr_init(r26997);
        mpfr_init(r26998);
}

double f_im(double x, double eps) {
        mpfr_set_d(r26993, x, MPFR_RNDN);
        mpfr_set_d(r26994, eps, MPFR_RNDN);
        mpfr_add(r26995, r26993, r26994, MPFR_RNDN);
        mpfr_tan(r26996, r26995, MPFR_RNDN);
        mpfr_tan(r26997, r26993, MPFR_RNDN);
        mpfr_sub(r26998, r26996, r26997, MPFR_RNDN);
        return mpfr_get_d(r26998, MPFR_RNDN);
}

static mpfr_t r26999, r27000, r27001, r27002, r27003, r27004, r27005, r27006, r27007, r27008, r27009, r27010, r27011, r27012, r27013, r27014, r27015, r27016, r27017, r27018, r27019, r27020, r27021, r27022;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r26999);
        mpfr_init(r27000);
        mpfr_init(r27001);
        mpfr_init(r27002);
        mpfr_init(r27003);
        mpfr_init_set_str(r27004, "1", 10, MPFR_RNDN);
        mpfr_init(r27005);
        mpfr_init(r27006);
        mpfr_init(r27007);
        mpfr_init(r27008);
        mpfr_init(r27009);
        mpfr_init(r27010);
        mpfr_init(r27011);
        mpfr_init(r27012);
        mpfr_init(r27013);
        mpfr_init(r27014);
        mpfr_init(r27015);
        mpfr_init(r27016);
        mpfr_init(r27017);
        mpfr_init(r27018);
        mpfr_init(r27019);
        mpfr_init(r27020);
        mpfr_init(r27021);
        mpfr_init(r27022);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r26999, x, MPFR_RNDN);
        mpfr_sin(r27000, r26999, MPFR_RNDN);
        mpfr_cos(r27001, r26999, MPFR_RNDN);
        mpfr_div(r27002, r27000, r27001, MPFR_RNDN);
        mpfr_mul(r27003, r27002, r27002, MPFR_RNDN);
        ;
        mpfr_set_d(r27005, eps, MPFR_RNDN);
        mpfr_sin(r27006, r27005, MPFR_RNDN);
        mpfr_div(r27007, r27006, r27001, MPFR_RNDN);
        mpfr_mul(r27008, r27007, r27007, MPFR_RNDN);
        mpfr_cos(r27009, r27005, MPFR_RNDN);
        mpfr_div(r27010, r27000, r27009, MPFR_RNDN);
        mpfr_mul(r27011, r27010, r27010, MPFR_RNDN);
        mpfr_mul(r27012, r27008, r27011, MPFR_RNDN);
        mpfr_sub(r27013, r27004, r27012, MPFR_RNDN);
        mpfr_div(r27014, r27003, r27013, MPFR_RNDN);
        mpfr_div(r27015, r27006, r27009, MPFR_RNDN);
        mpfr_div(r27016, r27015, r27013, MPFR_RNDN);
        mpfr_fma(r27017, r27014, r27015, r27016, MPFR_RNDN);
        mpfr_fma(r27018, r27015, r27015, r27004, MPFR_RNDN);
        mpfr_mul(r27019, r27018, r27002, MPFR_RNDN);
        mpfr_div(r27020, r27019, r27013, MPFR_RNDN);
        mpfr_sub(r27021, r27020, r27002, MPFR_RNDN);
        mpfr_add(r27022, r27017, r27021, MPFR_RNDN);
        return mpfr_get_d(r27022, MPFR_RNDN);
}

static mpfr_t r27023, r27024, r27025, r27026, r27027, r27028, r27029, r27030, r27031, r27032, r27033, r27034, r27035, r27036, r27037, r27038, r27039, r27040, r27041, r27042, r27043, r27044, r27045, r27046;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27023);
        mpfr_init(r27024);
        mpfr_init(r27025);
        mpfr_init(r27026);
        mpfr_init(r27027);
        mpfr_init_set_str(r27028, "1", 10, MPFR_RNDN);
        mpfr_init(r27029);
        mpfr_init(r27030);
        mpfr_init(r27031);
        mpfr_init(r27032);
        mpfr_init(r27033);
        mpfr_init(r27034);
        mpfr_init(r27035);
        mpfr_init(r27036);
        mpfr_init(r27037);
        mpfr_init(r27038);
        mpfr_init(r27039);
        mpfr_init(r27040);
        mpfr_init(r27041);
        mpfr_init(r27042);
        mpfr_init(r27043);
        mpfr_init(r27044);
        mpfr_init(r27045);
        mpfr_init(r27046);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r27023, x, MPFR_RNDN);
        mpfr_sin(r27024, r27023, MPFR_RNDN);
        mpfr_cos(r27025, r27023, MPFR_RNDN);
        mpfr_div(r27026, r27024, r27025, MPFR_RNDN);
        mpfr_mul(r27027, r27026, r27026, MPFR_RNDN);
        ;
        mpfr_set_d(r27029, eps, MPFR_RNDN);
        mpfr_sin(r27030, r27029, MPFR_RNDN);
        mpfr_div(r27031, r27030, r27025, MPFR_RNDN);
        mpfr_mul(r27032, r27031, r27031, MPFR_RNDN);
        mpfr_cos(r27033, r27029, MPFR_RNDN);
        mpfr_div(r27034, r27024, r27033, MPFR_RNDN);
        mpfr_mul(r27035, r27034, r27034, MPFR_RNDN);
        mpfr_mul(r27036, r27032, r27035, MPFR_RNDN);
        mpfr_sub(r27037, r27028, r27036, MPFR_RNDN);
        mpfr_div(r27038, r27027, r27037, MPFR_RNDN);
        mpfr_div(r27039, r27030, r27033, MPFR_RNDN);
        mpfr_div(r27040, r27039, r27037, MPFR_RNDN);
        mpfr_fma(r27041, r27038, r27039, r27040, MPFR_RNDN);
        mpfr_fma(r27042, r27039, r27039, r27028, MPFR_RNDN);
        mpfr_mul(r27043, r27042, r27026, MPFR_RNDN);
        mpfr_div(r27044, r27043, r27037, MPFR_RNDN);
        mpfr_sub(r27045, r27044, r27026, MPFR_RNDN);
        mpfr_add(r27046, r27041, r27045, MPFR_RNDN);
        return mpfr_get_d(r27046, MPFR_RNDN);
}