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

char *name = "Hyperbolic arc-(co)secant";

double f_if(float x) {
        float r9277 = 1;
        float r9278 = x;
        float r9279 = r9277 / r9278;
        float r9280 = r9278 * r9278;
        float r9281 = r9277 - r9280;
        float r9282 = sqrt(r9281);
        float r9283 = r9282 / r9278;
        float r9284 = r9279 + r9283;
        float r9285 = log(r9284);
        return r9285;
}

double f_id(double x) {
        double r9286 = 1;
        double r9287 = x;
        double r9288 = r9286 / r9287;
        double r9289 = r9287 * r9287;
        double r9290 = r9286 - r9289;
        double r9291 = sqrt(r9290);
        double r9292 = r9291 / r9287;
        double r9293 = r9288 + r9292;
        double r9294 = log(r9293);
        return r9294;
}


double f_of(float x) {
        float r9295 = 1;
        float r9296 = x;
        float r9297 = r9295 / r9296;
        float r9298 = r9296 * r9296;
        float r9299 = r9295 - r9298;
        float r9300 = sqrt(r9299);
        float r9301 = r9300 / r9296;
        float r9302 = r9297 + r9301;
        float r9303 = sqrt(r9302);
        float r9304 = log(r9303);
        float r9305 = r9304 + r9304;
        return r9305;
}

double f_od(double x) {
        double r9306 = 1;
        double r9307 = x;
        double r9308 = r9306 / r9307;
        double r9309 = r9307 * r9307;
        double r9310 = r9306 - r9309;
        double r9311 = sqrt(r9310);
        double r9312 = r9311 / r9307;
        double r9313 = r9308 + r9312;
        double r9314 = sqrt(r9313);
        double r9315 = log(r9314);
        double r9316 = r9315 + r9315;
        return r9316;
}

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 r9317, r9318, r9319, r9320, r9321, r9322, r9323, r9324, r9325;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9317, "1", 10, MPFR_RNDN);
        mpfr_init(r9318);
        mpfr_init(r9319);
        mpfr_init(r9320);
        mpfr_init(r9321);
        mpfr_init(r9322);
        mpfr_init(r9323);
        mpfr_init(r9324);
        mpfr_init(r9325);
}

double f_im(double x) {
        ;
        mpfr_set_d(r9318, x, MPFR_RNDN);
        mpfr_div(r9319, r9317, r9318, MPFR_RNDN);
        mpfr_mul(r9320, r9318, r9318, MPFR_RNDN);
        mpfr_sub(r9321, r9317, r9320, MPFR_RNDN);
        mpfr_sqrt(r9322, r9321, MPFR_RNDN);
        mpfr_div(r9323, r9322, r9318, MPFR_RNDN);
        mpfr_add(r9324, r9319, r9323, MPFR_RNDN);
        mpfr_log(r9325, r9324, MPFR_RNDN);
        return mpfr_get_d(r9325, MPFR_RNDN);
}

static mpfr_t r9326, r9327, r9328, r9329, r9330, r9331, r9332, r9333, r9334, r9335, r9336;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9326, "1", 10, MPFR_RNDN);
        mpfr_init(r9327);
        mpfr_init(r9328);
        mpfr_init(r9329);
        mpfr_init(r9330);
        mpfr_init(r9331);
        mpfr_init(r9332);
        mpfr_init(r9333);
        mpfr_init(r9334);
        mpfr_init(r9335);
        mpfr_init(r9336);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r9327, x, MPFR_RNDN);
        mpfr_div(r9328, r9326, r9327, MPFR_RNDN);
        mpfr_mul(r9329, r9327, r9327, MPFR_RNDN);
        mpfr_sub(r9330, r9326, r9329, MPFR_RNDN);
        mpfr_sqrt(r9331, r9330, MPFR_RNDN);
        mpfr_div(r9332, r9331, r9327, MPFR_RNDN);
        mpfr_add(r9333, r9328, r9332, MPFR_RNDN);
        mpfr_sqrt(r9334, r9333, MPFR_RNDN);
        mpfr_log(r9335, r9334, MPFR_RNDN);
        mpfr_add(r9336, r9335, r9335, MPFR_RNDN);
        return mpfr_get_d(r9336, MPFR_RNDN);
}

static mpfr_t r9337, r9338, r9339, r9340, r9341, r9342, r9343, r9344, r9345, r9346, r9347;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9337, "1", 10, MPFR_RNDN);
        mpfr_init(r9338);
        mpfr_init(r9339);
        mpfr_init(r9340);
        mpfr_init(r9341);
        mpfr_init(r9342);
        mpfr_init(r9343);
        mpfr_init(r9344);
        mpfr_init(r9345);
        mpfr_init(r9346);
        mpfr_init(r9347);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r9338, x, MPFR_RNDN);
        mpfr_div(r9339, r9337, r9338, MPFR_RNDN);
        mpfr_mul(r9340, r9338, r9338, MPFR_RNDN);
        mpfr_sub(r9341, r9337, r9340, MPFR_RNDN);
        mpfr_sqrt(r9342, r9341, MPFR_RNDN);
        mpfr_div(r9343, r9342, r9338, MPFR_RNDN);
        mpfr_add(r9344, r9339, r9343, MPFR_RNDN);
        mpfr_sqrt(r9345, r9344, MPFR_RNDN);
        mpfr_log(r9346, r9345, MPFR_RNDN);
        mpfr_add(r9347, r9346, r9346, MPFR_RNDN);
        return mpfr_get_d(r9347, MPFR_RNDN);
}