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