#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (20:1,3)"; double f_if(float v, float t) { float r24874 = 1; float r24875 = 5; float r24876 = v; float r24877 = r24876 * r24876; float r24878 = r24875 * r24877; float r24879 = r24874 - r24878; float r24880 = atan2(1.0, 0.0); float r24881 = t; float r24882 = r24880 * r24881; float r24883 = 2; float r24884 = 3; float r24885 = r24884 * r24877; float r24886 = r24874 - r24885; float r24887 = r24883 * r24886; float r24888 = sqrt(r24887); float r24889 = r24882 * r24888; float r24890 = r24874 - r24877; float r24891 = r24889 * r24890; float r24892 = r24879 / r24891; return r24892; } double f_id(double v, double t) { double r24893 = 1; double r24894 = 5; double r24895 = v; double r24896 = r24895 * r24895; double r24897 = r24894 * r24896; double r24898 = r24893 - r24897; double r24899 = atan2(1.0, 0.0); double r24900 = t; double r24901 = r24899 * r24900; double r24902 = 2; double r24903 = 3; double r24904 = r24903 * r24896; double r24905 = r24893 - r24904; double r24906 = r24902 * r24905; double r24907 = sqrt(r24906); double r24908 = r24901 * r24907; double r24909 = r24893 - r24896; double r24910 = r24908 * r24909; double r24911 = r24898 / r24910; return r24911; } double f_of(float v, float t) { float r24912 = 1; float r24913 = v; float r24914 = 5; float r24915 = r24913 * r24914; float r24916 = r24913 * r24915; float r24917 = r24912 - r24916; float r24918 = t; float r24919 = r24917 / r24918; float r24920 = atan2(1.0, 0.0); float r24921 = r24913 * r24913; float r24922 = r24921 * r24921; float r24923 = r24912 - r24922; float r24924 = r24920 * r24923; float r24925 = 2; float r24926 = 3; float r24927 = r24926 * r24921; float r24928 = r24912 - r24927; float r24929 = r24925 * r24928; float r24930 = sqrt(r24929); float r24931 = r24924 * r24930; float r24932 = r24919 / r24931; float r24933 = r24912 + r24921; float r24934 = r24932 * r24933; return r24934; } double f_od(double v, double t) { double r24935 = 1; double r24936 = v; double r24937 = 5; double r24938 = r24936 * r24937; double r24939 = r24936 * r24938; double r24940 = r24935 - r24939; double r24941 = t; double r24942 = r24940 / r24941; double r24943 = atan2(1.0, 0.0); double r24944 = r24936 * r24936; double r24945 = r24944 * r24944; double r24946 = r24935 - r24945; double r24947 = r24943 * r24946; double r24948 = 2; double r24949 = 3; double r24950 = r24949 * r24944; double r24951 = r24935 - r24950; double r24952 = r24948 * r24951; double r24953 = sqrt(r24952); double r24954 = r24947 * r24953; double r24955 = r24942 / r24954; double r24956 = r24935 + r24944; double r24957 = r24955 * r24956; return r24957; } 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 r24958, r24959, r24960, r24961, r24962, r24963, r24964, r24965, r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24958, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24959, "5", 10, MPFR_RNDN); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init(r24965); mpfr_init(r24966); mpfr_init_set_str(r24967, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24968, "3", 10, MPFR_RNDN); mpfr_init(r24969); mpfr_init(r24970); mpfr_init(r24971); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); } double f_im(double v, double t) { ; ; mpfr_set_d(r24960, v, MPFR_RNDN); mpfr_mul(r24961, r24960, r24960, MPFR_RNDN); mpfr_mul(r24962, r24959, r24961, MPFR_RNDN); mpfr_sub(r24963, r24958, r24962, MPFR_RNDN); mpfr_const_pi(r24964, MPFR_RNDN); mpfr_set_d(r24965, t, MPFR_RNDN); mpfr_mul(r24966, r24964, r24965, MPFR_RNDN); ; ; mpfr_mul(r24969, r24968, r24961, MPFR_RNDN); mpfr_sub(r24970, r24958, r24969, MPFR_RNDN); mpfr_mul(r24971, r24967, r24970, MPFR_RNDN); mpfr_sqrt(r24972, r24971, MPFR_RNDN); mpfr_mul(r24973, r24966, r24972, MPFR_RNDN); mpfr_sub(r24974, r24958, r24961, MPFR_RNDN); mpfr_mul(r24975, r24973, r24974, MPFR_RNDN); mpfr_div(r24976, r24963, r24975, MPFR_RNDN); return mpfr_get_d(r24976, MPFR_RNDN); } static mpfr_t r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987, r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997, r24998, r24999; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24977, "1", 10, MPFR_RNDN); mpfr_init(r24978); mpfr_init_set_str(r24979, "5", 10, MPFR_RNDN); mpfr_init(r24980); mpfr_init(r24981); mpfr_init(r24982); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); mpfr_init(r24988); mpfr_init(r24989); mpfr_init_set_str(r24990, "2", 10, MPFR_RNDN); mpfr_init_set_str(r24991, "3", 10, MPFR_RNDN); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); mpfr_init(r24998); mpfr_init(r24999); } double f_fm(double v, double t) { ; mpfr_set_d(r24978, v, MPFR_RNDN); ; mpfr_mul(r24980, r24978, r24979, MPFR_RNDN); mpfr_mul(r24981, r24978, r24980, MPFR_RNDN); mpfr_sub(r24982, r24977, r24981, MPFR_RNDN); mpfr_set_d(r24983, t, MPFR_RNDN); mpfr_div(r24984, r24982, r24983, MPFR_RNDN); mpfr_const_pi(r24985, MPFR_RNDN); mpfr_mul(r24986, r24978, r24978, MPFR_RNDN); mpfr_mul(r24987, r24986, r24986, MPFR_RNDN); mpfr_sub(r24988, r24977, r24987, MPFR_RNDN); mpfr_mul(r24989, r24985, r24988, MPFR_RNDN); ; ; mpfr_mul(r24992, r24991, r24986, MPFR_RNDN); mpfr_sub(r24993, r24977, r24992, MPFR_RNDN); mpfr_mul(r24994, r24990, r24993, MPFR_RNDN); mpfr_sqrt(r24995, r24994, MPFR_RNDN); mpfr_mul(r24996, r24989, r24995, MPFR_RNDN); mpfr_div(r24997, r24984, r24996, MPFR_RNDN); mpfr_add(r24998, r24977, r24986, MPFR_RNDN); mpfr_mul(r24999, r24997, r24998, MPFR_RNDN); return mpfr_get_d(r24999, MPFR_RNDN); } static mpfr_t r25000, r25001, r25002, r25003, r25004, r25005, 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(400); mpfr_init_set_str(r25000, "1", 10, MPFR_RNDN); mpfr_init(r25001); mpfr_init_set_str(r25002, "5", 10, MPFR_RNDN); 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); mpfr_init(r25012); mpfr_init_set_str(r25013, "2", 10, MPFR_RNDN); mpfr_init_set_str(r25014, "3", 10, MPFR_RNDN); 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, double t) { ; mpfr_set_d(r25001, v, MPFR_RNDN); ; mpfr_mul(r25003, r25001, r25002, MPFR_RNDN); mpfr_mul(r25004, r25001, r25003, MPFR_RNDN); mpfr_sub(r25005, r25000, r25004, MPFR_RNDN); mpfr_set_d(r25006, t, MPFR_RNDN); mpfr_div(r25007, r25005, r25006, MPFR_RNDN); mpfr_const_pi(r25008, MPFR_RNDN); mpfr_mul(r25009, r25001, r25001, MPFR_RNDN); mpfr_mul(r25010, r25009, r25009, MPFR_RNDN); mpfr_sub(r25011, r25000, r25010, MPFR_RNDN); mpfr_mul(r25012, r25008, r25011, MPFR_RNDN); ; ; mpfr_mul(r25015, r25014, r25009, MPFR_RNDN); mpfr_sub(r25016, r25000, r25015, MPFR_RNDN); mpfr_mul(r25017, r25013, r25016, MPFR_RNDN); mpfr_sqrt(r25018, r25017, MPFR_RNDN); mpfr_mul(r25019, r25012, r25018, MPFR_RNDN); mpfr_div(r25020, r25007, r25019, MPFR_RNDN); mpfr_add(r25021, r25000, r25009, MPFR_RNDN); mpfr_mul(r25022, r25020, r25021, MPFR_RNDN); return mpfr_get_d(r25022, MPFR_RNDN); }