#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r21509 = 2.0f; float r21510 = sqrt(r21509); float r21511 = 4.0f; float r21512 = r21510 / r21511; float r21513 = 1.0f; float r21514 = 3.0f; float r21515 = v; float r21516 = r21515 * r21515; float r21517 = r21514 * r21516; float r21518 = r21513 - r21517; float r21519 = sqrt(r21518); float r21520 = r21512 * r21519; float r21521 = r21513 - r21516; float r21522 = r21520 * r21521; return r21522; } double f_id(double v) { double r21523 = 2.0; double r21524 = sqrt(r21523); double r21525 = 4.0; double r21526 = r21524 / r21525; double r21527 = 1.0; double r21528 = 3.0; double r21529 = v; double r21530 = r21529 * r21529; double r21531 = r21528 * r21530; double r21532 = r21527 - r21531; double r21533 = sqrt(r21532); double r21534 = r21526 * r21533; double r21535 = r21527 - r21530; double r21536 = r21534 * r21535; return r21536; } double f_of(float v) { float r21537 = 2.0f; float r21538 = sqrt(r21537); float r21539 = 4.0f; float r21540 = r21538 / r21539; float r21541 = 1.0f; float r21542 = 3.0f; float r21543 = v; float r21544 = r21543 * r21543; float r21545 = r21542 * r21544; float r21546 = r21541 - r21545; float r21547 = sqrt(r21546); float r21548 = r21540 * r21547; float r21549 = log(r21548); float r21550 = exp(r21549); float r21551 = r21541 - r21544; float r21552 = r21550 * r21551; return r21552; } double f_od(double v) { double r21553 = 2.0; double r21554 = sqrt(r21553); double r21555 = 4.0; double r21556 = r21554 / r21555; double r21557 = 1.0; double r21558 = 3.0; double r21559 = v; double r21560 = r21559 * r21559; double r21561 = r21558 * r21560; double r21562 = r21557 - r21561; double r21563 = sqrt(r21562); double r21564 = r21556 * r21563; double r21565 = log(r21564); double r21566 = exp(r21565); double r21567 = r21557 - r21560; double r21568 = r21566 * r21567; return r21568; } 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 r21569, r21570, r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579, r21580, r21581, r21582; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21569, "2", 10, MPFR_RNDN); mpfr_init(r21570); mpfr_init_set_str(r21571, "4", 10, MPFR_RNDN); mpfr_init(r21572); mpfr_init_set_str(r21573, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21574, "3", 10, MPFR_RNDN); mpfr_init(r21575); mpfr_init(r21576); mpfr_init(r21577); mpfr_init(r21578); mpfr_init(r21579); mpfr_init(r21580); mpfr_init(r21581); mpfr_init(r21582); } double f_im(double v) { ; mpfr_sqrt(r21570, r21569, MPFR_RNDN); ; mpfr_div(r21572, r21570, r21571, MPFR_RNDN); ; ; mpfr_set_d(r21575, v, MPFR_RNDN); mpfr_sqr(r21576, r21575, MPFR_RNDN); mpfr_mul(r21577, r21574, r21576, MPFR_RNDN); mpfr_sub(r21578, r21573, r21577, MPFR_RNDN); mpfr_sqrt(r21579, r21578, MPFR_RNDN); mpfr_mul(r21580, r21572, r21579, MPFR_RNDN); mpfr_sub(r21581, r21573, r21576, MPFR_RNDN); mpfr_mul(r21582, r21580, r21581, MPFR_RNDN); return mpfr_get_d(r21582, MPFR_RNDN); } static mpfr_t r21583, r21584, r21585, r21586, r21587, r21588, r21589, r21590, r21591, r21592, r21593, r21594, r21595, r21596, r21597, r21598; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21583, "2", 10, MPFR_RNDN); mpfr_init(r21584); mpfr_init_set_str(r21585, "4", 10, MPFR_RNDN); mpfr_init(r21586); mpfr_init_set_str(r21587, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21588, "3", 10, MPFR_RNDN); mpfr_init(r21589); mpfr_init(r21590); mpfr_init(r21591); mpfr_init(r21592); mpfr_init(r21593); mpfr_init(r21594); mpfr_init(r21595); mpfr_init(r21596); mpfr_init(r21597); mpfr_init(r21598); } double f_fm(double v) { ; mpfr_sqrt(r21584, r21583, MPFR_RNDN); ; mpfr_div(r21586, r21584, r21585, MPFR_RNDN); ; ; mpfr_set_d(r21589, v, MPFR_RNDN); mpfr_sqr(r21590, r21589, MPFR_RNDN); mpfr_mul(r21591, r21588, r21590, MPFR_RNDN); mpfr_sub(r21592, r21587, r21591, MPFR_RNDN); mpfr_sqrt(r21593, r21592, MPFR_RNDN); mpfr_mul(r21594, r21586, r21593, MPFR_RNDN); mpfr_log(r21595, r21594, MPFR_RNDN); mpfr_exp(r21596, r21595, MPFR_RNDN); mpfr_sub(r21597, r21587, r21590, MPFR_RNDN); mpfr_mul(r21598, r21596, r21597, MPFR_RNDN); return mpfr_get_d(r21598, MPFR_RNDN); } static mpfr_t r21599, r21600, r21601, r21602, r21603, r21604, r21605, r21606, r21607, r21608, r21609, r21610, r21611, r21612, r21613, r21614; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21599, "2", 10, MPFR_RNDN); mpfr_init(r21600); mpfr_init_set_str(r21601, "4", 10, MPFR_RNDN); mpfr_init(r21602); mpfr_init_set_str(r21603, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21604, "3", 10, MPFR_RNDN); mpfr_init(r21605); mpfr_init(r21606); mpfr_init(r21607); mpfr_init(r21608); mpfr_init(r21609); mpfr_init(r21610); mpfr_init(r21611); mpfr_init(r21612); mpfr_init(r21613); mpfr_init(r21614); } double f_dm(double v) { ; mpfr_sqrt(r21600, r21599, MPFR_RNDN); ; mpfr_div(r21602, r21600, r21601, MPFR_RNDN); ; ; mpfr_set_d(r21605, v, MPFR_RNDN); mpfr_sqr(r21606, r21605, MPFR_RNDN); mpfr_mul(r21607, r21604, r21606, MPFR_RNDN); mpfr_sub(r21608, r21603, r21607, MPFR_RNDN); mpfr_sqrt(r21609, r21608, MPFR_RNDN); mpfr_mul(r21610, r21602, r21609, MPFR_RNDN); mpfr_log(r21611, r21610, MPFR_RNDN); mpfr_exp(r21612, r21611, MPFR_RNDN); mpfr_sub(r21613, r21603, r21606, MPFR_RNDN); mpfr_mul(r21614, r21612, r21613, MPFR_RNDN); return mpfr_get_d(r21614, MPFR_RNDN); }