#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 r24475 = eps;
        float r24476 = a;
        float r24477 = b;
        float r24478 = r24476 + r24477;
        float r24479 = r24478 * r24475;
        float r24480 = exp(r24479);
        float r24481 = 1;
        float r24482 = r24480 - r24481;
        float r24483 = r24475 * r24482;
        float r24484 = r24476 * r24475;
        float r24485 = exp(r24484);
        float r24486 = r24485 - r24481;
        float r24487 = r24477 * r24475;
        float r24488 = exp(r24487);
        float r24489 = r24488 - r24481;
        float r24490 = r24486 * r24489;
        float r24491 = r24483 / r24490;
        return r24491;
}

double f_id(double a, double b, double eps) {
        double r24492 = eps;
        double r24493 = a;
        double r24494 = b;
        double r24495 = r24493 + r24494;
        double r24496 = r24495 * r24492;
        double r24497 = exp(r24496);
        double r24498 = 1;
        double r24499 = r24497 - r24498;
        double r24500 = r24492 * r24499;
        double r24501 = r24493 * r24492;
        double r24502 = exp(r24501);
        double r24503 = r24502 - r24498;
        double r24504 = r24494 * r24492;
        double r24505 = exp(r24504);
        double r24506 = r24505 - r24498;
        double r24507 = r24503 * r24506;
        double r24508 = r24500 / r24507;
        return r24508;
}


double f_of(float a, float b, float eps) {
        float r24509 = 1;
        float r24510 = a;
        float r24511 = r24509 / r24510;
        float r24512 = b;
        float r24513 = r24509 / r24512;
        float r24514 = r24511 + r24513;
        float r24515 = -4.549790872176706e-158;
        bool r24516 = r24514 <= r24515;
        float r24517 = 2.201452965804001e-159;
        bool r24518 = r24514 <= r24517;
        float r24519 = !r24518;
        bool r24520 = r24516 || r24519;
        float r24521 = eps;
        float r24522 = r24510 + r24512;
        float r24523 = r24522 * r24521;
        float r24524 = exp(r24523);
        float r24525 = r24524 - r24509;
        float r24526 = r24521 * r24525;
        float r24527 = r24521 * r24512;
        float r24528 = exp(r24527);
        float r24529 = r24528 - r24509;
        float r24530 = exp(r24529);
        float r24531 = log(r24530);
        float r24532 = r24521 * r24510;
        float r24533 = exp(r24532);
        float r24534 = r24533 - r24509;
        float r24535 = r24531 * r24534;
        float r24536 = r24526 / r24535;
        float r24537 = r24520 ? r24514 : r24536;
        return r24537;
}

double f_od(double a, double b, double eps) {
        double r24538 = 1;
        double r24539 = a;
        double r24540 = r24538 / r24539;
        double r24541 = b;
        double r24542 = r24538 / r24541;
        double r24543 = r24540 + r24542;
        double r24544 = -4.549790872176706e-158;
        bool r24545 = r24543 <= r24544;
        double r24546 = 2.201452965804001e-159;
        bool r24547 = r24543 <= r24546;
        double r24548 = !r24547;
        bool r24549 = r24545 || r24548;
        double r24550 = eps;
        double r24551 = r24539 + r24541;
        double r24552 = r24551 * r24550;
        double r24553 = exp(r24552);
        double r24554 = r24553 - r24538;
        double r24555 = r24550 * r24554;
        double r24556 = r24550 * r24541;
        double r24557 = exp(r24556);
        double r24558 = r24557 - r24538;
        double r24559 = exp(r24558);
        double r24560 = log(r24559);
        double r24561 = r24550 * r24539;
        double r24562 = exp(r24561);
        double r24563 = r24562 - r24538;
        double r24564 = r24560 * r24563;
        double r24565 = r24555 / r24564;
        double r24566 = r24549 ? r24543 : r24565;
        return r24566;
}

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 r24567, r24568, r24569, r24570, r24571, r24572, r24573, r24574, r24575, r24576, r24577, r24578, r24579, r24580, r24581, r24582, r24583;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init(r24567);
        mpfr_init(r24568);
        mpfr_init(r24569);
        mpfr_init(r24570);
        mpfr_init(r24571);
        mpfr_init(r24572);
        mpfr_init_set_str(r24573, "1", 10, MPFR_RNDN);
        mpfr_init(r24574);
        mpfr_init(r24575);
        mpfr_init(r24576);
        mpfr_init(r24577);
        mpfr_init(r24578);
        mpfr_init(r24579);
        mpfr_init(r24580);
        mpfr_init(r24581);
        mpfr_init(r24582);
        mpfr_init(r24583);
}

double f_im(double a, double b, double eps) {
        mpfr_set_d(r24567, eps, MPFR_RNDN);
        mpfr_set_d(r24568, a, MPFR_RNDN);
        mpfr_set_d(r24569, b, MPFR_RNDN);
        mpfr_add(r24570, r24568, r24569, MPFR_RNDN);
        mpfr_mul(r24571, r24570, r24567, MPFR_RNDN);
        mpfr_exp(r24572, r24571, MPFR_RNDN);
        ;
        mpfr_sub(r24574, r24572, r24573, MPFR_RNDN);
        mpfr_mul(r24575, r24567, r24574, MPFR_RNDN);
        mpfr_mul(r24576, r24568, r24567, MPFR_RNDN);
        mpfr_exp(r24577, r24576, MPFR_RNDN);
        mpfr_sub(r24578, r24577, r24573, MPFR_RNDN);
        mpfr_mul(r24579, r24569, r24567, MPFR_RNDN);
        mpfr_exp(r24580, r24579, MPFR_RNDN);
        mpfr_sub(r24581, r24580, r24573, MPFR_RNDN);
        mpfr_mul(r24582, r24578, r24581, MPFR_RNDN);
        mpfr_div(r24583, r24575, r24582, MPFR_RNDN);
        return mpfr_get_d(r24583, MPFR_RNDN);
}

static mpfr_t r24584, r24585, r24586, r24587, r24588, r24589, r24590, r24591, r24592, r24593, r24594, r24595, r24596, r24597, r24598, r24599, r24600, r24601, r24602, r24603, r24604, r24605, r24606, r24607, r24608, r24609, r24610, r24611, r24612;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r24584, "1", 10, MPFR_RNDN);
        mpfr_init(r24585);
        mpfr_init(r24586);
        mpfr_init(r24587);
        mpfr_init(r24588);
        mpfr_init(r24589);
        mpfr_init_set_str(r24590, "-4.549790872176706e-158", 10, MPFR_RNDN);
        mpfr_init(r24591);
        mpfr_init_set_str(r24592, "2.201452965804001e-159", 10, MPFR_RNDN);
        mpfr_init(r24593);
        mpfr_init(r24594);
        mpfr_init(r24595);
        mpfr_init(r24596);
        mpfr_init(r24597);
        mpfr_init(r24598);
        mpfr_init(r24599);
        mpfr_init(r24600);
        mpfr_init(r24601);
        mpfr_init(r24602);
        mpfr_init(r24603);
        mpfr_init(r24604);
        mpfr_init(r24605);
        mpfr_init(r24606);
        mpfr_init(r24607);
        mpfr_init(r24608);
        mpfr_init(r24609);
        mpfr_init(r24610);
        mpfr_init(r24611);
        mpfr_init(r24612);
}

double f_fm(double a, double b, double eps) {
        ;
        mpfr_set_d(r24585, a, MPFR_RNDN);
        mpfr_div(r24586, r24584, r24585, MPFR_RNDN);
        mpfr_set_d(r24587, b, MPFR_RNDN);
        mpfr_div(r24588, r24584, r24587, MPFR_RNDN);
        mpfr_add(r24589, r24586, r24588, MPFR_RNDN);
        ;
        mpfr_set_si(r24591, mpfr_cmp(r24589, r24590) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24593, mpfr_cmp(r24589, r24592) <= 0, MPFR_RNDN);
        mpfr_set_si(r24594, !mpfr_get_si(r24593, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r24595, mpfr_get_si(r24591, MPFR_RNDN) || mpfr_get_si(r24594, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r24596, eps, MPFR_RNDN);
        mpfr_add(r24597, r24585, r24587, MPFR_RNDN);
        mpfr_mul(r24598, r24597, r24596, MPFR_RNDN);
        mpfr_exp(r24599, r24598, MPFR_RNDN);
        mpfr_sub(r24600, r24599, r24584, MPFR_RNDN);
        mpfr_mul(r24601, r24596, r24600, MPFR_RNDN);
        mpfr_mul(r24602, r24596, r24587, MPFR_RNDN);
        mpfr_exp(r24603, r24602, MPFR_RNDN);
        mpfr_sub(r24604, r24603, r24584, MPFR_RNDN);
        mpfr_exp(r24605, r24604, MPFR_RNDN);
        mpfr_log(r24606, r24605, MPFR_RNDN);
        mpfr_mul(r24607, r24596, r24585, MPFR_RNDN);
        mpfr_exp(r24608, r24607, MPFR_RNDN);
        mpfr_sub(r24609, r24608, r24584, MPFR_RNDN);
        mpfr_mul(r24610, r24606, r24609, MPFR_RNDN);
        mpfr_div(r24611, r24601, r24610, MPFR_RNDN);
        if (mpfr_get_si(r24595, MPFR_RNDN)) { mpfr_set(r24612, r24589, MPFR_RNDN); } else { mpfr_set(r24612, r24611, MPFR_RNDN); };
        return mpfr_get_d(r24612, MPFR_RNDN);
}

static mpfr_t r24613, r24614, r24615, r24616, r24617, r24618, r24619, r24620, r24621, r24622, r24623, r24624, r24625, r24626, r24627, r24628, r24629, r24630, r24631, r24632, r24633, r24634, r24635, r24636, r24637, r24638, r24639, r24640, r24641;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r24613, "1", 10, MPFR_RNDN);
        mpfr_init(r24614);
        mpfr_init(r24615);
        mpfr_init(r24616);
        mpfr_init(r24617);
        mpfr_init(r24618);
        mpfr_init_set_str(r24619, "-4.549790872176706e-158", 10, MPFR_RNDN);
        mpfr_init(r24620);
        mpfr_init_set_str(r24621, "2.201452965804001e-159", 10, MPFR_RNDN);
        mpfr_init(r24622);
        mpfr_init(r24623);
        mpfr_init(r24624);
        mpfr_init(r24625);
        mpfr_init(r24626);
        mpfr_init(r24627);
        mpfr_init(r24628);
        mpfr_init(r24629);
        mpfr_init(r24630);
        mpfr_init(r24631);
        mpfr_init(r24632);
        mpfr_init(r24633);
        mpfr_init(r24634);
        mpfr_init(r24635);
        mpfr_init(r24636);
        mpfr_init(r24637);
        mpfr_init(r24638);
        mpfr_init(r24639);
        mpfr_init(r24640);
        mpfr_init(r24641);
}

double f_dm(double a, double b, double eps) {
        ;
        mpfr_set_d(r24614, a, MPFR_RNDN);
        mpfr_div(r24615, r24613, r24614, MPFR_RNDN);
        mpfr_set_d(r24616, b, MPFR_RNDN);
        mpfr_div(r24617, r24613, r24616, MPFR_RNDN);
        mpfr_add(r24618, r24615, r24617, MPFR_RNDN);
        ;
        mpfr_set_si(r24620, mpfr_cmp(r24618, r24619) <= 0, MPFR_RNDN);
        ;
        mpfr_set_si(r24622, mpfr_cmp(r24618, r24621) <= 0, MPFR_RNDN);
        mpfr_set_si(r24623, !mpfr_get_si(r24622, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_si(r24624, mpfr_get_si(r24620, MPFR_RNDN) || mpfr_get_si(r24623, MPFR_RNDN), MPFR_RNDN);
        mpfr_set_d(r24625, eps, MPFR_RNDN);
        mpfr_add(r24626, r24614, r24616, MPFR_RNDN);
        mpfr_mul(r24627, r24626, r24625, MPFR_RNDN);
        mpfr_exp(r24628, r24627, MPFR_RNDN);
        mpfr_sub(r24629, r24628, r24613, MPFR_RNDN);
        mpfr_mul(r24630, r24625, r24629, MPFR_RNDN);
        mpfr_mul(r24631, r24625, r24616, MPFR_RNDN);
        mpfr_exp(r24632, r24631, MPFR_RNDN);
        mpfr_sub(r24633, r24632, r24613, MPFR_RNDN);
        mpfr_exp(r24634, r24633, MPFR_RNDN);
        mpfr_log(r24635, r24634, MPFR_RNDN);
        mpfr_mul(r24636, r24625, r24614, MPFR_RNDN);
        mpfr_exp(r24637, r24636, MPFR_RNDN);
        mpfr_sub(r24638, r24637, r24613, MPFR_RNDN);
        mpfr_mul(r24639, r24635, r24638, MPFR_RNDN);
        mpfr_div(r24640, r24630, r24639, MPFR_RNDN);
        if (mpfr_get_si(r24624, MPFR_RNDN)) { mpfr_set(r24641, r24618, MPFR_RNDN); } else { mpfr_set(r24641, r24640, MPFR_RNDN); };
        return mpfr_get_d(r24641, MPFR_RNDN);
}