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

char *name = "Octave 3.8, jcobi/4, as called";

double f_if(float i) {
        float r9753 = i;
        float r9754 = r9753 * r9753;
        float r9755 = r9754 * r9754;
        float r9756 = 2.0f;
        float r9757 = r9756 * r9753;
        float r9758 = r9757 * r9757;
        float r9759 = r9755 / r9758;
        float r9760 = 1.0f;
        float r9761 = r9758 - r9760;
        float r9762 = r9759 / r9761;
        return r9762;
}

double f_id(double i) {
        double r9763 = i;
        double r9764 = r9763 * r9763;
        double r9765 = r9764 * r9764;
        double r9766 = 2.0;
        double r9767 = r9766 * r9763;
        double r9768 = r9767 * r9767;
        double r9769 = r9765 / r9768;
        double r9770 = 1.0;
        double r9771 = r9768 - r9770;
        double r9772 = r9769 / r9771;
        return r9772;
}


double f_of(float i) {
        float r9773 = i;
        float r9774 = 221.95421906676754f;
        bool r9775 = r9773 <= r9774;
        float r9776 = 2.0f;
        float r9777 = r9773 / r9776;
        float r9778 = r9777 * r9777;
        float r9779 = r9776 * r9773;
        float r9780 = r9779 * r9779;
        float r9781 = 1.0f;
        float r9782 = r9780 - r9781;
        float r9783 = r9778 / r9782;
        float r9784 = 0.00390625f;
        float r9785 = 4.0f;
        float r9786 = pow(r9773, r9785);
        float r9787 = r9784 / r9786;
        float r9788 = 0.0625f;
        float r9789 = r9787 + r9788;
        float r9790 = 0.015625f;
        float r9791 = r9790 / r9773;
        float r9792 = r9791 / r9773;
        float r9793 = r9789 + r9792;
        float r9794 = r9775 ? r9783 : r9793;
        return r9794;
}

double f_od(double i) {
        double r9795 = i;
        double r9796 = 221.95421906676754;
        bool r9797 = r9795 <= r9796;
        double r9798 = 2.0;
        double r9799 = r9795 / r9798;
        double r9800 = r9799 * r9799;
        double r9801 = r9798 * r9795;
        double r9802 = r9801 * r9801;
        double r9803 = 1.0;
        double r9804 = r9802 - r9803;
        double r9805 = r9800 / r9804;
        double r9806 = 0.00390625;
        double r9807 = 4.0;
        double r9808 = pow(r9795, r9807);
        double r9809 = r9806 / r9808;
        double r9810 = 0.0625;
        double r9811 = r9809 + r9810;
        double r9812 = 0.015625;
        double r9813 = r9812 / r9795;
        double r9814 = r9813 / r9795;
        double r9815 = r9811 + r9814;
        double r9816 = r9797 ? r9805 : r9815;
        return r9816;
}

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 r9817, r9818, r9819, r9820, r9821, r9822, r9823, r9824, r9825, r9826;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9817);
        mpfr_init(r9818);
        mpfr_init(r9819);
        mpfr_init_set_str(r9820, "2", 10, MPFR_RNDN);
        mpfr_init(r9821);
        mpfr_init(r9822);
        mpfr_init(r9823);
        mpfr_init_set_str(r9824, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9825);
        mpfr_init(r9826);
}

double f_im(double i) {
        mpfr_set_d(r9817, i, MPFR_RNDN);
        mpfr_mul(r9818, r9817, r9817, MPFR_RNDN);
        mpfr_mul(r9819, r9818, r9818, MPFR_RNDN);
        ;
        mpfr_mul(r9821, r9820, r9817, MPFR_RNDN);
        mpfr_mul(r9822, r9821, r9821, MPFR_RNDN);
        mpfr_div(r9823, r9819, r9822, MPFR_RNDN);
        ;
        mpfr_sub(r9825, r9822, r9824, MPFR_RNDN);
        mpfr_div(r9826, r9823, r9825, MPFR_RNDN);
        return mpfr_get_d(r9826, MPFR_RNDN);
}

static mpfr_t r9827, r9828, r9829, r9830, r9831, r9832, r9833, r9834, r9835, r9836, r9837, r9838, r9839, r9840, r9841, r9842, r9843, r9844, r9845, r9846, r9847, r9848;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9827);
        mpfr_init_set_str(r9828, "221.95421906676754", 10, MPFR_RNDN);
        mpfr_init(r9829);
        mpfr_init_set_str(r9830, "2", 10, MPFR_RNDN);
        mpfr_init(r9831);
        mpfr_init(r9832);
        mpfr_init(r9833);
        mpfr_init(r9834);
        mpfr_init_set_str(r9835, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9836);
        mpfr_init(r9837);
        mpfr_init_set_str(r9838, "0.00390625", 10, MPFR_RNDN);
        mpfr_init_set_str(r9839, "4", 10, MPFR_RNDN);
        mpfr_init(r9840);
        mpfr_init(r9841);
        mpfr_init_set_str(r9842, "1/16", 10, MPFR_RNDN);
        mpfr_init(r9843);
        mpfr_init_set_str(r9844, "0.015625", 10, MPFR_RNDN);
        mpfr_init(r9845);
        mpfr_init(r9846);
        mpfr_init(r9847);
        mpfr_init(r9848);
}

double f_fm(double i) {
        mpfr_set_d(r9827, i, MPFR_RNDN);
        ;
        mpfr_set_si(r9829, mpfr_cmp(r9827, r9828) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9831, r9827, r9830, MPFR_RNDN);
        mpfr_mul(r9832, r9831, r9831, MPFR_RNDN);
        mpfr_mul(r9833, r9830, r9827, MPFR_RNDN);
        mpfr_mul(r9834, r9833, r9833, MPFR_RNDN);
        ;
        mpfr_sub(r9836, r9834, r9835, MPFR_RNDN);
        mpfr_div(r9837, r9832, r9836, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9840, r9827, r9839, MPFR_RNDN);
        mpfr_div(r9841, r9838, r9840, MPFR_RNDN);
        ;
        mpfr_add(r9843, r9841, r9842, MPFR_RNDN);
        ;
        mpfr_div(r9845, r9844, r9827, MPFR_RNDN);
        mpfr_div(r9846, r9845, r9827, MPFR_RNDN);
        mpfr_add(r9847, r9843, r9846, MPFR_RNDN);
        if (mpfr_get_si(r9829, MPFR_RNDN)) { mpfr_set(r9848, r9837, MPFR_RNDN); } else { mpfr_set(r9848, r9847, MPFR_RNDN); };
        return mpfr_get_d(r9848, MPFR_RNDN);
}

static mpfr_t r9849, r9850, r9851, r9852, r9853, r9854, r9855, r9856, r9857, r9858, r9859, r9860, r9861, r9862, r9863, r9864, r9865, r9866, r9867, r9868, r9869, r9870;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init(r9849);
        mpfr_init_set_str(r9850, "221.95421906676754", 10, MPFR_RNDN);
        mpfr_init(r9851);
        mpfr_init_set_str(r9852, "2", 10, MPFR_RNDN);
        mpfr_init(r9853);
        mpfr_init(r9854);
        mpfr_init(r9855);
        mpfr_init(r9856);
        mpfr_init_set_str(r9857, "1.0", 10, MPFR_RNDN);
        mpfr_init(r9858);
        mpfr_init(r9859);
        mpfr_init_set_str(r9860, "0.00390625", 10, MPFR_RNDN);
        mpfr_init_set_str(r9861, "4", 10, MPFR_RNDN);
        mpfr_init(r9862);
        mpfr_init(r9863);
        mpfr_init_set_str(r9864, "1/16", 10, MPFR_RNDN);
        mpfr_init(r9865);
        mpfr_init_set_str(r9866, "0.015625", 10, MPFR_RNDN);
        mpfr_init(r9867);
        mpfr_init(r9868);
        mpfr_init(r9869);
        mpfr_init(r9870);
}

double f_dm(double i) {
        mpfr_set_d(r9849, i, MPFR_RNDN);
        ;
        mpfr_set_si(r9851, mpfr_cmp(r9849, r9850) <= 0, MPFR_RNDN);
        ;
        mpfr_div(r9853, r9849, r9852, MPFR_RNDN);
        mpfr_mul(r9854, r9853, r9853, MPFR_RNDN);
        mpfr_mul(r9855, r9852, r9849, MPFR_RNDN);
        mpfr_mul(r9856, r9855, r9855, MPFR_RNDN);
        ;
        mpfr_sub(r9858, r9856, r9857, MPFR_RNDN);
        mpfr_div(r9859, r9854, r9858, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9862, r9849, r9861, MPFR_RNDN);
        mpfr_div(r9863, r9860, r9862, MPFR_RNDN);
        ;
        mpfr_add(r9865, r9863, r9864, MPFR_RNDN);
        ;
        mpfr_div(r9867, r9866, r9849, MPFR_RNDN);
        mpfr_div(r9868, r9867, r9849, MPFR_RNDN);
        mpfr_add(r9869, r9865, r9868, MPFR_RNDN);
        if (mpfr_get_si(r9851, MPFR_RNDN)) { mpfr_set(r9870, r9859, MPFR_RNDN); } else { mpfr_set(r9870, r9869, MPFR_RNDN); };
        return mpfr_get_d(r9870, MPFR_RNDN);
}