#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r21878 = b; float r21879 = -r21878; float r21880 = r21878 * r21878; float r21881 = 4; float r21882 = a; float r21883 = r21881 * r21882; float r21884 = c; float r21885 = r21883 * r21884; float r21886 = r21880 - r21885; float r21887 = sqrt(r21886); float r21888 = r21879 + r21887; float r21889 = 2; float r21890 = r21889 * r21882; float r21891 = r21888 / r21890; return r21891; } double f_id(double a, double b, double c) { double r21892 = b; double r21893 = -r21892; double r21894 = r21892 * r21892; double r21895 = 4; double r21896 = a; double r21897 = r21895 * r21896; double r21898 = c; double r21899 = r21897 * r21898; double r21900 = r21894 - r21899; double r21901 = sqrt(r21900); double r21902 = r21893 + r21901; double r21903 = 2; double r21904 = r21903 * r21896; double r21905 = r21902 / r21904; return r21905; } double f_of(float a, float b, float c) { float r21906 = c; float r21907 = a; float r21908 = r21906 * r21907; float r21909 = 4; float r21910 = -r21909; float r21911 = r21908 * r21910; float r21912 = r21909 * r21907; float r21913 = -r21906; float r21914 = b; float r21915 = r21914 * r21914; float r21916 = fma(r21912, r21913, r21915); float r21917 = sqrt(r21916); float r21918 = 3; float r21919 = pow(r21917, r21918); float r21920 = pow(r21914, r21918); float r21921 = r21919 + r21920; float r21922 = r21911 / r21921; float r21923 = 2; float r21924 = r21907 * r21923; float r21925 = -r21914; float r21926 = r21914 + r21914; float r21927 = r21926 * r21914; float r21928 = r21909 * r21908; float r21929 = r21927 - r21928; float r21930 = fma(r21925, r21917, r21929); float r21931 = r21924 / r21930; float r21932 = r21922 / r21931; return r21932; } double f_od(double a, double b, double c) { double r21933 = c; double r21934 = a; double r21935 = r21933 * r21934; double r21936 = 4; double r21937 = -r21936; double r21938 = r21935 * r21937; double r21939 = r21936 * r21934; double r21940 = -r21933; double r21941 = b; double r21942 = r21941 * r21941; double r21943 = fma(r21939, r21940, r21942); double r21944 = sqrt(r21943); double r21945 = 3; double r21946 = pow(r21944, r21945); double r21947 = pow(r21941, r21945); double r21948 = r21946 + r21947; double r21949 = r21938 / r21948; double r21950 = 2; double r21951 = r21934 * r21950; double r21952 = -r21941; double r21953 = r21941 + r21941; double r21954 = r21953 * r21941; double r21955 = r21936 * r21935; double r21956 = r21954 - r21955; double r21957 = fma(r21952, r21944, r21956); double r21958 = r21951 / r21957; double r21959 = r21949 / r21958; return r21959; } 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 r21960, r21961, r21962, r21963, r21964, r21965, r21966, r21967, r21968, r21969, r21970, r21971, r21972, r21973; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r21960); mpfr_init(r21961); mpfr_init(r21962); mpfr_init_set_str(r21963, "4", 10, MPFR_RNDN); mpfr_init(r21964); mpfr_init(r21965); mpfr_init(r21966); mpfr_init(r21967); mpfr_init(r21968); mpfr_init(r21969); mpfr_init(r21970); mpfr_init_set_str(r21971, "2", 10, MPFR_RNDN); mpfr_init(r21972); mpfr_init(r21973); } double f_im(double a, double b, double c) { mpfr_set_d(r21960, b, MPFR_RNDN); mpfr_neg(r21961, r21960, MPFR_RNDN); mpfr_mul(r21962, r21960, r21960, MPFR_RNDN); ; mpfr_set_d(r21964, a, MPFR_RNDN); mpfr_mul(r21965, r21963, r21964, MPFR_RNDN); mpfr_set_d(r21966, c, MPFR_RNDN); mpfr_mul(r21967, r21965, r21966, MPFR_RNDN); mpfr_sub(r21968, r21962, r21967, MPFR_RNDN); mpfr_sqrt(r21969, r21968, MPFR_RNDN); mpfr_add(r21970, r21961, r21969, MPFR_RNDN); ; mpfr_mul(r21972, r21971, r21964, MPFR_RNDN); mpfr_div(r21973, r21970, r21972, MPFR_RNDN); return mpfr_get_d(r21973, MPFR_RNDN); } static mpfr_t r21974, r21975, r21976, r21977, r21978, r21979, r21980, r21981, r21982, r21983, r21984, r21985, r21986, r21987, r21988, r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r21974); mpfr_init(r21975); mpfr_init(r21976); mpfr_init_set_str(r21977, "4", 10, MPFR_RNDN); mpfr_init(r21978); mpfr_init(r21979); mpfr_init(r21980); mpfr_init(r21981); mpfr_init(r21982); mpfr_init(r21983); mpfr_init(r21984); mpfr_init(r21985); mpfr_init_set_str(r21986, "3", 10, MPFR_RNDN); mpfr_init(r21987); mpfr_init(r21988); mpfr_init(r21989); mpfr_init(r21990); mpfr_init_set_str(r21991, "2", 10, MPFR_RNDN); mpfr_init(r21992); mpfr_init(r21993); mpfr_init(r21994); mpfr_init(r21995); mpfr_init(r21996); mpfr_init(r21997); mpfr_init(r21998); mpfr_init(r21999); mpfr_init(r22000); } double f_fm(double a, double b, double c) { mpfr_set_d(r21974, c, MPFR_RNDN); mpfr_set_d(r21975, a, MPFR_RNDN); mpfr_mul(r21976, r21974, r21975, MPFR_RNDN); ; mpfr_neg(r21978, r21977, MPFR_RNDN); mpfr_mul(r21979, r21976, r21978, MPFR_RNDN); mpfr_mul(r21980, r21977, r21975, MPFR_RNDN); mpfr_neg(r21981, r21974, MPFR_RNDN); mpfr_set_d(r21982, b, MPFR_RNDN); mpfr_mul(r21983, r21982, r21982, MPFR_RNDN); mpfr_fma(r21984, r21980, r21981, r21983, MPFR_RNDN); mpfr_sqrt(r21985, r21984, MPFR_RNDN); ; mpfr_pow(r21987, r21985, r21986, MPFR_RNDN); mpfr_pow(r21988, r21982, r21986, MPFR_RNDN); mpfr_add(r21989, r21987, r21988, MPFR_RNDN); mpfr_div(r21990, r21979, r21989, MPFR_RNDN); ; mpfr_mul(r21992, r21975, r21991, MPFR_RNDN); mpfr_neg(r21993, r21982, MPFR_RNDN); mpfr_add(r21994, r21982, r21982, MPFR_RNDN); mpfr_mul(r21995, r21994, r21982, MPFR_RNDN); mpfr_mul(r21996, r21977, r21976, MPFR_RNDN); mpfr_sub(r21997, r21995, r21996, MPFR_RNDN); mpfr_fma(r21998, r21993, r21985, r21997, MPFR_RNDN); mpfr_div(r21999, r21992, r21998, MPFR_RNDN); mpfr_div(r22000, r21990, r21999, MPFR_RNDN); return mpfr_get_d(r22000, MPFR_RNDN); } static mpfr_t r22001, r22002, r22003, r22004, r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22001); mpfr_init(r22002); mpfr_init(r22003); mpfr_init_set_str(r22004, "4", 10, MPFR_RNDN); mpfr_init(r22005); mpfr_init(r22006); mpfr_init(r22007); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init(r22012); mpfr_init_set_str(r22013, "3", 10, MPFR_RNDN); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init_set_str(r22018, "2", 10, MPFR_RNDN); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init(r22025); mpfr_init(r22026); mpfr_init(r22027); } double f_dm(double a, double b, double c) { mpfr_set_d(r22001, c, MPFR_RNDN); mpfr_set_d(r22002, a, MPFR_RNDN); mpfr_mul(r22003, r22001, r22002, MPFR_RNDN); ; mpfr_neg(r22005, r22004, MPFR_RNDN); mpfr_mul(r22006, r22003, r22005, MPFR_RNDN); mpfr_mul(r22007, r22004, r22002, MPFR_RNDN); mpfr_neg(r22008, r22001, MPFR_RNDN); mpfr_set_d(r22009, b, MPFR_RNDN); mpfr_mul(r22010, r22009, r22009, MPFR_RNDN); mpfr_fma(r22011, r22007, r22008, r22010, MPFR_RNDN); mpfr_sqrt(r22012, r22011, MPFR_RNDN); ; mpfr_pow(r22014, r22012, r22013, MPFR_RNDN); mpfr_pow(r22015, r22009, r22013, MPFR_RNDN); mpfr_add(r22016, r22014, r22015, MPFR_RNDN); mpfr_div(r22017, r22006, r22016, MPFR_RNDN); ; mpfr_mul(r22019, r22002, r22018, MPFR_RNDN); mpfr_neg(r22020, r22009, MPFR_RNDN); mpfr_add(r22021, r22009, r22009, MPFR_RNDN); mpfr_mul(r22022, r22021, r22009, MPFR_RNDN); mpfr_mul(r22023, r22004, r22003, MPFR_RNDN); mpfr_sub(r22024, r22022, r22023, MPFR_RNDN); mpfr_fma(r22025, r22020, r22012, r22024, MPFR_RNDN); mpfr_div(r22026, r22019, r22025, MPFR_RNDN); mpfr_div(r22027, r22017, r22026, MPFR_RNDN); return mpfr_get_d(r22027, MPFR_RNDN); }