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

char *name = "Kahan's exp quotient";

double f_if(float x) {
        float r23312 = x;
        float r23313 = exp(r23312);
        float r23314 = 1;
        float r23315 = r23313 - r23314;
        float r23316 = r23315 / r23312;
        return r23316;
}

double f_id(double x) {
        double r23317 = x;
        double r23318 = exp(r23317);
        double r23319 = 1;
        double r23320 = r23318 - r23319;
        double r23321 = r23320 / r23317;
        return r23321;
}


double f_of(float x) {
        float r23322 = x;
        float r23323 = exp(r23322);
        float r23324 = 1;
        float r23325 = r23323 - r23324;
        float r23326 = r23325 / r23322;
        float r23327 = 0.0;
        bool r23328 = r23326 <= r23327;
        float r23329 = 1/6;
        float r23330 = 2;
        float r23331 = pow(r23322, r23330);
        float r23332 = r23329 * r23331;
        float r23333 = 1/2;
        float r23334 = r23333 * r23322;
        float r23335 = r23324 + r23334;
        float r23336 = r23332 + r23335;
        float r23337 = 0.049221242662298555;
        bool r23338 = r23326 <= r23337;
        float r23339 = r23322 + r23322;
        float r23340 = exp(r23339);
        float r23341 = r23340 - r23324;
        float r23342 = r23323 + r23324;
        float r23343 = r23341 / r23342;
        float r23344 = r23343 / r23322;
        float r23345 = r23338 ? r23344 : r23336;
        float r23346 = r23328 ? r23336 : r23345;
        return r23346;
}

double f_od(double x) {
        double r23347 = x;
        double r23348 = exp(r23347);
        double r23349 = 1;
        double r23350 = r23348 - r23349;
        double r23351 = r23350 / r23347;
        double r23352 = 0.0;
        bool r23353 = r23351 <= r23352;
        double r23354 = 1/6;
        double r23355 = 2;
        double r23356 = pow(r23347, r23355);
        double r23357 = r23354 * r23356;
        double r23358 = 1/2;
        double r23359 = r23358 * r23347;
        double r23360 = r23349 + r23359;
        double r23361 = r23357 + r23360;
        double r23362 = 0.049221242662298555;
        bool r23363 = r23351 <= r23362;
        double r23364 = r23347 + r23347;
        double r23365 = exp(r23364);
        double r23366 = r23365 - r23349;
        double r23367 = r23348 + r23349;
        double r23368 = r23366 / r23367;
        double r23369 = r23368 / r23347;
        double r23370 = r23363 ? r23369 : r23361;
        double r23371 = r23353 ? r23361 : r23370;
        return r23371;
}

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 r23372, r23373, r23374, r23375, r23376;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r23372);
        mpfr_init(r23373);
        mpfr_init_set_str(r23374, "1", 10, MPFR_RNDN);
        mpfr_init(r23375);
        mpfr_init(r23376);
}

double f_im(double x) {
        mpfr_set_d(r23372, x, MPFR_RNDN);
        mpfr_exp(r23373, r23372, MPFR_RNDN);
        ;
        mpfr_sub(r23375, r23373, r23374, MPFR_RNDN);
        mpfr_div(r23376, r23375, r23372, MPFR_RNDN);
        return mpfr_get_d(r23376, MPFR_RNDN);
}

static mpfr_t r23377, r23378, r23379, r23380, r23381, r23382, r23383, r23384, r23385, r23386, r23387, r23388, r23389, r23390, r23391, r23392, r23393, r23394, r23395, r23396, r23397, r23398, r23399, r23400, r23401;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r23377);
        mpfr_init(r23378);
        mpfr_init_set_str(r23379, "1", 10, MPFR_RNDN);
        mpfr_init(r23380);
        mpfr_init(r23381);
        mpfr_init_set_str(r23382, "0.0", 10, MPFR_RNDN);
        mpfr_init(r23383);
        mpfr_init_set_str(r23384, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r23385, "2", 10, MPFR_RNDN);
        mpfr_init(r23386);
        mpfr_init(r23387);
        mpfr_init_set_str(r23388, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23389);
        mpfr_init(r23390);
        mpfr_init(r23391);
        mpfr_init_set_str(r23392, "0.049221242662298555", 10, MPFR_RNDN);
        mpfr_init(r23393);
        mpfr_init(r23394);
        mpfr_init(r23395);
        mpfr_init(r23396);
        mpfr_init(r23397);
        mpfr_init(r23398);
        mpfr_init(r23399);
        mpfr_init(r23400);
        mpfr_init(r23401);
}

double f_fm(double x) {
        mpfr_set_d(r23377, x, MPFR_RNDN);
        mpfr_exp(r23378, r23377, MPFR_RNDN);
        ;
        mpfr_sub(r23380, r23378, r23379, MPFR_RNDN);
        mpfr_div(r23381, r23380, r23377, MPFR_RNDN);
        ;
        mpfr_set_si(r23383, mpfr_cmp(r23381, r23382) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r23386, r23377, r23385, MPFR_RNDN);
        mpfr_mul(r23387, r23384, r23386, MPFR_RNDN);
        ;
        mpfr_mul(r23389, r23388, r23377, MPFR_RNDN);
        mpfr_add(r23390, r23379, r23389, MPFR_RNDN);
        mpfr_add(r23391, r23387, r23390, MPFR_RNDN);
        ;
        mpfr_set_si(r23393, mpfr_cmp(r23381, r23392) <= 0, MPFR_RNDN);
        mpfr_add(r23394, r23377, r23377, MPFR_RNDN);
        mpfr_exp(r23395, r23394, MPFR_RNDN);
        mpfr_sub(r23396, r23395, r23379, MPFR_RNDN);
        mpfr_add(r23397, r23378, r23379, MPFR_RNDN);
        mpfr_div(r23398, r23396, r23397, MPFR_RNDN);
        mpfr_div(r23399, r23398, r23377, MPFR_RNDN);
        if (mpfr_get_si(r23393, MPFR_RNDN)) { mpfr_set(r23400, r23399, MPFR_RNDN); } else { mpfr_set(r23400, r23391, MPFR_RNDN); };
        if (mpfr_get_si(r23383, MPFR_RNDN)) { mpfr_set(r23401, r23391, MPFR_RNDN); } else { mpfr_set(r23401, r23400, MPFR_RNDN); };
        return mpfr_get_d(r23401, MPFR_RNDN);
}

static mpfr_t r23402, r23403, r23404, r23405, r23406, r23407, r23408, r23409, r23410, r23411, r23412, r23413, r23414, r23415, r23416, r23417, r23418, r23419, r23420, r23421, r23422, r23423, r23424, r23425, r23426;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r23402);
        mpfr_init(r23403);
        mpfr_init_set_str(r23404, "1", 10, MPFR_RNDN);
        mpfr_init(r23405);
        mpfr_init(r23406);
        mpfr_init_set_str(r23407, "0.0", 10, MPFR_RNDN);
        mpfr_init(r23408);
        mpfr_init_set_str(r23409, "1/6", 10, MPFR_RNDN);
        mpfr_init_set_str(r23410, "2", 10, MPFR_RNDN);
        mpfr_init(r23411);
        mpfr_init(r23412);
        mpfr_init_set_str(r23413, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23414);
        mpfr_init(r23415);
        mpfr_init(r23416);
        mpfr_init_set_str(r23417, "0.049221242662298555", 10, MPFR_RNDN);
        mpfr_init(r23418);
        mpfr_init(r23419);
        mpfr_init(r23420);
        mpfr_init(r23421);
        mpfr_init(r23422);
        mpfr_init(r23423);
        mpfr_init(r23424);
        mpfr_init(r23425);
        mpfr_init(r23426);
}

double f_dm(double x) {
        mpfr_set_d(r23402, x, MPFR_RNDN);
        mpfr_exp(r23403, r23402, MPFR_RNDN);
        ;
        mpfr_sub(r23405, r23403, r23404, MPFR_RNDN);
        mpfr_div(r23406, r23405, r23402, MPFR_RNDN);
        ;
        mpfr_set_si(r23408, mpfr_cmp(r23406, r23407) <= 0, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r23411, r23402, r23410, MPFR_RNDN);
        mpfr_mul(r23412, r23409, r23411, MPFR_RNDN);
        ;
        mpfr_mul(r23414, r23413, r23402, MPFR_RNDN);
        mpfr_add(r23415, r23404, r23414, MPFR_RNDN);
        mpfr_add(r23416, r23412, r23415, MPFR_RNDN);
        ;
        mpfr_set_si(r23418, mpfr_cmp(r23406, r23417) <= 0, MPFR_RNDN);
        mpfr_add(r23419, r23402, r23402, MPFR_RNDN);
        mpfr_exp(r23420, r23419, MPFR_RNDN);
        mpfr_sub(r23421, r23420, r23404, MPFR_RNDN);
        mpfr_add(r23422, r23403, r23404, MPFR_RNDN);
        mpfr_div(r23423, r23421, r23422, MPFR_RNDN);
        mpfr_div(r23424, r23423, r23402, MPFR_RNDN);
        if (mpfr_get_si(r23418, MPFR_RNDN)) { mpfr_set(r23425, r23424, MPFR_RNDN); } else { mpfr_set(r23425, r23416, MPFR_RNDN); };
        if (mpfr_get_si(r23408, MPFR_RNDN)) { mpfr_set(r23426, r23416, MPFR_RNDN); } else { mpfr_set(r23426, r23425, MPFR_RNDN); };
        return mpfr_get_d(r23426, MPFR_RNDN);
}