#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 r27239 = x;
        float r27240 = r27239 * r27239;
        float r27241 = 1;
        float r27242 = r27240 + r27241;
        float r27243 = sqrt(r27242);
        float r27244 = r27239 + r27243;
        float r27245 = log(r27244);
        return r27245;
}

double f_id(double x) {
        double r27246 = x;
        double r27247 = r27246 * r27246;
        double r27248 = 1;
        double r27249 = r27247 + r27248;
        double r27250 = sqrt(r27249);
        double r27251 = r27246 + r27250;
        double r27252 = log(r27251);
        return r27252;
}


double f_of(float x) {
        float r27253 = x;
        float r27254 = -1.143928309884049;
        bool r27255 = r27253 <= r27254;
        float r27256 = 1/8;
        float r27257 = r27256 / r27253;
        float r27258 = r27253 * r27253;
        float r27259 = r27257 / r27258;
        float r27260 = 1/2;
        float r27261 = r27260 / r27253;
        float r27262 = r27259 - r27261;
        float r27263 = log(r27262);
        float r27264 = 0.9637119404562783;
        bool r27265 = r27253 <= r27264;
        float r27266 = 3/40;
        float r27267 = 5;
        float r27268 = pow(r27253, r27267);
        float r27269 = r27266 * r27268;
        float r27270 = r27269 + r27253;
        float r27271 = 1/6;
        float r27272 = 3;
        float r27273 = pow(r27253, r27272);
        float r27274 = r27271 * r27273;
        float r27275 = r27270 - r27274;
        float r27276 = r27253 + r27253;
        float r27277 = r27276 + r27261;
        float r27278 = r27277 - r27259;
        float r27279 = log(r27278);
        float r27280 = r27265 ? r27275 : r27279;
        float r27281 = r27255 ? r27263 : r27280;
        return r27281;
}

double f_od(double x) {
        double r27282 = x;
        double r27283 = -1.143928309884049;
        bool r27284 = r27282 <= r27283;
        double r27285 = 1/8;
        double r27286 = r27285 / r27282;
        double r27287 = r27282 * r27282;
        double r27288 = r27286 / r27287;
        double r27289 = 1/2;
        double r27290 = r27289 / r27282;
        double r27291 = r27288 - r27290;
        double r27292 = log(r27291);
        double r27293 = 0.9637119404562783;
        bool r27294 = r27282 <= r27293;
        double r27295 = 3/40;
        double r27296 = 5;
        double r27297 = pow(r27282, r27296);
        double r27298 = r27295 * r27297;
        double r27299 = r27298 + r27282;
        double r27300 = 1/6;
        double r27301 = 3;
        double r27302 = pow(r27282, r27301);
        double r27303 = r27300 * r27302;
        double r27304 = r27299 - r27303;
        double r27305 = r27282 + r27282;
        double r27306 = r27305 + r27290;
        double r27307 = r27306 - r27288;
        double r27308 = log(r27307);
        double r27309 = r27294 ? r27304 : r27308;
        double r27310 = r27284 ? r27292 : r27309;
        return r27310;
}

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 r27311, r27312, r27313, r27314, r27315, r27316, r27317;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27311);
        mpfr_init(r27312);
        mpfr_init_set_str(r27313, "1", 10, MPFR_RNDN);
        mpfr_init(r27314);
        mpfr_init(r27315);
        mpfr_init(r27316);
        mpfr_init(r27317);
}

double f_im(double x) {
        mpfr_set_d(r27311, x, MPFR_RNDN);
        mpfr_mul(r27312, r27311, r27311, MPFR_RNDN);
        ;
        mpfr_add(r27314, r27312, r27313, MPFR_RNDN);
        mpfr_sqrt(r27315, r27314, MPFR_RNDN);
        mpfr_add(r27316, r27311, r27315, MPFR_RNDN);
        mpfr_log(r27317, r27316, MPFR_RNDN);
        return mpfr_get_d(r27317, MPFR_RNDN);
}

static mpfr_t r27318, r27319, r27320, r27321, r27322, r27323, r27324, r27325, r27326, r27327, r27328, r27329, r27330, r27331, r27332, r27333, r27334, r27335, r27336, r27337, r27338, r27339, r27340, r27341, r27342, r27343, r27344, r27345, r27346;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27318);
        mpfr_init_set_str(r27319, "-1.143928309884049", 10, MPFR_RNDN);
        mpfr_init(r27320);
        mpfr_init_set_str(r27321, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27322);
        mpfr_init(r27323);
        mpfr_init(r27324);
        mpfr_init_set_str(r27325, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27326);
        mpfr_init(r27327);
        mpfr_init(r27328);
        mpfr_init_set_str(r27329, "0.9637119404562783", 10, MPFR_RNDN);
        mpfr_init(r27330);
        mpfr_init_set_str(r27331, "3/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r27332, "5", 10, MPFR_RNDN);
        mpfr_init(r27333);
        mpfr_init(r27334);
        mpfr_init(r27335);
        mpfr_init_set_str(r27336, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27337, "3", 10, MPFR_RNDN);
        mpfr_init(r27338);
        mpfr_init(r27339);
        mpfr_init(r27340);
        mpfr_init(r27341);
        mpfr_init(r27342);
        mpfr_init(r27343);
        mpfr_init(r27344);
        mpfr_init(r27345);
        mpfr_init(r27346);
}

double f_fm(double x) {
        mpfr_set_d(r27318, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27320, mpfr_cmp(r27318, r27319) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27322, r27321, r27318, MPFR_RNDN);
        mpfr_mul(r27323, r27318, r27318, MPFR_RNDN);
        mpfr_div(r27324, r27322, r27323, MPFR_RNDN);
        ;
        mpfr_div(r27326, r27325, r27318, MPFR_RNDN);
        mpfr_sub(r27327, r27324, r27326, MPFR_RNDN);
        mpfr_log(r27328, r27327, MPFR_RNDN);
        ;
        mpfr_set_si(r27330, mpfr_cmp(r27318, r27329) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27333, r27318, r27332, MPFR_RNDN);
        mpfr_mul(r27334, r27331, r27333, MPFR_RNDN);
        mpfr_add(r27335, r27334, r27318, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27338, r27318, r27337, MPFR_RNDN);
        mpfr_mul(r27339, r27336, r27338, MPFR_RNDN);
        mpfr_sub(r27340, r27335, r27339, MPFR_RNDN);
        mpfr_add(r27341, r27318, r27318, MPFR_RNDN);
        mpfr_add(r27342, r27341, r27326, MPFR_RNDN);
        mpfr_sub(r27343, r27342, r27324, MPFR_RNDN);
        mpfr_log(r27344, r27343, MPFR_RNDN);
        if (mpfr_get_si(r27330, MPFR_RNDN)) { mpfr_set(r27345, r27340, MPFR_RNDN); } else { mpfr_set(r27345, r27344, MPFR_RNDN); };
        if (mpfr_get_si(r27320, MPFR_RNDN)) { mpfr_set(r27346, r27328, MPFR_RNDN); } else { mpfr_set(r27346, r27345, MPFR_RNDN); };
        return mpfr_get_d(r27346, MPFR_RNDN);
}

static mpfr_t r27347, r27348, r27349, r27350, r27351, r27352, r27353, r27354, r27355, r27356, r27357, r27358, r27359, r27360, r27361, r27362, r27363, r27364, r27365, r27366, r27367, r27368, r27369, r27370, r27371, r27372, r27373, r27374, r27375;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27347);
        mpfr_init_set_str(r27348, "-1.143928309884049", 10, MPFR_RNDN);
        mpfr_init(r27349);
        mpfr_init_set_str(r27350, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27351);
        mpfr_init(r27352);
        mpfr_init(r27353);
        mpfr_init_set_str(r27354, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27355);
        mpfr_init(r27356);
        mpfr_init(r27357);
        mpfr_init_set_str(r27358, "0.9637119404562783", 10, MPFR_RNDN);
        mpfr_init(r27359);
        mpfr_init_set_str(r27360, "3/40", 10, MPFR_RNDN);
        mpfr_init_set_str(r27361, "5", 10, MPFR_RNDN);
        mpfr_init(r27362);
        mpfr_init(r27363);
        mpfr_init(r27364);
        mpfr_init_set_str(r27365, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27366, "3", 10, MPFR_RNDN);
        mpfr_init(r27367);
        mpfr_init(r27368);
        mpfr_init(r27369);
        mpfr_init(r27370);
        mpfr_init(r27371);
        mpfr_init(r27372);
        mpfr_init(r27373);
        mpfr_init(r27374);
        mpfr_init(r27375);
}

double f_dm(double x) {
        mpfr_set_d(r27347, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27349, mpfr_cmp(r27347, r27348) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27351, r27350, r27347, MPFR_RNDN);
        mpfr_mul(r27352, r27347, r27347, MPFR_RNDN);
        mpfr_div(r27353, r27351, r27352, MPFR_RNDN);
        ;
        mpfr_div(r27355, r27354, r27347, MPFR_RNDN);
        mpfr_sub(r27356, r27353, r27355, MPFR_RNDN);
        mpfr_log(r27357, r27356, MPFR_RNDN);
        ;
        mpfr_set_si(r27359, mpfr_cmp(r27347, r27358) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27362, r27347, r27361, MPFR_RNDN);
        mpfr_mul(r27363, r27360, r27362, MPFR_RNDN);
        mpfr_add(r27364, r27363, r27347, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27367, r27347, r27366, MPFR_RNDN);
        mpfr_mul(r27368, r27365, r27367, MPFR_RNDN);
        mpfr_sub(r27369, r27364, r27368, MPFR_RNDN);
        mpfr_add(r27370, r27347, r27347, MPFR_RNDN);
        mpfr_add(r27371, r27370, r27355, MPFR_RNDN);
        mpfr_sub(r27372, r27371, r27353, MPFR_RNDN);
        mpfr_log(r27373, r27372, MPFR_RNDN);
        if (mpfr_get_si(r27359, MPFR_RNDN)) { mpfr_set(r27374, r27369, MPFR_RNDN); } else { mpfr_set(r27374, r27373, MPFR_RNDN); };
        if (mpfr_get_si(r27349, MPFR_RNDN)) { mpfr_set(r27375, r27357, MPFR_RNDN); } else { mpfr_set(r27375, r27374, MPFR_RNDN); };
        return mpfr_get_d(r27375, MPFR_RNDN);
}