#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 r27563 = x;
        float r27564 = r27563 * r27563;
        float r27565 = 1;
        float r27566 = r27564 + r27565;
        float r27567 = sqrt(r27566);
        float r27568 = r27563 + r27567;
        float r27569 = log(r27568);
        return r27569;
}

double f_id(double x) {
        double r27570 = x;
        double r27571 = r27570 * r27570;
        double r27572 = 1;
        double r27573 = r27571 + r27572;
        double r27574 = sqrt(r27573);
        double r27575 = r27570 + r27574;
        double r27576 = log(r27575);
        return r27576;
}


double f_of(float x) {
        float r27577 = x;
        float r27578 = -1.0578922685264147;
        bool r27579 = r27577 <= r27578;
        float r27580 = 1/8;
        float r27581 = r27580 / r27577;
        float r27582 = r27577 * r27577;
        float r27583 = r27581 / r27582;
        float r27584 = 1/2;
        float r27585 = r27584 / r27577;
        float r27586 = r27583 - r27585;
        float r27587 = 1/16;
        float r27588 = 5;
        float r27589 = pow(r27577, r27588);
        float r27590 = r27587 / r27589;
        float r27591 = r27586 - r27590;
        float r27592 = log(r27591);
        float r27593 = 0.9635279732307327;
        bool r27594 = r27577 <= r27593;
        float r27595 = 3/40;
        float r27596 = r27595 * r27589;
        float r27597 = r27596 + r27577;
        float r27598 = 1/6;
        float r27599 = 3;
        float r27600 = pow(r27577, r27599);
        float r27601 = r27598 * r27600;
        float r27602 = r27597 - r27601;
        float r27603 = 2;
        float r27604 = r27577 * r27603;
        float r27605 = r27585 + r27604;
        float r27606 = r27605 - r27583;
        float r27607 = log(r27606);
        float r27608 = r27594 ? r27602 : r27607;
        float r27609 = r27579 ? r27592 : r27608;
        return r27609;
}

double f_od(double x) {
        double r27610 = x;
        double r27611 = -1.0578922685264147;
        bool r27612 = r27610 <= r27611;
        double r27613 = 1/8;
        double r27614 = r27613 / r27610;
        double r27615 = r27610 * r27610;
        double r27616 = r27614 / r27615;
        double r27617 = 1/2;
        double r27618 = r27617 / r27610;
        double r27619 = r27616 - r27618;
        double r27620 = 1/16;
        double r27621 = 5;
        double r27622 = pow(r27610, r27621);
        double r27623 = r27620 / r27622;
        double r27624 = r27619 - r27623;
        double r27625 = log(r27624);
        double r27626 = 0.9635279732307327;
        bool r27627 = r27610 <= r27626;
        double r27628 = 3/40;
        double r27629 = r27628 * r27622;
        double r27630 = r27629 + r27610;
        double r27631 = 1/6;
        double r27632 = 3;
        double r27633 = pow(r27610, r27632);
        double r27634 = r27631 * r27633;
        double r27635 = r27630 - r27634;
        double r27636 = 2;
        double r27637 = r27610 * r27636;
        double r27638 = r27618 + r27637;
        double r27639 = r27638 - r27616;
        double r27640 = log(r27639);
        double r27641 = r27627 ? r27635 : r27640;
        double r27642 = r27612 ? r27625 : r27641;
        return r27642;
}

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 r27643, r27644, r27645, r27646, r27647, r27648, r27649;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27643);
        mpfr_init(r27644);
        mpfr_init_set_str(r27645, "1", 10, MPFR_RNDN);
        mpfr_init(r27646);
        mpfr_init(r27647);
        mpfr_init(r27648);
        mpfr_init(r27649);
}

double f_im(double x) {
        mpfr_set_d(r27643, x, MPFR_RNDN);
        mpfr_mul(r27644, r27643, r27643, MPFR_RNDN);
        ;
        mpfr_add(r27646, r27644, r27645, MPFR_RNDN);
        mpfr_sqrt(r27647, r27646, MPFR_RNDN);
        mpfr_add(r27648, r27643, r27647, MPFR_RNDN);
        mpfr_log(r27649, r27648, MPFR_RNDN);
        return mpfr_get_d(r27649, MPFR_RNDN);
}

static mpfr_t r27650, r27651, r27652, r27653, r27654, r27655, r27656, r27657, r27658, r27659, r27660, r27661, r27662, r27663, r27664, r27665, r27666, r27667, r27668, r27669, r27670, r27671, r27672, r27673, r27674, r27675, r27676, r27677, r27678, r27679, r27680, r27681, r27682;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27650);
        mpfr_init_set_str(r27651, "-1.0578922685264147", 10, MPFR_RNDN);
        mpfr_init(r27652);
        mpfr_init_set_str(r27653, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27654);
        mpfr_init(r27655);
        mpfr_init(r27656);
        mpfr_init_set_str(r27657, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27658);
        mpfr_init(r27659);
        mpfr_init_set_str(r27660, "1/16", 10, MPFR_RNDN);
        mpfr_init_set_str(r27661, "5", 10, MPFR_RNDN);
        mpfr_init(r27662);
        mpfr_init(r27663);
        mpfr_init(r27664);
        mpfr_init(r27665);
        mpfr_init_set_str(r27666, "0.9635279732307327", 10, MPFR_RNDN);
        mpfr_init(r27667);
        mpfr_init_set_str(r27668, "3/40", 10, MPFR_RNDN);
        mpfr_init(r27669);
        mpfr_init(r27670);
        mpfr_init_set_str(r27671, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27672, "3", 10, MPFR_RNDN);
        mpfr_init(r27673);
        mpfr_init(r27674);
        mpfr_init(r27675);
        mpfr_init_set_str(r27676, "2", 10, MPFR_RNDN);
        mpfr_init(r27677);
        mpfr_init(r27678);
        mpfr_init(r27679);
        mpfr_init(r27680);
        mpfr_init(r27681);
        mpfr_init(r27682);
}

double f_fm(double x) {
        mpfr_set_d(r27650, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27652, mpfr_cmp(r27650, r27651) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27654, r27653, r27650, MPFR_RNDN);
        mpfr_mul(r27655, r27650, r27650, MPFR_RNDN);
        mpfr_div(r27656, r27654, r27655, MPFR_RNDN);
        ;
        mpfr_div(r27658, r27657, r27650, MPFR_RNDN);
        mpfr_sub(r27659, r27656, r27658, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27662, r27650, r27661, MPFR_RNDN);
        mpfr_div(r27663, r27660, r27662, MPFR_RNDN);
        mpfr_sub(r27664, r27659, r27663, MPFR_RNDN);
        mpfr_log(r27665, r27664, MPFR_RNDN);
        ;
        mpfr_set_si(r27667, mpfr_cmp(r27650, r27666) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r27669, r27668, r27662, MPFR_RNDN);
        mpfr_add(r27670, r27669, r27650, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27673, r27650, r27672, MPFR_RNDN);
        mpfr_mul(r27674, r27671, r27673, MPFR_RNDN);
        mpfr_sub(r27675, r27670, r27674, MPFR_RNDN);
        ;
        mpfr_mul(r27677, r27650, r27676, MPFR_RNDN);
        mpfr_add(r27678, r27658, r27677, MPFR_RNDN);
        mpfr_sub(r27679, r27678, r27656, MPFR_RNDN);
        mpfr_log(r27680, r27679, MPFR_RNDN);
        if (mpfr_get_si(r27667, MPFR_RNDN)) { mpfr_set(r27681, r27675, MPFR_RNDN); } else { mpfr_set(r27681, r27680, MPFR_RNDN); };
        if (mpfr_get_si(r27652, MPFR_RNDN)) { mpfr_set(r27682, r27665, MPFR_RNDN); } else { mpfr_set(r27682, r27681, MPFR_RNDN); };
        return mpfr_get_d(r27682, MPFR_RNDN);
}

static mpfr_t r27683, r27684, r27685, r27686, r27687, r27688, r27689, r27690, r27691, r27692, r27693, r27694, r27695, r27696, r27697, r27698, r27699, r27700, r27701, r27702, r27703, r27704, r27705, r27706, r27707, r27708, r27709, r27710, r27711, r27712, r27713, r27714, r27715;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init(r27683);
        mpfr_init_set_str(r27684, "-1.0578922685264147", 10, MPFR_RNDN);
        mpfr_init(r27685);
        mpfr_init_set_str(r27686, "1/8", 10, MPFR_RNDN);
        mpfr_init(r27687);
        mpfr_init(r27688);
        mpfr_init(r27689);
        mpfr_init_set_str(r27690, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27691);
        mpfr_init(r27692);
        mpfr_init_set_str(r27693, "1/16", 10, MPFR_RNDN);
        mpfr_init_set_str(r27694, "5", 10, MPFR_RNDN);
        mpfr_init(r27695);
        mpfr_init(r27696);
        mpfr_init(r27697);
        mpfr_init(r27698);
        mpfr_init_set_str(r27699, "0.9635279732307327", 10, MPFR_RNDN);
        mpfr_init(r27700);
        mpfr_init_set_str(r27701, "3/40", 10, MPFR_RNDN);
        mpfr_init(r27702);
        mpfr_init(r27703);
        mpfr_init_set_str(r27704, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r27705, "3", 10, MPFR_RNDN);
        mpfr_init(r27706);
        mpfr_init(r27707);
        mpfr_init(r27708);
        mpfr_init_set_str(r27709, "2", 10, MPFR_RNDN);
        mpfr_init(r27710);
        mpfr_init(r27711);
        mpfr_init(r27712);
        mpfr_init(r27713);
        mpfr_init(r27714);
        mpfr_init(r27715);
}

double f_dm(double x) {
        mpfr_set_d(r27683, x, MPFR_RNDN);
        ;
        mpfr_set_si(r27685, mpfr_cmp(r27683, r27684) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r27687, r27686, r27683, MPFR_RNDN);
        mpfr_mul(r27688, r27683, r27683, MPFR_RNDN);
        mpfr_div(r27689, r27687, r27688, MPFR_RNDN);
        ;
        mpfr_div(r27691, r27690, r27683, MPFR_RNDN);
        mpfr_sub(r27692, r27689, r27691, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27695, r27683, r27694, MPFR_RNDN);
        mpfr_div(r27696, r27693, r27695, MPFR_RNDN);
        mpfr_sub(r27697, r27692, r27696, MPFR_RNDN);
        mpfr_log(r27698, r27697, MPFR_RNDN);
        ;
        mpfr_set_si(r27700, mpfr_cmp(r27683, r27699) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r27702, r27701, r27695, MPFR_RNDN);
        mpfr_add(r27703, r27702, r27683, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27706, r27683, r27705, MPFR_RNDN);
        mpfr_mul(r27707, r27704, r27706, MPFR_RNDN);
        mpfr_sub(r27708, r27703, r27707, MPFR_RNDN);
        ;
        mpfr_mul(r27710, r27683, r27709, MPFR_RNDN);
        mpfr_add(r27711, r27691, r27710, MPFR_RNDN);
        mpfr_sub(r27712, r27711, r27689, MPFR_RNDN);
        mpfr_log(r27713, r27712, MPFR_RNDN);
        if (mpfr_get_si(r27700, MPFR_RNDN)) { mpfr_set(r27714, r27708, MPFR_RNDN); } else { mpfr_set(r27714, r27713, MPFR_RNDN); };
        if (mpfr_get_si(r27685, MPFR_RNDN)) { mpfr_set(r27715, r27698, MPFR_RNDN); } else { mpfr_set(r27715, r27714, MPFR_RNDN); };
        return mpfr_get_d(r27715, MPFR_RNDN);
}