#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 r22268 = eps;
        float r22269 = a;
        float r22270 = b;
        float r22271 = r22269 + r22270;
        float r22272 = r22271 * r22268;
        float r22273 = exp(r22272);
        float r22274 = 1;
        float r22275 = r22273 - r22274;
        float r22276 = r22268 * r22275;
        float r22277 = r22269 * r22268;
        float r22278 = exp(r22277);
        float r22279 = r22278 - r22274;
        float r22280 = r22270 * r22268;
        float r22281 = exp(r22280);
        float r22282 = r22281 - r22274;
        float r22283 = r22279 * r22282;
        float r22284 = r22276 / r22283;
        return r22284;
}

double f_id(double a, double b, double eps) {
        double r22285 = eps;
        double r22286 = a;
        double r22287 = b;
        double r22288 = r22286 + r22287;
        double r22289 = r22288 * r22285;
        double r22290 = exp(r22289);
        double r22291 = 1;
        double r22292 = r22290 - r22291;
        double r22293 = r22285 * r22292;
        double r22294 = r22286 * r22285;
        double r22295 = exp(r22294);
        double r22296 = r22295 - r22291;
        double r22297 = r22287 * r22285;
        double r22298 = exp(r22297);
        double r22299 = r22298 - r22291;
        double r22300 = r22296 * r22299;
        double r22301 = r22293 / r22300;
        return r22301;
}


double f_of(float a, float b, float eps) {
        float r22302 = 1;
        float r22303 = b;
        float r22304 = r22302 / r22303;
        float r22305 = cbrt(r22304);
        float r22306 = r22305 * r22305;
        float r22307 = a;
        float r22308 = r22302 / r22307;
        float r22309 = fma(r22306, r22305, r22308);
        float r22310 = -2.1394519520469999e-97;
        bool r22311 = r22309 <= r22310;
        float r22312 = r22304 + r22308;
        float r22313 = 1.1557750579452375e-09;
        bool r22314 = r22309 <= r22313;
        float r22315 = r22307 + r22303;
        float r22316 = eps;
        float r22317 = r22315 * r22316;
        float r22318 = expm1(r22317);
        float r22319 = r22316 * r22303;
        float r22320 = expm1(r22319);
        float r22321 = r22318 / r22320;
        float r22322 = r22307 * r22316;
        float r22323 = expm1(r22322);
        float r22324 = r22302 / r22316;
        float r22325 = r22323 * r22324;
        float r22326 = r22321 / r22325;
        float r22327 = r22314 ? r22326 : r22312;
        float r22328 = r22311 ? r22312 : r22327;
        return r22328;
}

double f_od(double a, double b, double eps) {
        double r22329 = 1;
        double r22330 = b;
        double r22331 = r22329 / r22330;
        double r22332 = cbrt(r22331);
        double r22333 = r22332 * r22332;
        double r22334 = a;
        double r22335 = r22329 / r22334;
        double r22336 = fma(r22333, r22332, r22335);
        double r22337 = -2.1394519520469999e-97;
        bool r22338 = r22336 <= r22337;
        double r22339 = r22331 + r22335;
        double r22340 = 1.1557750579452375e-09;
        bool r22341 = r22336 <= r22340;
        double r22342 = r22334 + r22330;
        double r22343 = eps;
        double r22344 = r22342 * r22343;
        double r22345 = expm1(r22344);
        double r22346 = r22343 * r22330;
        double r22347 = expm1(r22346);
        double r22348 = r22345 / r22347;
        double r22349 = r22334 * r22343;
        double r22350 = expm1(r22349);
        double r22351 = r22329 / r22343;
        double r22352 = r22350 * r22351;
        double r22353 = r22348 / r22352;
        double r22354 = r22341 ? r22353 : r22339;
        double r22355 = r22338 ? r22339 : r22354;
        return r22355;
}

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 r22356, r22357, r22358, r22359, r22360, r22361, r22362, r22363, r22364, r22365, r22366, r22367, r22368, r22369, r22370, r22371, r22372;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init(r22356);
        mpfr_init(r22357);
        mpfr_init(r22358);
        mpfr_init(r22359);
        mpfr_init(r22360);
        mpfr_init(r22361);
        mpfr_init_set_str(r22362, "1", 10, MPFR_RNDN);
        mpfr_init(r22363);
        mpfr_init(r22364);
        mpfr_init(r22365);
        mpfr_init(r22366);
        mpfr_init(r22367);
        mpfr_init(r22368);
        mpfr_init(r22369);
        mpfr_init(r22370);
        mpfr_init(r22371);
        mpfr_init(r22372);
}

double f_im(double a, double b, double eps) {
        mpfr_set_d(r22356, eps, MPFR_RNDN);
        mpfr_set_d(r22357, a, MPFR_RNDN);
        mpfr_set_d(r22358, b, MPFR_RNDN);
        mpfr_add(r22359, r22357, r22358, MPFR_RNDN);
        mpfr_mul(r22360, r22359, r22356, MPFR_RNDN);
        mpfr_exp(r22361, r22360, MPFR_RNDN);
        ;
        mpfr_sub(r22363, r22361, r22362, MPFR_RNDN);
        mpfr_mul(r22364, r22356, r22363, MPFR_RNDN);
        mpfr_mul(r22365, r22357, r22356, MPFR_RNDN);
        mpfr_exp(r22366, r22365, MPFR_RNDN);
        mpfr_sub(r22367, r22366, r22362, MPFR_RNDN);
        mpfr_mul(r22368, r22358, r22356, MPFR_RNDN);
        mpfr_exp(r22369, r22368, MPFR_RNDN);
        mpfr_sub(r22370, r22369, r22362, MPFR_RNDN);
        mpfr_mul(r22371, r22367, r22370, MPFR_RNDN);
        mpfr_div(r22372, r22364, r22371, MPFR_RNDN);
        return mpfr_get_d(r22372, MPFR_RNDN);
}

static mpfr_t r22373, r22374, r22375, r22376, r22377, r22378, r22379, r22380, r22381, r22382, r22383, r22384, r22385, r22386, r22387, r22388, r22389, r22390, r22391, r22392, r22393, r22394, r22395, r22396, r22397, r22398, r22399;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r22373, "1", 10, MPFR_RNDN);
        mpfr_init(r22374);
        mpfr_init(r22375);
        mpfr_init(r22376);
        mpfr_init(r22377);
        mpfr_init(r22378);
        mpfr_init(r22379);
        mpfr_init(r22380);
        mpfr_init_set_str(r22381, "-2.1394519520469999e-97", 10, MPFR_RNDN);
        mpfr_init(r22382);
        mpfr_init(r22383);
        mpfr_init_set_str(r22384, "1.1557750579452375e-09", 10, MPFR_RNDN);
        mpfr_init(r22385);
        mpfr_init(r22386);
        mpfr_init(r22387);
        mpfr_init(r22388);
        mpfr_init(r22389);
        mpfr_init(r22390);
        mpfr_init(r22391);
        mpfr_init(r22392);
        mpfr_init(r22393);
        mpfr_init(r22394);
        mpfr_init(r22395);
        mpfr_init(r22396);
        mpfr_init(r22397);
        mpfr_init(r22398);
        mpfr_init(r22399);
}

double f_fm(double a, double b, double eps) {
        ;
        mpfr_set_d(r22374, b, MPFR_RNDN);
        mpfr_div(r22375, r22373, r22374, MPFR_RNDN);
        mpfr_cbrt(r22376, r22375, MPFR_RNDN);
        mpfr_mul(r22377, r22376, r22376, MPFR_RNDN);
        mpfr_set_d(r22378, a, MPFR_RNDN);
        mpfr_div(r22379, r22373, r22378, MPFR_RNDN);
        mpfr_fma(r22380, r22377, r22376, r22379, MPFR_RNDN);
        ;
        mpfr_set_si(r22382, mpfr_cmp(r22380, r22381) <= 0, MPFR_RNDN);
        mpfr_add(r22383, r22375, r22379, MPFR_RNDN);
        ;
        mpfr_set_si(r22385, mpfr_cmp(r22380, r22384) <= 0, MPFR_RNDN);
        mpfr_add(r22386, r22378, r22374, MPFR_RNDN);
        mpfr_set_d(r22387, eps, MPFR_RNDN);
        mpfr_mul(r22388, r22386, r22387, MPFR_RNDN);
        mpfr_expm1(r22389, r22388, MPFR_RNDN);
        mpfr_mul(r22390, r22387, r22374, MPFR_RNDN);
        mpfr_expm1(r22391, r22390, MPFR_RNDN);
        mpfr_div(r22392, r22389, r22391, MPFR_RNDN);
        mpfr_mul(r22393, r22378, r22387, MPFR_RNDN);
        mpfr_expm1(r22394, r22393, MPFR_RNDN);
        mpfr_div(r22395, r22373, r22387, MPFR_RNDN);
        mpfr_mul(r22396, r22394, r22395, MPFR_RNDN);
        mpfr_div(r22397, r22392, r22396, MPFR_RNDN);
        if (mpfr_get_si(r22385, MPFR_RNDN)) { mpfr_set(r22398, r22397, MPFR_RNDN); } else { mpfr_set(r22398, r22383, MPFR_RNDN); };
        if (mpfr_get_si(r22382, MPFR_RNDN)) { mpfr_set(r22399, r22383, MPFR_RNDN); } else { mpfr_set(r22399, r22398, MPFR_RNDN); };
        return mpfr_get_d(r22399, MPFR_RNDN);
}

static mpfr_t r22400, r22401, r22402, r22403, r22404, r22405, r22406, r22407, r22408, r22409, r22410, r22411, r22412, r22413, r22414, r22415, r22416, r22417, r22418, r22419, r22420, r22421, r22422, r22423, r22424, r22425, r22426;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r22400, "1", 10, MPFR_RNDN);
        mpfr_init(r22401);
        mpfr_init(r22402);
        mpfr_init(r22403);
        mpfr_init(r22404);
        mpfr_init(r22405);
        mpfr_init(r22406);
        mpfr_init(r22407);
        mpfr_init_set_str(r22408, "-2.1394519520469999e-97", 10, MPFR_RNDN);
        mpfr_init(r22409);
        mpfr_init(r22410);
        mpfr_init_set_str(r22411, "1.1557750579452375e-09", 10, MPFR_RNDN);
        mpfr_init(r22412);
        mpfr_init(r22413);
        mpfr_init(r22414);
        mpfr_init(r22415);
        mpfr_init(r22416);
        mpfr_init(r22417);
        mpfr_init(r22418);
        mpfr_init(r22419);
        mpfr_init(r22420);
        mpfr_init(r22421);
        mpfr_init(r22422);
        mpfr_init(r22423);
        mpfr_init(r22424);
        mpfr_init(r22425);
        mpfr_init(r22426);
}

double f_dm(double a, double b, double eps) {
        ;
        mpfr_set_d(r22401, b, MPFR_RNDN);
        mpfr_div(r22402, r22400, r22401, MPFR_RNDN);
        mpfr_cbrt(r22403, r22402, MPFR_RNDN);
        mpfr_mul(r22404, r22403, r22403, MPFR_RNDN);
        mpfr_set_d(r22405, a, MPFR_RNDN);
        mpfr_div(r22406, r22400, r22405, MPFR_RNDN);
        mpfr_fma(r22407, r22404, r22403, r22406, MPFR_RNDN);
        ;
        mpfr_set_si(r22409, mpfr_cmp(r22407, r22408) <= 0, MPFR_RNDN);
        mpfr_add(r22410, r22402, r22406, MPFR_RNDN);
        ;
        mpfr_set_si(r22412, mpfr_cmp(r22407, r22411) <= 0, MPFR_RNDN);
        mpfr_add(r22413, r22405, r22401, MPFR_RNDN);
        mpfr_set_d(r22414, eps, MPFR_RNDN);
        mpfr_mul(r22415, r22413, r22414, MPFR_RNDN);
        mpfr_expm1(r22416, r22415, MPFR_RNDN);
        mpfr_mul(r22417, r22414, r22401, MPFR_RNDN);
        mpfr_expm1(r22418, r22417, MPFR_RNDN);
        mpfr_div(r22419, r22416, r22418, MPFR_RNDN);
        mpfr_mul(r22420, r22405, r22414, MPFR_RNDN);
        mpfr_expm1(r22421, r22420, MPFR_RNDN);
        mpfr_div(r22422, r22400, r22414, MPFR_RNDN);
        mpfr_mul(r22423, r22421, r22422, MPFR_RNDN);
        mpfr_div(r22424, r22419, r22423, MPFR_RNDN);
        if (mpfr_get_si(r22412, MPFR_RNDN)) { mpfr_set(r22425, r22424, MPFR_RNDN); } else { mpfr_set(r22425, r22410, MPFR_RNDN); };
        if (mpfr_get_si(r22409, MPFR_RNDN)) { mpfr_set(r22426, r22410, MPFR_RNDN); } else { mpfr_set(r22426, r22425, MPFR_RNDN); };
        return mpfr_get_d(r22426, MPFR_RNDN);
}