#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r24849 = a; float r24850 = k; float r24851 = m; float r24852 = pow(r24850, r24851); float r24853 = r24849 * r24852; float r24854 = 1; float r24855 = 10; float r24856 = r24855 * r24850; float r24857 = r24854 + r24856; float r24858 = r24850 * r24850; float r24859 = r24857 + r24858; float r24860 = r24853 / r24859; return r24860; } double f_id(double a, double k, double m) { double r24861 = a; double r24862 = k; double r24863 = m; double r24864 = pow(r24862, r24863); double r24865 = r24861 * r24864; double r24866 = 1; double r24867 = 10; double r24868 = r24867 * r24862; double r24869 = r24866 + r24868; double r24870 = r24862 * r24862; double r24871 = r24869 + r24870; double r24872 = r24865 / r24871; return r24872; } double f_of(float a, float k, float m) { float r24873 = k; float r24874 = 2.2029610971913975e+145; bool r24875 = r24873 <= r24874; float r24876 = a; float r24877 = 10; float r24878 = r24873 + r24877; float r24879 = r24873 * r24878; float r24880 = 1; float r24881 = r24879 + r24880; float r24882 = sqrt(r24881); float r24883 = r24876 / r24882; float r24884 = m; float r24885 = pow(r24873, r24884); float r24886 = r24877 + r24873; float r24887 = r24873 * r24886; float r24888 = r24880 + r24887; float r24889 = sqrt(r24888); float r24890 = r24885 / r24889; float r24891 = r24883 * r24890; float r24892 = -r24884; float r24893 = pow(r24873, r24892); float r24894 = r24893 / r24876; float r24895 = r24873 / r24876; float r24896 = r24895 / r24885; float r24897 = r24896 * r24886; float r24898 = r24894 + r24897; float r24899 = r24880 / r24898; float r24900 = r24875 ? r24891 : r24899; return r24900; } double f_od(double a, double k, double m) { double r24901 = k; double r24902 = 2.2029610971913975e+145; bool r24903 = r24901 <= r24902; double r24904 = a; double r24905 = 10; double r24906 = r24901 + r24905; double r24907 = r24901 * r24906; double r24908 = 1; double r24909 = r24907 + r24908; double r24910 = sqrt(r24909); double r24911 = r24904 / r24910; double r24912 = m; double r24913 = pow(r24901, r24912); double r24914 = r24905 + r24901; double r24915 = r24901 * r24914; double r24916 = r24908 + r24915; double r24917 = sqrt(r24916); double r24918 = r24913 / r24917; double r24919 = r24911 * r24918; double r24920 = -r24912; double r24921 = pow(r24901, r24920); double r24922 = r24921 / r24904; double r24923 = r24901 / r24904; double r24924 = r24923 / r24913; double r24925 = r24924 * r24914; double r24926 = r24922 + r24925; double r24927 = r24908 / r24926; double r24928 = r24903 ? r24919 : r24927; return r24928; } 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 r24929, r24930, r24931, r24932, r24933, r24934, r24935, r24936, r24937, r24938, r24939, r24940; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r24929); mpfr_init(r24930); mpfr_init(r24931); mpfr_init(r24932); mpfr_init(r24933); mpfr_init_set_str(r24934, "1", 10, MPFR_RNDN); mpfr_init_set_str(r24935, "10", 10, MPFR_RNDN); mpfr_init(r24936); mpfr_init(r24937); mpfr_init(r24938); mpfr_init(r24939); mpfr_init(r24940); } double f_im(double a, double k, double m) { mpfr_set_d(r24929, a, MPFR_RNDN); mpfr_set_d(r24930, k, MPFR_RNDN); mpfr_set_d(r24931, m, MPFR_RNDN); mpfr_pow(r24932, r24930, r24931, MPFR_RNDN); mpfr_mul(r24933, r24929, r24932, MPFR_RNDN); ; ; mpfr_mul(r24936, r24935, r24930, MPFR_RNDN); mpfr_add(r24937, r24934, r24936, MPFR_RNDN); mpfr_mul(r24938, r24930, r24930, MPFR_RNDN); mpfr_add(r24939, r24937, r24938, MPFR_RNDN); mpfr_div(r24940, r24933, r24939, MPFR_RNDN); return mpfr_get_d(r24940, MPFR_RNDN); } static mpfr_t r24941, r24942, r24943, r24944, r24945, r24946, r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964, r24965, r24966, r24967, r24968; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r24941); mpfr_init_set_str(r24942, "2.2029610971913975e+145", 10, MPFR_RNDN); mpfr_init(r24943); mpfr_init(r24944); mpfr_init_set_str(r24945, "10", 10, MPFR_RNDN); mpfr_init(r24946); mpfr_init(r24947); mpfr_init_set_str(r24948, "1", 10, MPFR_RNDN); mpfr_init(r24949); mpfr_init(r24950); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); mpfr_init(r24957); mpfr_init(r24958); mpfr_init(r24959); mpfr_init(r24960); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init(r24965); mpfr_init(r24966); mpfr_init(r24967); mpfr_init(r24968); } double f_fm(double a, double k, double m) { mpfr_set_d(r24941, k, MPFR_RNDN); ; mpfr_set_si(r24943, mpfr_cmp(r24941, r24942) <= 0, MPFR_RNDN); mpfr_set_d(r24944, a, MPFR_RNDN); ; mpfr_add(r24946, r24941, r24945, MPFR_RNDN); mpfr_mul(r24947, r24941, r24946, MPFR_RNDN); ; mpfr_add(r24949, r24947, r24948, MPFR_RNDN); mpfr_sqrt(r24950, r24949, MPFR_RNDN); mpfr_div(r24951, r24944, r24950, MPFR_RNDN); mpfr_set_d(r24952, m, MPFR_RNDN); mpfr_pow(r24953, r24941, r24952, MPFR_RNDN); mpfr_add(r24954, r24945, r24941, MPFR_RNDN); mpfr_mul(r24955, r24941, r24954, MPFR_RNDN); mpfr_add(r24956, r24948, r24955, MPFR_RNDN); mpfr_sqrt(r24957, r24956, MPFR_RNDN); mpfr_div(r24958, r24953, r24957, MPFR_RNDN); mpfr_mul(r24959, r24951, r24958, MPFR_RNDN); mpfr_neg(r24960, r24952, MPFR_RNDN); mpfr_pow(r24961, r24941, r24960, MPFR_RNDN); mpfr_div(r24962, r24961, r24944, MPFR_RNDN); mpfr_div(r24963, r24941, r24944, MPFR_RNDN); mpfr_div(r24964, r24963, r24953, MPFR_RNDN); mpfr_mul(r24965, r24964, r24954, MPFR_RNDN); mpfr_add(r24966, r24962, r24965, MPFR_RNDN); mpfr_div(r24967, r24948, r24966, MPFR_RNDN); if (mpfr_get_si(r24943, MPFR_RNDN)) { mpfr_set(r24968, r24959, MPFR_RNDN); } else { mpfr_set(r24968, r24967, MPFR_RNDN); }; return mpfr_get_d(r24968, MPFR_RNDN); } static mpfr_t r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981, r24982, r24983, r24984, r24985, r24986, r24987, r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24969); mpfr_init_set_str(r24970, "2.2029610971913975e+145", 10, MPFR_RNDN); mpfr_init(r24971); mpfr_init(r24972); mpfr_init_set_str(r24973, "10", 10, MPFR_RNDN); mpfr_init(r24974); mpfr_init(r24975); mpfr_init_set_str(r24976, "1", 10, MPFR_RNDN); 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); mpfr_init(r24988); mpfr_init(r24989); mpfr_init(r24990); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); } double f_dm(double a, double k, double m) { mpfr_set_d(r24969, k, MPFR_RNDN); ; mpfr_set_si(r24971, mpfr_cmp(r24969, r24970) <= 0, MPFR_RNDN); mpfr_set_d(r24972, a, MPFR_RNDN); ; mpfr_add(r24974, r24969, r24973, MPFR_RNDN); mpfr_mul(r24975, r24969, r24974, MPFR_RNDN); ; mpfr_add(r24977, r24975, r24976, MPFR_RNDN); mpfr_sqrt(r24978, r24977, MPFR_RNDN); mpfr_div(r24979, r24972, r24978, MPFR_RNDN); mpfr_set_d(r24980, m, MPFR_RNDN); mpfr_pow(r24981, r24969, r24980, MPFR_RNDN); mpfr_add(r24982, r24973, r24969, MPFR_RNDN); mpfr_mul(r24983, r24969, r24982, MPFR_RNDN); mpfr_add(r24984, r24976, r24983, MPFR_RNDN); mpfr_sqrt(r24985, r24984, MPFR_RNDN); mpfr_div(r24986, r24981, r24985, MPFR_RNDN); mpfr_mul(r24987, r24979, r24986, MPFR_RNDN); mpfr_neg(r24988, r24980, MPFR_RNDN); mpfr_pow(r24989, r24969, r24988, MPFR_RNDN); mpfr_div(r24990, r24989, r24972, MPFR_RNDN); mpfr_div(r24991, r24969, r24972, MPFR_RNDN); mpfr_div(r24992, r24991, r24981, MPFR_RNDN); mpfr_mul(r24993, r24992, r24982, MPFR_RNDN); mpfr_add(r24994, r24990, r24993, MPFR_RNDN); mpfr_div(r24995, r24976, r24994, MPFR_RNDN); if (mpfr_get_si(r24971, MPFR_RNDN)) { mpfr_set(r24996, r24987, MPFR_RNDN); } else { mpfr_set(r24996, r24995, MPFR_RNDN); }; return mpfr_get_d(r24996, MPFR_RNDN); }