#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r24526 = 4; float r24527 = 3; float r24528 = atan2(1.0, 0.0); float r24529 = r24527 * r24528; float r24530 = 1; float r24531 = v; float r24532 = r24531 * r24531; float r24533 = r24530 - r24532; float r24534 = r24529 * r24533; float r24535 = 2; float r24536 = 6; float r24537 = r24536 * r24532; float r24538 = r24535 - r24537; float r24539 = sqrt(r24538); float r24540 = r24534 * r24539; float r24541 = r24526 / r24540; return r24541; } double f_id(double v) { double r24542 = 4; double r24543 = 3; double r24544 = atan2(1.0, 0.0); double r24545 = r24543 * r24544; double r24546 = 1; double r24547 = v; double r24548 = r24547 * r24547; double r24549 = r24546 - r24548; double r24550 = r24545 * r24549; double r24551 = 2; double r24552 = 6; double r24553 = r24552 * r24548; double r24554 = r24551 - r24553; double r24555 = sqrt(r24554); double r24556 = r24550 * r24555; double r24557 = r24542 / r24556; return r24557; } double f_of(float v) { float r24558 = 4; float r24559 = atan2(1.0, 0.0); float r24560 = 3; float r24561 = r24559 * r24560; float r24562 = r24558 / r24561; float r24563 = 1; float r24564 = v; float r24565 = r24564 * r24564; float r24566 = r24563 - r24565; float r24567 = r24562 / r24566; float r24568 = 6; float r24569 = -r24568; float r24570 = 2; float r24571 = fma(r24569, r24565, r24570); float r24572 = sqrt(r24571); float r24573 = r24567 / r24572; float r24574 = pow(r24573, r24560); float r24575 = cbrt(r24574); return r24575; } double f_od(double v) { double r24576 = 4; double r24577 = atan2(1.0, 0.0); double r24578 = 3; double r24579 = r24577 * r24578; double r24580 = r24576 / r24579; double r24581 = 1; double r24582 = v; double r24583 = r24582 * r24582; double r24584 = r24581 - r24583; double r24585 = r24580 / r24584; double r24586 = 6; double r24587 = -r24586; double r24588 = 2; double r24589 = fma(r24587, r24583, r24588); double r24590 = sqrt(r24589); double r24591 = r24585 / r24590; double r24592 = pow(r24591, r24578); double r24593 = cbrt(r24592); return r24593; } 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 r24594, r24595, r24596, r24597, r24598, r24599, r24600, r24601, r24602, r24603, r24604, r24605, r24606, r24607, r24608, r24609; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r24594, "4", 10, MPFR_RNDN); mpfr_init_set_str(r24595, "3", 10, MPFR_RNDN); mpfr_init(r24596); mpfr_init(r24597); mpfr_init_set_str(r24598, "1", 10, MPFR_RNDN); mpfr_init(r24599); mpfr_init(r24600); mpfr_init(r24601); mpfr_init(r24602); mpfr_init_set_str(r24603, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24604, "6", 10, MPFR_RNDN); mpfr_init(r24605); mpfr_init(r24606); mpfr_init(r24607); mpfr_init(r24608); mpfr_init(r24609); } double f_im(double v) { ; ; mpfr_const_pi(r24596, MPFR_RNDN); mpfr_mul(r24597, r24595, r24596, MPFR_RNDN); ; mpfr_set_d(r24599, v, MPFR_RNDN); mpfr_mul(r24600, r24599, r24599, MPFR_RNDN); mpfr_sub(r24601, r24598, r24600, MPFR_RNDN); mpfr_mul(r24602, r24597, r24601, MPFR_RNDN); ; ; mpfr_mul(r24605, r24604, r24600, MPFR_RNDN); mpfr_sub(r24606, r24603, r24605, MPFR_RNDN); mpfr_sqrt(r24607, r24606, MPFR_RNDN); mpfr_mul(r24608, r24602, r24607, MPFR_RNDN); mpfr_div(r24609, r24594, r24608, MPFR_RNDN); return mpfr_get_d(r24609, MPFR_RNDN); } static mpfr_t r24610, r24611, r24612, r24613, r24614, r24615, r24616, r24617, r24618, r24619, r24620, r24621, r24622, r24623, r24624, r24625, r24626, r24627; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24610, "4", 10, MPFR_RNDN); mpfr_init(r24611); mpfr_init_set_str(r24612, "3", 10, MPFR_RNDN); mpfr_init(r24613); mpfr_init(r24614); mpfr_init_set_str(r24615, "1", 10, MPFR_RNDN); mpfr_init(r24616); mpfr_init(r24617); mpfr_init(r24618); mpfr_init(r24619); mpfr_init_set_str(r24620, "6", 10, MPFR_RNDN); mpfr_init(r24621); mpfr_init_set_str(r24622, "2", 10, MPFR_RNDN); mpfr_init(r24623); mpfr_init(r24624); mpfr_init(r24625); mpfr_init(r24626); mpfr_init(r24627); } double f_fm(double v) { ; mpfr_const_pi(r24611, MPFR_RNDN); ; mpfr_mul(r24613, r24611, r24612, MPFR_RNDN); mpfr_div(r24614, r24610, r24613, MPFR_RNDN); ; mpfr_set_d(r24616, v, MPFR_RNDN); mpfr_mul(r24617, r24616, r24616, MPFR_RNDN); mpfr_sub(r24618, r24615, r24617, MPFR_RNDN); mpfr_div(r24619, r24614, r24618, MPFR_RNDN); ; mpfr_neg(r24621, r24620, MPFR_RNDN); ; mpfr_fma(r24623, r24621, r24617, r24622, MPFR_RNDN); mpfr_sqrt(r24624, r24623, MPFR_RNDN); mpfr_div(r24625, r24619, r24624, MPFR_RNDN); mpfr_pow(r24626, r24625, r24612, MPFR_RNDN); mpfr_cbrt(r24627, r24626, MPFR_RNDN); return mpfr_get_d(r24627, MPFR_RNDN); } static mpfr_t r24628, r24629, r24630, r24631, r24632, r24633, r24634, r24635, r24636, r24637, r24638, r24639, r24640, r24641, r24642, r24643, r24644, r24645; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24628, "4", 10, MPFR_RNDN); mpfr_init(r24629); mpfr_init_set_str(r24630, "3", 10, MPFR_RNDN); mpfr_init(r24631); mpfr_init(r24632); mpfr_init_set_str(r24633, "1", 10, MPFR_RNDN); mpfr_init(r24634); mpfr_init(r24635); mpfr_init(r24636); mpfr_init(r24637); mpfr_init_set_str(r24638, "6", 10, MPFR_RNDN); mpfr_init(r24639); mpfr_init_set_str(r24640, "2", 10, MPFR_RNDN); mpfr_init(r24641); mpfr_init(r24642); mpfr_init(r24643); mpfr_init(r24644); mpfr_init(r24645); } double f_dm(double v) { ; mpfr_const_pi(r24629, MPFR_RNDN); ; mpfr_mul(r24631, r24629, r24630, MPFR_RNDN); mpfr_div(r24632, r24628, r24631, MPFR_RNDN); ; mpfr_set_d(r24634, v, MPFR_RNDN); mpfr_mul(r24635, r24634, r24634, MPFR_RNDN); mpfr_sub(r24636, r24633, r24635, MPFR_RNDN); mpfr_div(r24637, r24632, r24636, MPFR_RNDN); ; mpfr_neg(r24639, r24638, MPFR_RNDN); ; mpfr_fma(r24641, r24639, r24635, r24640, MPFR_RNDN); mpfr_sqrt(r24642, r24641, MPFR_RNDN); mpfr_div(r24643, r24637, r24642, MPFR_RNDN); mpfr_pow(r24644, r24643, r24630, MPFR_RNDN); mpfr_cbrt(r24645, r24644, MPFR_RNDN); return mpfr_get_d(r24645, MPFR_RNDN); }