#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r24880 = a; float r24881 = k; float r24882 = m; float r24883 = pow(r24881, r24882); float r24884 = r24880 * r24883; float r24885 = 1; float r24886 = 10; float r24887 = r24886 * r24881; float r24888 = r24885 + r24887; float r24889 = r24881 * r24881; float r24890 = r24888 + r24889; float r24891 = r24884 / r24890; return r24891; } double f_id(double a, double k, double m) { double r24892 = a; double r24893 = k; double r24894 = m; double r24895 = pow(r24893, r24894); double r24896 = r24892 * r24895; double r24897 = 1; double r24898 = 10; double r24899 = r24898 * r24893; double r24900 = r24897 + r24899; double r24901 = r24893 * r24893; double r24902 = r24900 + r24901; double r24903 = r24896 / r24902; return r24903; } double f_of(float a, float k, float m) { float r24904 = k; float r24905 = 3.6364048152768866e+66; bool r24906 = r24904 <= r24905; float r24907 = a; float r24908 = m; float r24909 = pow(r24904, r24908); float r24910 = r24907 * r24909; float r24911 = 1; float r24912 = 10; float r24913 = r24912 * r24904; float r24914 = r24911 + r24913; float r24915 = r24904 * r24904; float r24916 = r24914 + r24915; float r24917 = r24910 / r24916; float r24918 = -r24908; float r24919 = pow(r24904, r24918); float r24920 = r24919 / r24907; float r24921 = r24904 / r24907; float r24922 = r24921 / r24909; float r24923 = r24912 + r24904; float r24924 = r24922 * r24923; float r24925 = r24920 + r24924; float r24926 = r24911 / r24925; float r24927 = r24906 ? r24917 : r24926; return r24927; } double f_od(double a, double k, double m) { double r24928 = k; double r24929 = 3.6364048152768866e+66; bool r24930 = r24928 <= r24929; double r24931 = a; double r24932 = m; double r24933 = pow(r24928, r24932); double r24934 = r24931 * r24933; double r24935 = 1; double r24936 = 10; double r24937 = r24936 * r24928; double r24938 = r24935 + r24937; double r24939 = r24928 * r24928; double r24940 = r24938 + r24939; double r24941 = r24934 / r24940; double r24942 = -r24932; double r24943 = pow(r24928, r24942); double r24944 = r24943 / r24931; double r24945 = r24928 / r24931; double r24946 = r24945 / r24933; double r24947 = r24936 + r24928; double r24948 = r24946 * r24947; double r24949 = r24944 + r24948; double r24950 = r24935 / r24949; double r24951 = r24930 ? r24941 : r24950; return r24951; } 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 r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961, r24962, r24963; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); mpfr_init_set_str(r24957, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24958, "10", 10, MPFR_RNDN); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); } double f_im(double a, double k, double m) { mpfr_set_d(r24952, a, MPFR_RNDN); mpfr_set_d(r24953, k, MPFR_RNDN); mpfr_set_d(r24954, m, MPFR_RNDN); mpfr_pow(r24955, r24953, r24954, MPFR_RNDN); mpfr_mul(r24956, r24952, r24955, MPFR_RNDN); ; ; mpfr_mul(r24959, r24958, r24953, MPFR_RNDN); mpfr_add(r24960, r24957, r24959, MPFR_RNDN); mpfr_mul(r24961, r24953, r24953, MPFR_RNDN); mpfr_add(r24962, r24960, r24961, MPFR_RNDN); mpfr_div(r24963, r24956, r24962, MPFR_RNDN); return mpfr_get_d(r24963, MPFR_RNDN); } static mpfr_t r24964, r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r24964); mpfr_init_set_str(r24965, "3.6364048152768866e+66", 10, MPFR_RNDN); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init(r24970); mpfr_init_set_str(r24971, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24972, "10", 10, MPFR_RNDN); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); } double f_fm(double a, double k, double m) { mpfr_set_d(r24964, k, MPFR_RNDN); ; mpfr_set_si(r24966, mpfr_cmp(r24964, r24965) <= 0, MPFR_RNDN); mpfr_set_d(r24967, a, MPFR_RNDN); mpfr_set_d(r24968, m, MPFR_RNDN); mpfr_pow(r24969, r24964, r24968, MPFR_RNDN); mpfr_mul(r24970, r24967, r24969, MPFR_RNDN); ; ; mpfr_mul(r24973, r24972, r24964, MPFR_RNDN); mpfr_add(r24974, r24971, r24973, MPFR_RNDN); mpfr_mul(r24975, r24964, r24964, MPFR_RNDN); mpfr_add(r24976, r24974, r24975, MPFR_RNDN); mpfr_div(r24977, r24970, r24976, MPFR_RNDN); mpfr_neg(r24978, r24968, MPFR_RNDN); mpfr_pow(r24979, r24964, r24978, MPFR_RNDN); mpfr_div(r24980, r24979, r24967, MPFR_RNDN); mpfr_div(r24981, r24964, r24967, MPFR_RNDN); mpfr_div(r24982, r24981, r24969, MPFR_RNDN); mpfr_add(r24983, r24972, r24964, MPFR_RNDN); mpfr_mul(r24984, r24982, r24983, MPFR_RNDN); mpfr_add(r24985, r24980, r24984, MPFR_RNDN); mpfr_div(r24986, r24971, r24985, MPFR_RNDN); if (mpfr_get_si(r24966, MPFR_RNDN)) { mpfr_set(r24987, r24977, MPFR_RNDN); } else { mpfr_set(r24987, r24986, MPFR_RNDN); }; return mpfr_get_d(r24987, MPFR_RNDN); } static mpfr_t r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r24988); mpfr_init_set_str(r24989, "3.6364048152768866e+66", 10, MPFR_RNDN); mpfr_init(r24990); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init_set_str(r24995, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24996, "10", 10, MPFR_RNDN); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); mpfr_init(r25003); mpfr_init(r25004); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); } double f_dm(double a, double k, double m) { mpfr_set_d(r24988, k, MPFR_RNDN); ; mpfr_set_si(r24990, mpfr_cmp(r24988, r24989) <= 0, MPFR_RNDN); mpfr_set_d(r24991, a, MPFR_RNDN); mpfr_set_d(r24992, m, MPFR_RNDN); mpfr_pow(r24993, r24988, r24992, MPFR_RNDN); mpfr_mul(r24994, r24991, r24993, MPFR_RNDN); ; ; mpfr_mul(r24997, r24996, r24988, MPFR_RNDN); mpfr_add(r24998, r24995, r24997, MPFR_RNDN); mpfr_mul(r24999, r24988, r24988, MPFR_RNDN); mpfr_add(r25000, r24998, r24999, MPFR_RNDN); mpfr_div(r25001, r24994, r25000, MPFR_RNDN); mpfr_neg(r25002, r24992, MPFR_RNDN); mpfr_pow(r25003, r24988, r25002, MPFR_RNDN); mpfr_div(r25004, r25003, r24991, MPFR_RNDN); mpfr_div(r25005, r24988, r24991, MPFR_RNDN); mpfr_div(r25006, r25005, r24993, MPFR_RNDN); mpfr_add(r25007, r24996, r24988, MPFR_RNDN); mpfr_mul(r25008, r25006, r25007, MPFR_RNDN); mpfr_add(r25009, r25004, r25008, MPFR_RNDN); mpfr_div(r25010, r24995, r25009, MPFR_RNDN); if (mpfr_get_si(r24990, MPFR_RNDN)) { mpfr_set(r25011, r25001, MPFR_RNDN); } else { mpfr_set(r25011, r25010, MPFR_RNDN); }; return mpfr_get_d(r25011, MPFR_RNDN); }