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

char *name = "expq3 (problem 3.4.2)";

double f_if(float a, float b, float eps) {
        float r24375 = eps;
        float r24376 = a;
        float r24377 = b;
        float r24378 = r24376 + r24377;
        float r24379 = r24378 * r24375;
        float r24380 = exp(r24379);
        float r24381 = 1;
        float r24382 = r24380 - r24381;
        float r24383 = r24375 * r24382;
        float r24384 = r24376 * r24375;
        float r24385 = exp(r24384);
        float r24386 = r24385 - r24381;
        float r24387 = r24377 * r24375;
        float r24388 = exp(r24387);
        float r24389 = r24388 - r24381;
        float r24390 = r24386 * r24389;
        float r24391 = r24383 / r24390;
        return r24391;
}

double f_id(double a, double b, double eps) {
        double r24392 = eps;
        double r24393 = a;
        double r24394 = b;
        double r24395 = r24393 + r24394;
        double r24396 = r24395 * r24392;
        double r24397 = exp(r24396);
        double r24398 = 1;
        double r24399 = r24397 - r24398;
        double r24400 = r24392 * r24399;
        double r24401 = r24393 * r24392;
        double r24402 = exp(r24401);
        double r24403 = r24402 - r24398;
        double r24404 = r24394 * r24392;
        double r24405 = exp(r24404);
        double r24406 = r24405 - r24398;
        double r24407 = r24403 * r24406;
        double r24408 = r24400 / r24407;
        return r24408;
}


double f_of(float a, float b, float eps) {
        float r24409 = 1;
        float r24410 = a;
        float r24411 = r24409 / r24410;
        float r24412 = b;
        float r24413 = r24409 / r24412;
        float r24414 = r24411 + r24413;
        float r24415 = -4.635212546899196e-100;
        bool r24416 = r24414 <= r24415;
        float r24417 = 1.7540289770605972e-197;
        bool r24418 = r24414 <= r24417;
        float r24419 = !r24418;
        bool r24420 = r24416 || r24419;
        float r24421 = eps;
        float r24422 = r24410 + r24412;
        float r24423 = r24422 * r24421;
        float r24424 = exp(r24423);
        float r24425 = r24424 - r24409;
        float r24426 = r24421 * r24425;
        float r24427 = r24421 * r24412;
        float r24428 = exp(r24427);
        float r24429 = r24428 - r24409;
        float r24430 = r24421 * r24410;
        float r24431 = exp(r24430);
        float r24432 = r24431 - r24409;
        float r24433 = r24429 * r24432;
        float r24434 = r24426 / r24433;
        float r24435 = r24420 ? r24414 : r24434;
        return r24435;
}

double f_od(double a, double b, double eps) {
        double r24436 = 1;
        double r24437 = a;
        double r24438 = r24436 / r24437;
        double r24439 = b;
        double r24440 = r24436 / r24439;
        double r24441 = r24438 + r24440;
        double r24442 = -4.635212546899196e-100;
        bool r24443 = r24441 <= r24442;
        double r24444 = 1.7540289770605972e-197;
        bool r24445 = r24441 <= r24444;
        double r24446 = !r24445;
        bool r24447 = r24443 || r24446;
        double r24448 = eps;
        double r24449 = r24437 + r24439;
        double r24450 = r24449 * r24448;
        double r24451 = exp(r24450);
        double r24452 = r24451 - r24436;
        double r24453 = r24448 * r24452;
        double r24454 = r24448 * r24439;
        double r24455 = exp(r24454);
        double r24456 = r24455 - r24436;
        double r24457 = r24448 * r24437;
        double r24458 = exp(r24457);
        double r24459 = r24458 - r24436;
        double r24460 = r24456 * r24459;
        double r24461 = r24453 / r24460;
        double r24462 = r24447 ? r24441 : r24461;
        return r24462;
}

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 r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471, r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r24463);
        mpfr_init(r24464);
        mpfr_init(r24465);
        mpfr_init(r24466);
        mpfr_init(r24467);
        mpfr_init(r24468);
        mpfr_init_set_str(r24469, "1", 10, MPFR_RNDN);
        mpfr_init(r24470);
        mpfr_init(r24471);
        mpfr_init(r24472);
        mpfr_init(r24473);
        mpfr_init(r24474);
        mpfr_init(r24475);
        mpfr_init(r24476);
        mpfr_init(r24477);
        mpfr_init(r24478);
        mpfr_init(r24479);
}

double f_im(double a, double b, double eps) {
        mpfr_set_d(r24463, eps, MPFR_RNDN);
        mpfr_set_d(r24464, a, MPFR_RNDN);
        mpfr_set_d(r24465, b, MPFR_RNDN);
        mpfr_add(r24466, r24464, r24465, MPFR_RNDN);
        mpfr_mul(r24467, r24466, r24463, MPFR_RNDN);
        mpfr_exp(r24468, r24467, MPFR_RNDN);
        ;
        mpfr_sub(r24470, r24468, r24469, MPFR_RNDN);
        mpfr_mul(r24471, r24463, r24470, MPFR_RNDN);
        mpfr_mul(r24472, r24464, r24463, MPFR_RNDN);
        mpfr_exp(r24473, r24472, MPFR_RNDN);
        mpfr_sub(r24474, r24473, r24469, MPFR_RNDN);
        mpfr_mul(r24475, r24465, r24463, MPFR_RNDN);
        mpfr_exp(r24476, r24475, MPFR_RNDN);
        mpfr_sub(r24477, r24476, r24469, MPFR_RNDN);
        mpfr_mul(r24478, r24474, r24477, MPFR_RNDN);
        mpfr_div(r24479, r24471, r24478, MPFR_RNDN);
        return mpfr_get_d(r24479, MPFR_RNDN);
}

static mpfr_t r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487, r24488, r24489, r24490, r24491, r24492, r24493, r24494, r24495, r24496, r24497, r24498, r24499, r24500, r24501, r24502, r24503, r24504, r24505, r24506;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r24480, "1", 10, MPFR_RNDN);
        mpfr_init(r24481);
        mpfr_init(r24482);
        mpfr_init(r24483);
        mpfr_init(r24484);
        mpfr_init(r24485);
        mpfr_init_set_str(r24486, "-4.635212546899196e-100", 10, MPFR_RNDN);
        mpfr_init(r24487);
        mpfr_init_set_str(r24488, "1.7540289770605972e-197", 10, MPFR_RNDN);
        mpfr_init(r24489);
        mpfr_init(r24490);
        mpfr_init(r24491);
        mpfr_init(r24492);
        mpfr_init(r24493);
        mpfr_init(r24494);
        mpfr_init(r24495);
        mpfr_init(r24496);
        mpfr_init(r24497);
        mpfr_init(r24498);
        mpfr_init(r24499);
        mpfr_init(r24500);
        mpfr_init(r24501);
        mpfr_init(r24502);
        mpfr_init(r24503);
        mpfr_init(r24504);
        mpfr_init(r24505);
        mpfr_init(r24506);
}

double f_fm(double a, double b, double eps) {
        ;
        mpfr_set_d(r24481, a, MPFR_RNDN);
        mpfr_div(r24482, r24480, r24481, MPFR_RNDN);
        mpfr_set_d(r24483, b, MPFR_RNDN);
        mpfr_div(r24484, r24480, r24483, MPFR_RNDN);
        mpfr_add(r24485, r24482, r24484, MPFR_RNDN);
        ;
        mpfr_set_si(r24487, mpfr_cmp(r24485, r24486) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24489, mpfr_cmp(r24485, r24488) <= 0, MPFR_RNDN);
        mpfr_set_si(r24490, !mpfr_get_si(r24489, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r24491, mpfr_get_si(r24487, MPFR_RNDN) || mpfr_get_si(r24490, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r24492, eps, MPFR_RNDN);
        mpfr_add(r24493, r24481, r24483, MPFR_RNDN);
        mpfr_mul(r24494, r24493, r24492, MPFR_RNDN);
        mpfr_exp(r24495, r24494, MPFR_RNDN);
        mpfr_sub(r24496, r24495, r24480, MPFR_RNDN);
        mpfr_mul(r24497, r24492, r24496, MPFR_RNDN);
        mpfr_mul(r24498, r24492, r24483, MPFR_RNDN);
        mpfr_exp(r24499, r24498, MPFR_RNDN);
        mpfr_sub(r24500, r24499, r24480, MPFR_RNDN);
        mpfr_mul(r24501, r24492, r24481, MPFR_RNDN);
        mpfr_exp(r24502, r24501, MPFR_RNDN);
        mpfr_sub(r24503, r24502, r24480, MPFR_RNDN);
        mpfr_mul(r24504, r24500, r24503, MPFR_RNDN);
        mpfr_div(r24505, r24497, r24504, MPFR_RNDN);
        if (mpfr_get_si(r24491, MPFR_RNDN)) { mpfr_set(r24506, r24485, MPFR_RNDN); } else { mpfr_set(r24506, r24505, MPFR_RNDN); };
        return mpfr_get_d(r24506, MPFR_RNDN);
}

static mpfr_t r24507, r24508, r24509, r24510, r24511, r24512, r24513, r24514, r24515, r24516, r24517, r24518, r24519, r24520, r24521, r24522, r24523, r24524, r24525, r24526, r24527, r24528, r24529, r24530, r24531, r24532, r24533;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r24507, "1", 10, MPFR_RNDN);
        mpfr_init(r24508);
        mpfr_init(r24509);
        mpfr_init(r24510);
        mpfr_init(r24511);
        mpfr_init(r24512);
        mpfr_init_set_str(r24513, "-4.635212546899196e-100", 10, MPFR_RNDN);
        mpfr_init(r24514);
        mpfr_init_set_str(r24515, "1.7540289770605972e-197", 10, MPFR_RNDN);
        mpfr_init(r24516);
        mpfr_init(r24517);
        mpfr_init(r24518);
        mpfr_init(r24519);
        mpfr_init(r24520);
        mpfr_init(r24521);
        mpfr_init(r24522);
        mpfr_init(r24523);
        mpfr_init(r24524);
        mpfr_init(r24525);
        mpfr_init(r24526);
        mpfr_init(r24527);
        mpfr_init(r24528);
        mpfr_init(r24529);
        mpfr_init(r24530);
        mpfr_init(r24531);
        mpfr_init(r24532);
        mpfr_init(r24533);
}

double f_dm(double a, double b, double eps) {
        ;
        mpfr_set_d(r24508, a, MPFR_RNDN);
        mpfr_div(r24509, r24507, r24508, MPFR_RNDN);
        mpfr_set_d(r24510, b, MPFR_RNDN);
        mpfr_div(r24511, r24507, r24510, MPFR_RNDN);
        mpfr_add(r24512, r24509, r24511, MPFR_RNDN);
        ;
        mpfr_set_si(r24514, mpfr_cmp(r24512, r24513) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24516, mpfr_cmp(r24512, r24515) <= 0, MPFR_RNDN);
        mpfr_set_si(r24517, !mpfr_get_si(r24516, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r24518, mpfr_get_si(r24514, MPFR_RNDN) || mpfr_get_si(r24517, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r24519, eps, MPFR_RNDN);
        mpfr_add(r24520, r24508, r24510, MPFR_RNDN);
        mpfr_mul(r24521, r24520, r24519, MPFR_RNDN);
        mpfr_exp(r24522, r24521, MPFR_RNDN);
        mpfr_sub(r24523, r24522, r24507, MPFR_RNDN);
        mpfr_mul(r24524, r24519, r24523, MPFR_RNDN);
        mpfr_mul(r24525, r24519, r24510, MPFR_RNDN);
        mpfr_exp(r24526, r24525, MPFR_RNDN);
        mpfr_sub(r24527, r24526, r24507, MPFR_RNDN);
        mpfr_mul(r24528, r24519, r24508, MPFR_RNDN);
        mpfr_exp(r24529, r24528, MPFR_RNDN);
        mpfr_sub(r24530, r24529, r24507, MPFR_RNDN);
        mpfr_mul(r24531, r24527, r24530, MPFR_RNDN);
        mpfr_div(r24532, r24524, r24531, MPFR_RNDN);
        if (mpfr_get_si(r24518, MPFR_RNDN)) { mpfr_set(r24533, r24512, MPFR_RNDN); } else { mpfr_set(r24533, r24532, MPFR_RNDN); };
        return mpfr_get_d(r24533, MPFR_RNDN);
}