#include #include #include #include #include char *name = "Harley's example"; double f_if(float c_p, float c_n, float t, float s) { float r11498 = 1; float r11499 = s; float r11500 = -r11499; float r11501 = exp(r11500); float r11502 = r11498 + r11501; float r11503 = r11498 / r11502; float r11504 = c_p; float r11505 = pow(r11503, r11504); float r11506 = r11498 - r11503; float r11507 = c_n; float r11508 = pow(r11506, r11507); float r11509 = r11505 * r11508; float r11510 = t; float r11511 = -r11510; float r11512 = exp(r11511); float r11513 = r11498 + r11512; float r11514 = r11498 / r11513; float r11515 = pow(r11514, r11504); float r11516 = r11498 - r11514; float r11517 = pow(r11516, r11507); float r11518 = r11515 * r11517; float r11519 = r11509 / r11518; return r11519; } double f_id(double c_p, double c_n, double t, double s) { double r11520 = 1; double r11521 = s; double r11522 = -r11521; double r11523 = exp(r11522); double r11524 = r11520 + r11523; double r11525 = r11520 / r11524; double r11526 = c_p; double r11527 = pow(r11525, r11526); double r11528 = r11520 - r11525; double r11529 = c_n; double r11530 = pow(r11528, r11529); double r11531 = r11527 * r11530; double r11532 = t; double r11533 = -r11532; double r11534 = exp(r11533); double r11535 = r11520 + r11534; double r11536 = r11520 / r11535; double r11537 = pow(r11536, r11526); double r11538 = r11520 - r11536; double r11539 = pow(r11538, r11529); double r11540 = r11537 * r11539; double r11541 = r11531 / r11540; return r11541; } double f_of(float c_p, float c_n, float t, float s) { float r11542 = 1; float r11543 = -r11542; float r11544 = s; float r11545 = -r11544; float r11546 = exp(r11545); float r11547 = r11546 + r11542; float r11548 = r11543 / r11547; float r11549 = log1p(r11548); float r11550 = t; float r11551 = -r11550; float r11552 = exp(r11551); float r11553 = r11552 + r11542; float r11554 = r11543 / r11553; float r11555 = log1p(r11554); float r11556 = r11549 - r11555; float r11557 = c_n; float r11558 = log1p(r11552); float r11559 = log1p(r11546); float r11560 = r11558 - r11559; float r11561 = c_p; float r11562 = r11560 * r11561; float r11563 = fma(r11556, r11557, r11562); float r11564 = exp(r11563); return r11564; } double f_od(double c_p, double c_n, double t, double s) { double r11565 = 1; double r11566 = -r11565; double r11567 = s; double r11568 = -r11567; double r11569 = exp(r11568); double r11570 = r11569 + r11565; double r11571 = r11566 / r11570; double r11572 = log1p(r11571); double r11573 = t; double r11574 = -r11573; double r11575 = exp(r11574); double r11576 = r11575 + r11565; double r11577 = r11566 / r11576; double r11578 = log1p(r11577); double r11579 = r11572 - r11578; double r11580 = c_n; double r11581 = log1p(r11575); double r11582 = log1p(r11569); double r11583 = r11581 - r11582; double r11584 = c_p; double r11585 = r11583 * r11584; double r11586 = fma(r11579, r11580, r11585); double r11587 = exp(r11586); return r11587; } 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 r11588, r11589, r11590, r11591, r11592, r11593, r11594, r11595, r11596, r11597, r11598, r11599, r11600, r11601, r11602, r11603, r11604, r11605, r11606, r11607, r11608, r11609; void setup_mpfr_f_im() { mpfr_set_default_prec(1872); mpfr_init_set_str(r11588, "1", 10, MPFR_RNDN); mpfr_init(r11589); mpfr_init(r11590); mpfr_init(r11591); mpfr_init(r11592); mpfr_init(r11593); mpfr_init(r11594); mpfr_init(r11595); mpfr_init(r11596); mpfr_init(r11597); mpfr_init(r11598); mpfr_init(r11599); mpfr_init(r11600); mpfr_init(r11601); mpfr_init(r11602); mpfr_init(r11603); mpfr_init(r11604); mpfr_init(r11605); mpfr_init(r11606); mpfr_init(r11607); mpfr_init(r11608); mpfr_init(r11609); } double f_im(double c_p, double c_n, double t, double s) { ; mpfr_set_d(r11589, s, MPFR_RNDN); mpfr_neg(r11590, r11589, MPFR_RNDN); mpfr_exp(r11591, r11590, MPFR_RNDN); mpfr_add(r11592, r11588, r11591, MPFR_RNDN); mpfr_div(r11593, r11588, r11592, MPFR_RNDN); mpfr_set_d(r11594, c_p, MPFR_RNDN); mpfr_pow(r11595, r11593, r11594, MPFR_RNDN); mpfr_sub(r11596, r11588, r11593, MPFR_RNDN); mpfr_set_d(r11597, c_n, MPFR_RNDN); mpfr_pow(r11598, r11596, r11597, MPFR_RNDN); mpfr_mul(r11599, r11595, r11598, MPFR_RNDN); mpfr_set_d(r11600, t, MPFR_RNDN); mpfr_neg(r11601, r11600, MPFR_RNDN); mpfr_exp(r11602, r11601, MPFR_RNDN); mpfr_add(r11603, r11588, r11602, MPFR_RNDN); mpfr_div(r11604, r11588, r11603, MPFR_RNDN); mpfr_pow(r11605, r11604, r11594, MPFR_RNDN); mpfr_sub(r11606, r11588, r11604, MPFR_RNDN); mpfr_pow(r11607, r11606, r11597, MPFR_RNDN); mpfr_mul(r11608, r11605, r11607, MPFR_RNDN); mpfr_div(r11609, r11599, r11608, MPFR_RNDN); return mpfr_get_d(r11609, MPFR_RNDN); } static mpfr_t r11610, r11611, r11612, r11613, r11614, r11615, r11616, r11617, r11618, r11619, r11620, r11621, r11622, r11623, r11624, r11625, r11626, r11627, r11628, r11629, r11630, r11631, r11632; void setup_mpfr_f_fm() { mpfr_set_default_prec(1872); mpfr_init_set_str(r11610, "1", 10, MPFR_RNDN); mpfr_init(r11611); mpfr_init(r11612); mpfr_init(r11613); mpfr_init(r11614); mpfr_init(r11615); mpfr_init(r11616); mpfr_init(r11617); mpfr_init(r11618); mpfr_init(r11619); mpfr_init(r11620); mpfr_init(r11621); mpfr_init(r11622); mpfr_init(r11623); mpfr_init(r11624); mpfr_init(r11625); mpfr_init(r11626); mpfr_init(r11627); mpfr_init(r11628); mpfr_init(r11629); mpfr_init(r11630); mpfr_init(r11631); mpfr_init(r11632); } double f_fm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r11611, r11610, MPFR_RNDN); mpfr_set_d(r11612, s, MPFR_RNDN); mpfr_neg(r11613, r11612, MPFR_RNDN); mpfr_exp(r11614, r11613, MPFR_RNDN); mpfr_add(r11615, r11614, r11610, MPFR_RNDN); mpfr_div(r11616, r11611, r11615, MPFR_RNDN); mpfr_log1p(r11617, r11616, MPFR_RNDN); mpfr_set_d(r11618, t, MPFR_RNDN); mpfr_neg(r11619, r11618, MPFR_RNDN); mpfr_exp(r11620, r11619, MPFR_RNDN); mpfr_add(r11621, r11620, r11610, MPFR_RNDN); mpfr_div(r11622, r11611, r11621, MPFR_RNDN); mpfr_log1p(r11623, r11622, MPFR_RNDN); mpfr_sub(r11624, r11617, r11623, MPFR_RNDN); mpfr_set_d(r11625, c_n, MPFR_RNDN); mpfr_log1p(r11626, r11620, MPFR_RNDN); mpfr_log1p(r11627, r11614, MPFR_RNDN); mpfr_sub(r11628, r11626, r11627, MPFR_RNDN); mpfr_set_d(r11629, c_p, MPFR_RNDN); mpfr_mul(r11630, r11628, r11629, MPFR_RNDN); mpfr_fma(r11631, r11624, r11625, r11630, MPFR_RNDN); mpfr_exp(r11632, r11631, MPFR_RNDN); return mpfr_get_d(r11632, MPFR_RNDN); } static mpfr_t r11633, r11634, r11635, r11636, r11637, r11638, r11639, r11640, r11641, r11642, r11643, r11644, r11645, r11646, r11647, r11648, r11649, r11650, r11651, r11652, r11653, r11654, r11655; void setup_mpfr_f_dm() { mpfr_set_default_prec(1872); mpfr_init_set_str(r11633, "1", 10, MPFR_RNDN); mpfr_init(r11634); mpfr_init(r11635); mpfr_init(r11636); mpfr_init(r11637); mpfr_init(r11638); mpfr_init(r11639); mpfr_init(r11640); mpfr_init(r11641); mpfr_init(r11642); mpfr_init(r11643); mpfr_init(r11644); mpfr_init(r11645); mpfr_init(r11646); mpfr_init(r11647); mpfr_init(r11648); mpfr_init(r11649); mpfr_init(r11650); mpfr_init(r11651); mpfr_init(r11652); mpfr_init(r11653); mpfr_init(r11654); mpfr_init(r11655); } double f_dm(double c_p, double c_n, double t, double s) { ; mpfr_neg(r11634, r11633, MPFR_RNDN); mpfr_set_d(r11635, s, MPFR_RNDN); mpfr_neg(r11636, r11635, MPFR_RNDN); mpfr_exp(r11637, r11636, MPFR_RNDN); mpfr_add(r11638, r11637, r11633, MPFR_RNDN); mpfr_div(r11639, r11634, r11638, MPFR_RNDN); mpfr_log1p(r11640, r11639, MPFR_RNDN); mpfr_set_d(r11641, t, MPFR_RNDN); mpfr_neg(r11642, r11641, MPFR_RNDN); mpfr_exp(r11643, r11642, MPFR_RNDN); mpfr_add(r11644, r11643, r11633, MPFR_RNDN); mpfr_div(r11645, r11634, r11644, MPFR_RNDN); mpfr_log1p(r11646, r11645, MPFR_RNDN); mpfr_sub(r11647, r11640, r11646, MPFR_RNDN); mpfr_set_d(r11648, c_n, MPFR_RNDN); mpfr_log1p(r11649, r11643, MPFR_RNDN); mpfr_log1p(r11650, r11637, MPFR_RNDN); mpfr_sub(r11651, r11649, r11650, MPFR_RNDN); mpfr_set_d(r11652, c_p, MPFR_RNDN); mpfr_mul(r11653, r11651, r11652, MPFR_RNDN); mpfr_fma(r11654, r11647, r11648, r11653, MPFR_RNDN); mpfr_exp(r11655, r11654, MPFR_RNDN); return mpfr_get_d(r11655, MPFR_RNDN); }