#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r24718 = a; float r24719 = k; float r24720 = m; float r24721 = pow(r24719, r24720); float r24722 = r24718 * r24721; float r24723 = 1; float r24724 = 10; float r24725 = r24724 * r24719; float r24726 = r24723 + r24725; float r24727 = r24719 * r24719; float r24728 = r24726 + r24727; float r24729 = r24722 / r24728; return r24729; } double f_id(double a, double k, double m) { double r24730 = a; double r24731 = k; double r24732 = m; double r24733 = pow(r24731, r24732); double r24734 = r24730 * r24733; double r24735 = 1; double r24736 = 10; double r24737 = r24736 * r24731; double r24738 = r24735 + r24737; double r24739 = r24731 * r24731; double r24740 = r24738 + r24739; double r24741 = r24734 / r24740; return r24741; } double f_of(float a, float k, float m) { float r24742 = k; float r24743 = 6.493800852448035e+60; bool r24744 = r24742 <= r24743; float r24745 = m; float r24746 = pow(r24742, r24745); float r24747 = a; float r24748 = r24746 * r24747; float r24749 = 1; float r24750 = 10; float r24751 = r24750 + r24742; float r24752 = r24751 * r24742; float r24753 = r24749 + r24752; float r24754 = r24749 / r24753; float r24755 = r24748 * r24754; float r24756 = -r24745; float r24757 = pow(r24742, r24756); float r24758 = r24757 / r24747; float r24759 = r24742 / r24747; float r24760 = r24759 / r24746; float r24761 = r24760 * r24751; float r24762 = r24758 + r24761; float r24763 = r24749 / r24762; float r24764 = r24744 ? r24755 : r24763; return r24764; } double f_od(double a, double k, double m) { double r24765 = k; double r24766 = 6.493800852448035e+60; bool r24767 = r24765 <= r24766; double r24768 = m; double r24769 = pow(r24765, r24768); double r24770 = a; double r24771 = r24769 * r24770; double r24772 = 1; double r24773 = 10; double r24774 = r24773 + r24765; double r24775 = r24774 * r24765; double r24776 = r24772 + r24775; double r24777 = r24772 / r24776; double r24778 = r24771 * r24777; double r24779 = -r24768; double r24780 = pow(r24765, r24779); double r24781 = r24780 / r24770; double r24782 = r24765 / r24770; double r24783 = r24782 / r24769; double r24784 = r24783 * r24774; double r24785 = r24781 + r24784; double r24786 = r24772 / r24785; double r24787 = r24767 ? r24778 : r24786; return r24787; } 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 r24788, r24789, r24790, r24791, r24792, r24793, r24794, r24795, r24796, r24797, r24798, r24799; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); mpfr_init_set_str(r24793, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24794, "10", 10, MPFR_RNDN); mpfr_init(r24795); mpfr_init(r24796); mpfr_init(r24797); mpfr_init(r24798); mpfr_init(r24799); } double f_im(double a, double k, double m) { mpfr_set_d(r24788, a, MPFR_RNDN); mpfr_set_d(r24789, k, MPFR_RNDN); mpfr_set_d(r24790, m, MPFR_RNDN); mpfr_pow(r24791, r24789, r24790, MPFR_RNDN); mpfr_mul(r24792, r24788, r24791, MPFR_RNDN); ; ; mpfr_mul(r24795, r24794, r24789, MPFR_RNDN); mpfr_add(r24796, r24793, r24795, MPFR_RNDN); mpfr_mul(r24797, r24789, r24789, MPFR_RNDN); mpfr_add(r24798, r24796, r24797, MPFR_RNDN); mpfr_div(r24799, r24792, r24798, MPFR_RNDN); return mpfr_get_d(r24799, MPFR_RNDN); } static mpfr_t r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808, r24809, r24810, r24811, r24812, r24813, r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r24800); mpfr_init_set_str(r24801, "6.493800852448035e+60", 10, MPFR_RNDN); mpfr_init(r24802); mpfr_init(r24803); mpfr_init(r24804); mpfr_init(r24805); mpfr_init(r24806); mpfr_init_set_str(r24807, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24808, "10", 10, MPFR_RNDN); mpfr_init(r24809); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); mpfr_init(r24813); mpfr_init(r24814); mpfr_init(r24815); mpfr_init(r24816); mpfr_init(r24817); mpfr_init(r24818); mpfr_init(r24819); mpfr_init(r24820); mpfr_init(r24821); mpfr_init(r24822); } double f_fm(double a, double k, double m) { mpfr_set_d(r24800, k, MPFR_RNDN); ; mpfr_set_si(r24802, mpfr_cmp(r24800, r24801) <= 0, MPFR_RNDN); mpfr_set_d(r24803, m, MPFR_RNDN); mpfr_pow(r24804, r24800, r24803, MPFR_RNDN); mpfr_set_d(r24805, a, MPFR_RNDN); mpfr_mul(r24806, r24804, r24805, MPFR_RNDN); ; ; mpfr_add(r24809, r24808, r24800, MPFR_RNDN); mpfr_mul(r24810, r24809, r24800, MPFR_RNDN); mpfr_add(r24811, r24807, r24810, MPFR_RNDN); mpfr_div(r24812, r24807, r24811, MPFR_RNDN); mpfr_mul(r24813, r24806, r24812, MPFR_RNDN); mpfr_neg(r24814, r24803, MPFR_RNDN); mpfr_pow(r24815, r24800, r24814, MPFR_RNDN); mpfr_div(r24816, r24815, r24805, MPFR_RNDN); mpfr_div(r24817, r24800, r24805, MPFR_RNDN); mpfr_div(r24818, r24817, r24804, MPFR_RNDN); mpfr_mul(r24819, r24818, r24809, MPFR_RNDN); mpfr_add(r24820, r24816, r24819, MPFR_RNDN); mpfr_div(r24821, r24807, r24820, MPFR_RNDN); if (mpfr_get_si(r24802, MPFR_RNDN)) { mpfr_set(r24822, r24813, MPFR_RNDN); } else { mpfr_set(r24822, r24821, MPFR_RNDN); }; return mpfr_get_d(r24822, MPFR_RNDN); } static mpfr_t r24823, r24824, r24825, r24826, r24827, r24828, r24829, r24830, r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838, r24839, r24840, r24841, r24842, r24843, r24844, r24845; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r24823); mpfr_init_set_str(r24824, "6.493800852448035e+60", 10, MPFR_RNDN); mpfr_init(r24825); mpfr_init(r24826); mpfr_init(r24827); mpfr_init(r24828); mpfr_init(r24829); mpfr_init_set_str(r24830, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24831, "10", 10, MPFR_RNDN); mpfr_init(r24832); mpfr_init(r24833); mpfr_init(r24834); mpfr_init(r24835); mpfr_init(r24836); mpfr_init(r24837); mpfr_init(r24838); mpfr_init(r24839); mpfr_init(r24840); mpfr_init(r24841); mpfr_init(r24842); mpfr_init(r24843); mpfr_init(r24844); mpfr_init(r24845); } double f_dm(double a, double k, double m) { mpfr_set_d(r24823, k, MPFR_RNDN); ; mpfr_set_si(r24825, mpfr_cmp(r24823, r24824) <= 0, MPFR_RNDN); mpfr_set_d(r24826, m, MPFR_RNDN); mpfr_pow(r24827, r24823, r24826, MPFR_RNDN); mpfr_set_d(r24828, a, MPFR_RNDN); mpfr_mul(r24829, r24827, r24828, MPFR_RNDN); ; ; mpfr_add(r24832, r24831, r24823, MPFR_RNDN); mpfr_mul(r24833, r24832, r24823, MPFR_RNDN); mpfr_add(r24834, r24830, r24833, MPFR_RNDN); mpfr_div(r24835, r24830, r24834, MPFR_RNDN); mpfr_mul(r24836, r24829, r24835, MPFR_RNDN); mpfr_neg(r24837, r24826, MPFR_RNDN); mpfr_pow(r24838, r24823, r24837, MPFR_RNDN); mpfr_div(r24839, r24838, r24828, MPFR_RNDN); mpfr_div(r24840, r24823, r24828, MPFR_RNDN); mpfr_div(r24841, r24840, r24827, MPFR_RNDN); mpfr_mul(r24842, r24841, r24832, MPFR_RNDN); mpfr_add(r24843, r24839, r24842, MPFR_RNDN); mpfr_div(r24844, r24830, r24843, MPFR_RNDN); if (mpfr_get_si(r24825, MPFR_RNDN)) { mpfr_set(r24845, r24836, MPFR_RNDN); } else { mpfr_set(r24845, r24844, MPFR_RNDN); }; return mpfr_get_d(r24845, MPFR_RNDN); }