#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r24669 = a; float r24670 = k; float r24671 = m; float r24672 = pow(r24670, r24671); float r24673 = r24669 * r24672; float r24674 = 1; float r24675 = 10; float r24676 = r24675 * r24670; float r24677 = r24674 + r24676; float r24678 = r24670 * r24670; float r24679 = r24677 + r24678; float r24680 = r24673 / r24679; return r24680; } double f_id(double a, double k, double m) { double r24681 = a; double r24682 = k; double r24683 = m; double r24684 = pow(r24682, r24683); double r24685 = r24681 * r24684; double r24686 = 1; double r24687 = 10; double r24688 = r24687 * r24682; double r24689 = r24686 + r24688; double r24690 = r24682 * r24682; double r24691 = r24689 + r24690; double r24692 = r24685 / r24691; return r24692; } double f_of(float a, float k, float m) { float r24693 = k; float r24694 = 3.6364048152768866e+66; bool r24695 = r24693 <= r24694; float r24696 = a; float r24697 = m; float r24698 = pow(r24693, r24697); float r24699 = r24696 * r24698; float r24700 = 1; float r24701 = 10; float r24702 = r24701 * r24693; float r24703 = r24700 + r24702; float r24704 = r24693 * r24693; float r24705 = r24703 + r24704; float r24706 = r24699 / r24705; float r24707 = r24693 / r24696; float r24708 = -r24697; float r24709 = pow(r24693, r24708); float r24710 = r24693 + r24701; float r24711 = r24709 * r24710; float r24712 = r24709 / r24696; float r24713 = fma(r24707, r24711, r24712); float r24714 = r24700 / r24713; float r24715 = r24695 ? r24706 : r24714; return r24715; } double f_od(double a, double k, double m) { double r24716 = k; double r24717 = 3.6364048152768866e+66; bool r24718 = r24716 <= r24717; double r24719 = a; double r24720 = m; double r24721 = pow(r24716, r24720); double r24722 = r24719 * r24721; double r24723 = 1; double r24724 = 10; double r24725 = r24724 * r24716; double r24726 = r24723 + r24725; double r24727 = r24716 * r24716; double r24728 = r24726 + r24727; double r24729 = r24722 / r24728; double r24730 = r24716 / r24719; double r24731 = -r24720; double r24732 = pow(r24716, r24731); double r24733 = r24716 + r24724; double r24734 = r24732 * r24733; double r24735 = r24732 / r24719; double r24736 = fma(r24730, r24734, r24735); double r24737 = r24723 / r24736; double r24738 = r24718 ? r24729 : r24737; return r24738; } 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 r24739, r24740, r24741, r24742, r24743, r24744, r24745, r24746, r24747, r24748, r24749, r24750; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24739); mpfr_init(r24740); mpfr_init(r24741); mpfr_init(r24742); mpfr_init(r24743); mpfr_init_set_str(r24744, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24745, "10", 10, MPFR_RNDN); mpfr_init(r24746); mpfr_init(r24747); mpfr_init(r24748); mpfr_init(r24749); mpfr_init(r24750); } double f_im(double a, double k, double m) { mpfr_set_d(r24739, a, MPFR_RNDN); mpfr_set_d(r24740, k, MPFR_RNDN); mpfr_set_d(r24741, m, MPFR_RNDN); mpfr_pow(r24742, r24740, r24741, MPFR_RNDN); mpfr_mul(r24743, r24739, r24742, MPFR_RNDN); ; ; mpfr_mul(r24746, r24745, r24740, MPFR_RNDN); mpfr_add(r24747, r24744, r24746, MPFR_RNDN); mpfr_mul(r24748, r24740, r24740, MPFR_RNDN); mpfr_add(r24749, r24747, r24748, MPFR_RNDN); mpfr_div(r24750, r24743, r24749, MPFR_RNDN); return mpfr_get_d(r24750, MPFR_RNDN); } static mpfr_t r24751, r24752, r24753, r24754, r24755, r24756, r24757, r24758, r24759, r24760, r24761, r24762, r24763, r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24751); mpfr_init_set_str(r24752, "3.6364048152768866e+66", 10, MPFR_RNDN); mpfr_init(r24753); mpfr_init(r24754); mpfr_init(r24755); mpfr_init(r24756); mpfr_init(r24757); mpfr_init_set_str(r24758, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24759, "10", 10, MPFR_RNDN); mpfr_init(r24760); mpfr_init(r24761); mpfr_init(r24762); mpfr_init(r24763); mpfr_init(r24764); mpfr_init(r24765); mpfr_init(r24766); mpfr_init(r24767); mpfr_init(r24768); mpfr_init(r24769); mpfr_init(r24770); mpfr_init(r24771); mpfr_init(r24772); mpfr_init(r24773); } double f_fm(double a, double k, double m) { mpfr_set_d(r24751, k, MPFR_RNDN); ; mpfr_set_si(r24753, mpfr_cmp(r24751, r24752) <= 0, MPFR_RNDN); mpfr_set_d(r24754, a, MPFR_RNDN); mpfr_set_d(r24755, m, MPFR_RNDN); mpfr_pow(r24756, r24751, r24755, MPFR_RNDN); mpfr_mul(r24757, r24754, r24756, MPFR_RNDN); ; ; mpfr_mul(r24760, r24759, r24751, MPFR_RNDN); mpfr_add(r24761, r24758, r24760, MPFR_RNDN); mpfr_mul(r24762, r24751, r24751, MPFR_RNDN); mpfr_add(r24763, r24761, r24762, MPFR_RNDN); mpfr_div(r24764, r24757, r24763, MPFR_RNDN); mpfr_div(r24765, r24751, r24754, MPFR_RNDN); mpfr_neg(r24766, r24755, MPFR_RNDN); mpfr_pow(r24767, r24751, r24766, MPFR_RNDN); mpfr_add(r24768, r24751, r24759, MPFR_RNDN); mpfr_mul(r24769, r24767, r24768, MPFR_RNDN); mpfr_div(r24770, r24767, r24754, MPFR_RNDN); mpfr_fma(r24771, r24765, r24769, r24770, MPFR_RNDN); mpfr_div(r24772, r24758, r24771, MPFR_RNDN); if (mpfr_get_si(r24753, MPFR_RNDN)) { mpfr_set(r24773, r24764, MPFR_RNDN); } else { mpfr_set(r24773, r24772, MPFR_RNDN); }; return mpfr_get_d(r24773, MPFR_RNDN); } static mpfr_t r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792, r24793, r24794, r24795, r24796; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24774); mpfr_init_set_str(r24775, "3.6364048152768866e+66", 10, MPFR_RNDN); mpfr_init(r24776); mpfr_init(r24777); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init_set_str(r24781, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24782, "10", 10, MPFR_RNDN); mpfr_init(r24783); mpfr_init(r24784); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); mpfr_init(r24793); mpfr_init(r24794); mpfr_init(r24795); mpfr_init(r24796); } double f_dm(double a, double k, double m) { mpfr_set_d(r24774, k, MPFR_RNDN); ; mpfr_set_si(r24776, mpfr_cmp(r24774, r24775) <= 0, MPFR_RNDN); mpfr_set_d(r24777, a, MPFR_RNDN); mpfr_set_d(r24778, m, MPFR_RNDN); mpfr_pow(r24779, r24774, r24778, MPFR_RNDN); mpfr_mul(r24780, r24777, r24779, MPFR_RNDN); ; ; mpfr_mul(r24783, r24782, r24774, MPFR_RNDN); mpfr_add(r24784, r24781, r24783, MPFR_RNDN); mpfr_mul(r24785, r24774, r24774, MPFR_RNDN); mpfr_add(r24786, r24784, r24785, MPFR_RNDN); mpfr_div(r24787, r24780, r24786, MPFR_RNDN); mpfr_div(r24788, r24774, r24777, MPFR_RNDN); mpfr_neg(r24789, r24778, MPFR_RNDN); mpfr_pow(r24790, r24774, r24789, MPFR_RNDN); mpfr_add(r24791, r24774, r24782, MPFR_RNDN); mpfr_mul(r24792, r24790, r24791, MPFR_RNDN); mpfr_div(r24793, r24790, r24777, MPFR_RNDN); mpfr_fma(r24794, r24788, r24792, r24793, MPFR_RNDN); mpfr_div(r24795, r24781, r24794, MPFR_RNDN); if (mpfr_get_si(r24776, MPFR_RNDN)) { mpfr_set(r24796, r24787, MPFR_RNDN); } else { mpfr_set(r24796, r24795, MPFR_RNDN); }; return mpfr_get_d(r24796, MPFR_RNDN); }