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

char *name = "Hyperbolic arc-cosine";

double f_if(float x) {
        float r9149 = x;
        float r9150 = r9149 * r9149;
        float r9151 = 1;
        float r9152 = r9150 - r9151;
        float r9153 = sqrt(r9152);
        float r9154 = r9149 + r9153;
        float r9155 = log(r9154);
        return r9155;
}

double f_id(double x) {
        double r9156 = x;
        double r9157 = r9156 * r9156;
        double r9158 = 1;
        double r9159 = r9157 - r9158;
        double r9160 = sqrt(r9159);
        double r9161 = r9156 + r9160;
        double r9162 = log(r9161);
        return r9162;
}


double f_of(float x) {
        float r9163 = x;
        float r9164 = 1;
        float r9165 = r9163 + r9164;
        float r9166 = sqrt(r9165);
        float r9167 = r9163 - r9164;
        float r9168 = sqrt(r9167);
        float r9169 = r9166 * r9168;
        float r9170 = r9163 + r9169;
        float r9171 = log(r9170);
        return r9171;
}

double f_od(double x) {
        double r9172 = x;
        double r9173 = 1;
        double r9174 = r9172 + r9173;
        double r9175 = sqrt(r9174);
        double r9176 = r9172 - r9173;
        double r9177 = sqrt(r9176);
        double r9178 = r9175 * r9177;
        double r9179 = r9172 + r9178;
        double r9180 = log(r9179);
        return r9180;
}

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 r9181, r9182, r9183, r9184, r9185, r9186, r9187;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init_set_str(r9183, "1", 10, MPFR_RNDN);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
}

double f_im(double x) {
        mpfr_set_d(r9181, x, MPFR_RNDN);
        mpfr_mul(r9182, r9181, r9181, MPFR_RNDN);
        ;
        mpfr_sub(r9184, r9182, r9183, MPFR_RNDN);
        mpfr_sqrt(r9185, r9184, MPFR_RNDN);
        mpfr_add(r9186, r9181, r9185, MPFR_RNDN);
        mpfr_log(r9187, r9186, MPFR_RNDN);
        return mpfr_get_d(r9187, MPFR_RNDN);
}

static mpfr_t r9188, r9189, r9190, r9191, r9192, r9193, r9194, r9195, r9196;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9188);
        mpfr_init_set_str(r9189, "1", 10, MPFR_RNDN);
        mpfr_init(r9190);
        mpfr_init(r9191);
        mpfr_init(r9192);
        mpfr_init(r9193);
        mpfr_init(r9194);
        mpfr_init(r9195);
        mpfr_init(r9196);
}

double f_fm(double x) {
        mpfr_set_d(r9188, x, MPFR_RNDN);
        ;
        mpfr_add(r9190, r9188, r9189, MPFR_RNDN);
        mpfr_sqrt(r9191, r9190, MPFR_RNDN);
        mpfr_sub(r9192, r9188, r9189, MPFR_RNDN);
        mpfr_sqrt(r9193, r9192, MPFR_RNDN);
        mpfr_mul(r9194, r9191, r9193, MPFR_RNDN);
        mpfr_add(r9195, r9188, r9194, MPFR_RNDN);
        mpfr_log(r9196, r9195, MPFR_RNDN);
        return mpfr_get_d(r9196, MPFR_RNDN);
}

static mpfr_t r9197, r9198, r9199, r9200, r9201, r9202, r9203, r9204, r9205;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9197);
        mpfr_init_set_str(r9198, "1", 10, MPFR_RNDN);
        mpfr_init(r9199);
        mpfr_init(r9200);
        mpfr_init(r9201);
        mpfr_init(r9202);
        mpfr_init(r9203);
        mpfr_init(r9204);
        mpfr_init(r9205);
}

double f_dm(double x) {
        mpfr_set_d(r9197, x, MPFR_RNDN);
        ;
        mpfr_add(r9199, r9197, r9198, MPFR_RNDN);
        mpfr_sqrt(r9200, r9199, MPFR_RNDN);
        mpfr_sub(r9201, r9197, r9198, MPFR_RNDN);
        mpfr_sqrt(r9202, r9201, MPFR_RNDN);
        mpfr_mul(r9203, r9200, r9202, MPFR_RNDN);
        mpfr_add(r9204, r9197, r9203, MPFR_RNDN);
        mpfr_log(r9205, r9204, MPFR_RNDN);
        return mpfr_get_d(r9205, MPFR_RNDN);
}