#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r27482 = 2; float r27483 = sqrt(r27482); float r27484 = 4; float r27485 = r27483 / r27484; float r27486 = 1; float r27487 = 3; float r27488 = v; float r27489 = r27488 * r27488; float r27490 = r27487 * r27489; float r27491 = r27486 - r27490; float r27492 = sqrt(r27491); float r27493 = r27485 * r27492; float r27494 = r27486 - r27489; float r27495 = r27493 * r27494; return r27495; } double f_id(double v) { double r27496 = 2; double r27497 = sqrt(r27496); double r27498 = 4; double r27499 = r27497 / r27498; double r27500 = 1; double r27501 = 3; double r27502 = v; double r27503 = r27502 * r27502; double r27504 = r27501 * r27503; double r27505 = r27500 - r27504; double r27506 = sqrt(r27505); double r27507 = r27499 * r27506; double r27508 = r27500 - r27503; double r27509 = r27507 * r27508; return r27509; } double f_of(float v) { float r27510 = 2; float r27511 = sqrt(r27510); float r27512 = 4; float r27513 = r27511 / r27512; float r27514 = log(r27513); float r27515 = 1; float r27516 = 3; float r27517 = v; float r27518 = r27517 * r27517; float r27519 = r27516 * r27518; float r27520 = r27515 - r27519; float r27521 = sqrt(r27520); float r27522 = log(r27521); float r27523 = r27514 + r27522; float r27524 = r27515 - r27518; float r27525 = log(r27524); float r27526 = r27523 + r27525; float r27527 = exp(r27526); return r27527; } double f_od(double v) { double r27528 = 2; double r27529 = sqrt(r27528); double r27530 = 4; double r27531 = r27529 / r27530; double r27532 = log(r27531); double r27533 = 1; double r27534 = 3; double r27535 = v; double r27536 = r27535 * r27535; double r27537 = r27534 * r27536; double r27538 = r27533 - r27537; double r27539 = sqrt(r27538); double r27540 = log(r27539); double r27541 = r27532 + r27540; double r27542 = r27533 - r27536; double r27543 = log(r27542); double r27544 = r27541 + r27543; double r27545 = exp(r27544); return r27545; } 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 r27546, r27547, r27548, r27549, r27550, r27551, r27552, r27553, r27554, r27555, r27556, r27557, r27558, r27559; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27546, "2", 10, MPFR_RNDN); mpfr_init(r27547); mpfr_init_set_str(r27548, "4", 10, MPFR_RNDN); mpfr_init(r27549); mpfr_init_set_str(r27550, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27551, "3", 10, MPFR_RNDN); mpfr_init(r27552); mpfr_init(r27553); mpfr_init(r27554); mpfr_init(r27555); mpfr_init(r27556); mpfr_init(r27557); mpfr_init(r27558); mpfr_init(r27559); } double f_im(double v) { ; mpfr_sqrt(r27547, r27546, MPFR_RNDN); ; mpfr_div(r27549, r27547, r27548, MPFR_RNDN); ; ; mpfr_set_d(r27552, v, MPFR_RNDN); mpfr_mul(r27553, r27552, r27552, MPFR_RNDN); mpfr_mul(r27554, r27551, r27553, MPFR_RNDN); mpfr_sub(r27555, r27550, r27554, MPFR_RNDN); mpfr_sqrt(r27556, r27555, MPFR_RNDN); mpfr_mul(r27557, r27549, r27556, MPFR_RNDN); mpfr_sub(r27558, r27550, r27553, MPFR_RNDN); mpfr_mul(r27559, r27557, r27558, MPFR_RNDN); return mpfr_get_d(r27559, MPFR_RNDN); } static mpfr_t r27560, r27561, r27562, r27563, r27564, r27565, r27566, r27567, r27568, r27569, r27570, r27571, r27572, r27573, r27574, r27575, r27576, r27577; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27560, "2", 10, MPFR_RNDN); mpfr_init(r27561); mpfr_init_set_str(r27562, "4", 10, MPFR_RNDN); mpfr_init(r27563); mpfr_init(r27564); mpfr_init_set_str(r27565, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27566, "3", 10, MPFR_RNDN); mpfr_init(r27567); mpfr_init(r27568); mpfr_init(r27569); mpfr_init(r27570); mpfr_init(r27571); mpfr_init(r27572); mpfr_init(r27573); mpfr_init(r27574); mpfr_init(r27575); mpfr_init(r27576); mpfr_init(r27577); } double f_fm(double v) { ; mpfr_sqrt(r27561, r27560, MPFR_RNDN); ; mpfr_div(r27563, r27561, r27562, MPFR_RNDN); mpfr_log(r27564, r27563, MPFR_RNDN); ; ; mpfr_set_d(r27567, v, MPFR_RNDN); mpfr_mul(r27568, r27567, r27567, MPFR_RNDN); mpfr_mul(r27569, r27566, r27568, MPFR_RNDN); mpfr_sub(r27570, r27565, r27569, MPFR_RNDN); mpfr_sqrt(r27571, r27570, MPFR_RNDN); mpfr_log(r27572, r27571, MPFR_RNDN); mpfr_add(r27573, r27564, r27572, MPFR_RNDN); mpfr_sub(r27574, r27565, r27568, MPFR_RNDN); mpfr_log(r27575, r27574, MPFR_RNDN); mpfr_add(r27576, r27573, r27575, MPFR_RNDN); mpfr_exp(r27577, r27576, MPFR_RNDN); return mpfr_get_d(r27577, MPFR_RNDN); } static mpfr_t r27578, r27579, r27580, r27581, r27582, r27583, r27584, r27585, r27586, r27587, r27588, r27589, r27590, r27591, r27592, r27593, r27594, r27595; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27578, "2", 10, MPFR_RNDN); mpfr_init(r27579); mpfr_init_set_str(r27580, "4", 10, MPFR_RNDN); mpfr_init(r27581); mpfr_init(r27582); mpfr_init_set_str(r27583, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27584, "3", 10, MPFR_RNDN); mpfr_init(r27585); mpfr_init(r27586); mpfr_init(r27587); mpfr_init(r27588); mpfr_init(r27589); mpfr_init(r27590); mpfr_init(r27591); mpfr_init(r27592); mpfr_init(r27593); mpfr_init(r27594); mpfr_init(r27595); } double f_dm(double v) { ; mpfr_sqrt(r27579, r27578, MPFR_RNDN); ; mpfr_div(r27581, r27579, r27580, MPFR_RNDN); mpfr_log(r27582, r27581, MPFR_RNDN); ; ; mpfr_set_d(r27585, v, MPFR_RNDN); mpfr_mul(r27586, r27585, r27585, MPFR_RNDN); mpfr_mul(r27587, r27584, r27586, MPFR_RNDN); mpfr_sub(r27588, r27583, r27587, MPFR_RNDN); mpfr_sqrt(r27589, r27588, MPFR_RNDN); mpfr_log(r27590, r27589, MPFR_RNDN); mpfr_add(r27591, r27582, r27590, MPFR_RNDN); mpfr_sub(r27592, r27583, r27586, MPFR_RNDN); mpfr_log(r27593, r27592, MPFR_RNDN); mpfr_add(r27594, r27591, r27593, MPFR_RNDN); mpfr_exp(r27595, r27594, MPFR_RNDN); return mpfr_get_d(r27595, MPFR_RNDN); }