#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 r20801 = x;
        float r20802 = eps;
        float r20803 = r20801 + r20802;
        float r20804 = tan(r20803);
        float r20805 = tan(r20801);
        float r20806 = r20804 - r20805;
        return r20806;
}

double f_id(double x, double eps) {
        double r20807 = x;
        double r20808 = eps;
        double r20809 = r20807 + r20808;
        double r20810 = tan(r20809);
        double r20811 = tan(r20807);
        double r20812 = r20810 - r20811;
        return r20812;
}


double f_of(float x, float eps) {
        float r20813 = eps;
        float r20814 = -2.1454937753948827e-23f;
        bool r20815 = r20813 <= r20814;
        float r20816 = x;
        float r20817 = tan(r20816);
        float r20818 = tan(r20813);
        float r20819 = r20817 + r20818;
        float r20820 = 1.0f;
        float r20821 = 3.0f;
        float r20822 = pow(r20817, r20821);
        float r20823 = pow(r20818, r20821);
        float r20824 = r20822 * r20823;
        float r20825 = cbrt(r20824);
        float r20826 = r20820 - r20825;
        float r20827 = r20819 / r20826;
        float r20828 = r20827 - r20817;
        float r20829 = 6.683176783612011e-65f;
        bool r20830 = r20813 <= r20829;
        float r20831 = pow(r20813, r20821);
        float r20832 = 2.0f;
        float r20833 = pow(r20816, r20832);
        float r20834 = r20831 * r20833;
        float r20835 = pow(r20813, r20832);
        float r20836 = r20835 * r20816;
        float r20837 = r20834 + r20836;
        float r20838 = r20813 + r20837;
        float r20839 = r20830 ? r20838 : r20828;
        float r20840 = r20815 ? r20828 : r20839;
        return r20840;
}

double f_od(double x, double eps) {
        double r20841 = eps;
        double r20842 = -2.1454937753948827e-23;
        bool r20843 = r20841 <= r20842;
        double r20844 = x;
        double r20845 = tan(r20844);
        double r20846 = tan(r20841);
        double r20847 = r20845 + r20846;
        double r20848 = 1.0;
        double r20849 = 3.0;
        double r20850 = pow(r20845, r20849);
        double r20851 = pow(r20846, r20849);
        double r20852 = r20850 * r20851;
        double r20853 = cbrt(r20852);
        double r20854 = r20848 - r20853;
        double r20855 = r20847 / r20854;
        double r20856 = r20855 - r20845;
        double r20857 = 6.683176783612011e-65;
        bool r20858 = r20841 <= r20857;
        double r20859 = pow(r20841, r20849);
        double r20860 = 2.0;
        double r20861 = pow(r20844, r20860);
        double r20862 = r20859 * r20861;
        double r20863 = pow(r20841, r20860);
        double r20864 = r20863 * r20844;
        double r20865 = r20862 + r20864;
        double r20866 = r20841 + r20865;
        double r20867 = r20858 ? r20866 : r20856;
        double r20868 = r20843 ? r20856 : r20867;
        return r20868;
}

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 r20869, r20870, r20871, r20872, r20873, r20874;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20869);
        mpfr_init(r20870);
        mpfr_init(r20871);
        mpfr_init(r20872);
        mpfr_init(r20873);
        mpfr_init(r20874);
}

double f_im(double x, double eps) {
        mpfr_set_d(r20869, x, MPFR_RNDN);
        mpfr_set_d(r20870, eps, MPFR_RNDN);
        mpfr_add(r20871, r20869, r20870, MPFR_RNDN);
        mpfr_tan(r20872, r20871, MPFR_RNDN);
        mpfr_tan(r20873, r20869, MPFR_RNDN);
        mpfr_sub(r20874, r20872, r20873, MPFR_RNDN);
        return mpfr_get_d(r20874, MPFR_RNDN);
}

static mpfr_t r20875, r20876, r20877, r20878, r20879, r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20875);
        mpfr_init_set_str(r20876, "-2.1454937753948827e-23", 10, MPFR_RNDN);
        mpfr_init(r20877);
        mpfr_init(r20878);
        mpfr_init(r20879);
        mpfr_init(r20880);
        mpfr_init(r20881);
        mpfr_init_set_str(r20882, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20883, "3", 10, MPFR_RNDN);
        mpfr_init(r20884);
        mpfr_init(r20885);
        mpfr_init(r20886);
        mpfr_init(r20887);
        mpfr_init(r20888);
        mpfr_init(r20889);
        mpfr_init(r20890);
        mpfr_init_set_str(r20891, "6.683176783612011e-65", 10, MPFR_RNDN);
        mpfr_init(r20892);
        mpfr_init(r20893);
        mpfr_init_set_str(r20894, "2", 10, MPFR_RNDN);
        mpfr_init(r20895);
        mpfr_init(r20896);
        mpfr_init(r20897);
        mpfr_init(r20898);
        mpfr_init(r20899);
        mpfr_init(r20900);
        mpfr_init(r20901);
        mpfr_init(r20902);
}

double f_fm(double x, double eps) {
        mpfr_set_d(r20875, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20877, mpfr_cmp(r20875, r20876) <= 0, MPFR_RNDN);
        mpfr_set_d(r20878, x, MPFR_RNDN);
        mpfr_tan(r20879, r20878, MPFR_RNDN);
        mpfr_tan(r20880, r20875, MPFR_RNDN);
        mpfr_add(r20881, r20879, r20880, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20884, r20879, r20883, MPFR_RNDN);
        mpfr_pow(r20885, r20880, r20883, MPFR_RNDN);
        mpfr_mul(r20886, r20884, r20885, MPFR_RNDN);
        mpfr_cbrt(r20887, r20886, MPFR_RNDN);
        mpfr_sub(r20888, r20882, r20887, MPFR_RNDN);
        mpfr_div(r20889, r20881, r20888, MPFR_RNDN);
        mpfr_sub(r20890, r20889, r20879, MPFR_RNDN);
        ;
        mpfr_set_si(r20892, mpfr_cmp(r20875, r20891) <= 0, MPFR_RNDN);
        mpfr_pow(r20893, r20875, r20883, MPFR_RNDN);
        ;
        mpfr_pow(r20895, r20878, r20894, MPFR_RNDN);
        mpfr_mul(r20896, r20893, r20895, MPFR_RNDN);
        mpfr_pow(r20897, r20875, r20894, MPFR_RNDN);
        mpfr_mul(r20898, r20897, r20878, MPFR_RNDN);
        mpfr_add(r20899, r20896, r20898, MPFR_RNDN);
        mpfr_add(r20900, r20875, r20899, MPFR_RNDN);
        if (mpfr_get_si(r20892, MPFR_RNDN)) { mpfr_set(r20901, r20900, MPFR_RNDN); } else { mpfr_set(r20901, r20890, MPFR_RNDN); };
        if (mpfr_get_si(r20877, MPFR_RNDN)) { mpfr_set(r20902, r20890, MPFR_RNDN); } else { mpfr_set(r20902, r20901, MPFR_RNDN); };
        return mpfr_get_d(r20902, MPFR_RNDN);
}

static mpfr_t r20903, r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922, r20923, r20924, r20925, r20926, r20927, r20928, r20929, r20930;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2192);
        mpfr_init(r20903);
        mpfr_init_set_str(r20904, "-2.1454937753948827e-23", 10, MPFR_RNDN);
        mpfr_init(r20905);
        mpfr_init(r20906);
        mpfr_init(r20907);
        mpfr_init(r20908);
        mpfr_init(r20909);
        mpfr_init_set_str(r20910, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r20911, "3", 10, MPFR_RNDN);
        mpfr_init(r20912);
        mpfr_init(r20913);
        mpfr_init(r20914);
        mpfr_init(r20915);
        mpfr_init(r20916);
        mpfr_init(r20917);
        mpfr_init(r20918);
        mpfr_init_set_str(r20919, "6.683176783612011e-65", 10, MPFR_RNDN);
        mpfr_init(r20920);
        mpfr_init(r20921);
        mpfr_init_set_str(r20922, "2", 10, MPFR_RNDN);
        mpfr_init(r20923);
        mpfr_init(r20924);
        mpfr_init(r20925);
        mpfr_init(r20926);
        mpfr_init(r20927);
        mpfr_init(r20928);
        mpfr_init(r20929);
        mpfr_init(r20930);
}

double f_dm(double x, double eps) {
        mpfr_set_d(r20903, eps, MPFR_RNDN);
        ;
        mpfr_set_si(r20905, mpfr_cmp(r20903, r20904) <= 0, MPFR_RNDN);
        mpfr_set_d(r20906, x, MPFR_RNDN);
        mpfr_tan(r20907, r20906, MPFR_RNDN);
        mpfr_tan(r20908, r20903, MPFR_RNDN);
        mpfr_add(r20909, r20907, r20908, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20912, r20907, r20911, MPFR_RNDN);
        mpfr_pow(r20913, r20908, r20911, MPFR_RNDN);
        mpfr_mul(r20914, r20912, r20913, MPFR_RNDN);
        mpfr_cbrt(r20915, r20914, MPFR_RNDN);
        mpfr_sub(r20916, r20910, r20915, MPFR_RNDN);
        mpfr_div(r20917, r20909, r20916, MPFR_RNDN);
        mpfr_sub(r20918, r20917, r20907, MPFR_RNDN);
        ;
        mpfr_set_si(r20920, mpfr_cmp(r20903, r20919) <= 0, MPFR_RNDN);
        mpfr_pow(r20921, r20903, r20911, MPFR_RNDN);
        ;
        mpfr_pow(r20923, r20906, r20922, MPFR_RNDN);
        mpfr_mul(r20924, r20921, r20923, MPFR_RNDN);
        mpfr_pow(r20925, r20903, r20922, MPFR_RNDN);
        mpfr_mul(r20926, r20925, r20906, MPFR_RNDN);
        mpfr_add(r20927, r20924, r20926, MPFR_RNDN);
        mpfr_add(r20928, r20903, r20927, MPFR_RNDN);
        if (mpfr_get_si(r20920, MPFR_RNDN)) { mpfr_set(r20929, r20928, MPFR_RNDN); } else { mpfr_set(r20929, r20918, MPFR_RNDN); };
        if (mpfr_get_si(r20905, MPFR_RNDN)) { mpfr_set(r20930, r20918, MPFR_RNDN); } else { mpfr_set(r20930, r20929, MPFR_RNDN); };
        return mpfr_get_d(r20930, MPFR_RNDN);
}

