#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r26659 = 4; float r26660 = 3; float r26661 = atan2(1.0, 0.0); float r26662 = r26660 * r26661; float r26663 = 1; float r26664 = v; float r26665 = r26664 * r26664; float r26666 = r26663 - r26665; float r26667 = r26662 * r26666; float r26668 = 2; float r26669 = 6; float r26670 = r26669 * r26665; float r26671 = r26668 - r26670; float r26672 = sqrt(r26671); float r26673 = r26667 * r26672; float r26674 = r26659 / r26673; return r26674; } double f_id(double v) { double r26675 = 4; double r26676 = 3; double r26677 = atan2(1.0, 0.0); double r26678 = r26676 * r26677; double r26679 = 1; double r26680 = v; double r26681 = r26680 * r26680; double r26682 = r26679 - r26681; double r26683 = r26678 * r26682; double r26684 = 2; double r26685 = 6; double r26686 = r26685 * r26681; double r26687 = r26684 - r26686; double r26688 = sqrt(r26687); double r26689 = r26683 * r26688; double r26690 = r26675 / r26689; return r26690; } double f_of(float v) { float r26691 = 4; float r26692 = atan2(1.0, 0.0); float r26693 = 3; float r26694 = r26692 * r26693; float r26695 = r26691 / r26694; float r26696 = 1; float r26697 = v; float r26698 = r26697 * r26697; float r26699 = r26696 - r26698; float r26700 = r26695 / r26699; float r26701 = 6; float r26702 = -r26701; float r26703 = 2; float r26704 = fma(r26702, r26698, r26703); float r26705 = sqrt(r26704); float r26706 = r26700 / r26705; float r26707 = pow(r26706, r26693); float r26708 = cbrt(r26707); return r26708; } double f_od(double v) { double r26709 = 4; double r26710 = atan2(1.0, 0.0); double r26711 = 3; double r26712 = r26710 * r26711; double r26713 = r26709 / r26712; double r26714 = 1; double r26715 = v; double r26716 = r26715 * r26715; double r26717 = r26714 - r26716; double r26718 = r26713 / r26717; double r26719 = 6; double r26720 = -r26719; double r26721 = 2; double r26722 = fma(r26720, r26716, r26721); double r26723 = sqrt(r26722); double r26724 = r26718 / r26723; double r26725 = pow(r26724, r26711); double r26726 = cbrt(r26725); return r26726; } 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 r26727, r26728, r26729, r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r26727, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26728, "3", 10, MPFR_RNDN); mpfr_init(r26729); mpfr_init(r26730); mpfr_init_set_str(r26731, "1", 10, MPFR_RNDN); mpfr_init(r26732); mpfr_init(r26733); mpfr_init(r26734); mpfr_init(r26735); mpfr_init_set_str(r26736, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26737, "6", 10, MPFR_RNDN); mpfr_init(r26738); mpfr_init(r26739); mpfr_init(r26740); mpfr_init(r26741); mpfr_init(r26742); } double f_im(double v) { ; ; mpfr_const_pi(r26729, MPFR_RNDN); mpfr_mul(r26730, r26728, r26729, MPFR_RNDN); ; mpfr_set_d(r26732, v, MPFR_RNDN); mpfr_mul(r26733, r26732, r26732, MPFR_RNDN); mpfr_sub(r26734, r26731, r26733, MPFR_RNDN); mpfr_mul(r26735, r26730, r26734, MPFR_RNDN); ; ; mpfr_mul(r26738, r26737, r26733, MPFR_RNDN); mpfr_sub(r26739, r26736, r26738, MPFR_RNDN); mpfr_sqrt(r26740, r26739, MPFR_RNDN); mpfr_mul(r26741, r26735, r26740, MPFR_RNDN); mpfr_div(r26742, r26727, r26741, MPFR_RNDN); return mpfr_get_d(r26742, MPFR_RNDN); } static mpfr_t r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757, r26758, r26759, r26760; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26743, "4", 10, MPFR_RNDN); mpfr_init(r26744); mpfr_init_set_str(r26745, "3", 10, MPFR_RNDN); mpfr_init(r26746); mpfr_init(r26747); mpfr_init_set_str(r26748, "1", 10, MPFR_RNDN); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init_set_str(r26753, "6", 10, MPFR_RNDN); mpfr_init(r26754); mpfr_init_set_str(r26755, "2", 10, MPFR_RNDN); mpfr_init(r26756); mpfr_init(r26757); mpfr_init(r26758); mpfr_init(r26759); mpfr_init(r26760); } double f_fm(double v) { ; mpfr_const_pi(r26744, MPFR_RNDN); ; mpfr_mul(r26746, r26744, r26745, MPFR_RNDN); mpfr_div(r26747, r26743, r26746, MPFR_RNDN); ; mpfr_set_d(r26749, v, MPFR_RNDN); mpfr_mul(r26750, r26749, r26749, MPFR_RNDN); mpfr_sub(r26751, r26748, r26750, MPFR_RNDN); mpfr_div(r26752, r26747, r26751, MPFR_RNDN); ; mpfr_neg(r26754, r26753, MPFR_RNDN); ; mpfr_fma(r26756, r26754, r26750, r26755, MPFR_RNDN); mpfr_sqrt(r26757, r26756, MPFR_RNDN); mpfr_div(r26758, r26752, r26757, MPFR_RNDN); mpfr_pow(r26759, r26758, r26745, MPFR_RNDN); mpfr_cbrt(r26760, r26759, MPFR_RNDN); return mpfr_get_d(r26760, MPFR_RNDN); } static mpfr_t r26761, r26762, r26763, r26764, r26765, r26766, r26767, r26768, r26769, r26770, r26771, r26772, r26773, r26774, r26775, r26776, r26777, r26778; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26761, "4", 10, MPFR_RNDN); mpfr_init(r26762); mpfr_init_set_str(r26763, "3", 10, MPFR_RNDN); mpfr_init(r26764); mpfr_init(r26765); mpfr_init_set_str(r26766, "1", 10, MPFR_RNDN); mpfr_init(r26767); mpfr_init(r26768); mpfr_init(r26769); mpfr_init(r26770); mpfr_init_set_str(r26771, "6", 10, MPFR_RNDN); mpfr_init(r26772); mpfr_init_set_str(r26773, "2", 10, MPFR_RNDN); mpfr_init(r26774); mpfr_init(r26775); mpfr_init(r26776); mpfr_init(r26777); mpfr_init(r26778); } double f_dm(double v) { ; mpfr_const_pi(r26762, MPFR_RNDN); ; mpfr_mul(r26764, r26762, r26763, MPFR_RNDN); mpfr_div(r26765, r26761, r26764, MPFR_RNDN); ; mpfr_set_d(r26767, v, MPFR_RNDN); mpfr_mul(r26768, r26767, r26767, MPFR_RNDN); mpfr_sub(r26769, r26766, r26768, MPFR_RNDN); mpfr_div(r26770, r26765, r26769, MPFR_RNDN); ; mpfr_neg(r26772, r26771, MPFR_RNDN); ; mpfr_fma(r26774, r26772, r26768, r26773, MPFR_RNDN); mpfr_sqrt(r26775, r26774, MPFR_RNDN); mpfr_div(r26776, r26770, r26775, MPFR_RNDN); mpfr_pow(r26777, r26776, r26763, MPFR_RNDN); mpfr_cbrt(r26778, r26777, MPFR_RNDN); return mpfr_get_d(r26778, MPFR_RNDN); }