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

char *name = "NMSE problem 3.4.1";

double f_if(float x) {
        float r15143 = 1.0f;
        float r15144 = x;
        float r15145 = cos(r15144);
        float r15146 = r15143 - r15145;
        float r15147 = r15144 * r15144;
        float r15148 = r15146 / r15147;
        return r15148;
}

double f_id(double x) {
        double r15149 = 1.0;
        double r15150 = x;
        double r15151 = cos(r15150);
        double r15152 = r15149 - r15151;
        double r15153 = r15150 * r15150;
        double r15154 = r15152 / r15153;
        return r15154;
}


double f_of(float x) {
        float r15155 = x;
        float r15156 = -0.3497152626514435f;
        bool r15157 = r15155 <= r15156;
        float r15158 = 1.0f;
        float r15159 = cos(r15155);
        float r15160 = r15158 - r15159;
        float r15161 = sqrt(r15160);
        float r15162 = r15161 / r15155;
        float r15163 = r15162 * r15162;
        float r15164 = 0.44031599164009094f;
        bool r15165 = r15155 <= r15164;
        float r15166 = 0.0013888889225199819f;
        float r15167 = 4.0f;
        float r15168 = pow(r15155, r15167);
        float r15169 = r15166 * r15168;
        float r15170 = 0.5f;
        float r15171 = r15169 + r15170;
        float r15172 = 0.0416666679084301f;
        float r15173 = r15155 * r15155;
        float r15174 = r15172 * r15173;
        float r15175 = r15171 - r15174;
        float r15176 = r15165 ? r15175 : r15163;
        float r15177 = r15157 ? r15163 : r15176;
        return r15177;
}

double f_od(double x) {
        double r15178 = x;
        double r15179 = -0.3497152626514435;
        bool r15180 = r15178 <= r15179;
        double r15181 = 1.0;
        double r15182 = cos(r15178);
        double r15183 = r15181 - r15182;
        double r15184 = sqrt(r15183);
        double r15185 = r15184 / r15178;
        double r15186 = r15185 * r15185;
        double r15187 = 0.44031599164009094;
        bool r15188 = r15178 <= r15187;
        double r15189 = 0.0013888889225199819;
        double r15190 = 4.0;
        double r15191 = pow(r15178, r15190);
        double r15192 = r15189 * r15191;
        double r15193 = 0.5;
        double r15194 = r15192 + r15193;
        double r15195 = 0.0416666679084301;
        double r15196 = r15178 * r15178;
        double r15197 = r15195 * r15196;
        double r15198 = r15194 - r15197;
        double r15199 = r15188 ? r15198 : r15186;
        double r15200 = r15180 ? r15186 : r15199;
        return r15200;
}

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 r15201, r15202, r15203, r15204, r15205, r15206;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r15201, "1", 10, MPFR_RNDN);
        mpfr_init(r15202);
        mpfr_init(r15203);
        mpfr_init(r15204);
        mpfr_init(r15205);
        mpfr_init(r15206);
}

double f_im(double x) {
        ;
        mpfr_set_d(r15202, x, MPFR_RNDN);
        mpfr_cos(r15203, r15202, MPFR_RNDN);
        mpfr_sub(r15204, r15201, r15203, MPFR_RNDN);
        mpfr_sqr(r15205, r15202, MPFR_RNDN);
        mpfr_div(r15206, r15204, r15205, MPFR_RNDN);
        return mpfr_get_d(r15206, MPFR_RNDN);
}

static mpfr_t r15207, r15208, r15209, r15210, r15211, r15212, r15213, r15214, r15215, r15216, r15217, r15218, r15219, r15220, r15221, r15222, r15223, r15224, r15225, r15226, r15227, r15228, r15229;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r15207);
        mpfr_init_set_str(r15208, "-0.34971526f0", 10, MPFR_RNDN);
        mpfr_init(r15209);
        mpfr_init_set_str(r15210, "1", 10, MPFR_RNDN);
        mpfr_init(r15211);
        mpfr_init(r15212);
        mpfr_init(r15213);
        mpfr_init(r15214);
        mpfr_init(r15215);
        mpfr_init_set_str(r15216, "0.440316f0", 10, MPFR_RNDN);
        mpfr_init(r15217);
        mpfr_init_set_str(r15218, "1/720", 10, MPFR_RNDN);
        mpfr_init_set_str(r15219, "4", 10, MPFR_RNDN);
        mpfr_init(r15220);
        mpfr_init(r15221);
        mpfr_init_set_str(r15222, "1/2", 10, MPFR_RNDN);
        mpfr_init(r15223);
        mpfr_init_set_str(r15224, "1/24", 10, MPFR_RNDN);
        mpfr_init(r15225);
        mpfr_init(r15226);
        mpfr_init(r15227);
        mpfr_init(r15228);
        mpfr_init(r15229);
}

double f_fm(double x) {
        mpfr_set_d(r15207, x, MPFR_RNDN);
        ;
        mpfr_set_si(r15209, mpfr_cmp(r15207, r15208) <= 0, MPFR_RNDN);
        ;
        mpfr_cos(r15211, r15207, MPFR_RNDN);
        mpfr_sub(r15212, r15210, r15211, MPFR_RNDN);
        mpfr_sqrt(r15213, r15212, MPFR_RNDN);
        mpfr_div(r15214, r15213, r15207, MPFR_RNDN);
        mpfr_sqr(r15215, r15214, MPFR_RNDN);
        ;
        mpfr_set_si(r15217, mpfr_cmp(r15207, r15216) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r15220, r15207, r15219, MPFR_RNDN);
        mpfr_mul(r15221, r15218, r15220, MPFR_RNDN);
        ;
        mpfr_add(r15223, r15221, r15222, MPFR_RNDN);
        ;
        mpfr_sqr(r15225, r15207, MPFR_RNDN);
        mpfr_mul(r15226, r15224, r15225, MPFR_RNDN);
        mpfr_sub(r15227, r15223, r15226, MPFR_RNDN);
        if (mpfr_get_si(r15217, MPFR_RNDN)) { mpfr_set(r15228, r15227, MPFR_RNDN); } else { mpfr_set(r15228, r15215, MPFR_RNDN); };
        if (mpfr_get_si(r15209, MPFR_RNDN)) { mpfr_set(r15229, r15215, MPFR_RNDN); } else { mpfr_set(r15229, r15228, MPFR_RNDN); };
        return mpfr_get_d(r15229, MPFR_RNDN);
}

static mpfr_t r15230, r15231, r15232, r15233, r15234, r15235, r15236, r15237, r15238, r15239, r15240, r15241, r15242, r15243, r15244, r15245, r15246, r15247, r15248, r15249, r15250, r15251, r15252;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r15230);
        mpfr_init_set_str(r15231, "-0.34971526f0", 10, MPFR_RNDN);
        mpfr_init(r15232);
        mpfr_init_set_str(r15233, "1", 10, MPFR_RNDN);
        mpfr_init(r15234);
        mpfr_init(r15235);
        mpfr_init(r15236);
        mpfr_init(r15237);
        mpfr_init(r15238);
        mpfr_init_set_str(r15239, "0.440316f0", 10, MPFR_RNDN);
        mpfr_init(r15240);
        mpfr_init_set_str(r15241, "1/720", 10, MPFR_RNDN);
        mpfr_init_set_str(r15242, "4", 10, MPFR_RNDN);
        mpfr_init(r15243);
        mpfr_init(r15244);
        mpfr_init_set_str(r15245, "1/2", 10, MPFR_RNDN);
        mpfr_init(r15246);
        mpfr_init_set_str(r15247, "1/24", 10, MPFR_RNDN);
        mpfr_init(r15248);
        mpfr_init(r15249);
        mpfr_init(r15250);
        mpfr_init(r15251);
        mpfr_init(r15252);
}

double f_dm(double x) {
        mpfr_set_d(r15230, x, MPFR_RNDN);
        ;
        mpfr_set_si(r15232, mpfr_cmp(r15230, r15231) <= 0, MPFR_RNDN);
        ;
        mpfr_cos(r15234, r15230, MPFR_RNDN);
        mpfr_sub(r15235, r15233, r15234, MPFR_RNDN);
        mpfr_sqrt(r15236, r15235, MPFR_RNDN);
        mpfr_div(r15237, r15236, r15230, MPFR_RNDN);
        mpfr_sqr(r15238, r15237, MPFR_RNDN);
        ;
        mpfr_set_si(r15240, mpfr_cmp(r15230, r15239) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r15243, r15230, r15242, MPFR_RNDN);
        mpfr_mul(r15244, r15241, r15243, MPFR_RNDN);
        ;
        mpfr_add(r15246, r15244, r15245, MPFR_RNDN);
        ;
        mpfr_sqr(r15248, r15230, MPFR_RNDN);
        mpfr_mul(r15249, r15247, r15248, MPFR_RNDN);
        mpfr_sub(r15250, r15246, r15249, MPFR_RNDN);
        if (mpfr_get_si(r15240, MPFR_RNDN)) { mpfr_set(r15251, r15250, MPFR_RNDN); } else { mpfr_set(r15251, r15238, MPFR_RNDN); };
        if (mpfr_get_si(r15232, MPFR_RNDN)) { mpfr_set(r15252, r15238, MPFR_RNDN); } else { mpfr_set(r15252, r15251, MPFR_RNDN); };
        return mpfr_get_d(r15252, MPFR_RNDN);
}