#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r22421 = 2; float r22422 = sqrt(r22421); float r22423 = 4; float r22424 = r22422 / r22423; float r22425 = 1; float r22426 = 3; float r22427 = v; float r22428 = r22427 * r22427; float r22429 = r22426 * r22428; float r22430 = r22425 - r22429; float r22431 = sqrt(r22430); float r22432 = r22424 * r22431; float r22433 = r22425 - r22428; float r22434 = r22432 * r22433; return r22434; } double f_id(double v) { double r22435 = 2; double r22436 = sqrt(r22435); double r22437 = 4; double r22438 = r22436 / r22437; double r22439 = 1; double r22440 = 3; double r22441 = v; double r22442 = r22441 * r22441; double r22443 = r22440 * r22442; double r22444 = r22439 - r22443; double r22445 = sqrt(r22444); double r22446 = r22438 * r22445; double r22447 = r22439 - r22442; double r22448 = r22446 * r22447; return r22448; } double f_of(float v) { float r22449 = 2; float r22450 = sqrt(r22449); float r22451 = 4; float r22452 = r22450 / r22451; float r22453 = 1; float r22454 = 3; float r22455 = v; float r22456 = r22455 * r22455; float r22457 = r22454 * r22456; float r22458 = r22453 - r22457; float r22459 = sqrt(r22458); float r22460 = r22452 * r22459; float r22461 = exp(r22460); float r22462 = log(r22461); float r22463 = r22453 - r22456; float r22464 = r22462 * r22463; return r22464; } double f_od(double v) { double r22465 = 2; double r22466 = sqrt(r22465); double r22467 = 4; double r22468 = r22466 / r22467; double r22469 = 1; double r22470 = 3; double r22471 = v; double r22472 = r22471 * r22471; double r22473 = r22470 * r22472; double r22474 = r22469 - r22473; double r22475 = sqrt(r22474); double r22476 = r22468 * r22475; double r22477 = exp(r22476); double r22478 = log(r22477); double r22479 = r22469 - r22472; double r22480 = r22478 * r22479; return r22480; } 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 r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490, r22491, r22492, r22493, r22494; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22481, "2", 10, MPFR_RNDN); mpfr_init(r22482); mpfr_init_set_str(r22483, "4", 10, MPFR_RNDN); mpfr_init(r22484); mpfr_init_set_str(r22485, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22486, "3", 10, MPFR_RNDN); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); mpfr_init(r22491); mpfr_init(r22492); mpfr_init(r22493); mpfr_init(r22494); } double f_im(double v) { ; mpfr_sqrt(r22482, r22481, MPFR_RNDN); ; mpfr_div(r22484, r22482, r22483, MPFR_RNDN); ; ; mpfr_set_d(r22487, v, MPFR_RNDN); mpfr_mul(r22488, r22487, r22487, MPFR_RNDN); mpfr_mul(r22489, r22486, r22488, MPFR_RNDN); mpfr_sub(r22490, r22485, r22489, MPFR_RNDN); mpfr_sqrt(r22491, r22490, MPFR_RNDN); mpfr_mul(r22492, r22484, r22491, MPFR_RNDN); mpfr_sub(r22493, r22485, r22488, MPFR_RNDN); mpfr_mul(r22494, r22492, r22493, MPFR_RNDN); return mpfr_get_d(r22494, MPFR_RNDN); } static mpfr_t r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22495, "2", 10, MPFR_RNDN); mpfr_init(r22496); mpfr_init_set_str(r22497, "4", 10, MPFR_RNDN); mpfr_init(r22498); mpfr_init_set_str(r22499, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22500, "3", 10, MPFR_RNDN); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init(r22507); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); } double f_fm(double v) { ; mpfr_sqrt(r22496, r22495, MPFR_RNDN); ; mpfr_div(r22498, r22496, r22497, MPFR_RNDN); ; ; mpfr_set_d(r22501, v, MPFR_RNDN); mpfr_mul(r22502, r22501, r22501, MPFR_RNDN); mpfr_mul(r22503, r22500, r22502, MPFR_RNDN); mpfr_sub(r22504, r22499, r22503, MPFR_RNDN); mpfr_sqrt(r22505, r22504, MPFR_RNDN); mpfr_mul(r22506, r22498, r22505, MPFR_RNDN); mpfr_exp(r22507, r22506, MPFR_RNDN); mpfr_log(r22508, r22507, MPFR_RNDN); mpfr_sub(r22509, r22499, r22502, MPFR_RNDN); mpfr_mul(r22510, r22508, r22509, MPFR_RNDN); return mpfr_get_d(r22510, MPFR_RNDN); } static mpfr_t r22511, r22512, r22513, r22514, r22515, r22516, r22517, r22518, r22519, r22520, r22521, r22522, r22523, r22524, r22525, r22526; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22511, "2", 10, MPFR_RNDN); mpfr_init(r22512); mpfr_init_set_str(r22513, "4", 10, MPFR_RNDN); mpfr_init(r22514); mpfr_init_set_str(r22515, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22516, "3", 10, MPFR_RNDN); mpfr_init(r22517); mpfr_init(r22518); mpfr_init(r22519); mpfr_init(r22520); mpfr_init(r22521); mpfr_init(r22522); mpfr_init(r22523); mpfr_init(r22524); mpfr_init(r22525); mpfr_init(r22526); } double f_dm(double v) { ; mpfr_sqrt(r22512, r22511, MPFR_RNDN); ; mpfr_div(r22514, r22512, r22513, MPFR_RNDN); ; ; mpfr_set_d(r22517, v, MPFR_RNDN); mpfr_mul(r22518, r22517, r22517, MPFR_RNDN); mpfr_mul(r22519, r22516, r22518, MPFR_RNDN); mpfr_sub(r22520, r22515, r22519, MPFR_RNDN); mpfr_sqrt(r22521, r22520, MPFR_RNDN); mpfr_mul(r22522, r22514, r22521, MPFR_RNDN); mpfr_exp(r22523, r22522, MPFR_RNDN); mpfr_log(r22524, r22523, MPFR_RNDN); mpfr_sub(r22525, r22515, r22518, MPFR_RNDN); mpfr_mul(r22526, r22524, r22525, MPFR_RNDN); return mpfr_get_d(r22526, MPFR_RNDN); }