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

char *name = "Hyperbolic arcsine";

double f_if(float x) {
        float r27323 = x;
        float r27324 = r27323 * r27323;
        float r27325 = 1;
        float r27326 = r27324 + r27325;
        float r27327 = sqrt(r27326);
        float r27328 = r27323 + r27327;
        float r27329 = log(r27328);
        return r27329;
}

double f_id(double x) {
        double r27330 = x;
        double r27331 = r27330 * r27330;
        double r27332 = 1;
        double r27333 = r27331 + r27332;
        double r27334 = sqrt(r27333);
        double r27335 = r27330 + r27334;
        double r27336 = log(r27335);
        return r27336;
}


double f_of(float x) {
        float r27337 = x;
        float r27338 = -1.143928309884049;
        bool r27339 = r27337 <= r27338;
        float r27340 = 1/8;
        float r27341 = r27340 / r27337;
        float r27342 = r27337 * r27337;
        float r27343 = r27341 / r27342;
        float r27344 = 1/2;
        float r27345 = r27344 / r27337;
        float r27346 = r27343 - r27345;
        float r27347 = log(r27346);
        float r27348 = 0.9637119404562783;
        bool r27349 = r27337 <= r27348;
        float r27350 = 3/40;
        float r27351 = 5;
        float r27352 = pow(r27337, r27351);
        float r27353 = r27350 * r27352;
        float r27354 = r27353 + r27337;
        float r27355 = 1/6;
        float r27356 = 3;
        float r27357 = pow(r27337, r27356);
        float r27358 = r27355 * r27357;
        float r27359 = r27354 - r27358;
        float r27360 = r27337 + r27337;
        float r27361 = r27360 + r27345;
        float r27362 = r27361 - r27343;
        float r27363 = log(r27362);
        float r27364 = r27349 ? r27359 : r27363;
        float r27365 = r27339 ? r27347 : r27364;
        return r27365;
}

double f_od(double x) {
        double r27366 = x;
        double r27367 = -1.143928309884049;
        bool r27368 = r27366 <= r27367;
        double r27369 = 1/8;
        double r27370 = r27369 / r27366;
        double r27371 = r27366 * r27366;
        double r27372 = r27370 / r27371;
        double r27373 = 1/2;
        double r27374 = r27373 / r27366;
        double r27375 = r27372 - r27374;
        double r27376 = log(r27375);
        double r27377 = 0.9637119404562783;
        bool r27378 = r27366 <= r27377;
        double r27379 = 3/40;
        double r27380 = 5;
        double r27381 = pow(r27366, r27380);
        double r27382 = r27379 * r27381;
        double r27383 = r27382 + r27366;
        double r27384 = 1/6;
        double r27385 = 3;
        double r27386 = pow(r27366, r27385);
        double r27387 = r27384 * r27386;
        double r27388 = r27383 - r27387;
        double r27389 = r27366 + r27366;
        double r27390 = r27389 + r27374;
        double r27391 = r27390 - r27372;
        double r27392 = log(r27391);
        double r27393 = r27378 ? r27388 : r27392;
        double r27394 = r27368 ? r27376 : r27393;
        return r27394;
}

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 r27395, r27396, r27397, r27398, r27399, r27400, r27401;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27395);
        mpfr_init(r27396);
        mpfr_init_set_str(r27397, "1", 10, MPFR_RNDN);
        mpfr_init(r27398);
        mpfr_init(r27399);
        mpfr_init(r27400);
        mpfr_init(r27401);
}

double f_im(double x) {
        mpfr_set_d(r27395, x, MPFR_RNDN);
        mpfr_mul(r27396, r27395, r27395, MPFR_RNDN);
        ;
        mpfr_add(r27398, r27396, r27397, MPFR_RNDN);
        mpfr_sqrt(r27399, r27398, MPFR_RNDN);
        mpfr_add(r27400, r27395, r27399, MPFR_RNDN);
        mpfr_log(r27401, r27400, MPFR_RNDN);
        return mpfr_get_d(r27401, MPFR_RNDN);
}

static mpfr_t r27402, r27403, r27404, r27405, r27406, r27407, r27408, r27409, r27410, r27411, r27412, r27413, r27414, r27415, r27416, r27417, r27418, r27419, r27420, r27421, r27422, r27423, r27424, r27425, r27426, r27427, r27428, r27429, r27430;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27402);
        mpfr_init_set_str(r27403, "-1.143928309884049", 10, MPFR_RNDN);
        mpfr_init(r27404);
        mpfr_init_set_str(r27405, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27406);
        mpfr_init(r27407);
        mpfr_init(r27408);
        mpfr_init_set_str(r27409, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27410);
        mpfr_init(r27411);
        mpfr_init(r27412);
        mpfr_init_set_str(r27413, "0.9637119404562783", 10, MPFR_RNDN);
        mpfr_init(r27414);
        mpfr_init_set_str(r27415, "3/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r27416, "5", 10, MPFR_RNDN);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init_set_str(r27420, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27421, "3", 10, MPFR_RNDN);
        mpfr_init(r27422);
        mpfr_init(r27423);
        mpfr_init(r27424);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init(r27427);
        mpfr_init(r27428);
        mpfr_init(r27429);
        mpfr_init(r27430);
}

double f_fm(double x) {
        mpfr_set_d(r27402, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27404, mpfr_cmp(r27402, r27403) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27406, r27405, r27402, MPFR_RNDN);
        mpfr_mul(r27407, r27402, r27402, MPFR_RNDN);
        mpfr_div(r27408, r27406, r27407, MPFR_RNDN);
        ;
        mpfr_div(r27410, r27409, r27402, MPFR_RNDN);
        mpfr_sub(r27411, r27408, r27410, MPFR_RNDN);
        mpfr_log(r27412, r27411, MPFR_RNDN);
        ;
        mpfr_set_si(r27414, mpfr_cmp(r27402, r27413) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27417, r27402, r27416, MPFR_RNDN);
        mpfr_mul(r27418, r27415, r27417, MPFR_RNDN);
        mpfr_add(r27419, r27418, r27402, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27422, r27402, r27421, MPFR_RNDN);
        mpfr_mul(r27423, r27420, r27422, MPFR_RNDN);
        mpfr_sub(r27424, r27419, r27423, MPFR_RNDN);
        mpfr_add(r27425, r27402, r27402, MPFR_RNDN);
        mpfr_add(r27426, r27425, r27410, MPFR_RNDN);
        mpfr_sub(r27427, r27426, r27408, MPFR_RNDN);
        mpfr_log(r27428, r27427, MPFR_RNDN);
        if (mpfr_get_si(r27414, MPFR_RNDN)) { mpfr_set(r27429, r27424, MPFR_RNDN); } else { mpfr_set(r27429, r27428, MPFR_RNDN); };
        if (mpfr_get_si(r27404, MPFR_RNDN)) { mpfr_set(r27430, r27412, MPFR_RNDN); } else { mpfr_set(r27430, r27429, MPFR_RNDN); };
        return mpfr_get_d(r27430, MPFR_RNDN);
}

static mpfr_t r27431, r27432, r27433, r27434, r27435, r27436, r27437, r27438, r27439, r27440, r27441, r27442, r27443, r27444, r27445, r27446, r27447, r27448, r27449, r27450, r27451, r27452, r27453, r27454, r27455, r27456, r27457, r27458, r27459;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27431);
        mpfr_init_set_str(r27432, "-1.143928309884049", 10, MPFR_RNDN);
        mpfr_init(r27433);
        mpfr_init_set_str(r27434, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27435);
        mpfr_init(r27436);
        mpfr_init(r27437);
        mpfr_init_set_str(r27438, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27439);
        mpfr_init(r27440);
        mpfr_init(r27441);
        mpfr_init_set_str(r27442, "0.9637119404562783", 10, MPFR_RNDN);
        mpfr_init(r27443);
        mpfr_init_set_str(r27444, "3/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r27445, "5", 10, MPFR_RNDN);
        mpfr_init(r27446);
        mpfr_init(r27447);
        mpfr_init(r27448);
        mpfr_init_set_str(r27449, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27450, "3", 10, MPFR_RNDN);
        mpfr_init(r27451);
        mpfr_init(r27452);
        mpfr_init(r27453);
        mpfr_init(r27454);
        mpfr_init(r27455);
        mpfr_init(r27456);
        mpfr_init(r27457);
        mpfr_init(r27458);
        mpfr_init(r27459);
}

double f_dm(double x) {
        mpfr_set_d(r27431, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27433, mpfr_cmp(r27431, r27432) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27435, r27434, r27431, MPFR_RNDN);
        mpfr_mul(r27436, r27431, r27431, MPFR_RNDN);
        mpfr_div(r27437, r27435, r27436, MPFR_RNDN);
        ;
        mpfr_div(r27439, r27438, r27431, MPFR_RNDN);
        mpfr_sub(r27440, r27437, r27439, MPFR_RNDN);
        mpfr_log(r27441, r27440, MPFR_RNDN);
        ;
        mpfr_set_si(r27443, mpfr_cmp(r27431, r27442) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27446, r27431, r27445, MPFR_RNDN);
        mpfr_mul(r27447, r27444, r27446, MPFR_RNDN);
        mpfr_add(r27448, r27447, r27431, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27451, r27431, r27450, MPFR_RNDN);
        mpfr_mul(r27452, r27449, r27451, MPFR_RNDN);
        mpfr_sub(r27453, r27448, r27452, MPFR_RNDN);
        mpfr_add(r27454, r27431, r27431, MPFR_RNDN);
        mpfr_add(r27455, r27454, r27439, MPFR_RNDN);
        mpfr_sub(r27456, r27455, r27437, MPFR_RNDN);
        mpfr_log(r27457, r27456, MPFR_RNDN);
        if (mpfr_get_si(r27443, MPFR_RNDN)) { mpfr_set(r27458, r27453, MPFR_RNDN); } else { mpfr_set(r27458, r27457, MPFR_RNDN); };
        if (mpfr_get_si(r27433, MPFR_RNDN)) { mpfr_set(r27459, r27441, MPFR_RNDN); } else { mpfr_set(r27459, r27458, MPFR_RNDN); };
        return mpfr_get_d(r27459, MPFR_RNDN);
}