#include #include #include #include #include char *name = "expq3 (problem 3.4.2)"; double f_if(float a, float b, float eps) { float r26387 = eps; float r26388 = a; float r26389 = b; float r26390 = r26388 + r26389; float r26391 = r26390 * r26387; float r26392 = exp(r26391); float r26393 = 1; float r26394 = r26392 - r26393; float r26395 = r26387 * r26394; float r26396 = r26388 * r26387; float r26397 = exp(r26396); float r26398 = r26397 - r26393; float r26399 = r26389 * r26387; float r26400 = exp(r26399); float r26401 = r26400 - r26393; float r26402 = r26398 * r26401; float r26403 = r26395 / r26402; return r26403; } double f_id(double a, double b, double eps) { double r26404 = eps; double r26405 = a; double r26406 = b; double r26407 = r26405 + r26406; double r26408 = r26407 * r26404; double r26409 = exp(r26408); double r26410 = 1; double r26411 = r26409 - r26410; double r26412 = r26404 * r26411; double r26413 = r26405 * r26404; double r26414 = exp(r26413); double r26415 = r26414 - r26410; double r26416 = r26406 * r26404; double r26417 = exp(r26416); double r26418 = r26417 - r26410; double r26419 = r26415 * r26418; double r26420 = r26412 / r26419; return r26420; } double f_of(float a, float b, float eps) { float r26421 = eps; float r26422 = a; float r26423 = b; float r26424 = r26422 + r26423; float r26425 = r26424 * r26421; float r26426 = exp(r26425); float r26427 = 1; float r26428 = r26426 - r26427; float r26429 = r26421 * r26428; float r26430 = r26422 * r26421; float r26431 = exp(r26430); float r26432 = r26431 - r26427; float r26433 = r26423 * r26421; float r26434 = exp(r26433); float r26435 = r26434 - r26427; float r26436 = r26432 * r26435; float r26437 = r26429 / r26436; float r26438 = -8.40239431222467e-53; bool r26439 = r26437 <= r26438; float r26440 = r26427 / r26423; float r26441 = r26427 / r26422; float r26442 = r26440 + r26441; float r26443 = 1.7821478557518553e+308; bool r26444 = r26437 <= r26443; float r26445 = r26421 * r26424; float r26446 = exp(r26445); float r26447 = -r26421; float r26448 = fma(r26421, r26446, r26447); float r26449 = r26421 * r26423; float r26450 = expm1(r26449); float r26451 = expm1(r26430); float r26452 = r26450 * r26451; float r26453 = r26448 / r26452; float r26454 = r26444 ? r26453 : r26442; float r26455 = r26439 ? r26442 : r26454; return r26455; } double f_od(double a, double b, double eps) { double r26456 = eps; double r26457 = a; double r26458 = b; double r26459 = r26457 + r26458; double r26460 = r26459 * r26456; double r26461 = exp(r26460); double r26462 = 1; double r26463 = r26461 - r26462; double r26464 = r26456 * r26463; double r26465 = r26457 * r26456; double r26466 = exp(r26465); double r26467 = r26466 - r26462; double r26468 = r26458 * r26456; double r26469 = exp(r26468); double r26470 = r26469 - r26462; double r26471 = r26467 * r26470; double r26472 = r26464 / r26471; double r26473 = -8.40239431222467e-53; bool r26474 = r26472 <= r26473; double r26475 = r26462 / r26458; double r26476 = r26462 / r26457; double r26477 = r26475 + r26476; double r26478 = 1.7821478557518553e+308; bool r26479 = r26472 <= r26478; double r26480 = r26456 * r26459; double r26481 = exp(r26480); double r26482 = -r26456; double r26483 = fma(r26456, r26481, r26482); double r26484 = r26456 * r26458; double r26485 = expm1(r26484); double r26486 = expm1(r26465); double r26487 = r26485 * r26486; double r26488 = r26483 / r26487; double r26489 = r26479 ? r26488 : r26477; double r26490 = r26474 ? r26477 : r26489; return r26490; } 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 r26491, r26492, r26493, r26494, r26495, r26496, r26497, r26498, r26499, r26500, r26501, r26502, r26503, r26504, r26505, r26506, r26507; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r26491); mpfr_init(r26492); mpfr_init(r26493); mpfr_init(r26494); mpfr_init(r26495); mpfr_init(r26496); mpfr_init_set_str(r26497, "1", 10, MPFR_RNDN); mpfr_init(r26498); mpfr_init(r26499); mpfr_init(r26500); mpfr_init(r26501); mpfr_init(r26502); mpfr_init(r26503); mpfr_init(r26504); mpfr_init(r26505); mpfr_init(r26506); mpfr_init(r26507); } double f_im(double a, double b, double eps) { mpfr_set_d(r26491, eps, MPFR_RNDN); mpfr_set_d(r26492, a, MPFR_RNDN); mpfr_set_d(r26493, b, MPFR_RNDN); mpfr_add(r26494, r26492, r26493, MPFR_RNDN); mpfr_mul(r26495, r26494, r26491, MPFR_RNDN); mpfr_exp(r26496, r26495, MPFR_RNDN); ; mpfr_sub(r26498, r26496, r26497, MPFR_RNDN); mpfr_mul(r26499, r26491, r26498, MPFR_RNDN); mpfr_mul(r26500, r26492, r26491, MPFR_RNDN); mpfr_exp(r26501, r26500, MPFR_RNDN); mpfr_sub(r26502, r26501, r26497, MPFR_RNDN); mpfr_mul(r26503, r26493, r26491, MPFR_RNDN); mpfr_exp(r26504, r26503, MPFR_RNDN); mpfr_sub(r26505, r26504, r26497, MPFR_RNDN); mpfr_mul(r26506, r26502, r26505, MPFR_RNDN); mpfr_div(r26507, r26499, r26506, MPFR_RNDN); return mpfr_get_d(r26507, MPFR_RNDN); } static mpfr_t r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519, r26520, r26521, r26522, r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r26508); mpfr_init(r26509); mpfr_init(r26510); mpfr_init(r26511); mpfr_init(r26512); mpfr_init(r26513); mpfr_init_set_str(r26514, "1", 10, MPFR_RNDN); mpfr_init(r26515); mpfr_init(r26516); mpfr_init(r26517); mpfr_init(r26518); mpfr_init(r26519); mpfr_init(r26520); mpfr_init(r26521); mpfr_init(r26522); mpfr_init(r26523); mpfr_init(r26524); mpfr_init_set_str(r26525, "-8.40239431222467e-53", 10, MPFR_RNDN); mpfr_init(r26526); mpfr_init(r26527); mpfr_init(r26528); mpfr_init(r26529); mpfr_init_set_str(r26530, "1.7821478557518553e+308", 10, MPFR_RNDN); mpfr_init(r26531); mpfr_init(r26532); mpfr_init(r26533); mpfr_init(r26534); mpfr_init(r26535); mpfr_init(r26536); mpfr_init(r26537); mpfr_init(r26538); mpfr_init(r26539); mpfr_init(r26540); mpfr_init(r26541); mpfr_init(r26542); } double f_fm(double a, double b, double eps) { mpfr_set_d(r26508, eps, MPFR_RNDN); mpfr_set_d(r26509, a, MPFR_RNDN); mpfr_set_d(r26510, b, MPFR_RNDN); mpfr_add(r26511, r26509, r26510, MPFR_RNDN); mpfr_mul(r26512, r26511, r26508, MPFR_RNDN); mpfr_exp(r26513, r26512, MPFR_RNDN); ; mpfr_sub(r26515, r26513, r26514, MPFR_RNDN); mpfr_mul(r26516, r26508, r26515, MPFR_RNDN); mpfr_mul(r26517, r26509, r26508, MPFR_RNDN); mpfr_exp(r26518, r26517, MPFR_RNDN); mpfr_sub(r26519, r26518, r26514, MPFR_RNDN); mpfr_mul(r26520, r26510, r26508, MPFR_RNDN); mpfr_exp(r26521, r26520, MPFR_RNDN); mpfr_sub(r26522, r26521, r26514, MPFR_RNDN); mpfr_mul(r26523, r26519, r26522, MPFR_RNDN); mpfr_div(r26524, r26516, r26523, MPFR_RNDN); ; mpfr_set_si(r26526, mpfr_cmp(r26524, r26525) <= 0, MPFR_RNDN); mpfr_div(r26527, r26514, r26510, MPFR_RNDN); mpfr_div(r26528, r26514, r26509, MPFR_RNDN); mpfr_add(r26529, r26527, r26528, MPFR_RNDN); ; mpfr_set_si(r26531, mpfr_cmp(r26524, r26530) <= 0, MPFR_RNDN); mpfr_mul(r26532, r26508, r26511, MPFR_RNDN); mpfr_exp(r26533, r26532, MPFR_RNDN); mpfr_neg(r26534, r26508, MPFR_RNDN); mpfr_fma(r26535, r26508, r26533, r26534, MPFR_RNDN); mpfr_mul(r26536, r26508, r26510, MPFR_RNDN); mpfr_expm1(r26537, r26536, MPFR_RNDN); mpfr_expm1(r26538, r26517, MPFR_RNDN); mpfr_mul(r26539, r26537, r26538, MPFR_RNDN); mpfr_div(r26540, r26535, r26539, MPFR_RNDN); if (mpfr_get_si(r26531, MPFR_RNDN)) { mpfr_set(r26541, r26540, MPFR_RNDN); } else { mpfr_set(r26541, r26529, MPFR_RNDN); }; if (mpfr_get_si(r26526, MPFR_RNDN)) { mpfr_set(r26542, r26529, MPFR_RNDN); } else { mpfr_set(r26542, r26541, MPFR_RNDN); }; return mpfr_get_d(r26542, MPFR_RNDN); } static mpfr_t r26543, r26544, r26545, r26546, r26547, r26548, r26549, r26550, r26551, r26552, r26553, r26554, r26555, r26556, r26557, r26558, r26559, r26560, r26561, r26562, r26563, r26564, r26565, r26566, r26567, r26568, r26569, r26570, r26571, r26572, r26573, r26574, r26575, r26576, r26577; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r26543); mpfr_init(r26544); mpfr_init(r26545); mpfr_init(r26546); mpfr_init(r26547); mpfr_init(r26548); mpfr_init_set_str(r26549, "1", 10, MPFR_RNDN); mpfr_init(r26550); mpfr_init(r26551); mpfr_init(r26552); mpfr_init(r26553); mpfr_init(r26554); mpfr_init(r26555); mpfr_init(r26556); mpfr_init(r26557); mpfr_init(r26558); mpfr_init(r26559); mpfr_init_set_str(r26560, "-8.40239431222467e-53", 10, MPFR_RNDN); mpfr_init(r26561); mpfr_init(r26562); mpfr_init(r26563); mpfr_init(r26564); mpfr_init_set_str(r26565, "1.7821478557518553e+308", 10, MPFR_RNDN); mpfr_init(r26566); mpfr_init(r26567); mpfr_init(r26568); mpfr_init(r26569); mpfr_init(r26570); mpfr_init(r26571); mpfr_init(r26572); mpfr_init(r26573); mpfr_init(r26574); mpfr_init(r26575); mpfr_init(r26576); mpfr_init(r26577); } double f_dm(double a, double b, double eps) { mpfr_set_d(r26543, eps, MPFR_RNDN); mpfr_set_d(r26544, a, MPFR_RNDN); mpfr_set_d(r26545, b, MPFR_RNDN); mpfr_add(r26546, r26544, r26545, MPFR_RNDN); mpfr_mul(r26547, r26546, r26543, MPFR_RNDN); mpfr_exp(r26548, r26547, MPFR_RNDN); ; mpfr_sub(r26550, r26548, r26549, MPFR_RNDN); mpfr_mul(r26551, r26543, r26550, MPFR_RNDN); mpfr_mul(r26552, r26544, r26543, MPFR_RNDN); mpfr_exp(r26553, r26552, MPFR_RNDN); mpfr_sub(r26554, r26553, r26549, MPFR_RNDN); mpfr_mul(r26555, r26545, r26543, MPFR_RNDN); mpfr_exp(r26556, r26555, MPFR_RNDN); mpfr_sub(r26557, r26556, r26549, MPFR_RNDN); mpfr_mul(r26558, r26554, r26557, MPFR_RNDN); mpfr_div(r26559, r26551, r26558, MPFR_RNDN); ; mpfr_set_si(r26561, mpfr_cmp(r26559, r26560) <= 0, MPFR_RNDN); mpfr_div(r26562, r26549, r26545, MPFR_RNDN); mpfr_div(r26563, r26549, r26544, MPFR_RNDN); mpfr_add(r26564, r26562, r26563, MPFR_RNDN); ; mpfr_set_si(r26566, mpfr_cmp(r26559, r26565) <= 0, MPFR_RNDN); mpfr_mul(r26567, r26543, r26546, MPFR_RNDN); mpfr_exp(r26568, r26567, MPFR_RNDN); mpfr_neg(r26569, r26543, MPFR_RNDN); mpfr_fma(r26570, r26543, r26568, r26569, MPFR_RNDN); mpfr_mul(r26571, r26543, r26545, MPFR_RNDN); mpfr_expm1(r26572, r26571, MPFR_RNDN); mpfr_expm1(r26573, r26552, MPFR_RNDN); mpfr_mul(r26574, r26572, r26573, MPFR_RNDN); mpfr_div(r26575, r26570, r26574, MPFR_RNDN); if (mpfr_get_si(r26566, MPFR_RNDN)) { mpfr_set(r26576, r26575, MPFR_RNDN); } else { mpfr_set(r26576, r26564, MPFR_RNDN); }; if (mpfr_get_si(r26561, MPFR_RNDN)) { mpfr_set(r26577, r26564, MPFR_RNDN); } else { mpfr_set(r26577, r26576, MPFR_RNDN); }; return mpfr_get_d(r26577, MPFR_RNDN); }