#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 r19994 = 1.0f;
        float r19995 = x;
        float r19996 = r19994 / r19995;
        float r19997 = tan(r19995);
        float r19998 = r19994 / r19997;
        float r19999 = r19996 - r19998;
        return r19999;
}

double f_id(double x) {
        double r20000 = 1.0;
        double r20001 = x;
        double r20002 = r20000 / r20001;
        double r20003 = tan(r20001);
        double r20004 = r20000 / r20003;
        double r20005 = r20002 - r20004;
        return r20005;
}


double f_of(float x) {
        float r20006 = 0.022222222222222223f;
        float r20007 = x;
        float r20008 = 3.0f;
        float r20009 = pow(r20007, r20008);
        float r20010 = r20006 * r20009;
        float r20011 = 0.0021164021164021165f;
        float r20012 = 5.0f;
        float r20013 = pow(r20007, r20012);
        float r20014 = r20011 * r20013;
        float r20015 = 0.3333333333333333f;
        float r20016 = r20015 * r20007;
        float r20017 = r20014 + r20016;
        float r20018 = r20010 + r20017;
        return r20018;
}

double f_od(double x) {
        double r20019 = 0.022222222222222223;
        double r20020 = x;
        double r20021 = 3.0;
        double r20022 = pow(r20020, r20021);
        double r20023 = r20019 * r20022;
        double r20024 = 0.0021164021164021165;
        double r20025 = 5.0;
        double r20026 = pow(r20020, r20025);
        double r20027 = r20024 * r20026;
        double r20028 = 0.3333333333333333;
        double r20029 = r20028 * r20020;
        double r20030 = r20027 + r20029;
        double r20031 = r20023 + r20030;
        return r20031;
}

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 r20032, r20033, r20034, r20035, r20036, r20037;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20032, "1", 10, MPFR_RNDN);
        mpfr_init(r20033);
        mpfr_init(r20034);
        mpfr_init(r20035);
        mpfr_init(r20036);
        mpfr_init(r20037);
}

double f_im(double x) {
        ;
        mpfr_set_d(r20033, x, MPFR_RNDN);
        mpfr_div(r20034, r20032, r20033, MPFR_RNDN);
        mpfr_tan(r20035, r20033, MPFR_RNDN);
        mpfr_div(r20036, r20032, r20035, MPFR_RNDN);
        mpfr_sub(r20037, r20034, r20036, MPFR_RNDN);
        return mpfr_get_d(r20037, MPFR_RNDN);
}

static mpfr_t r20038, r20039, r20040, r20041, r20042, r20043, r20044, r20045, r20046, r20047, r20048, r20049, r20050;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20038, "1/45", 10, MPFR_RNDN);
        mpfr_init(r20039);
        mpfr_init_set_str(r20040, "3", 10, MPFR_RNDN);
        mpfr_init(r20041);
        mpfr_init(r20042);
        mpfr_init_set_str(r20043, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r20044, "5", 10, MPFR_RNDN);
        mpfr_init(r20045);
        mpfr_init(r20046);
        mpfr_init_set_str(r20047, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20048);
        mpfr_init(r20049);
        mpfr_init(r20050);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r20039, x, MPFR_RNDN);
        ;
        mpfr_pow(r20041, r20039, r20040, MPFR_RNDN);
        mpfr_mul(r20042, r20038, r20041, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20045, r20039, r20044, MPFR_RNDN);
        mpfr_mul(r20046, r20043, r20045, MPFR_RNDN);
        ;
        mpfr_mul(r20048, r20047, r20039, MPFR_RNDN);
        mpfr_add(r20049, r20046, r20048, MPFR_RNDN);
        mpfr_add(r20050, r20042, r20049, MPFR_RNDN);
        return mpfr_get_d(r20050, MPFR_RNDN);
}

static mpfr_t r20051, r20052, r20053, r20054, r20055, r20056, r20057, r20058, r20059, r20060, r20061, r20062, r20063;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20051, "1/45", 10, MPFR_RNDN);
        mpfr_init(r20052);
        mpfr_init_set_str(r20053, "3", 10, MPFR_RNDN);
        mpfr_init(r20054);
        mpfr_init(r20055);
        mpfr_init_set_str(r20056, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r20057, "5", 10, MPFR_RNDN);
        mpfr_init(r20058);
        mpfr_init(r20059);
        mpfr_init_set_str(r20060, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20061);
        mpfr_init(r20062);
        mpfr_init(r20063);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r20052, x, MPFR_RNDN);
        ;
        mpfr_pow(r20054, r20052, r20053, MPFR_RNDN);
        mpfr_mul(r20055, r20051, r20054, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20058, r20052, r20057, MPFR_RNDN);
        mpfr_mul(r20059, r20056, r20058, MPFR_RNDN);
        ;
        mpfr_mul(r20061, r20060, r20052, MPFR_RNDN);
        mpfr_add(r20062, r20059, r20061, MPFR_RNDN);
        mpfr_add(r20063, r20055, r20062, MPFR_RNDN);
        return mpfr_get_d(r20063, MPFR_RNDN);
}