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

char *name = "cos2 (problem 3.4.1)";

double f_if(float x) {
        float r18894 = 1.0f;
        float r18895 = x;
        float r18896 = cos(r18895);
        float r18897 = r18894 - r18896;
        float r18898 = r18895 * r18895;
        float r18899 = r18897 / r18898;
        return r18899;
}

double f_id(double x) {
        double r18900 = 1.0;
        double r18901 = x;
        double r18902 = cos(r18901);
        double r18903 = r18900 - r18902;
        double r18904 = r18901 * r18901;
        double r18905 = r18903 / r18904;
        return r18905;
}


double f_of(float x) {
        float r18906 = x;
        float r18907 = sin(r18906);
        float r18908 = r18907 / r18906;
        float r18909 = 1.0f;
        float r18910 = cos(r18906);
        float r18911 = r18909 + r18910;
        float r18912 = r18907 / r18911;
        float r18913 = r18912 / r18906;
        float r18914 = r18908 * r18913;
        return r18914;
}

double f_od(double x) {
        double r18915 = x;
        double r18916 = sin(r18915);
        double r18917 = r18916 / r18915;
        double r18918 = 1.0;
        double r18919 = cos(r18915);
        double r18920 = r18918 + r18919;
        double r18921 = r18916 / r18920;
        double r18922 = r18921 / r18915;
        double r18923 = r18917 * r18922;
        return r18923;
}

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 r18924, r18925, r18926, r18927, r18928, r18929;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r18924, "1", 10, MPFR_RNDN);
        mpfr_init(r18925);
        mpfr_init(r18926);
        mpfr_init(r18927);
        mpfr_init(r18928);
        mpfr_init(r18929);
}

double f_im(double x) {
        ;
        mpfr_set_d(r18925, x, MPFR_RNDN);
        mpfr_cos(r18926, r18925, MPFR_RNDN);
        mpfr_sub(r18927, r18924, r18926, MPFR_RNDN);
        mpfr_mul(r18928, r18925, r18925, MPFR_RNDN);
        mpfr_div(r18929, r18927, r18928, MPFR_RNDN);
        return mpfr_get_d(r18929, MPFR_RNDN);
}

static mpfr_t r18930, r18931, r18932, r18933, r18934, r18935, r18936, r18937, r18938;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r18930);
        mpfr_init(r18931);
        mpfr_init(r18932);
        mpfr_init_set_str(r18933, "1", 10, MPFR_RNDN);
        mpfr_init(r18934);
        mpfr_init(r18935);
        mpfr_init(r18936);
        mpfr_init(r18937);
        mpfr_init(r18938);
}

double f_fm(double x) {
        mpfr_set_d(r18930, x, MPFR_RNDN);
        mpfr_sin(r18931, r18930, MPFR_RNDN);
        mpfr_div(r18932, r18931, r18930, MPFR_RNDN);
        ;
        mpfr_cos(r18934, r18930, MPFR_RNDN);
        mpfr_add(r18935, r18933, r18934, MPFR_RNDN);
        mpfr_div(r18936, r18931, r18935, MPFR_RNDN);
        mpfr_div(r18937, r18936, r18930, MPFR_RNDN);
        mpfr_mul(r18938, r18932, r18937, MPFR_RNDN);
        return mpfr_get_d(r18938, MPFR_RNDN);
}

static mpfr_t r18939, r18940, r18941, r18942, r18943, r18944, r18945, r18946, r18947;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init(r18939);
        mpfr_init(r18940);
        mpfr_init(r18941);
        mpfr_init_set_str(r18942, "1", 10, MPFR_RNDN);
        mpfr_init(r18943);
        mpfr_init(r18944);
        mpfr_init(r18945);
        mpfr_init(r18946);
        mpfr_init(r18947);
}

double f_dm(double x) {
        mpfr_set_d(r18939, x, MPFR_RNDN);
        mpfr_sin(r18940, r18939, MPFR_RNDN);
        mpfr_div(r18941, r18940, r18939, MPFR_RNDN);
        ;
        mpfr_cos(r18943, r18939, MPFR_RNDN);
        mpfr_add(r18944, r18942, r18943, MPFR_RNDN);
        mpfr_div(r18945, r18940, r18944, MPFR_RNDN);
        mpfr_div(r18946, r18945, r18939, MPFR_RNDN);
        mpfr_mul(r18947, r18941, r18946, MPFR_RNDN);
        return mpfr_get_d(r18947, MPFR_RNDN);
}