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