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