#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix B, 2"; double f_if(float v) { float r24913 = 2; float r24914 = sqrt(r24913); float r24915 = 4; float r24916 = r24914 / r24915; float r24917 = 1; float r24918 = 3; float r24919 = v; float r24920 = r24919 * r24919; float r24921 = r24918 * r24920; float r24922 = r24917 - r24921; float r24923 = sqrt(r24922); float r24924 = r24916 * r24923; float r24925 = r24917 - r24920; float r24926 = r24924 * r24925; return r24926; } double f_id(double v) { double r24927 = 2; double r24928 = sqrt(r24927); double r24929 = 4; double r24930 = r24928 / r24929; double r24931 = 1; double r24932 = 3; double r24933 = v; double r24934 = r24933 * r24933; double r24935 = r24932 * r24934; double r24936 = r24931 - r24935; double r24937 = sqrt(r24936); double r24938 = r24930 * r24937; double r24939 = r24931 - r24934; double r24940 = r24938 * r24939; return r24940; } double f_of(float v) { float r24941 = 2; float r24942 = sqrt(r24941); float r24943 = 4; float r24944 = r24942 / r24943; float r24945 = log(r24944); float r24946 = 1; float r24947 = 3; float r24948 = v; float r24949 = r24948 * r24948; float r24950 = r24947 * r24949; float r24951 = r24946 - r24950; float r24952 = sqrt(r24951); float r24953 = log(r24952); float r24954 = r24945 + r24953; float r24955 = exp(r24954); float r24956 = r24946 - r24949; float r24957 = r24955 * r24956; return r24957; } double f_od(double v) { double r24958 = 2; double r24959 = sqrt(r24958); double r24960 = 4; double r24961 = r24959 / r24960; double r24962 = log(r24961); double r24963 = 1; double r24964 = 3; double r24965 = v; double r24966 = r24965 * r24965; double r24967 = r24964 * r24966; double r24968 = r24963 - r24967; double r24969 = sqrt(r24968); double r24970 = log(r24969); double r24971 = r24962 + r24970; double r24972 = exp(r24971); double r24973 = r24963 - r24966; double r24974 = r24972 * r24973; return r24974; } 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 r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987, r24988; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r24975, "2", 10, MPFR_RNDN); mpfr_init(r24976); mpfr_init_set_str(r24977, "4", 10, MPFR_RNDN); mpfr_init(r24978); mpfr_init_set_str(r24979, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24980, "3", 10, MPFR_RNDN); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); mpfr_init(r24988); } double f_im(double v) { ; mpfr_sqrt(r24976, r24975, MPFR_RNDN); ; mpfr_div(r24978, r24976, r24977, MPFR_RNDN); ; ; mpfr_set_d(r24981, v, MPFR_RNDN); mpfr_mul(r24982, r24981, r24981, MPFR_RNDN); mpfr_mul(r24983, r24980, r24982, MPFR_RNDN); mpfr_sub(r24984, r24979, r24983, MPFR_RNDN); mpfr_sqrt(r24985, r24984, MPFR_RNDN); mpfr_mul(r24986, r24978, r24985, MPFR_RNDN); mpfr_sub(r24987, r24979, r24982, MPFR_RNDN); mpfr_mul(r24988, r24986, r24987, MPFR_RNDN); return mpfr_get_d(r24988, MPFR_RNDN); } static mpfr_t r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r24989, "2", 10, MPFR_RNDN); mpfr_init(r24990); mpfr_init_set_str(r24991, "4", 10, MPFR_RNDN); mpfr_init(r24992); mpfr_init(r24993); mpfr_init_set_str(r24994, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24995, "3", 10, MPFR_RNDN); mpfr_init(r24996); 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); } double f_fm(double v) { ; mpfr_sqrt(r24990, r24989, MPFR_RNDN); ; mpfr_div(r24992, r24990, r24991, MPFR_RNDN); mpfr_log(r24993, r24992, MPFR_RNDN); ; ; mpfr_set_d(r24996, v, MPFR_RNDN); mpfr_mul(r24997, r24996, r24996, MPFR_RNDN); mpfr_mul(r24998, r24995, r24997, MPFR_RNDN); mpfr_sub(r24999, r24994, r24998, MPFR_RNDN); mpfr_sqrt(r25000, r24999, MPFR_RNDN); mpfr_log(r25001, r25000, MPFR_RNDN); mpfr_add(r25002, r24993, r25001, MPFR_RNDN); mpfr_exp(r25003, r25002, MPFR_RNDN); mpfr_sub(r25004, r24994, r24997, MPFR_RNDN); mpfr_mul(r25005, r25003, r25004, MPFR_RNDN); return mpfr_get_d(r25005, MPFR_RNDN); } static mpfr_t r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013, r25014, r25015, r25016, r25017, r25018, r25019, r25020, r25021, r25022; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r25006, "2", 10, MPFR_RNDN); mpfr_init(r25007); mpfr_init_set_str(r25008, "4", 10, MPFR_RNDN); mpfr_init(r25009); mpfr_init(r25010); mpfr_init_set_str(r25011, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25012, "3", 10, MPFR_RNDN); mpfr_init(r25013); mpfr_init(r25014); mpfr_init(r25015); mpfr_init(r25016); mpfr_init(r25017); mpfr_init(r25018); mpfr_init(r25019); mpfr_init(r25020); mpfr_init(r25021); mpfr_init(r25022); } double f_dm(double v) { ; mpfr_sqrt(r25007, r25006, MPFR_RNDN); ; mpfr_div(r25009, r25007, r25008, MPFR_RNDN); mpfr_log(r25010, r25009, MPFR_RNDN); ; ; mpfr_set_d(r25013, v, MPFR_RNDN); mpfr_mul(r25014, r25013, r25013, MPFR_RNDN); mpfr_mul(r25015, r25012, r25014, MPFR_RNDN); mpfr_sub(r25016, r25011, r25015, MPFR_RNDN); mpfr_sqrt(r25017, r25016, MPFR_RNDN); mpfr_log(r25018, r25017, MPFR_RNDN); mpfr_add(r25019, r25010, r25018, MPFR_RNDN); mpfr_exp(r25020, r25019, MPFR_RNDN); mpfr_sub(r25021, r25011, r25014, MPFR_RNDN); mpfr_mul(r25022, r25020, r25021, MPFR_RNDN); return mpfr_get_d(r25022, MPFR_RNDN); }