#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 r31241 = 1;
        float r31242 = x;
        float r31243 = r31241 / r31242;
        float r31244 = r31242 * r31242;
        float r31245 = r31241 - r31244;
        float r31246 = sqrt(r31245);
        float r31247 = r31246 / r31242;
        float r31248 = r31243 + r31247;
        float r31249 = log(r31248);
        return r31249;
}

double f_id(double x) {
        double r31250 = 1;
        double r31251 = x;
        double r31252 = r31250 / r31251;
        double r31253 = r31251 * r31251;
        double r31254 = r31250 - r31253;
        double r31255 = sqrt(r31254);
        double r31256 = r31255 / r31251;
        double r31257 = r31252 + r31256;
        double r31258 = log(r31257);
        return r31258;
}


double f_of(float x) {
        float r31259 = 1;
        float r31260 = x;
        float r31261 = r31259 / r31260;
        float r31262 = r31260 * r31260;
        float r31263 = r31259 - r31262;
        float r31264 = sqrt(r31263);
        float r31265 = r31264 / r31260;
        float r31266 = r31261 + r31265;
        float r31267 = log(r31266);
        return r31267;
}

double f_od(double x) {
        double r31268 = 1;
        double r31269 = x;
        double r31270 = r31268 / r31269;
        double r31271 = r31269 * r31269;
        double r31272 = r31268 - r31271;
        double r31273 = sqrt(r31272);
        double r31274 = r31273 / r31269;
        double r31275 = r31270 + r31274;
        double r31276 = log(r31275);
        return r31276;
}

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 r31277, r31278, r31279, r31280, r31281, r31282, r31283, r31284, r31285;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r31277, "1", 10, MPFR_RNDN);
        mpfr_init(r31278);
        mpfr_init(r31279);
        mpfr_init(r31280);
        mpfr_init(r31281);
        mpfr_init(r31282);
        mpfr_init(r31283);
        mpfr_init(r31284);
        mpfr_init(r31285);
}

double f_im(double x) {
        ;
        mpfr_set_d(r31278, x, MPFR_RNDN);
        mpfr_div(r31279, r31277, r31278, MPFR_RNDN);
        mpfr_mul(r31280, r31278, r31278, MPFR_RNDN);
        mpfr_sub(r31281, r31277, r31280, MPFR_RNDN);
        mpfr_sqrt(r31282, r31281, MPFR_RNDN);
        mpfr_div(r31283, r31282, r31278, MPFR_RNDN);
        mpfr_add(r31284, r31279, r31283, MPFR_RNDN);
        mpfr_log(r31285, r31284, MPFR_RNDN);
        return mpfr_get_d(r31285, MPFR_RNDN);
}

static mpfr_t r31286, r31287, r31288, r31289, r31290, r31291, r31292, r31293, r31294;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r31286, "1", 10, MPFR_RNDN);
        mpfr_init(r31287);
        mpfr_init(r31288);
        mpfr_init(r31289);
        mpfr_init(r31290);
        mpfr_init(r31291);
        mpfr_init(r31292);
        mpfr_init(r31293);
        mpfr_init(r31294);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r31287, x, MPFR_RNDN);
        mpfr_div(r31288, r31286, r31287, MPFR_RNDN);
        mpfr_mul(r31289, r31287, r31287, MPFR_RNDN);
        mpfr_sub(r31290, r31286, r31289, MPFR_RNDN);
        mpfr_sqrt(r31291, r31290, MPFR_RNDN);
        mpfr_div(r31292, r31291, r31287, MPFR_RNDN);
        mpfr_add(r31293, r31288, r31292, MPFR_RNDN);
        mpfr_log(r31294, r31293, MPFR_RNDN);
        return mpfr_get_d(r31294, MPFR_RNDN);
}

static mpfr_t r31295, r31296, r31297, r31298, r31299, r31300, r31301, r31302, r31303;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r31295, "1", 10, MPFR_RNDN);
        mpfr_init(r31296);
        mpfr_init(r31297);
        mpfr_init(r31298);
        mpfr_init(r31299);
        mpfr_init(r31300);
        mpfr_init(r31301);
        mpfr_init(r31302);
        mpfr_init(r31303);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r31296, x, MPFR_RNDN);
        mpfr_div(r31297, r31295, r31296, MPFR_RNDN);
        mpfr_mul(r31298, r31296, r31296, MPFR_RNDN);
        mpfr_sub(r31299, r31295, r31298, MPFR_RNDN);
        mpfr_sqrt(r31300, r31299, MPFR_RNDN);
        mpfr_div(r31301, r31300, r31296, MPFR_RNDN);
        mpfr_add(r31302, r31297, r31301, MPFR_RNDN);
        mpfr_log(r31303, r31302, MPFR_RNDN);
        return mpfr_get_d(r31303, MPFR_RNDN);
}