#include #include #include #include #include char *name = "Falkner and Boettcher, Appendix A"; double f_if(float a, float k, float m) { float r10887 = a; float r10888 = k; float r10889 = m; float r10890 = pow(r10888, r10889); float r10891 = r10887 * r10890; float r10892 = 1.0f; float r10893 = 10.0f; float r10894 = r10893 * r10888; float r10895 = r10892 + r10894; float r10896 = r10888 * r10888; float r10897 = r10895 + r10896; float r10898 = r10891 / r10897; return r10898; } double f_id(double a, double k, double m) { double r10899 = a; double r10900 = k; double r10901 = m; double r10902 = pow(r10900, r10901); double r10903 = r10899 * r10902; double r10904 = 1.0; double r10905 = 10.0; double r10906 = r10905 * r10900; double r10907 = r10904 + r10906; double r10908 = r10900 * r10900; double r10909 = r10907 + r10908; double r10910 = r10903 / r10909; return r10910; } double f_of(float a, float k, float m) { float r10911 = k; float r10912 = 2.439311431159999e+130f; bool r10913 = r10911 <= r10912; float r10914 = a; float r10915 = m; float r10916 = pow(r10911, r10915); float r10917 = r10914 * r10916; float r10918 = r10911 * r10911; float r10919 = 10.0f; float r10920 = r10911 * r10919; float r10921 = 1.0f; float r10922 = r10920 + r10921; float r10923 = r10918 + r10922; float r10924 = r10917 / r10923; float r10925 = exp(r10915); float r10926 = log(r10911); float r10927 = pow(r10925, r10926); float r10928 = 4.0f; float r10929 = pow(r10911, r10928); float r10930 = r10914 / r10929; float r10931 = r10927 * r10930; float r10932 = 99.0f; float r10933 = r10914 / r10911; float r10934 = r10919 / r10911; float r10935 = r10934 * r10933; float r10936 = r10933 - r10935; float r10937 = r10927 / r10911; float r10938 = r10936 * r10937; float r10939 = fma(r10931, r10932, r10938); float r10940 = r10913 ? r10924 : r10939; return r10940; } double f_od(double a, double k, double m) { double r10941 = k; double r10942 = 2.439311431159999e+130; bool r10943 = r10941 <= r10942; double r10944 = a; double r10945 = m; double r10946 = pow(r10941, r10945); double r10947 = r10944 * r10946; double r10948 = r10941 * r10941; double r10949 = 10.0; double r10950 = r10941 * r10949; double r10951 = 1.0; double r10952 = r10950 + r10951; double r10953 = r10948 + r10952; double r10954 = r10947 / r10953; double r10955 = exp(r10945); double r10956 = log(r10941); double r10957 = pow(r10955, r10956); double r10958 = 4.0; double r10959 = pow(r10941, r10958); double r10960 = r10944 / r10959; double r10961 = r10957 * r10960; double r10962 = 99.0; double r10963 = r10944 / r10941; double r10964 = r10949 / r10941; double r10965 = r10964 * r10963; double r10966 = r10963 - r10965; double r10967 = r10957 / r10941; double r10968 = r10966 * r10967; double r10969 = fma(r10961, r10962, r10968); double r10970 = r10943 ? r10954 : r10969; return r10970; } 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 r10971, r10972, r10973, r10974, r10975, r10976, r10977, r10978, r10979, r10980, r10981, r10982; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r10971); mpfr_init(r10972); mpfr_init(r10973); mpfr_init(r10974); mpfr_init(r10975); mpfr_init_set_str(r10976, "1", 10, MPFR_RNDN); mpfr_init_set_str(r10977, "10", 10, MPFR_RNDN); mpfr_init(r10978); mpfr_init(r10979); mpfr_init(r10980); mpfr_init(r10981); mpfr_init(r10982); } double f_im(double a, double k, double m) { mpfr_set_d(r10971, a, MPFR_RNDN); mpfr_set_d(r10972, k, MPFR_RNDN); mpfr_set_d(r10973, m, MPFR_RNDN); mpfr_pow(r10974, r10972, r10973, MPFR_RNDN); mpfr_mul(r10975, r10971, r10974, MPFR_RNDN); ; ; mpfr_mul(r10978, r10977, r10972, MPFR_RNDN); mpfr_add(r10979, r10976, r10978, MPFR_RNDN); mpfr_mul(r10980, r10972, r10972, MPFR_RNDN); mpfr_add(r10981, r10979, r10980, MPFR_RNDN); mpfr_div(r10982, r10975, r10981, MPFR_RNDN); return mpfr_get_d(r10982, MPFR_RNDN); } static mpfr_t r10983, r10984, r10985, r10986, r10987, r10988, r10989, r10990, r10991, r10992, r10993, r10994, r10995, r10996, r10997, r10998, r10999, r11000, r11001, r11002, r11003, r11004, r11005, r11006, r11007, r11008, r11009, r11010, r11011, r11012; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r10983); mpfr_init_set_str(r10984, "2.439311431159999e+130", 10, MPFR_RNDN); mpfr_init(r10985); mpfr_init(r10986); mpfr_init(r10987); mpfr_init(r10988); mpfr_init(r10989); mpfr_init(r10990); mpfr_init_set_str(r10991, "10", 10, MPFR_RNDN); mpfr_init(r10992); mpfr_init_set_str(r10993, "1", 10, MPFR_RNDN); mpfr_init(r10994); mpfr_init(r10995); mpfr_init(r10996); mpfr_init(r10997); mpfr_init(r10998); mpfr_init(r10999); mpfr_init_set_str(r11000, "4", 10, MPFR_RNDN); mpfr_init(r11001); mpfr_init(r11002); mpfr_init(r11003); mpfr_init_set_str(r11004, "99", 10, MPFR_RNDN); mpfr_init(r11005); mpfr_init(r11006); mpfr_init(r11007); mpfr_init(r11008); mpfr_init(r11009); mpfr_init(r11010); mpfr_init(r11011); mpfr_init(r11012); } double f_fm(double a, double k, double m) { mpfr_set_d(r10983, k, MPFR_RNDN); ; mpfr_set_si(r10985, mpfr_cmp(r10983, r10984) <= 0, MPFR_RNDN); mpfr_set_d(r10986, a, MPFR_RNDN); mpfr_set_d(r10987, m, MPFR_RNDN); mpfr_pow(r10988, r10983, r10987, MPFR_RNDN); mpfr_mul(r10989, r10986, r10988, MPFR_RNDN); mpfr_mul(r10990, r10983, r10983, MPFR_RNDN); ; mpfr_mul(r10992, r10983, r10991, MPFR_RNDN); ; mpfr_add(r10994, r10992, r10993, MPFR_RNDN); mpfr_add(r10995, r10990, r10994, MPFR_RNDN); mpfr_div(r10996, r10989, r10995, MPFR_RNDN); mpfr_exp(r10997, r10987, MPFR_RNDN); mpfr_log(r10998, r10983, MPFR_RNDN); mpfr_pow(r10999, r10997, r10998, MPFR_RNDN); ; mpfr_pow(r11001, r10983, r11000, MPFR_RNDN); mpfr_div(r11002, r10986, r11001, MPFR_RNDN); mpfr_mul(r11003, r10999, r11002, MPFR_RNDN); ; mpfr_div(r11005, r10986, r10983, MPFR_RNDN); mpfr_div(r11006, r10991, r10983, MPFR_RNDN); mpfr_mul(r11007, r11006, r11005, MPFR_RNDN); mpfr_sub(r11008, r11005, r11007, MPFR_RNDN); mpfr_div(r11009, r10999, r10983, MPFR_RNDN); mpfr_mul(r11010, r11008, r11009, MPFR_RNDN); mpfr_fma(r11011, r11003, r11004, r11010, MPFR_RNDN); if (mpfr_get_si(r10985, MPFR_RNDN)) { mpfr_set(r11012, r10996, MPFR_RNDN); } else { mpfr_set(r11012, r11011, MPFR_RNDN); }; return mpfr_get_d(r11012, MPFR_RNDN); } static mpfr_t r11013, r11014, r11015, r11016, r11017, r11018, r11019, r11020, r11021, r11022, r11023, r11024, r11025, r11026, r11027, r11028, r11029, r11030, r11031, r11032, r11033, r11034, r11035, r11036, r11037, r11038, r11039, r11040, r11041, r11042; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r11013); mpfr_init_set_str(r11014, "2.439311431159999e+130", 10, MPFR_RNDN); mpfr_init(r11015); mpfr_init(r11016); mpfr_init(r11017); mpfr_init(r11018); mpfr_init(r11019); mpfr_init(r11020); mpfr_init_set_str(r11021, "10", 10, MPFR_RNDN); mpfr_init(r11022); mpfr_init_set_str(r11023, "1", 10, MPFR_RNDN); mpfr_init(r11024); mpfr_init(r11025); mpfr_init(r11026); mpfr_init(r11027); mpfr_init(r11028); mpfr_init(r11029); mpfr_init_set_str(r11030, "4", 10, MPFR_RNDN); mpfr_init(r11031); mpfr_init(r11032); mpfr_init(r11033); mpfr_init_set_str(r11034, "99", 10, MPFR_RNDN); mpfr_init(r11035); mpfr_init(r11036); mpfr_init(r11037); mpfr_init(r11038); mpfr_init(r11039); mpfr_init(r11040); mpfr_init(r11041); mpfr_init(r11042); } double f_dm(double a, double k, double m) { mpfr_set_d(r11013, k, MPFR_RNDN); ; mpfr_set_si(r11015, mpfr_cmp(r11013, r11014) <= 0, MPFR_RNDN); mpfr_set_d(r11016, a, MPFR_RNDN); mpfr_set_d(r11017, m, MPFR_RNDN); mpfr_pow(r11018, r11013, r11017, MPFR_RNDN); mpfr_mul(r11019, r11016, r11018, MPFR_RNDN); mpfr_mul(r11020, r11013, r11013, MPFR_RNDN); ; mpfr_mul(r11022, r11013, r11021, MPFR_RNDN); ; mpfr_add(r11024, r11022, r11023, MPFR_RNDN); mpfr_add(r11025, r11020, r11024, MPFR_RNDN); mpfr_div(r11026, r11019, r11025, MPFR_RNDN); mpfr_exp(r11027, r11017, MPFR_RNDN); mpfr_log(r11028, r11013, MPFR_RNDN); mpfr_pow(r11029, r11027, r11028, MPFR_RNDN); ; mpfr_pow(r11031, r11013, r11030, MPFR_RNDN); mpfr_div(r11032, r11016, r11031, MPFR_RNDN); mpfr_mul(r11033, r11029, r11032, MPFR_RNDN); ; mpfr_div(r11035, r11016, r11013, MPFR_RNDN); mpfr_div(r11036, r11021, r11013, MPFR_RNDN); mpfr_mul(r11037, r11036, r11035, MPFR_RNDN); mpfr_sub(r11038, r11035, r11037, MPFR_RNDN); mpfr_div(r11039, r11029, r11013, MPFR_RNDN); mpfr_mul(r11040, r11038, r11039, MPFR_RNDN); mpfr_fma(r11041, r11033, r11034, r11040, MPFR_RNDN); if (mpfr_get_si(r11015, MPFR_RNDN)) { mpfr_set(r11042, r11026, MPFR_RNDN); } else { mpfr_set(r11042, r11041, MPFR_RNDN); }; return mpfr_get_d(r11042, MPFR_RNDN); }