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

char *name = "Asymptote C";

double f_if(float x) {
        float r26675 = x;
        float r26676 = 1;
        float r26677 = r26675 + r26676;
        float r26678 = r26675 / r26677;
        float r26679 = r26675 - r26676;
        float r26680 = r26677 / r26679;
        float r26681 = r26678 - r26680;
        return r26681;
}

double f_id(double x) {
        double r26682 = x;
        double r26683 = 1;
        double r26684 = r26682 + r26683;
        double r26685 = r26682 / r26684;
        double r26686 = r26682 - r26683;
        double r26687 = r26684 / r26686;
        double r26688 = r26685 - r26687;
        return r26688;
}


double f_of(float x) {
        float r26689 = 1;
        float r26690 = 3;
        float r26691 = x;
        float r26692 = r26690 / r26691;
        float r26693 = r26689 + r26692;
        float r26694 = -r26689;
        float r26695 = r26691 * r26691;
        float r26696 = r26694 / r26695;
        float r26697 = -r26690;
        float r26698 = r26697 / r26691;
        float r26699 = fma(r26693, r26696, r26698);
        float r26700 = pow(r26699, r26690);
        float r26701 = cbrt(r26700);
        float r26702 = -0.00010071291504276669;
        bool r26703 = r26701 <= r26702;
        float r26704 = r26689 + r26691;
        float r26705 = r26689 / r26704;
        float r26706 = r26691 - r26689;
        float r26707 = r26704 / r26706;
        float r26708 = -r26707;
        float r26709 = fma(r26691, r26705, r26708);
        float r26710 = 2.8259313836317426e-05;
        bool r26711 = r26701 <= r26710;
        float r26712 = r26694 / r26706;
        float r26713 = r26691 / r26704;
        float r26714 = fma(r26712, r26704, r26713);
        float r26715 = r26711 ? r26699 : r26714;
        float r26716 = r26703 ? r26709 : r26715;
        return r26716;
}

double f_od(double x) {
        double r26717 = 1;
        double r26718 = 3;
        double r26719 = x;
        double r26720 = r26718 / r26719;
        double r26721 = r26717 + r26720;
        double r26722 = -r26717;
        double r26723 = r26719 * r26719;
        double r26724 = r26722 / r26723;
        double r26725 = -r26718;
        double r26726 = r26725 / r26719;
        double r26727 = fma(r26721, r26724, r26726);
        double r26728 = pow(r26727, r26718);
        double r26729 = cbrt(r26728);
        double r26730 = -0.00010071291504276669;
        bool r26731 = r26729 <= r26730;
        double r26732 = r26717 + r26719;
        double r26733 = r26717 / r26732;
        double r26734 = r26719 - r26717;
        double r26735 = r26732 / r26734;
        double r26736 = -r26735;
        double r26737 = fma(r26719, r26733, r26736);
        double r26738 = 2.8259313836317426e-05;
        bool r26739 = r26729 <= r26738;
        double r26740 = r26722 / r26734;
        double r26741 = r26719 / r26732;
        double r26742 = fma(r26740, r26732, r26741);
        double r26743 = r26739 ? r26727 : r26742;
        double r26744 = r26731 ? r26737 : r26743;
        return r26744;
}

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 r26745, r26746, r26747, r26748, r26749, r26750, r26751;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r26745);
        mpfr_init_set_str(r26746, "1", 10, MPFR_RNDN);
        mpfr_init(r26747);
        mpfr_init(r26748);
        mpfr_init(r26749);
        mpfr_init(r26750);
        mpfr_init(r26751);
}

double f_im(double x) {
        mpfr_set_d(r26745, x, MPFR_RNDN);
        ;
        mpfr_add(r26747, r26745, r26746, MPFR_RNDN);
        mpfr_div(r26748, r26745, r26747, MPFR_RNDN);
        mpfr_sub(r26749, r26745, r26746, MPFR_RNDN);
        mpfr_div(r26750, r26747, r26749, MPFR_RNDN);
        mpfr_sub(r26751, r26748, r26750, MPFR_RNDN);
        return mpfr_get_d(r26751, MPFR_RNDN);
}

static mpfr_t r26752, r26753, r26754, r26755, r26756, r26757, r26758, r26759, r26760, r26761, r26762, r26763, r26764, r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778, r26779;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26752, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26753, "3", 10, MPFR_RNDN);
        mpfr_init(r26754);
        mpfr_init(r26755);
        mpfr_init(r26756);
        mpfr_init(r26757);
        mpfr_init(r26758);
        mpfr_init(r26759);
        mpfr_init(r26760);
        mpfr_init(r26761);
        mpfr_init(r26762);
        mpfr_init(r26763);
        mpfr_init(r26764);
        mpfr_init_set_str(r26765, "-0.00010071291504276669", 10, MPFR_RNDN);
        mpfr_init(r26766);
        mpfr_init(r26767);
        mpfr_init(r26768);
        mpfr_init(r26769);
        mpfr_init(r26770);
        mpfr_init(r26771);
        mpfr_init(r26772);
        mpfr_init_set_str(r26773, "2.8259313836317426e-05", 10, MPFR_RNDN);
        mpfr_init(r26774);
        mpfr_init(r26775);
        mpfr_init(r26776);
        mpfr_init(r26777);
        mpfr_init(r26778);
        mpfr_init(r26779);
}

double f_fm(double x) {
        ;
        ;
        mpfr_set_d(r26754, x, MPFR_RNDN);
        mpfr_div(r26755, r26753, r26754, MPFR_RNDN);
        mpfr_add(r26756, r26752, r26755, MPFR_RNDN);
        mpfr_neg(r26757, r26752, MPFR_RNDN);
        mpfr_mul(r26758, r26754, r26754, MPFR_RNDN);
        mpfr_div(r26759, r26757, r26758, MPFR_RNDN);
        mpfr_neg(r26760, r26753, MPFR_RNDN);
        mpfr_div(r26761, r26760, r26754, MPFR_RNDN);
        mpfr_fma(r26762, r26756, r26759, r26761, MPFR_RNDN);
        mpfr_pow(r26763, r26762, r26753, MPFR_RNDN);
        mpfr_cbrt(r26764, r26763, MPFR_RNDN);
        ;
        mpfr_set_si(r26766, mpfr_cmp(r26764, r26765) <= 0, MPFR_RNDN);
        mpfr_add(r26767, r26752, r26754, MPFR_RNDN);
        mpfr_div(r26768, r26752, r26767, MPFR_RNDN);
        mpfr_sub(r26769, r26754, r26752, MPFR_RNDN);
        mpfr_div(r26770, r26767, r26769, MPFR_RNDN);
        mpfr_neg(r26771, r26770, MPFR_RNDN);
        mpfr_fma(r26772, r26754, r26768, r26771, MPFR_RNDN);
        ;
        mpfr_set_si(r26774, mpfr_cmp(r26764, r26773) <= 0, MPFR_RNDN);
        mpfr_div(r26775, r26757, r26769, MPFR_RNDN);
        mpfr_div(r26776, r26754, r26767, MPFR_RNDN);
        mpfr_fma(r26777, r26775, r26767, r26776, MPFR_RNDN);
        if (mpfr_get_si(r26774, MPFR_RNDN)) { mpfr_set(r26778, r26762, MPFR_RNDN); } else { mpfr_set(r26778, r26777, MPFR_RNDN); };
        if (mpfr_get_si(r26766, MPFR_RNDN)) { mpfr_set(r26779, r26772, MPFR_RNDN); } else { mpfr_set(r26779, r26778, MPFR_RNDN); };
        return mpfr_get_d(r26779, MPFR_RNDN);
}

static mpfr_t r26780, r26781, r26782, r26783, r26784, r26785, r26786, r26787, r26788, r26789, r26790, r26791, r26792, r26793, r26794, r26795, r26796, r26797, r26798, r26799, r26800, r26801, r26802, r26803, r26804, r26805, r26806, r26807;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r26780, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r26781, "3", 10, MPFR_RNDN);
        mpfr_init(r26782);
        mpfr_init(r26783);
        mpfr_init(r26784);
        mpfr_init(r26785);
        mpfr_init(r26786);
        mpfr_init(r26787);
        mpfr_init(r26788);
        mpfr_init(r26789);
        mpfr_init(r26790);
        mpfr_init(r26791);
        mpfr_init(r26792);
        mpfr_init_set_str(r26793, "-0.00010071291504276669", 10, MPFR_RNDN);
        mpfr_init(r26794);
        mpfr_init(r26795);
        mpfr_init(r26796);
        mpfr_init(r26797);
        mpfr_init(r26798);
        mpfr_init(r26799);
        mpfr_init(r26800);
        mpfr_init_set_str(r26801, "2.8259313836317426e-05", 10, MPFR_RNDN);
        mpfr_init(r26802);
        mpfr_init(r26803);
        mpfr_init(r26804);
        mpfr_init(r26805);
        mpfr_init(r26806);
        mpfr_init(r26807);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r26782, x, MPFR_RNDN);
        mpfr_div(r26783, r26781, r26782, MPFR_RNDN);
        mpfr_add(r26784, r26780, r26783, MPFR_RNDN);
        mpfr_neg(r26785, r26780, MPFR_RNDN);
        mpfr_mul(r26786, r26782, r26782, MPFR_RNDN);
        mpfr_div(r26787, r26785, r26786, MPFR_RNDN);
        mpfr_neg(r26788, r26781, MPFR_RNDN);
        mpfr_div(r26789, r26788, r26782, MPFR_RNDN);
        mpfr_fma(r26790, r26784, r26787, r26789, MPFR_RNDN);
        mpfr_pow(r26791, r26790, r26781, MPFR_RNDN);
        mpfr_cbrt(r26792, r26791, MPFR_RNDN);
        ;
        mpfr_set_si(r26794, mpfr_cmp(r26792, r26793) <= 0, MPFR_RNDN);
        mpfr_add(r26795, r26780, r26782, MPFR_RNDN);
        mpfr_div(r26796, r26780, r26795, MPFR_RNDN);
        mpfr_sub(r26797, r26782, r26780, MPFR_RNDN);
        mpfr_div(r26798, r26795, r26797, MPFR_RNDN);
        mpfr_neg(r26799, r26798, MPFR_RNDN);
        mpfr_fma(r26800, r26782, r26796, r26799, MPFR_RNDN);
        ;
        mpfr_set_si(r26802, mpfr_cmp(r26792, r26801) <= 0, MPFR_RNDN);
        mpfr_div(r26803, r26785, r26797, MPFR_RNDN);
        mpfr_div(r26804, r26782, r26795, MPFR_RNDN);
        mpfr_fma(r26805, r26803, r26795, r26804, MPFR_RNDN);
        if (mpfr_get_si(r26802, MPFR_RNDN)) { mpfr_set(r26806, r26790, MPFR_RNDN); } else { mpfr_set(r26806, r26805, MPFR_RNDN); };
        if (mpfr_get_si(r26794, MPFR_RNDN)) { mpfr_set(r26807, r26800, MPFR_RNDN); } else { mpfr_set(r26807, r26806, MPFR_RNDN); };
        return mpfr_get_d(r26807, MPFR_RNDN);
}