#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r24937 = a; float r24938 = k; float r24939 = m; float r24940 = pow(r24938, r24939); float r24941 = r24937 * r24940; float r24942 = 1; float r24943 = 10; float r24944 = r24943 * r24938; float r24945 = r24942 + r24944; float r24946 = r24938 * r24938; float r24947 = r24945 + r24946; float r24948 = r24941 / r24947; return r24948; } double f_id(double a, double k, double m) { double r24949 = a; double r24950 = k; double r24951 = m; double r24952 = pow(r24950, r24951); double r24953 = r24949 * r24952; double r24954 = 1; double r24955 = 10; double r24956 = r24955 * r24950; double r24957 = r24954 + r24956; double r24958 = r24950 * r24950; double r24959 = r24957 + r24958; double r24960 = r24953 / r24959; return r24960; } double f_of(float a, float k, float m) { float r24961 = k; float r24962 = 3.8534943108273086e+128; bool r24963 = r24961 <= r24962; float r24964 = m; float r24965 = pow(r24961, r24964); float r24966 = a; float r24967 = r24965 * r24966; float r24968 = 1; float r24969 = 10; float r24970 = r24969 + r24961; float r24971 = r24970 * r24961; float r24972 = r24968 + r24971; float r24973 = r24967 / r24972; float r24974 = -r24964; float r24975 = pow(r24961, r24974); float r24976 = r24975 / r24966; float r24977 = r24961 / r24966; float r24978 = r24977 / r24965; float r24979 = r24978 * r24970; float r24980 = r24976 + r24979; float r24981 = r24968 / r24980; float r24982 = r24963 ? r24973 : r24981; return r24982; } double f_od(double a, double k, double m) { double r24983 = k; double r24984 = 3.8534943108273086e+128; bool r24985 = r24983 <= r24984; double r24986 = m; double r24987 = pow(r24983, r24986); double r24988 = a; double r24989 = r24987 * r24988; double r24990 = 1; double r24991 = 10; double r24992 = r24991 + r24983; double r24993 = r24992 * r24983; double r24994 = r24990 + r24993; double r24995 = r24989 / r24994; double r24996 = -r24986; double r24997 = pow(r24983, r24996); double r24998 = r24997 / r24988; double r24999 = r24983 / r24988; double r25000 = r24999 / r24987; double r25001 = r25000 * r24992; double r25002 = r24998 + r25001; double r25003 = r24990 / r25002; double r25004 = r24985 ? r24995 : r25003; return r25004; } 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 r25005, r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init_set_str(r25010, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25011, "10", 10, MPFR_RNDN); mpfr_init(r25012); mpfr_init(r25013); mpfr_init(r25014); mpfr_init(r25015); mpfr_init(r25016); } double f_im(double a, double k, double m) { mpfr_set_d(r25005, a, MPFR_RNDN); mpfr_set_d(r25006, k, MPFR_RNDN); mpfr_set_d(r25007, m, MPFR_RNDN); mpfr_pow(r25008, r25006, r25007, MPFR_RNDN); mpfr_mul(r25009, r25005, r25008, MPFR_RNDN); ; ; mpfr_mul(r25012, r25011, r25006, MPFR_RNDN); mpfr_add(r25013, r25010, r25012, MPFR_RNDN); mpfr_mul(r25014, r25006, r25006, MPFR_RNDN); mpfr_add(r25015, r25013, r25014, MPFR_RNDN); mpfr_div(r25016, r25009, r25015, MPFR_RNDN); return mpfr_get_d(r25016, MPFR_RNDN); } static mpfr_t r25017, r25018, r25019, r25020, r25021, r25022, r25023, r25024, r25025, r25026, r25027, r25028, r25029, r25030, r25031, r25032, r25033, r25034, r25035, r25036, r25037, r25038; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25017); mpfr_init_set_str(r25018, "3.8534943108273086e+128", 10, MPFR_RNDN); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); mpfr_init(r25023); mpfr_init_set_str(r25024, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25025, "10", 10, MPFR_RNDN); mpfr_init(r25026); mpfr_init(r25027); mpfr_init(r25028); mpfr_init(r25029); mpfr_init(r25030); mpfr_init(r25031); mpfr_init(r25032); mpfr_init(r25033); mpfr_init(r25034); mpfr_init(r25035); mpfr_init(r25036); mpfr_init(r25037); mpfr_init(r25038); } double f_fm(double a, double k, double m) { mpfr_set_d(r25017, k, MPFR_RNDN); ; mpfr_set_si(r25019, mpfr_cmp(r25017, r25018) <= 0, MPFR_RNDN); mpfr_set_d(r25020, m, MPFR_RNDN); mpfr_pow(r25021, r25017, r25020, MPFR_RNDN); mpfr_set_d(r25022, a, MPFR_RNDN); mpfr_mul(r25023, r25021, r25022, MPFR_RNDN); ; ; mpfr_add(r25026, r25025, r25017, MPFR_RNDN); mpfr_mul(r25027, r25026, r25017, MPFR_RNDN); mpfr_add(r25028, r25024, r25027, MPFR_RNDN); mpfr_div(r25029, r25023, r25028, MPFR_RNDN); mpfr_neg(r25030, r25020, MPFR_RNDN); mpfr_pow(r25031, r25017, r25030, MPFR_RNDN); mpfr_div(r25032, r25031, r25022, MPFR_RNDN); mpfr_div(r25033, r25017, r25022, MPFR_RNDN); mpfr_div(r25034, r25033, r25021, MPFR_RNDN); mpfr_mul(r25035, r25034, r25026, MPFR_RNDN); mpfr_add(r25036, r25032, r25035, MPFR_RNDN); mpfr_div(r25037, r25024, r25036, MPFR_RNDN); if (mpfr_get_si(r25019, MPFR_RNDN)) { mpfr_set(r25038, r25029, MPFR_RNDN); } else { mpfr_set(r25038, r25037, MPFR_RNDN); }; return mpfr_get_d(r25038, MPFR_RNDN); } static mpfr_t r25039, r25040, r25041, r25042, r25043, r25044, r25045, r25046, r25047, r25048, r25049, r25050, r25051, r25052, r25053, r25054, r25055, r25056, r25057, r25058, r25059, r25060; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25039); mpfr_init_set_str(r25040, "3.8534943108273086e+128", 10, MPFR_RNDN); mpfr_init(r25041); mpfr_init(r25042); mpfr_init(r25043); mpfr_init(r25044); mpfr_init(r25045); mpfr_init_set_str(r25046, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25047, "10", 10, MPFR_RNDN); mpfr_init(r25048); mpfr_init(r25049); mpfr_init(r25050); mpfr_init(r25051); mpfr_init(r25052); mpfr_init(r25053); mpfr_init(r25054); mpfr_init(r25055); mpfr_init(r25056); mpfr_init(r25057); mpfr_init(r25058); mpfr_init(r25059); mpfr_init(r25060); } double f_dm(double a, double k, double m) { mpfr_set_d(r25039, k, MPFR_RNDN); ; mpfr_set_si(r25041, mpfr_cmp(r25039, r25040) <= 0, MPFR_RNDN); mpfr_set_d(r25042, m, MPFR_RNDN); mpfr_pow(r25043, r25039, r25042, MPFR_RNDN); mpfr_set_d(r25044, a, MPFR_RNDN); mpfr_mul(r25045, r25043, r25044, MPFR_RNDN); ; ; mpfr_add(r25048, r25047, r25039, MPFR_RNDN); mpfr_mul(r25049, r25048, r25039, MPFR_RNDN); mpfr_add(r25050, r25046, r25049, MPFR_RNDN); mpfr_div(r25051, r25045, r25050, MPFR_RNDN); mpfr_neg(r25052, r25042, MPFR_RNDN); mpfr_pow(r25053, r25039, r25052, MPFR_RNDN); mpfr_div(r25054, r25053, r25044, MPFR_RNDN); mpfr_div(r25055, r25039, r25044, MPFR_RNDN); mpfr_div(r25056, r25055, r25043, MPFR_RNDN); mpfr_mul(r25057, r25056, r25048, MPFR_RNDN); mpfr_add(r25058, r25054, r25057, MPFR_RNDN); mpfr_div(r25059, r25046, r25058, MPFR_RNDN); if (mpfr_get_si(r25041, MPFR_RNDN)) { mpfr_set(r25060, r25051, MPFR_RNDN); } else { mpfr_set(r25060, r25059, MPFR_RNDN); }; return mpfr_get_d(r25060, MPFR_RNDN); }