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

char *name = "NMSE problem 3.4.5";

double f_if(float x) {
        float r14863 = x;
        float r14864 = sin(r14863);
        float r14865 = r14863 - r14864;
        float r14866 = tan(r14863);
        float r14867 = r14863 - r14866;
        float r14868 = r14865 / r14867;
        return r14868;
}

double f_id(double x) {
        double r14869 = x;
        double r14870 = sin(r14869);
        double r14871 = r14869 - r14870;
        double r14872 = tan(r14869);
        double r14873 = r14869 - r14872;
        double r14874 = r14871 / r14873;
        return r14874;
}


double f_of(float x) {
        float r14875 = x;
        float r14876 = -3.2784791869173944e-06f;
        bool r14877 = r14875 <= r14876;
        float r14878 = sin(r14875);
        float r14879 = r14875 - r14878;
        float r14880 = tan(r14875);
        float r14881 = r14875 - r14880;
        float r14882 = r14879 / r14881;
        float r14883 = 501837.3637557024f;
        bool r14884 = r14875 <= r14883;
        float r14885 = 0.225f;
        float r14886 = r14875 * r14875;
        float r14887 = r14885 * r14886;
        float r14888 = 0.009642857142857142f;
        float r14889 = 4.0f;
        float r14890 = pow(r14875, r14889);
        float r14891 = r14888 * r14890;
        float r14892 = 0.5f;
        float r14893 = r14891 + r14892;
        float r14894 = r14887 - r14893;
        float r14895 = r14884 ? r14894 : r14882;
        float r14896 = r14877 ? r14882 : r14895;
        return r14896;
}

double f_od(double x) {
        double r14897 = x;
        double r14898 = -3.2784791869173944e-06;
        bool r14899 = r14897 <= r14898;
        double r14900 = sin(r14897);
        double r14901 = r14897 - r14900;
        double r14902 = tan(r14897);
        double r14903 = r14897 - r14902;
        double r14904 = r14901 / r14903;
        double r14905 = 501837.3637557024;
        bool r14906 = r14897 <= r14905;
        double r14907 = 0.225;
        double r14908 = r14897 * r14897;
        double r14909 = r14907 * r14908;
        double r14910 = 0.009642857142857142;
        double r14911 = 4.0;
        double r14912 = pow(r14897, r14911);
        double r14913 = r14910 * r14912;
        double r14914 = 0.5;
        double r14915 = r14913 + r14914;
        double r14916 = r14909 - r14915;
        double r14917 = r14906 ? r14916 : r14904;
        double r14918 = r14899 ? r14904 : r14917;
        return r14918;
}

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 r14919, r14920, r14921, r14922, r14923, r14924;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r14919);
        mpfr_init(r14920);
        mpfr_init(r14921);
        mpfr_init(r14922);
        mpfr_init(r14923);
        mpfr_init(r14924);
}

double f_im(double x) {
        mpfr_set_d(r14919, x, MPFR_RNDN);
        mpfr_sin(r14920, r14919, MPFR_RNDN);
        mpfr_sub(r14921, r14919, r14920, MPFR_RNDN);
        mpfr_tan(r14922, r14919, MPFR_RNDN);
        mpfr_sub(r14923, r14919, r14922, MPFR_RNDN);
        mpfr_div(r14924, r14921, r14923, MPFR_RNDN);
        return mpfr_get_d(r14924, MPFR_RNDN);
}

static mpfr_t r14925, r14926, r14927, r14928, r14929, r14930, r14931, r14932, r14933, r14934, r14935, r14936, r14937, r14938, r14939, r14940, r14941, r14942, r14943, r14944, r14945, r14946;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r14925);
        mpfr_init_set_str(r14926, "-3.2784791869173944e-06", 10, MPFR_RNDN);
        mpfr_init(r14927);
        mpfr_init(r14928);
        mpfr_init(r14929);
        mpfr_init(r14930);
        mpfr_init(r14931);
        mpfr_init(r14932);
        mpfr_init_set_str(r14933, "501837.3637557024", 10, MPFR_RNDN);
        mpfr_init(r14934);
        mpfr_init_set_str(r14935, "9/40", 10, MPFR_RNDN);
        mpfr_init(r14936);
        mpfr_init(r14937);
        mpfr_init_set_str(r14938, "27/2800", 10, MPFR_RNDN);
        mpfr_init_set_str(r14939, "4", 10, MPFR_RNDN);
        mpfr_init(r14940);
        mpfr_init(r14941);
        mpfr_init_set_str(r14942, "1/2", 10, MPFR_RNDN);
        mpfr_init(r14943);
        mpfr_init(r14944);
        mpfr_init(r14945);
        mpfr_init(r14946);
}

double f_fm(double x) {
        mpfr_set_d(r14925, x, MPFR_RNDN);
        ;
        mpfr_set_si(r14927, mpfr_cmp(r14925, r14926) <= 0, MPFR_RNDN);
        mpfr_sin(r14928, r14925, MPFR_RNDN);
        mpfr_sub(r14929, r14925, r14928, MPFR_RNDN);
        mpfr_tan(r14930, r14925, MPFR_RNDN);
        mpfr_sub(r14931, r14925, r14930, MPFR_RNDN);
        mpfr_div(r14932, r14929, r14931, MPFR_RNDN);
        ;
        mpfr_set_si(r14934, mpfr_cmp(r14925, r14933) <= 0, MPFR_RNDN);
        ;
        mpfr_sqr(r14936, r14925, MPFR_RNDN);
        mpfr_mul(r14937, r14935, r14936, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r14940, r14925, r14939, MPFR_RNDN);
        mpfr_mul(r14941, r14938, r14940, MPFR_RNDN);
        ;
        mpfr_add(r14943, r14941, r14942, MPFR_RNDN);
        mpfr_sub(r14944, r14937, r14943, MPFR_RNDN);
        if (mpfr_get_si(r14934, MPFR_RNDN)) { mpfr_set(r14945, r14944, MPFR_RNDN); } else { mpfr_set(r14945, r14932, MPFR_RNDN); };
        if (mpfr_get_si(r14927, MPFR_RNDN)) { mpfr_set(r14946, r14932, MPFR_RNDN); } else { mpfr_set(r14946, r14945, MPFR_RNDN); };
        return mpfr_get_d(r14946, MPFR_RNDN);
}

static mpfr_t r14947, r14948, r14949, r14950, r14951, r14952, r14953, r14954, r14955, r14956, r14957, r14958, r14959, r14960, r14961, r14962, r14963, r14964, r14965, r14966, r14967, r14968;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r14947);
        mpfr_init_set_str(r14948, "-3.2784791869173944e-06", 10, MPFR_RNDN);
        mpfr_init(r14949);
        mpfr_init(r14950);
        mpfr_init(r14951);
        mpfr_init(r14952);
        mpfr_init(r14953);
        mpfr_init(r14954);
        mpfr_init_set_str(r14955, "501837.3637557024", 10, MPFR_RNDN);
        mpfr_init(r14956);
        mpfr_init_set_str(r14957, "9/40", 10, MPFR_RNDN);
        mpfr_init(r14958);
        mpfr_init(r14959);
        mpfr_init_set_str(r14960, "27/2800", 10, MPFR_RNDN);
        mpfr_init_set_str(r14961, "4", 10, MPFR_RNDN);
        mpfr_init(r14962);
        mpfr_init(r14963);
        mpfr_init_set_str(r14964, "1/2", 10, MPFR_RNDN);
        mpfr_init(r14965);
        mpfr_init(r14966);
        mpfr_init(r14967);
        mpfr_init(r14968);
}

double f_dm(double x) {
        mpfr_set_d(r14947, x, MPFR_RNDN);
        ;
        mpfr_set_si(r14949, mpfr_cmp(r14947, r14948) <= 0, MPFR_RNDN);
        mpfr_sin(r14950, r14947, MPFR_RNDN);
        mpfr_sub(r14951, r14947, r14950, MPFR_RNDN);
        mpfr_tan(r14952, r14947, MPFR_RNDN);
        mpfr_sub(r14953, r14947, r14952, MPFR_RNDN);
        mpfr_div(r14954, r14951, r14953, MPFR_RNDN);
        ;
        mpfr_set_si(r14956, mpfr_cmp(r14947, r14955) <= 0, MPFR_RNDN);
        ;
        mpfr_sqr(r14958, r14947, MPFR_RNDN);
        mpfr_mul(r14959, r14957, r14958, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r14962, r14947, r14961, MPFR_RNDN);
        mpfr_mul(r14963, r14960, r14962, MPFR_RNDN);
        ;
        mpfr_add(r14965, r14963, r14964, MPFR_RNDN);
        mpfr_sub(r14966, r14959, r14965, MPFR_RNDN);
        if (mpfr_get_si(r14956, MPFR_RNDN)) { mpfr_set(r14967, r14966, MPFR_RNDN); } else { mpfr_set(r14967, r14954, MPFR_RNDN); };
        if (mpfr_get_si(r14949, MPFR_RNDN)) { mpfr_set(r14968, r14954, MPFR_RNDN); } else { mpfr_set(r14968, r14967, MPFR_RNDN); };
        return mpfr_get_d(r14968, MPFR_RNDN);
}