#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 r22230 = 2;
        float r22231 = 1;
        float r22232 = -2;
        float r22233 = x;
        float r22234 = r22232 * r22233;
        float r22235 = exp(r22234);
        float r22236 = r22231 + r22235;
        float r22237 = r22230 / r22236;
        float r22238 = r22237 - r22231;
        return r22238;
}

double f_id(double x, double __attribute__((unused)) y) {
        double r22239 = 2;
        double r22240 = 1;
        double r22241 = -2;
        double r22242 = x;
        double r22243 = r22241 * r22242;
        double r22244 = exp(r22243);
        double r22245 = r22240 + r22244;
        double r22246 = r22239 / r22245;
        double r22247 = r22246 - r22240;
        return r22247;
}


double f_of(float x, float __attribute__((unused)) y) {
        float r22248 = 2;
        float r22249 = 1;
        float r22250 = -2;
        float r22251 = x;
        float r22252 = r22250 * r22251;
        float r22253 = exp(r22252);
        float r22254 = r22249 + r22253;
        float r22255 = r22248 / r22254;
        float r22256 = r22255 - r22249;
        float r22257 = -2.7755575615628914e-17;
        bool r22258 = r22256 <= r22257;
        float r22259 = 8.353037248554074e-07;
        bool r22260 = r22256 <= r22259;
        float r22261 = 2/15;
        float r22262 = 5;
        float r22263 = pow(r22251, r22262);
        float r22264 = r22261 * r22263;
        float r22265 = r22264 + r22251;
        float r22266 = 1/3;
        float r22267 = 3;
        float r22268 = pow(r22251, r22267);
        float r22269 = r22266 * r22268;
        float r22270 = r22265 - r22269;
        float r22271 = r22260 ? r22270 : r22256;
        float r22272 = r22258 ? r22256 : r22271;
        return r22272;
}

double f_od(double x, double __attribute__((unused)) y) {
        double r22273 = 2;
        double r22274 = 1;
        double r22275 = -2;
        double r22276 = x;
        double r22277 = r22275 * r22276;
        double r22278 = exp(r22277);
        double r22279 = r22274 + r22278;
        double r22280 = r22273 / r22279;
        double r22281 = r22280 - r22274;
        double r22282 = -2.7755575615628914e-17;
        bool r22283 = r22281 <= r22282;
        double r22284 = 8.353037248554074e-07;
        bool r22285 = r22281 <= r22284;
        double r22286 = 2/15;
        double r22287 = 5;
        double r22288 = pow(r22276, r22287);
        double r22289 = r22286 * r22288;
        double r22290 = r22289 + r22276;
        double r22291 = 1/3;
        double r22292 = 3;
        double r22293 = pow(r22276, r22292);
        double r22294 = r22291 * r22293;
        double r22295 = r22290 - r22294;
        double r22296 = r22285 ? r22295 : r22281;
        double r22297 = r22283 ? r22281 : r22296;
        return r22297;
}

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 r22298, r22299, r22300, r22301, r22302, r22303, r22304, r22305, r22306;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22298, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22299, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22300, "-2", 10, MPFR_RNDN);
        mpfr_init(r22301);
        mpfr_init(r22302);
        mpfr_init(r22303);
        mpfr_init(r22304);
        mpfr_init(r22305);
        mpfr_init(r22306);
}

double f_im(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22301, x, MPFR_RNDN);
        mpfr_mul(r22302, r22300, r22301, MPFR_RNDN);
        mpfr_exp(r22303, r22302, MPFR_RNDN);
        mpfr_add(r22304, r22299, r22303, MPFR_RNDN);
        mpfr_div(r22305, r22298, r22304, MPFR_RNDN);
        mpfr_sub(r22306, r22305, r22299, MPFR_RNDN);
        return mpfr_get_d(r22306, MPFR_RNDN);
}

static mpfr_t r22307, r22308, r22309, r22310, r22311, r22312, r22313, r22314, r22315, r22316, r22317, r22318, r22319, r22320, r22321, r22322, r22323, r22324, r22325, r22326, r22327, r22328, r22329, r22330, r22331;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22307, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22308, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22309, "-2", 10, MPFR_RNDN);
        mpfr_init(r22310);
        mpfr_init(r22311);
        mpfr_init(r22312);
        mpfr_init(r22313);
        mpfr_init(r22314);
        mpfr_init(r22315);
        mpfr_init_set_str(r22316, "-2.7755575615628914e-17", 10, MPFR_RNDN);
        mpfr_init(r22317);
        mpfr_init_set_str(r22318, "8.353037248554074e-07", 10, MPFR_RNDN);
        mpfr_init(r22319);
        mpfr_init_set_str(r22320, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22321, "5", 10, MPFR_RNDN);
        mpfr_init(r22322);
        mpfr_init(r22323);
        mpfr_init(r22324);
        mpfr_init_set_str(r22325, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r22326, "3", 10, MPFR_RNDN);
        mpfr_init(r22327);
        mpfr_init(r22328);
        mpfr_init(r22329);
        mpfr_init(r22330);
        mpfr_init(r22331);
}

double f_fm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22310, x, MPFR_RNDN);
        mpfr_mul(r22311, r22309, r22310, MPFR_RNDN);
        mpfr_exp(r22312, r22311, MPFR_RNDN);
        mpfr_add(r22313, r22308, r22312, MPFR_RNDN);
        mpfr_div(r22314, r22307, r22313, MPFR_RNDN);
        mpfr_sub(r22315, r22314, r22308, MPFR_RNDN);
        ;
        mpfr_set_si(r22317, mpfr_cmp(r22315, r22316) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r22319, mpfr_cmp(r22315, r22318) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22322, r22310, r22321, MPFR_RNDN);
        mpfr_mul(r22323, r22320, r22322, MPFR_RNDN);
        mpfr_add(r22324, r22323, r22310, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22327, r22310, r22326, MPFR_RNDN);
        mpfr_mul(r22328, r22325, r22327, MPFR_RNDN);
        mpfr_sub(r22329, r22324, r22328, MPFR_RNDN);
        if (mpfr_get_si(r22319, MPFR_RNDN)) { mpfr_set(r22330, r22329, MPFR_RNDN); } else { mpfr_set(r22330, r22315, MPFR_RNDN); };
        if (mpfr_get_si(r22317, MPFR_RNDN)) { mpfr_set(r22331, r22315, MPFR_RNDN); } else { mpfr_set(r22331, r22330, MPFR_RNDN); };
        return mpfr_get_d(r22331, MPFR_RNDN);
}

static mpfr_t r22332, r22333, r22334, r22335, r22336, r22337, r22338, r22339, r22340, r22341, r22342, r22343, r22344, r22345, r22346, r22347, r22348, r22349, r22350, r22351, r22352, r22353, r22354, r22355, r22356;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r22332, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22333, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22334, "-2", 10, MPFR_RNDN);
        mpfr_init(r22335);
        mpfr_init(r22336);
        mpfr_init(r22337);
        mpfr_init(r22338);
        mpfr_init(r22339);
        mpfr_init(r22340);
        mpfr_init_set_str(r22341, "-2.7755575615628914e-17", 10, MPFR_RNDN);
        mpfr_init(r22342);
        mpfr_init_set_str(r22343, "8.353037248554074e-07", 10, MPFR_RNDN);
        mpfr_init(r22344);
        mpfr_init_set_str(r22345, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22346, "5", 10, MPFR_RNDN);
        mpfr_init(r22347);
        mpfr_init(r22348);
        mpfr_init(r22349);
        mpfr_init_set_str(r22350, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r22351, "3", 10, MPFR_RNDN);
        mpfr_init(r22352);
        mpfr_init(r22353);
        mpfr_init(r22354);
        mpfr_init(r22355);
        mpfr_init(r22356);
}

double f_dm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22335, x, MPFR_RNDN);
        mpfr_mul(r22336, r22334, r22335, MPFR_RNDN);
        mpfr_exp(r22337, r22336, MPFR_RNDN);
        mpfr_add(r22338, r22333, r22337, MPFR_RNDN);
        mpfr_div(r22339, r22332, r22338, MPFR_RNDN);
        mpfr_sub(r22340, r22339, r22333, MPFR_RNDN);
        ;
        mpfr_set_si(r22342, mpfr_cmp(r22340, r22341) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r22344, mpfr_cmp(r22340, r22343) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22347, r22335, r22346, MPFR_RNDN);
        mpfr_mul(r22348, r22345, r22347, MPFR_RNDN);
        mpfr_add(r22349, r22348, r22335, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22352, r22335, r22351, MPFR_RNDN);
        mpfr_mul(r22353, r22350, r22352, MPFR_RNDN);
        mpfr_sub(r22354, r22349, r22353, MPFR_RNDN);
        if (mpfr_get_si(r22344, MPFR_RNDN)) { mpfr_set(r22355, r22354, MPFR_RNDN); } else { mpfr_set(r22355, r22340, MPFR_RNDN); };
        if (mpfr_get_si(r22342, MPFR_RNDN)) { mpfr_set(r22356, r22340, MPFR_RNDN); } else { mpfr_set(r22356, r22355, MPFR_RNDN); };
        return mpfr_get_d(r22356, MPFR_RNDN);
}

