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

char *name = "arccos";

double f_if(float x) {
        float r7786 = 2;
        float r7787 = 1;
        float r7788 = x;
        float r7789 = r7787 - r7788;
        float r7790 = r7787 + r7788;
        float r7791 = r7789 / r7790;
        float r7792 = sqrt(r7791);
        float r7793 = atan(r7792);
        float r7794 = r7786 * r7793;
        return r7794;
}

double f_id(double x) {
        double r7795 = 2;
        double r7796 = 1;
        double r7797 = x;
        double r7798 = r7796 - r7797;
        double r7799 = r7796 + r7797;
        double r7800 = r7798 / r7799;
        double r7801 = sqrt(r7800);
        double r7802 = atan(r7801);
        double r7803 = r7795 * r7802;
        return r7803;
}


double f_of(float x) {
        float r7804 = 2;
        float r7805 = 1;
        float r7806 = x;
        float r7807 = r7805 + r7806;
        float r7808 = r7805 / r7807;
        float r7809 = r7806 / r7807;
        float r7810 = r7808 - r7809;
        float r7811 = sqrt(r7810);
        float r7812 = atan(r7811);
        float r7813 = r7804 * r7812;
        return r7813;
}

double f_od(double x) {
        double r7814 = 2;
        double r7815 = 1;
        double r7816 = x;
        double r7817 = r7815 + r7816;
        double r7818 = r7815 / r7817;
        double r7819 = r7816 / r7817;
        double r7820 = r7818 - r7819;
        double r7821 = sqrt(r7820);
        double r7822 = atan(r7821);
        double r7823 = r7814 * r7822;
        return r7823;
}

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 r7824, r7825, r7826, r7827, r7828, r7829, r7830, r7831, r7832;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7824, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r7825, "1", 10, MPFR_RNDN);
        mpfr_init(r7826);
        mpfr_init(r7827);
        mpfr_init(r7828);
        mpfr_init(r7829);
        mpfr_init(r7830);
        mpfr_init(r7831);
        mpfr_init(r7832);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r7826, x, MPFR_RNDN);
        mpfr_sub(r7827, r7825, r7826, MPFR_RNDN);
        mpfr_add(r7828, r7825, r7826, MPFR_RNDN);
        mpfr_div(r7829, r7827, r7828, MPFR_RNDN);
        mpfr_sqrt(r7830, r7829, MPFR_RNDN);
        mpfr_atan(r7831, r7830, MPFR_RNDN);
        mpfr_mul(r7832, r7824, r7831, MPFR_RNDN);
        return mpfr_get_d(r7832, MPFR_RNDN);
}

static mpfr_t r7833, r7834, r7835, r7836, r7837, r7838, r7839, r7840, r7841, r7842;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7833, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r7834, "1", 10, MPFR_RNDN);
        mpfr_init(r7835);
        mpfr_init(r7836);
        mpfr_init(r7837);
        mpfr_init(r7838);
        mpfr_init(r7839);
        mpfr_init(r7840);
        mpfr_init(r7841);
        mpfr_init(r7842);
}

double f_fm(double x) {
        ;
        ;
        mpfr_set_d(r7835, x, MPFR_RNDN);
        mpfr_add(r7836, r7834, r7835, MPFR_RNDN);
        mpfr_div(r7837, r7834, r7836, MPFR_RNDN);
        mpfr_div(r7838, r7835, r7836, MPFR_RNDN);
        mpfr_sub(r7839, r7837, r7838, MPFR_RNDN);
        mpfr_sqrt(r7840, r7839, MPFR_RNDN);
        mpfr_atan(r7841, r7840, MPFR_RNDN);
        mpfr_mul(r7842, r7833, r7841, MPFR_RNDN);
        return mpfr_get_d(r7842, MPFR_RNDN);
}

static mpfr_t r7843, r7844, r7845, r7846, r7847, r7848, r7849, r7850, r7851, r7852;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r7843, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r7844, "1", 10, MPFR_RNDN);
        mpfr_init(r7845);
        mpfr_init(r7846);
        mpfr_init(r7847);
        mpfr_init(r7848);
        mpfr_init(r7849);
        mpfr_init(r7850);
        mpfr_init(r7851);
        mpfr_init(r7852);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r7845, x, MPFR_RNDN);
        mpfr_add(r7846, r7844, r7845, MPFR_RNDN);
        mpfr_div(r7847, r7844, r7846, MPFR_RNDN);
        mpfr_div(r7848, r7845, r7846, MPFR_RNDN);
        mpfr_sub(r7849, r7847, r7848, MPFR_RNDN);
        mpfr_sqrt(r7850, r7849, MPFR_RNDN);
        mpfr_atan(r7851, r7850, MPFR_RNDN);
        mpfr_mul(r7852, r7843, r7851, MPFR_RNDN);
        return mpfr_get_d(r7852, MPFR_RNDN);
}