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

char *name = "Logistic function from Lakshay Garg";

double f_if(float x, float __attribute__((unused)) y) {
        float r8716 = 2;
        float r8717 = 1;
        float r8718 = -2;
        float r8719 = x;
        float r8720 = r8718 * r8719;
        float r8721 = exp(r8720);
        float r8722 = r8717 + r8721;
        float r8723 = r8716 / r8722;
        float r8724 = r8723 - r8717;
        return r8724;
}

double f_id(double x, double __attribute__((unused)) y) {
        double r8725 = 2;
        double r8726 = 1;
        double r8727 = -2;
        double r8728 = x;
        double r8729 = r8727 * r8728;
        double r8730 = exp(r8729);
        double r8731 = r8726 + r8730;
        double r8732 = r8725 / r8731;
        double r8733 = r8732 - r8726;
        return r8733;
}


double f_of(float x, float __attribute__((unused)) y) {
        float r8734 = 2;
        float r8735 = 1;
        float r8736 = -2;
        float r8737 = x;
        float r8738 = r8736 * r8737;
        float r8739 = exp(r8738);
        float r8740 = r8735 + r8739;
        float r8741 = r8734 / r8740;
        float r8742 = r8741 - r8735;
        float r8743 = -1.6888443106222433e-07;
        bool r8744 = r8742 <= r8743;
        float r8745 = log(r8734);
        float r8746 = log1p(r8739);
        float r8747 = r8745 - r8746;
        float r8748 = exp(r8747);
        float r8749 = 3;
        float r8750 = pow(r8748, r8749);
        float r8751 = pow(r8735, r8749);
        float r8752 = r8750 - r8751;
        float r8753 = exp(r8746);
        float r8754 = r8734 / r8753;
        float r8755 = r8735 + r8754;
        float r8756 = fma(r8754, r8754, r8755);
        float r8757 = r8752 / r8756;
        float r8758 = 0.0032654260743848235;
        bool r8759 = r8742 <= r8758;
        float r8760 = 2/15;
        float r8761 = 5;
        float r8762 = pow(r8737, r8761);
        float r8763 = r8760 * r8762;
        float r8764 = r8737 + r8763;
        float r8765 = 1/3;
        float r8766 = pow(r8737, r8749);
        float r8767 = r8765 * r8766;
        float r8768 = r8764 - r8767;
        float r8769 = expm1(r8747);
        float r8770 = log(r8769);
        float r8771 = exp(r8770);
        float r8772 = r8759 ? r8768 : r8771;
        float r8773 = r8744 ? r8757 : r8772;
        return r8773;
}

double f_od(double x, double __attribute__((unused)) y) {
        double r8774 = 2;
        double r8775 = 1;
        double r8776 = -2;
        double r8777 = x;
        double r8778 = r8776 * r8777;
        double r8779 = exp(r8778);
        double r8780 = r8775 + r8779;
        double r8781 = r8774 / r8780;
        double r8782 = r8781 - r8775;
        double r8783 = -1.6888443106222433e-07;
        bool r8784 = r8782 <= r8783;
        double r8785 = log(r8774);
        double r8786 = log1p(r8779);
        double r8787 = r8785 - r8786;
        double r8788 = exp(r8787);
        double r8789 = 3;
        double r8790 = pow(r8788, r8789);
        double r8791 = pow(r8775, r8789);
        double r8792 = r8790 - r8791;
        double r8793 = exp(r8786);
        double r8794 = r8774 / r8793;
        double r8795 = r8775 + r8794;
        double r8796 = fma(r8794, r8794, r8795);
        double r8797 = r8792 / r8796;
        double r8798 = 0.0032654260743848235;
        bool r8799 = r8782 <= r8798;
        double r8800 = 2/15;
        double r8801 = 5;
        double r8802 = pow(r8777, r8801);
        double r8803 = r8800 * r8802;
        double r8804 = r8777 + r8803;
        double r8805 = 1/3;
        double r8806 = pow(r8777, r8789);
        double r8807 = r8805 * r8806;
        double r8808 = r8804 - r8807;
        double r8809 = expm1(r8787);
        double r8810 = log(r8809);
        double r8811 = exp(r8810);
        double r8812 = r8799 ? r8808 : r8811;
        double r8813 = r8784 ? r8797 : r8812;
        return r8813;
}

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 r8814, r8815, r8816, r8817, r8818, r8819, r8820, r8821, r8822;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r8814, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r8815, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8816, "-2", 10, MPFR_RNDN);
        mpfr_init(r8817);
        mpfr_init(r8818);
        mpfr_init(r8819);
        mpfr_init(r8820);
        mpfr_init(r8821);
        mpfr_init(r8822);
}

double f_im(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r8817, x, MPFR_RNDN);
        mpfr_mul(r8818, r8816, r8817, MPFR_RNDN);
        mpfr_exp(r8819, r8818, MPFR_RNDN);
        mpfr_add(r8820, r8815, r8819, MPFR_RNDN);
        mpfr_div(r8821, r8814, r8820, MPFR_RNDN);
        mpfr_sub(r8822, r8821, r8815, MPFR_RNDN);
        return mpfr_get_d(r8822, MPFR_RNDN);
}

static mpfr_t r8823, r8824, r8825, r8826, r8827, r8828, r8829, r8830, r8831, r8832, r8833, r8834, r8835, r8836, r8837, r8838, r8839, r8840, r8841, r8842, r8843, r8844, r8845, r8846, r8847, r8848, r8849, r8850, r8851, r8852, r8853, r8854, r8855, r8856, r8857, r8858, r8859, r8860, r8861, r8862;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r8823, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r8824, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8825, "-2", 10, MPFR_RNDN);
        mpfr_init(r8826);
        mpfr_init(r8827);
        mpfr_init(r8828);
        mpfr_init(r8829);
        mpfr_init(r8830);
        mpfr_init(r8831);
        mpfr_init_set_str(r8832, "-1.6888443106222433e-07", 10, MPFR_RNDN);
        mpfr_init(r8833);
        mpfr_init(r8834);
        mpfr_init(r8835);
        mpfr_init(r8836);
        mpfr_init(r8837);
        mpfr_init_set_str(r8838, "3", 10, MPFR_RNDN);
        mpfr_init(r8839);
        mpfr_init(r8840);
        mpfr_init(r8841);
        mpfr_init(r8842);
        mpfr_init(r8843);
        mpfr_init(r8844);
        mpfr_init(r8845);
        mpfr_init(r8846);
        mpfr_init_set_str(r8847, "0.0032654260743848235", 10, MPFR_RNDN);
        mpfr_init(r8848);
        mpfr_init_set_str(r8849, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r8850, "5", 10, MPFR_RNDN);
        mpfr_init(r8851);
        mpfr_init(r8852);
        mpfr_init(r8853);
        mpfr_init_set_str(r8854, "1/3", 10, MPFR_RNDN);
        mpfr_init(r8855);
        mpfr_init(r8856);
        mpfr_init(r8857);
        mpfr_init(r8858);
        mpfr_init(r8859);
        mpfr_init(r8860);
        mpfr_init(r8861);
        mpfr_init(r8862);
}

double f_fm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r8826, x, MPFR_RNDN);
        mpfr_mul(r8827, r8825, r8826, MPFR_RNDN);
        mpfr_exp(r8828, r8827, MPFR_RNDN);
        mpfr_add(r8829, r8824, r8828, MPFR_RNDN);
        mpfr_div(r8830, r8823, r8829, MPFR_RNDN);
        mpfr_sub(r8831, r8830, r8824, MPFR_RNDN);
        ;
        mpfr_set_si(r8833, mpfr_cmp(r8831, r8832) <= 0, MPFR_RNDN);
        mpfr_log(r8834, r8823, MPFR_RNDN);
        mpfr_log1p(r8835, r8828, MPFR_RNDN);
        mpfr_sub(r8836, r8834, r8835, MPFR_RNDN);
        mpfr_exp(r8837, r8836, MPFR_RNDN);
        ;
        mpfr_pow(r8839, r8837, r8838, MPFR_RNDN);
        mpfr_pow(r8840, r8824, r8838, MPFR_RNDN);
        mpfr_sub(r8841, r8839, r8840, MPFR_RNDN);
        mpfr_exp(r8842, r8835, MPFR_RNDN);
        mpfr_div(r8843, r8823, r8842, MPFR_RNDN);
        mpfr_add(r8844, r8824, r8843, MPFR_RNDN);
        mpfr_fma(r8845, r8843, r8843, r8844, MPFR_RNDN);
        mpfr_div(r8846, r8841, r8845, MPFR_RNDN);
        ;
        mpfr_set_si(r8848, mpfr_cmp(r8831, r8847) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8851, r8826, r8850, MPFR_RNDN);
        mpfr_mul(r8852, r8849, r8851, MPFR_RNDN);
        mpfr_add(r8853, r8826, r8852, MPFR_RNDN);
        ;
        mpfr_pow(r8855, r8826, r8838, MPFR_RNDN);
        mpfr_mul(r8856, r8854, r8855, MPFR_RNDN);
        mpfr_sub(r8857, r8853, r8856, MPFR_RNDN);
        mpfr_expm1(r8858, r8836, MPFR_RNDN);
        mpfr_log(r8859, r8858, MPFR_RNDN);
        mpfr_exp(r8860, r8859, MPFR_RNDN);
        if (mpfr_get_si(r8848, MPFR_RNDN)) { mpfr_set(r8861, r8857, MPFR_RNDN); } else { mpfr_set(r8861, r8860, MPFR_RNDN); };
        if (mpfr_get_si(r8833, MPFR_RNDN)) { mpfr_set(r8862, r8846, MPFR_RNDN); } else { mpfr_set(r8862, r8861, MPFR_RNDN); };
        return mpfr_get_d(r8862, MPFR_RNDN);
}

static mpfr_t r8863, r8864, r8865, r8866, r8867, r8868, r8869, r8870, r8871, r8872, r8873, r8874, r8875, r8876, r8877, r8878, r8879, r8880, r8881, r8882, r8883, r8884, r8885, r8886, r8887, r8888, r8889, r8890, r8891, r8892, r8893, r8894, r8895, r8896, r8897, r8898, r8899, r8900, r8901, r8902;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r8863, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r8864, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r8865, "-2", 10, MPFR_RNDN);
        mpfr_init(r8866);
        mpfr_init(r8867);
        mpfr_init(r8868);
        mpfr_init(r8869);
        mpfr_init(r8870);
        mpfr_init(r8871);
        mpfr_init_set_str(r8872, "-1.6888443106222433e-07", 10, MPFR_RNDN);
        mpfr_init(r8873);
        mpfr_init(r8874);
        mpfr_init(r8875);
        mpfr_init(r8876);
        mpfr_init(r8877);
        mpfr_init_set_str(r8878, "3", 10, MPFR_RNDN);
        mpfr_init(r8879);
        mpfr_init(r8880);
        mpfr_init(r8881);
        mpfr_init(r8882);
        mpfr_init(r8883);
        mpfr_init(r8884);
        mpfr_init(r8885);
        mpfr_init(r8886);
        mpfr_init_set_str(r8887, "0.0032654260743848235", 10, MPFR_RNDN);
        mpfr_init(r8888);
        mpfr_init_set_str(r8889, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r8890, "5", 10, MPFR_RNDN);
        mpfr_init(r8891);
        mpfr_init(r8892);
        mpfr_init(r8893);
        mpfr_init_set_str(r8894, "1/3", 10, MPFR_RNDN);
        mpfr_init(r8895);
        mpfr_init(r8896);
        mpfr_init(r8897);
        mpfr_init(r8898);
        mpfr_init(r8899);
        mpfr_init(r8900);
        mpfr_init(r8901);
        mpfr_init(r8902);
}

double f_dm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r8866, x, MPFR_RNDN);
        mpfr_mul(r8867, r8865, r8866, MPFR_RNDN);
        mpfr_exp(r8868, r8867, MPFR_RNDN);
        mpfr_add(r8869, r8864, r8868, MPFR_RNDN);
        mpfr_div(r8870, r8863, r8869, MPFR_RNDN);
        mpfr_sub(r8871, r8870, r8864, MPFR_RNDN);
        ;
        mpfr_set_si(r8873, mpfr_cmp(r8871, r8872) <= 0, MPFR_RNDN);
        mpfr_log(r8874, r8863, MPFR_RNDN);
        mpfr_log1p(r8875, r8868, MPFR_RNDN);
        mpfr_sub(r8876, r8874, r8875, MPFR_RNDN);
        mpfr_exp(r8877, r8876, MPFR_RNDN);
        ;
        mpfr_pow(r8879, r8877, r8878, MPFR_RNDN);
        mpfr_pow(r8880, r8864, r8878, MPFR_RNDN);
        mpfr_sub(r8881, r8879, r8880, MPFR_RNDN);
        mpfr_exp(r8882, r8875, MPFR_RNDN);
        mpfr_div(r8883, r8863, r8882, MPFR_RNDN);
        mpfr_add(r8884, r8864, r8883, MPFR_RNDN);
        mpfr_fma(r8885, r8883, r8883, r8884, MPFR_RNDN);
        mpfr_div(r8886, r8881, r8885, MPFR_RNDN);
        ;
        mpfr_set_si(r8888, mpfr_cmp(r8871, r8887) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8891, r8866, r8890, MPFR_RNDN);
        mpfr_mul(r8892, r8889, r8891, MPFR_RNDN);
        mpfr_add(r8893, r8866, r8892, MPFR_RNDN);
        ;
        mpfr_pow(r8895, r8866, r8878, MPFR_RNDN);
        mpfr_mul(r8896, r8894, r8895, MPFR_RNDN);
        mpfr_sub(r8897, r8893, r8896, MPFR_RNDN);
        mpfr_expm1(r8898, r8876, MPFR_RNDN);
        mpfr_log(r8899, r8898, MPFR_RNDN);
        mpfr_exp(r8900, r8899, MPFR_RNDN);
        if (mpfr_get_si(r8888, MPFR_RNDN)) { mpfr_set(r8901, r8897, MPFR_RNDN); } else { mpfr_set(r8901, r8900, MPFR_RNDN); };
        if (mpfr_get_si(r8873, MPFR_RNDN)) { mpfr_set(r8902, r8886, MPFR_RNDN); } else { mpfr_set(r8902, r8901, MPFR_RNDN); };
        return mpfr_get_d(r8902, MPFR_RNDN);
}