#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 r22321 = 2;
        float r22322 = 1;
        float r22323 = -2;
        float r22324 = x;
        float r22325 = r22323 * r22324;
        float r22326 = exp(r22325);
        float r22327 = r22322 + r22326;
        float r22328 = r22321 / r22327;
        float r22329 = r22328 - r22322;
        return r22329;
}

double f_id(double x, double __attribute__((unused)) y) {
        double r22330 = 2;
        double r22331 = 1;
        double r22332 = -2;
        double r22333 = x;
        double r22334 = r22332 * r22333;
        double r22335 = exp(r22334);
        double r22336 = r22331 + r22335;
        double r22337 = r22330 / r22336;
        double r22338 = r22337 - r22331;
        return r22338;
}


double f_of(float x, float __attribute__((unused)) y) {
        float r22339 = 2;
        float r22340 = 1;
        float r22341 = -2;
        float r22342 = x;
        float r22343 = r22341 * r22342;
        float r22344 = exp(r22343);
        float r22345 = r22340 + r22344;
        float r22346 = r22339 / r22345;
        float r22347 = r22346 - r22340;
        float r22348 = -0.0021405943016779886;
        bool r22349 = r22347 <= r22348;
        float r22350 = 3;
        float r22351 = pow(r22346, r22350);
        float r22352 = pow(r22340, r22350);
        float r22353 = r22351 - r22352;
        float r22354 = r22346 + r22340;
        float r22355 = r22346 * r22346;
        float r22356 = r22354 + r22355;
        float r22357 = r22353 / r22356;
        float r22358 = 0.06250821168851108;
        bool r22359 = r22347 <= r22358;
        float r22360 = 2/15;
        float r22361 = 5;
        float r22362 = pow(r22342, r22361);
        float r22363 = r22360 * r22362;
        float r22364 = r22363 + r22342;
        float r22365 = 1/3;
        float r22366 = pow(r22342, r22350);
        float r22367 = r22365 * r22366;
        float r22368 = r22364 - r22367;
        float r22369 = r22359 ? r22368 : r22357;
        float r22370 = r22349 ? r22357 : r22369;
        return r22370;
}

double f_od(double x, double __attribute__((unused)) y) {
        double r22371 = 2;
        double r22372 = 1;
        double r22373 = -2;
        double r22374 = x;
        double r22375 = r22373 * r22374;
        double r22376 = exp(r22375);
        double r22377 = r22372 + r22376;
        double r22378 = r22371 / r22377;
        double r22379 = r22378 - r22372;
        double r22380 = -0.0021405943016779886;
        bool r22381 = r22379 <= r22380;
        double r22382 = 3;
        double r22383 = pow(r22378, r22382);
        double r22384 = pow(r22372, r22382);
        double r22385 = r22383 - r22384;
        double r22386 = r22378 + r22372;
        double r22387 = r22378 * r22378;
        double r22388 = r22386 + r22387;
        double r22389 = r22385 / r22388;
        double r22390 = 0.06250821168851108;
        bool r22391 = r22379 <= r22390;
        double r22392 = 2/15;
        double r22393 = 5;
        double r22394 = pow(r22374, r22393);
        double r22395 = r22392 * r22394;
        double r22396 = r22395 + r22374;
        double r22397 = 1/3;
        double r22398 = pow(r22374, r22382);
        double r22399 = r22397 * r22398;
        double r22400 = r22396 - r22399;
        double r22401 = r22391 ? r22400 : r22389;
        double r22402 = r22381 ? r22389 : r22401;
        return r22402;
}

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 r22403, r22404, r22405, r22406, r22407, r22408, r22409, r22410, r22411;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22403, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22404, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22405, "-2", 10, MPFR_RNDN);
        mpfr_init(r22406);
        mpfr_init(r22407);
        mpfr_init(r22408);
        mpfr_init(r22409);
        mpfr_init(r22410);
        mpfr_init(r22411);
}

double f_im(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22406, x, MPFR_RNDN);
        mpfr_mul(r22407, r22405, r22406, MPFR_RNDN);
        mpfr_exp(r22408, r22407, MPFR_RNDN);
        mpfr_add(r22409, r22404, r22408, MPFR_RNDN);
        mpfr_div(r22410, r22403, r22409, MPFR_RNDN);
        mpfr_sub(r22411, r22410, r22404, MPFR_RNDN);
        return mpfr_get_d(r22411, MPFR_RNDN);
}

static mpfr_t r22412, r22413, r22414, r22415, r22416, r22417, r22418, r22419, r22420, r22421, r22422, r22423, r22424, r22425, r22426, r22427, r22428, r22429, r22430, r22431, r22432, r22433, r22434, r22435, r22436, r22437, r22438, r22439, r22440, r22441, r22442, r22443;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22412, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22413, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22414, "-2", 10, MPFR_RNDN);
        mpfr_init(r22415);
        mpfr_init(r22416);
        mpfr_init(r22417);
        mpfr_init(r22418);
        mpfr_init(r22419);
        mpfr_init(r22420);
        mpfr_init_set_str(r22421, "-0.0021405943016779886", 10, MPFR_RNDN);
        mpfr_init(r22422);
        mpfr_init_set_str(r22423, "3", 10, MPFR_RNDN);
        mpfr_init(r22424);
        mpfr_init(r22425);
        mpfr_init(r22426);
        mpfr_init(r22427);
        mpfr_init(r22428);
        mpfr_init(r22429);
        mpfr_init(r22430);
        mpfr_init_set_str(r22431, "0.06250821168851108", 10, MPFR_RNDN);
        mpfr_init(r22432);
        mpfr_init_set_str(r22433, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22434, "5", 10, MPFR_RNDN);
        mpfr_init(r22435);
        mpfr_init(r22436);
        mpfr_init(r22437);
        mpfr_init_set_str(r22438, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22439);
        mpfr_init(r22440);
        mpfr_init(r22441);
        mpfr_init(r22442);
        mpfr_init(r22443);
}

double f_fm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22415, x, MPFR_RNDN);
        mpfr_mul(r22416, r22414, r22415, MPFR_RNDN);
        mpfr_exp(r22417, r22416, MPFR_RNDN);
        mpfr_add(r22418, r22413, r22417, MPFR_RNDN);
        mpfr_div(r22419, r22412, r22418, MPFR_RNDN);
        mpfr_sub(r22420, r22419, r22413, MPFR_RNDN);
        ;
        mpfr_set_si(r22422, mpfr_cmp(r22420, r22421) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r22424, r22419, r22423, MPFR_RNDN);
        mpfr_pow(r22425, r22413, r22423, MPFR_RNDN);
        mpfr_sub(r22426, r22424, r22425, MPFR_RNDN);
        mpfr_add(r22427, r22419, r22413, MPFR_RNDN);
        mpfr_mul(r22428, r22419, r22419, MPFR_RNDN);
        mpfr_add(r22429, r22427, r22428, MPFR_RNDN);
        mpfr_div(r22430, r22426, r22429, MPFR_RNDN);
        ;
        mpfr_set_si(r22432, mpfr_cmp(r22420, r22431) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22435, r22415, r22434, MPFR_RNDN);
        mpfr_mul(r22436, r22433, r22435, MPFR_RNDN);
        mpfr_add(r22437, r22436, r22415, MPFR_RNDN);
        ;
        mpfr_pow(r22439, r22415, r22423, MPFR_RNDN);
        mpfr_mul(r22440, r22438, r22439, MPFR_RNDN);
        mpfr_sub(r22441, r22437, r22440, MPFR_RNDN);
        if (mpfr_get_si(r22432, MPFR_RNDN)) { mpfr_set(r22442, r22441, MPFR_RNDN); } else { mpfr_set(r22442, r22430, MPFR_RNDN); };
        if (mpfr_get_si(r22422, MPFR_RNDN)) { mpfr_set(r22443, r22430, MPFR_RNDN); } else { mpfr_set(r22443, r22442, MPFR_RNDN); };
        return mpfr_get_d(r22443, MPFR_RNDN);
}

static mpfr_t r22444, r22445, r22446, r22447, r22448, r22449, r22450, r22451, r22452, r22453, r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461, r22462, r22463, r22464, r22465, r22466, r22467, r22468, r22469, r22470, r22471, r22472, r22473, r22474, r22475;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22444, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22445, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22446, "-2", 10, MPFR_RNDN);
        mpfr_init(r22447);
        mpfr_init(r22448);
        mpfr_init(r22449);
        mpfr_init(r22450);
        mpfr_init(r22451);
        mpfr_init(r22452);
        mpfr_init_set_str(r22453, "-0.0021405943016779886", 10, MPFR_RNDN);
        mpfr_init(r22454);
        mpfr_init_set_str(r22455, "3", 10, MPFR_RNDN);
        mpfr_init(r22456);
        mpfr_init(r22457);
        mpfr_init(r22458);
        mpfr_init(r22459);
        mpfr_init(r22460);
        mpfr_init(r22461);
        mpfr_init(r22462);
        mpfr_init_set_str(r22463, "0.06250821168851108", 10, MPFR_RNDN);
        mpfr_init(r22464);
        mpfr_init_set_str(r22465, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22466, "5", 10, MPFR_RNDN);
        mpfr_init(r22467);
        mpfr_init(r22468);
        mpfr_init(r22469);
        mpfr_init_set_str(r22470, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22471);
        mpfr_init(r22472);
        mpfr_init(r22473);
        mpfr_init(r22474);
        mpfr_init(r22475);
}

double f_dm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22447, x, MPFR_RNDN);
        mpfr_mul(r22448, r22446, r22447, MPFR_RNDN);
        mpfr_exp(r22449, r22448, MPFR_RNDN);
        mpfr_add(r22450, r22445, r22449, MPFR_RNDN);
        mpfr_div(r22451, r22444, r22450, MPFR_RNDN);
        mpfr_sub(r22452, r22451, r22445, MPFR_RNDN);
        ;
        mpfr_set_si(r22454, mpfr_cmp(r22452, r22453) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r22456, r22451, r22455, MPFR_RNDN);
        mpfr_pow(r22457, r22445, r22455, MPFR_RNDN);
        mpfr_sub(r22458, r22456, r22457, MPFR_RNDN);
        mpfr_add(r22459, r22451, r22445, MPFR_RNDN);
        mpfr_mul(r22460, r22451, r22451, MPFR_RNDN);
        mpfr_add(r22461, r22459, r22460, MPFR_RNDN);
        mpfr_div(r22462, r22458, r22461, MPFR_RNDN);
        ;
        mpfr_set_si(r22464, mpfr_cmp(r22452, r22463) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22467, r22447, r22466, MPFR_RNDN);
        mpfr_mul(r22468, r22465, r22467, MPFR_RNDN);
        mpfr_add(r22469, r22468, r22447, MPFR_RNDN);
        ;
        mpfr_pow(r22471, r22447, r22455, MPFR_RNDN);
        mpfr_mul(r22472, r22470, r22471, MPFR_RNDN);
        mpfr_sub(r22473, r22469, r22472, MPFR_RNDN);
        if (mpfr_get_si(r22464, MPFR_RNDN)) { mpfr_set(r22474, r22473, MPFR_RNDN); } else { mpfr_set(r22474, r22462, MPFR_RNDN); };
        if (mpfr_get_si(r22454, MPFR_RNDN)) { mpfr_set(r22475, r22462, MPFR_RNDN); } else { mpfr_set(r22475, r22474, MPFR_RNDN); };
        return mpfr_get_d(r22475, MPFR_RNDN);
}