#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 r22357 = 2;
        float r22358 = 1;
        float r22359 = -2;
        float r22360 = x;
        float r22361 = r22359 * r22360;
        float r22362 = exp(r22361);
        float r22363 = r22358 + r22362;
        float r22364 = r22357 / r22363;
        float r22365 = r22364 - r22358;
        return r22365;
}

double f_id(double x, double __attribute__((unused)) y) {
        double r22366 = 2;
        double r22367 = 1;
        double r22368 = -2;
        double r22369 = x;
        double r22370 = r22368 * r22369;
        double r22371 = exp(r22370);
        double r22372 = r22367 + r22371;
        double r22373 = r22366 / r22372;
        double r22374 = r22373 - r22367;
        return r22374;
}


double f_of(float x, float __attribute__((unused)) y) {
        float r22375 = 2;
        float r22376 = 1;
        float r22377 = -2;
        float r22378 = x;
        float r22379 = r22377 * r22378;
        float r22380 = exp(r22379);
        float r22381 = r22376 + r22380;
        float r22382 = r22375 / r22381;
        float r22383 = r22382 - r22376;
        float r22384 = -0.0010361760913228334;
        bool r22385 = r22383 <= r22384;
        float r22386 = 3;
        float r22387 = pow(r22382, r22386);
        float r22388 = pow(r22376, r22386);
        float r22389 = r22387 - r22388;
        float r22390 = r22382 + r22376;
        float r22391 = r22382 * r22382;
        float r22392 = r22390 + r22391;
        float r22393 = r22389 / r22392;
        float r22394 = 0.0039062543158841865;
        bool r22395 = r22383 <= r22394;
        float r22396 = 2/15;
        float r22397 = 5;
        float r22398 = pow(r22378, r22397);
        float r22399 = r22396 * r22398;
        float r22400 = r22399 + r22378;
        float r22401 = 1/3;
        float r22402 = pow(r22378, r22386);
        float r22403 = r22401 * r22402;
        float r22404 = r22400 - r22403;
        float r22405 = r22395 ? r22404 : r22393;
        float r22406 = r22385 ? r22393 : r22405;
        return r22406;
}

double f_od(double x, double __attribute__((unused)) y) {
        double r22407 = 2;
        double r22408 = 1;
        double r22409 = -2;
        double r22410 = x;
        double r22411 = r22409 * r22410;
        double r22412 = exp(r22411);
        double r22413 = r22408 + r22412;
        double r22414 = r22407 / r22413;
        double r22415 = r22414 - r22408;
        double r22416 = -0.0010361760913228334;
        bool r22417 = r22415 <= r22416;
        double r22418 = 3;
        double r22419 = pow(r22414, r22418);
        double r22420 = pow(r22408, r22418);
        double r22421 = r22419 - r22420;
        double r22422 = r22414 + r22408;
        double r22423 = r22414 * r22414;
        double r22424 = r22422 + r22423;
        double r22425 = r22421 / r22424;
        double r22426 = 0.0039062543158841865;
        bool r22427 = r22415 <= r22426;
        double r22428 = 2/15;
        double r22429 = 5;
        double r22430 = pow(r22410, r22429);
        double r22431 = r22428 * r22430;
        double r22432 = r22431 + r22410;
        double r22433 = 1/3;
        double r22434 = pow(r22410, r22418);
        double r22435 = r22433 * r22434;
        double r22436 = r22432 - r22435;
        double r22437 = r22427 ? r22436 : r22425;
        double r22438 = r22417 ? r22425 : r22437;
        return r22438;
}

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 r22439, r22440, r22441, r22442, r22443, r22444, r22445, r22446, r22447;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22439, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22440, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22441, "-2", 10, MPFR_RNDN);
        mpfr_init(r22442);
        mpfr_init(r22443);
        mpfr_init(r22444);
        mpfr_init(r22445);
        mpfr_init(r22446);
        mpfr_init(r22447);
}

double f_im(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22442, x, MPFR_RNDN);
        mpfr_mul(r22443, r22441, r22442, MPFR_RNDN);
        mpfr_exp(r22444, r22443, MPFR_RNDN);
        mpfr_add(r22445, r22440, r22444, MPFR_RNDN);
        mpfr_div(r22446, r22439, r22445, MPFR_RNDN);
        mpfr_sub(r22447, r22446, r22440, MPFR_RNDN);
        return mpfr_get_d(r22447, MPFR_RNDN);
}

static mpfr_t 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, r22476, r22477, r22478, r22479;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22448, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22449, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22450, "-2", 10, MPFR_RNDN);
        mpfr_init(r22451);
        mpfr_init(r22452);
        mpfr_init(r22453);
        mpfr_init(r22454);
        mpfr_init(r22455);
        mpfr_init(r22456);
        mpfr_init_set_str(r22457, "-0.0010361760913228334", 10, MPFR_RNDN);
        mpfr_init(r22458);
        mpfr_init_set_str(r22459, "3", 10, MPFR_RNDN);
        mpfr_init(r22460);
        mpfr_init(r22461);
        mpfr_init(r22462);
        mpfr_init(r22463);
        mpfr_init(r22464);
        mpfr_init(r22465);
        mpfr_init(r22466);
        mpfr_init_set_str(r22467, "0.0039062543158841865", 10, MPFR_RNDN);
        mpfr_init(r22468);
        mpfr_init_set_str(r22469, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22470, "5", 10, MPFR_RNDN);
        mpfr_init(r22471);
        mpfr_init(r22472);
        mpfr_init(r22473);
        mpfr_init_set_str(r22474, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22475);
        mpfr_init(r22476);
        mpfr_init(r22477);
        mpfr_init(r22478);
        mpfr_init(r22479);
}

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

static mpfr_t r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510, r22511;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r22480, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r22481, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22482, "-2", 10, MPFR_RNDN);
        mpfr_init(r22483);
        mpfr_init(r22484);
        mpfr_init(r22485);
        mpfr_init(r22486);
        mpfr_init(r22487);
        mpfr_init(r22488);
        mpfr_init_set_str(r22489, "-0.0010361760913228334", 10, MPFR_RNDN);
        mpfr_init(r22490);
        mpfr_init_set_str(r22491, "3", 10, MPFR_RNDN);
        mpfr_init(r22492);
        mpfr_init(r22493);
        mpfr_init(r22494);
        mpfr_init(r22495);
        mpfr_init(r22496);
        mpfr_init(r22497);
        mpfr_init(r22498);
        mpfr_init_set_str(r22499, "0.0039062543158841865", 10, MPFR_RNDN);
        mpfr_init(r22500);
        mpfr_init_set_str(r22501, "2/15", 10, MPFR_RNDN);
        mpfr_init_set_str(r22502, "5", 10, MPFR_RNDN);
        mpfr_init(r22503);
        mpfr_init(r22504);
        mpfr_init(r22505);
        mpfr_init_set_str(r22506, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22507);
        mpfr_init(r22508);
        mpfr_init(r22509);
        mpfr_init(r22510);
        mpfr_init(r22511);
}

double f_dm(double x, double __attribute__((unused)) y) {
        ;
        ;
        ;
        mpfr_set_d(r22483, x, MPFR_RNDN);
        mpfr_mul(r22484, r22482, r22483, MPFR_RNDN);
        mpfr_exp(r22485, r22484, MPFR_RNDN);
        mpfr_add(r22486, r22481, r22485, MPFR_RNDN);
        mpfr_div(r22487, r22480, r22486, MPFR_RNDN);
        mpfr_sub(r22488, r22487, r22481, MPFR_RNDN);
        ;
        mpfr_set_si(r22490, mpfr_cmp(r22488, r22489) <= 0, MPFR_RNDN);
        ;
        mpfr_pow(r22492, r22487, r22491, MPFR_RNDN);
        mpfr_pow(r22493, r22481, r22491, MPFR_RNDN);
        mpfr_sub(r22494, r22492, r22493, MPFR_RNDN);
        mpfr_add(r22495, r22487, r22481, MPFR_RNDN);
        mpfr_mul(r22496, r22487, r22487, MPFR_RNDN);
        mpfr_add(r22497, r22495, r22496, MPFR_RNDN);
        mpfr_div(r22498, r22494, r22497, MPFR_RNDN);
        ;
        mpfr_set_si(r22500, mpfr_cmp(r22488, r22499) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22503, r22483, r22502, MPFR_RNDN);
        mpfr_mul(r22504, r22501, r22503, MPFR_RNDN);
        mpfr_add(r22505, r22504, r22483, MPFR_RNDN);
        ;
        mpfr_pow(r22507, r22483, r22491, MPFR_RNDN);
        mpfr_mul(r22508, r22506, r22507, MPFR_RNDN);
        mpfr_sub(r22509, r22505, r22508, MPFR_RNDN);
        if (mpfr_get_si(r22500, MPFR_RNDN)) { mpfr_set(r22510, r22509, MPFR_RNDN); } else { mpfr_set(r22510, r22498, MPFR_RNDN); };
        if (mpfr_get_si(r22490, MPFR_RNDN)) { mpfr_set(r22511, r22498, MPFR_RNDN); } else { mpfr_set(r22511, r22510, MPFR_RNDN); };
        return mpfr_get_d(r22511, MPFR_RNDN);
}