#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 r9129 = x;
        float r9130 = r9129 * r9129;
        float r9131 = 1;
        float r9132 = r9130 - r9131;
        float r9133 = sqrt(r9132);
        float r9134 = r9129 + r9133;
        float r9135 = log(r9134);
        return r9135;
}

double f_id(double x) {
        double r9136 = x;
        double r9137 = r9136 * r9136;
        double r9138 = 1;
        double r9139 = r9137 - r9138;
        double r9140 = sqrt(r9139);
        double r9141 = r9136 + r9140;
        double r9142 = log(r9141);
        return r9142;
}


double f_of(float x) {
        float r9143 = x;
        float r9144 = 1;
        float r9145 = r9143 + r9144;
        float r9146 = sqrt(r9145);
        float r9147 = r9143 - r9144;
        float r9148 = sqrt(r9147);
        float r9149 = r9146 * r9148;
        float r9150 = r9143 + r9149;
        float r9151 = log(r9150);
        return r9151;
}

double f_od(double x) {
        double r9152 = x;
        double r9153 = 1;
        double r9154 = r9152 + r9153;
        double r9155 = sqrt(r9154);
        double r9156 = r9152 - r9153;
        double r9157 = sqrt(r9156);
        double r9158 = r9155 * r9157;
        double r9159 = r9152 + r9158;
        double r9160 = log(r9159);
        return r9160;
}

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 r9161, r9162, r9163, r9164, r9165, r9166, r9167;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9161);
        mpfr_init(r9162);
        mpfr_init_set_str(r9163, "1", 10, MPFR_RNDN);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init(r9167);
}

double f_im(double x) {
        mpfr_set_d(r9161, x, MPFR_RNDN);
        mpfr_mul(r9162, r9161, r9161, MPFR_RNDN);
        ;
        mpfr_sub(r9164, r9162, r9163, MPFR_RNDN);
        mpfr_sqrt(r9165, r9164, MPFR_RNDN);
        mpfr_add(r9166, r9161, r9165, MPFR_RNDN);
        mpfr_log(r9167, r9166, MPFR_RNDN);
        return mpfr_get_d(r9167, MPFR_RNDN);
}

static mpfr_t r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9168);
        mpfr_init_set_str(r9169, "1", 10, MPFR_RNDN);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
}

double f_fm(double x) {
        mpfr_set_d(r9168, x, MPFR_RNDN);
        ;
        mpfr_add(r9170, r9168, r9169, MPFR_RNDN);
        mpfr_sqrt(r9171, r9170, MPFR_RNDN);
        mpfr_sub(r9172, r9168, r9169, MPFR_RNDN);
        mpfr_sqrt(r9173, r9172, MPFR_RNDN);
        mpfr_mul(r9174, r9171, r9173, MPFR_RNDN);
        mpfr_add(r9175, r9168, r9174, MPFR_RNDN);
        mpfr_log(r9176, r9175, MPFR_RNDN);
        return mpfr_get_d(r9176, MPFR_RNDN);
}

static mpfr_t r9177, r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r9177);
        mpfr_init_set_str(r9178, "1", 10, MPFR_RNDN);
        mpfr_init(r9179);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init(r9183);
        mpfr_init(r9184);
        mpfr_init(r9185);
}

double f_dm(double x) {
        mpfr_set_d(r9177, x, MPFR_RNDN);
        ;
        mpfr_add(r9179, r9177, r9178, MPFR_RNDN);
        mpfr_sqrt(r9180, r9179, MPFR_RNDN);
        mpfr_sub(r9181, r9177, r9178, MPFR_RNDN);
        mpfr_sqrt(r9182, r9181, MPFR_RNDN);
        mpfr_mul(r9183, r9180, r9182, MPFR_RNDN);
        mpfr_add(r9184, r9177, r9183, MPFR_RNDN);
        mpfr_log(r9185, r9184, MPFR_RNDN);
        return mpfr_get_d(r9185, MPFR_RNDN);
}