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

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

double f_if(float x) {
        float r31238 = 1;
        float r31239 = 2;
        float r31240 = r31238 / r31239;
        float r31241 = x;
        float r31242 = r31238 + r31241;
        float r31243 = r31238 - r31241;
        float r31244 = r31242 / r31243;
        float r31245 = log(r31244);
        float r31246 = r31240 * r31245;
        return r31246;
}

double f_id(double x) {
        double r31247 = 1;
        double r31248 = 2;
        double r31249 = r31247 / r31248;
        double r31250 = x;
        double r31251 = r31247 + r31250;
        double r31252 = r31247 - r31250;
        double r31253 = r31251 / r31252;
        double r31254 = log(r31253);
        double r31255 = r31249 * r31254;
        return r31255;
}


double f_of(float x) {
        float r31256 = 2/3;
        float r31257 = x;
        float r31258 = 3;
        float r31259 = pow(r31257, r31258);
        float r31260 = r31256 * r31259;
        float r31261 = 2/5;
        float r31262 = 5;
        float r31263 = pow(r31257, r31262);
        float r31264 = r31261 * r31263;
        float r31265 = 2;
        float r31266 = r31265 * r31257;
        float r31267 = r31264 + r31266;
        float r31268 = r31260 + r31267;
        float r31269 = r31268 / r31265;
        return r31269;
}

double f_od(double x) {
        double r31270 = 2/3;
        double r31271 = x;
        double r31272 = 3;
        double r31273 = pow(r31271, r31272);
        double r31274 = r31270 * r31273;
        double r31275 = 2/5;
        double r31276 = 5;
        double r31277 = pow(r31271, r31276);
        double r31278 = r31275 * r31277;
        double r31279 = 2;
        double r31280 = r31279 * r31271;
        double r31281 = r31278 + r31280;
        double r31282 = r31274 + r31281;
        double r31283 = r31282 / r31279;
        return r31283;
}

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 r31284, r31285, r31286, r31287, r31288, r31289, r31290, r31291, r31292;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r31284, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r31285, "2", 10, MPFR_RNDN);
        mpfr_init(r31286);
        mpfr_init(r31287);
        mpfr_init(r31288);
        mpfr_init(r31289);
        mpfr_init(r31290);
        mpfr_init(r31291);
        mpfr_init(r31292);
}

double f_im(double x) {
        ;
        ;
        mpfr_div(r31286, r31284, r31285, MPFR_RNDN);
        mpfr_set_d(r31287, x, MPFR_RNDN);
        mpfr_add(r31288, r31284, r31287, MPFR_RNDN);
        mpfr_sub(r31289, r31284, r31287, MPFR_RNDN);
        mpfr_div(r31290, r31288, r31289, MPFR_RNDN);
        mpfr_log(r31291, r31290, MPFR_RNDN);
        mpfr_mul(r31292, r31286, r31291, MPFR_RNDN);
        return mpfr_get_d(r31292, MPFR_RNDN);
}

static mpfr_t r31293, r31294, r31295, r31296, r31297, r31298, r31299, r31300, r31301, r31302, r31303, r31304, r31305, r31306;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r31293, "2/3", 10, MPFR_RNDN);
        mpfr_init(r31294);
        mpfr_init_set_str(r31295, "3", 10, MPFR_RNDN);
        mpfr_init(r31296);
        mpfr_init(r31297);
        mpfr_init_set_str(r31298, "2/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r31299, "5", 10, MPFR_RNDN);
        mpfr_init(r31300);
        mpfr_init(r31301);
        mpfr_init_set_str(r31302, "2", 10, MPFR_RNDN);
        mpfr_init(r31303);
        mpfr_init(r31304);
        mpfr_init(r31305);
        mpfr_init(r31306);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r31294, x, MPFR_RNDN);
        ;
        mpfr_pow(r31296, r31294, r31295, MPFR_RNDN);
        mpfr_mul(r31297, r31293, r31296, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r31300, r31294, r31299, MPFR_RNDN);
        mpfr_mul(r31301, r31298, r31300, MPFR_RNDN);
        ;
        mpfr_mul(r31303, r31302, r31294, MPFR_RNDN);
        mpfr_add(r31304, r31301, r31303, MPFR_RNDN);
        mpfr_add(r31305, r31297, r31304, MPFR_RNDN);
        mpfr_div(r31306, r31305, r31302, MPFR_RNDN);
        return mpfr_get_d(r31306, MPFR_RNDN);
}

static mpfr_t r31307, r31308, r31309, r31310, r31311, r31312, r31313, r31314, r31315, r31316, r31317, r31318, r31319, r31320;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r31307, "2/3", 10, MPFR_RNDN);
        mpfr_init(r31308);
        mpfr_init_set_str(r31309, "3", 10, MPFR_RNDN);
        mpfr_init(r31310);
        mpfr_init(r31311);
        mpfr_init_set_str(r31312, "2/5", 10, MPFR_RNDN);
        mpfr_init_set_str(r31313, "5", 10, MPFR_RNDN);
        mpfr_init(r31314);
        mpfr_init(r31315);
        mpfr_init_set_str(r31316, "2", 10, MPFR_RNDN);
        mpfr_init(r31317);
        mpfr_init(r31318);
        mpfr_init(r31319);
        mpfr_init(r31320);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r31308, x, MPFR_RNDN);
        ;
        mpfr_pow(r31310, r31308, r31309, MPFR_RNDN);
        mpfr_mul(r31311, r31307, r31310, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r31314, r31308, r31313, MPFR_RNDN);
        mpfr_mul(r31315, r31312, r31314, MPFR_RNDN);
        ;
        mpfr_mul(r31317, r31316, r31308, MPFR_RNDN);
        mpfr_add(r31318, r31315, r31317, MPFR_RNDN);
        mpfr_add(r31319, r31311, r31318, MPFR_RNDN);
        mpfr_div(r31320, r31319, r31316, MPFR_RNDN);
        return mpfr_get_d(r31320, MPFR_RNDN);
}