#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 r25306 = 1;
        float r25307 = x;
        float r25308 = r25306 / r25307;
        float r25309 = r25307 * r25307;
        float r25310 = r25306 - r25309;
        float r25311 = sqrt(r25310);
        float r25312 = r25311 / r25307;
        float r25313 = r25308 + r25312;
        float r25314 = log(r25313);
        return r25314;
}

double f_id(double x) {
        double r25315 = 1;
        double r25316 = x;
        double r25317 = r25315 / r25316;
        double r25318 = r25316 * r25316;
        double r25319 = r25315 - r25318;
        double r25320 = sqrt(r25319);
        double r25321 = r25320 / r25316;
        double r25322 = r25317 + r25321;
        double r25323 = log(r25322);
        return r25323;
}


double f_of(float x) {
        float r25324 = 1;
        float r25325 = x;
        float r25326 = r25324 / r25325;
        float r25327 = r25325 * r25325;
        float r25328 = r25324 - r25327;
        float r25329 = sqrt(r25328);
        float r25330 = r25329 / r25325;
        float r25331 = r25326 + r25330;
        float r25332 = log(r25331);
        return r25332;
}

double f_od(double x) {
        double r25333 = 1;
        double r25334 = x;
        double r25335 = r25333 / r25334;
        double r25336 = r25334 * r25334;
        double r25337 = r25333 - r25336;
        double r25338 = sqrt(r25337);
        double r25339 = r25338 / r25334;
        double r25340 = r25335 + r25339;
        double r25341 = log(r25340);
        return r25341;
}

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 r25342, r25343, r25344, r25345, r25346, r25347, r25348, r25349, r25350;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25342, "1", 10, MPFR_RNDN);
        mpfr_init(r25343);
        mpfr_init(r25344);
        mpfr_init(r25345);
        mpfr_init(r25346);
        mpfr_init(r25347);
        mpfr_init(r25348);
        mpfr_init(r25349);
        mpfr_init(r25350);
}

double f_im(double x) {
        ;
        mpfr_set_d(r25343, x, MPFR_RNDN);
        mpfr_div(r25344, r25342, r25343, MPFR_RNDN);
        mpfr_mul(r25345, r25343, r25343, MPFR_RNDN);
        mpfr_sub(r25346, r25342, r25345, MPFR_RNDN);
        mpfr_sqrt(r25347, r25346, MPFR_RNDN);
        mpfr_div(r25348, r25347, r25343, MPFR_RNDN);
        mpfr_add(r25349, r25344, r25348, MPFR_RNDN);
        mpfr_log(r25350, r25349, MPFR_RNDN);
        return mpfr_get_d(r25350, MPFR_RNDN);
}

static mpfr_t r25351, r25352, r25353, r25354, r25355, r25356, r25357, r25358, r25359;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25351, "1", 10, MPFR_RNDN);
        mpfr_init(r25352);
        mpfr_init(r25353);
        mpfr_init(r25354);
        mpfr_init(r25355);
        mpfr_init(r25356);
        mpfr_init(r25357);
        mpfr_init(r25358);
        mpfr_init(r25359);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r25352, x, MPFR_RNDN);
        mpfr_div(r25353, r25351, r25352, MPFR_RNDN);
        mpfr_mul(r25354, r25352, r25352, MPFR_RNDN);
        mpfr_sub(r25355, r25351, r25354, MPFR_RNDN);
        mpfr_sqrt(r25356, r25355, MPFR_RNDN);
        mpfr_div(r25357, r25356, r25352, MPFR_RNDN);
        mpfr_add(r25358, r25353, r25357, MPFR_RNDN);
        mpfr_log(r25359, r25358, MPFR_RNDN);
        return mpfr_get_d(r25359, MPFR_RNDN);
}

static mpfr_t r25360, r25361, r25362, r25363, r25364, r25365, r25366, r25367, r25368;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r25360, "1", 10, MPFR_RNDN);
        mpfr_init(r25361);
        mpfr_init(r25362);
        mpfr_init(r25363);
        mpfr_init(r25364);
        mpfr_init(r25365);
        mpfr_init(r25366);
        mpfr_init(r25367);
        mpfr_init(r25368);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r25361, x, MPFR_RNDN);
        mpfr_div(r25362, r25360, r25361, MPFR_RNDN);
        mpfr_mul(r25363, r25361, r25361, MPFR_RNDN);
        mpfr_sub(r25364, r25360, r25363, MPFR_RNDN);
        mpfr_sqrt(r25365, r25364, MPFR_RNDN);
        mpfr_div(r25366, r25365, r25361, MPFR_RNDN);
        mpfr_add(r25367, r25362, r25366, MPFR_RNDN);
        mpfr_log(r25368, r25367, MPFR_RNDN);
        return mpfr_get_d(r25368, MPFR_RNDN);
}