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

char *name = "invcot (example 3.9)";

double f_if(float x) {
        float r10056 = 1.0f;
        float r10057 = x;
        float r10058 = r10056 / r10057;
        float r10059 = tan(r10057);
        float r10060 = r10056 / r10059;
        float r10061 = r10058 - r10060;
        return r10061;
}

double f_id(double x) {
        double r10062 = 1.0;
        double r10063 = x;
        double r10064 = r10062 / r10063;
        double r10065 = tan(r10063);
        double r10066 = r10062 / r10065;
        double r10067 = r10064 - r10066;
        return r10067;
}


double f_of(float x) {
        float r10068 = 0.022222222222222223f;
        float r10069 = x;
        float r10070 = r10069 * r10069;
        float r10071 = 0.3333333333333333f;
        float r10072 = fma(r10068, r10070, r10071);
        float r10073 = 0.0021164021164021165f;
        float r10074 = 5.0f;
        float r10075 = pow(r10069, r10074);
        float r10076 = r10073 * r10075;
        float r10077 = fma(r10072, r10069, r10076);
        return r10077;
}

double f_od(double x) {
        double r10078 = 0.022222222222222223;
        double r10079 = x;
        double r10080 = r10079 * r10079;
        double r10081 = 0.3333333333333333;
        double r10082 = fma(r10078, r10080, r10081);
        double r10083 = 0.0021164021164021165;
        double r10084 = 5.0;
        double r10085 = pow(r10079, r10084);
        double r10086 = r10083 * r10085;
        double r10087 = fma(r10082, r10079, r10086);
        return r10087;
}

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 r10088, r10089, r10090, r10091, r10092, r10093;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10088, "1", 10, MPFR_RNDN);
        mpfr_init(r10089);
        mpfr_init(r10090);
        mpfr_init(r10091);
        mpfr_init(r10092);
        mpfr_init(r10093);
}

double f_im(double x) {
        ;
        mpfr_set_d(r10089, x, MPFR_RNDN);
        mpfr_div(r10090, r10088, r10089, MPFR_RNDN);
        mpfr_tan(r10091, r10089, MPFR_RNDN);
        mpfr_div(r10092, r10088, r10091, MPFR_RNDN);
        mpfr_sub(r10093, r10090, r10092, MPFR_RNDN);
        return mpfr_get_d(r10093, MPFR_RNDN);
}

static mpfr_t r10094, r10095, r10096, r10097, r10098, r10099, r10100, r10101, r10102, r10103;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10094, "1/45", 10, MPFR_RNDN);
        mpfr_init(r10095);
        mpfr_init(r10096);
        mpfr_init_set_str(r10097, "1/3", 10, MPFR_RNDN);
        mpfr_init(r10098);
        mpfr_init_set_str(r10099, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r10100, "5", 10, MPFR_RNDN);
        mpfr_init(r10101);
        mpfr_init(r10102);
        mpfr_init(r10103);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10095, x, MPFR_RNDN);
        mpfr_mul(r10096, r10095, r10095, MPFR_RNDN);
        ;
        mpfr_fma(r10098, r10094, r10096, r10097, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10101, r10095, r10100, MPFR_RNDN);
        mpfr_mul(r10102, r10099, r10101, MPFR_RNDN);
        mpfr_fma(r10103, r10098, r10095, r10102, MPFR_RNDN);
        return mpfr_get_d(r10103, MPFR_RNDN);
}

static mpfr_t r10104, r10105, r10106, r10107, r10108, r10109, r10110, r10111, r10112, r10113;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10104, "1/45", 10, MPFR_RNDN);
        mpfr_init(r10105);
        mpfr_init(r10106);
        mpfr_init_set_str(r10107, "1/3", 10, MPFR_RNDN);
        mpfr_init(r10108);
        mpfr_init_set_str(r10109, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r10110, "5", 10, MPFR_RNDN);
        mpfr_init(r10111);
        mpfr_init(r10112);
        mpfr_init(r10113);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10105, x, MPFR_RNDN);
        mpfr_mul(r10106, r10105, r10105, MPFR_RNDN);
        ;
        mpfr_fma(r10108, r10104, r10106, r10107, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10111, r10105, r10110, MPFR_RNDN);
        mpfr_mul(r10112, r10109, r10111, MPFR_RNDN);
        mpfr_fma(r10113, r10108, r10105, r10112, MPFR_RNDN);
        return mpfr_get_d(r10113, MPFR_RNDN);
}