#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 r27461 = x;
        float r27462 = r27461 * r27461;
        float r27463 = 1;
        float r27464 = r27462 + r27463;
        float r27465 = sqrt(r27464);
        float r27466 = r27461 + r27465;
        float r27467 = log(r27466);
        return r27467;
}

double f_id(double x) {
        double r27468 = x;
        double r27469 = r27468 * r27468;
        double r27470 = 1;
        double r27471 = r27469 + r27470;
        double r27472 = sqrt(r27471);
        double r27473 = r27468 + r27472;
        double r27474 = log(r27473);
        return r27474;
}


double f_of(float x) {
        float r27475 = x;
        float r27476 = -1.0916759748291127;
        bool r27477 = r27475 <= r27476;
        float r27478 = 1;
        float r27479 = r27478 / r27475;
        float r27480 = 1/8;
        float r27481 = r27480 / r27475;
        float r27482 = 1/2;
        float r27483 = -r27482;
        float r27484 = fma(r27479, r27481, r27483);
        float r27485 = 1/16;
        float r27486 = -r27485;
        float r27487 = 5;
        float r27488 = pow(r27475, r27487);
        float r27489 = r27486 / r27488;
        float r27490 = fma(r27479, r27484, r27489);
        float r27491 = log(r27490);
        float r27492 = 0.007483018826987035;
        bool r27493 = r27475 <= r27492;
        float r27494 = 3/40;
        float r27495 = r27494 * r27488;
        float r27496 = r27495 + r27475;
        float r27497 = 1/6;
        float r27498 = 3;
        float r27499 = pow(r27475, r27498);
        float r27500 = r27497 * r27499;
        float r27501 = r27496 - r27500;
        float r27502 = hypot(r27478, r27475);
        float r27503 = sqrt(r27502);
        float r27504 = fma(r27503, r27503, r27475);
        float r27505 = log(r27504);
        float r27506 = r27493 ? r27501 : r27505;
        float r27507 = r27477 ? r27491 : r27506;
        return r27507;
}

double f_od(double x) {
        double r27508 = x;
        double r27509 = -1.0916759748291127;
        bool r27510 = r27508 <= r27509;
        double r27511 = 1;
        double r27512 = r27511 / r27508;
        double r27513 = 1/8;
        double r27514 = r27513 / r27508;
        double r27515 = 1/2;
        double r27516 = -r27515;
        double r27517 = fma(r27512, r27514, r27516);
        double r27518 = 1/16;
        double r27519 = -r27518;
        double r27520 = 5;
        double r27521 = pow(r27508, r27520);
        double r27522 = r27519 / r27521;
        double r27523 = fma(r27512, r27517, r27522);
        double r27524 = log(r27523);
        double r27525 = 0.007483018826987035;
        bool r27526 = r27508 <= r27525;
        double r27527 = 3/40;
        double r27528 = r27527 * r27521;
        double r27529 = r27528 + r27508;
        double r27530 = 1/6;
        double r27531 = 3;
        double r27532 = pow(r27508, r27531);
        double r27533 = r27530 * r27532;
        double r27534 = r27529 - r27533;
        double r27535 = hypot(r27511, r27508);
        double r27536 = sqrt(r27535);
        double r27537 = fma(r27536, r27536, r27508);
        double r27538 = log(r27537);
        double r27539 = r27526 ? r27534 : r27538;
        double r27540 = r27510 ? r27524 : r27539;
        return r27540;
}

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 r27541, r27542, r27543, r27544, r27545, r27546, r27547;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27541);
        mpfr_init(r27542);
        mpfr_init_set_str(r27543, "1", 10, MPFR_RNDN);
        mpfr_init(r27544);
        mpfr_init(r27545);
        mpfr_init(r27546);
        mpfr_init(r27547);
}

double f_im(double x) {
        mpfr_set_d(r27541, x, MPFR_RNDN);
        mpfr_mul(r27542, r27541, r27541, MPFR_RNDN);
        ;
        mpfr_add(r27544, r27542, r27543, MPFR_RNDN);
        mpfr_sqrt(r27545, r27544, MPFR_RNDN);
        mpfr_add(r27546, r27541, r27545, MPFR_RNDN);
        mpfr_log(r27547, r27546, MPFR_RNDN);
        return mpfr_get_d(r27547, MPFR_RNDN);
}

static mpfr_t r27548, r27549, r27550, r27551, r27552, r27553, r27554, r27555, r27556, r27557, r27558, r27559, r27560, r27561, r27562, r27563, r27564, r27565, r27566, r27567, r27568, r27569, r27570, r27571, r27572, r27573, r27574, r27575, r27576, r27577, r27578, r27579, r27580;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27548);
        mpfr_init_set_str(r27549, "-1.0916759748291127", 10, MPFR_RNDN);
        mpfr_init(r27550);
        mpfr_init_set_str(r27551, "1", 10, MPFR_RNDN);
        mpfr_init(r27552);
        mpfr_init_set_str(r27553, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27554);
        mpfr_init_set_str(r27555, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27556);
        mpfr_init(r27557);
        mpfr_init_set_str(r27558, "1/16", 10, MPFR_RNDN);
        mpfr_init(r27559);
        mpfr_init_set_str(r27560, "5", 10, MPFR_RNDN);
        mpfr_init(r27561);
        mpfr_init(r27562);
        mpfr_init(r27563);
        mpfr_init(r27564);
        mpfr_init_set_str(r27565, "0.007483018826987035", 10, MPFR_RNDN);
        mpfr_init(r27566);
        mpfr_init_set_str(r27567, "3/40", 10, MPFR_RNDN);
        mpfr_init(r27568);
        mpfr_init(r27569);
        mpfr_init_set_str(r27570, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27571, "3", 10, MPFR_RNDN);
        mpfr_init(r27572);
        mpfr_init(r27573);
        mpfr_init(r27574);
        mpfr_init(r27575);
        mpfr_init(r27576);
        mpfr_init(r27577);
        mpfr_init(r27578);
        mpfr_init(r27579);
        mpfr_init(r27580);
}

double f_fm(double x) {
        mpfr_set_d(r27548, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27550, mpfr_cmp(r27548, r27549) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27552, r27551, r27548, MPFR_RNDN);
        ;
        mpfr_div(r27554, r27553, r27548, MPFR_RNDN);
        ;
        mpfr_neg(r27556, r27555, MPFR_RNDN);
        mpfr_fma(r27557, r27552, r27554, r27556, MPFR_RNDN);
        ;
        mpfr_neg(r27559, r27558, MPFR_RNDN);
        ;
        mpfr_pow(r27561, r27548, r27560, MPFR_RNDN);
        mpfr_div(r27562, r27559, r27561, MPFR_RNDN);
        mpfr_fma(r27563, r27552, r27557, r27562, MPFR_RNDN);
        mpfr_log(r27564, r27563, MPFR_RNDN);
        ;
        mpfr_set_si(r27566, mpfr_cmp(r27548, r27565) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r27568, r27567, r27561, MPFR_RNDN);
        mpfr_add(r27569, r27568, r27548, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27572, r27548, r27571, MPFR_RNDN);
        mpfr_mul(r27573, r27570, r27572, MPFR_RNDN);
        mpfr_sub(r27574, r27569, r27573, MPFR_RNDN);
        mpfr_hypot(r27575, r27551, r27548, MPFR_RNDN);
        mpfr_sqrt(r27576, r27575, MPFR_RNDN);
        mpfr_fma(r27577, r27576, r27576, r27548, MPFR_RNDN);
        mpfr_log(r27578, r27577, MPFR_RNDN);
        if (mpfr_get_si(r27566, MPFR_RNDN)) { mpfr_set(r27579, r27574, MPFR_RNDN); } else { mpfr_set(r27579, r27578, MPFR_RNDN); };
        if (mpfr_get_si(r27550, MPFR_RNDN)) { mpfr_set(r27580, r27564, MPFR_RNDN); } else { mpfr_set(r27580, r27579, MPFR_RNDN); };
        return mpfr_get_d(r27580, MPFR_RNDN);
}

static mpfr_t r27581, r27582, r27583, r27584, r27585, r27586, r27587, r27588, r27589, r27590, r27591, r27592, r27593, r27594, r27595, r27596, r27597, r27598, r27599, r27600, r27601, r27602, r27603, r27604, r27605, r27606, r27607, r27608, r27609, r27610, r27611, r27612, r27613;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27581);
        mpfr_init_set_str(r27582, "-1.0916759748291127", 10, MPFR_RNDN);
        mpfr_init(r27583);
        mpfr_init_set_str(r27584, "1", 10, MPFR_RNDN);
        mpfr_init(r27585);
        mpfr_init_set_str(r27586, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27587);
        mpfr_init_set_str(r27588, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27589);
        mpfr_init(r27590);
        mpfr_init_set_str(r27591, "1/16", 10, MPFR_RNDN);
        mpfr_init(r27592);
        mpfr_init_set_str(r27593, "5", 10, MPFR_RNDN);
        mpfr_init(r27594);
        mpfr_init(r27595);
        mpfr_init(r27596);
        mpfr_init(r27597);
        mpfr_init_set_str(r27598, "0.007483018826987035", 10, MPFR_RNDN);
        mpfr_init(r27599);
        mpfr_init_set_str(r27600, "3/40", 10, MPFR_RNDN);
        mpfr_init(r27601);
        mpfr_init(r27602);
        mpfr_init_set_str(r27603, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27604, "3", 10, MPFR_RNDN);
        mpfr_init(r27605);
        mpfr_init(r27606);
        mpfr_init(r27607);
        mpfr_init(r27608);
        mpfr_init(r27609);
        mpfr_init(r27610);
        mpfr_init(r27611);
        mpfr_init(r27612);
        mpfr_init(r27613);
}

double f_dm(double x) {
        mpfr_set_d(r27581, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27583, mpfr_cmp(r27581, r27582) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27585, r27584, r27581, MPFR_RNDN);
        ;
        mpfr_div(r27587, r27586, r27581, MPFR_RNDN);
        ;
        mpfr_neg(r27589, r27588, MPFR_RNDN);
        mpfr_fma(r27590, r27585, r27587, r27589, MPFR_RNDN);
        ;
        mpfr_neg(r27592, r27591, MPFR_RNDN);
        ;
        mpfr_pow(r27594, r27581, r27593, MPFR_RNDN);
        mpfr_div(r27595, r27592, r27594, MPFR_RNDN);
        mpfr_fma(r27596, r27585, r27590, r27595, MPFR_RNDN);
        mpfr_log(r27597, r27596, MPFR_RNDN);
        ;
        mpfr_set_si(r27599, mpfr_cmp(r27581, r27598) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r27601, r27600, r27594, MPFR_RNDN);
        mpfr_add(r27602, r27601, r27581, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27605, r27581, r27604, MPFR_RNDN);
        mpfr_mul(r27606, r27603, r27605, MPFR_RNDN);
        mpfr_sub(r27607, r27602, r27606, MPFR_RNDN);
        mpfr_hypot(r27608, r27584, r27581, MPFR_RNDN);
        mpfr_sqrt(r27609, r27608, MPFR_RNDN);
        mpfr_fma(r27610, r27609, r27609, r27581, MPFR_RNDN);
        mpfr_log(r27611, r27610, MPFR_RNDN);
        if (mpfr_get_si(r27599, MPFR_RNDN)) { mpfr_set(r27612, r27607, MPFR_RNDN); } else { mpfr_set(r27612, r27611, MPFR_RNDN); };
        if (mpfr_get_si(r27583, MPFR_RNDN)) { mpfr_set(r27613, r27597, MPFR_RNDN); } else { mpfr_set(r27613, r27612, MPFR_RNDN); };
        return mpfr_get_d(r27613, MPFR_RNDN);
}