#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r21915 = b; float r21916 = -r21915; float r21917 = r21915 * r21915; float r21918 = 4; float r21919 = a; float r21920 = r21918 * r21919; float r21921 = c; float r21922 = r21920 * r21921; float r21923 = r21917 - r21922; float r21924 = sqrt(r21923); float r21925 = r21916 + r21924; float r21926 = 2; float r21927 = r21926 * r21919; float r21928 = r21925 / r21927; return r21928; } double f_id(double a, double b, double c) { double r21929 = b; double r21930 = -r21929; double r21931 = r21929 * r21929; double r21932 = 4; double r21933 = a; double r21934 = r21932 * r21933; double r21935 = c; double r21936 = r21934 * r21935; double r21937 = r21931 - r21936; double r21938 = sqrt(r21937); double r21939 = r21930 + r21938; double r21940 = 2; double r21941 = r21940 * r21933; double r21942 = r21939 / r21941; return r21942; } double f_of(float a, float b, float c) { float r21943 = 1; float r21944 = 2; float r21945 = r21943 / r21944; float r21946 = c; float r21947 = r21946 / r21943; float r21948 = 4; float r21949 = r21947 * r21948; float r21950 = b; float r21951 = -r21950; float r21952 = r21950 * r21950; float r21953 = 3; float r21954 = pow(r21952, r21953); float r21955 = a; float r21956 = r21948 * r21955; float r21957 = r21946 * r21956; float r21958 = pow(r21957, r21953); float r21959 = r21954 - r21958; float r21960 = r21946 * r21948; float r21961 = r21960 * r21955; float r21962 = r21952 + r21961; float r21963 = r21962 * r21961; float r21964 = r21952 * r21952; float r21965 = r21963 + r21964; float r21966 = r21959 / r21965; float r21967 = sqrt(r21966); float r21968 = r21951 - r21967; float r21969 = r21949 / r21968; float r21970 = r21945 * r21969; return r21970; } double f_od(double a, double b, double c) { double r21971 = 1; double r21972 = 2; double r21973 = r21971 / r21972; double r21974 = c; double r21975 = r21974 / r21971; double r21976 = 4; double r21977 = r21975 * r21976; double r21978 = b; double r21979 = -r21978; double r21980 = r21978 * r21978; double r21981 = 3; double r21982 = pow(r21980, r21981); double r21983 = a; double r21984 = r21976 * r21983; double r21985 = r21974 * r21984; double r21986 = pow(r21985, r21981); double r21987 = r21982 - r21986; double r21988 = r21974 * r21976; double r21989 = r21988 * r21983; double r21990 = r21980 + r21989; double r21991 = r21990 * r21989; double r21992 = r21980 * r21980; double r21993 = r21991 + r21992; double r21994 = r21987 / r21993; double r21995 = sqrt(r21994); double r21996 = r21979 - r21995; double r21997 = r21977 / r21996; double r21998 = r21973 * r21997; return r21998; } 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 r21999, r22000, r22001, r22002, r22003, r22004, r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21999); mpfr_init(r22000); mpfr_init(r22001); mpfr_init_set_str(r22002, "4", 10, MPFR_RNDN); mpfr_init(r22003); mpfr_init(r22004); mpfr_init(r22005); mpfr_init(r22006); mpfr_init(r22007); mpfr_init(r22008); mpfr_init(r22009); mpfr_init_set_str(r22010, "2", 10, MPFR_RNDN); mpfr_init(r22011); mpfr_init(r22012); } double f_im(double a, double b, double c) { mpfr_set_d(r21999, b, MPFR_RNDN); mpfr_neg(r22000, r21999, MPFR_RNDN); mpfr_mul(r22001, r21999, r21999, MPFR_RNDN); ; mpfr_set_d(r22003, a, MPFR_RNDN); mpfr_mul(r22004, r22002, r22003, MPFR_RNDN); mpfr_set_d(r22005, c, MPFR_RNDN); mpfr_mul(r22006, r22004, r22005, MPFR_RNDN); mpfr_sub(r22007, r22001, r22006, MPFR_RNDN); mpfr_sqrt(r22008, r22007, MPFR_RNDN); mpfr_add(r22009, r22000, r22008, MPFR_RNDN); ; mpfr_mul(r22011, r22010, r22003, MPFR_RNDN); mpfr_div(r22012, r22009, r22011, MPFR_RNDN); return mpfr_get_d(r22012, MPFR_RNDN); } static mpfr_t r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039, r22040; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r22013, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22014, "2", 10, MPFR_RNDN); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init_set_str(r22018, "4", 10, MPFR_RNDN); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); mpfr_init(r22022); mpfr_init_set_str(r22023, "3", 10, MPFR_RNDN); mpfr_init(r22024); mpfr_init(r22025); mpfr_init(r22026); mpfr_init(r22027); mpfr_init(r22028); mpfr_init(r22029); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); mpfr_init(r22040); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r22015, r22013, r22014, MPFR_RNDN); mpfr_set_d(r22016, c, MPFR_RNDN); mpfr_div(r22017, r22016, r22013, MPFR_RNDN); ; mpfr_mul(r22019, r22017, r22018, MPFR_RNDN); mpfr_set_d(r22020, b, MPFR_RNDN); mpfr_neg(r22021, r22020, MPFR_RNDN); mpfr_mul(r22022, r22020, r22020, MPFR_RNDN); ; mpfr_pow(r22024, r22022, r22023, MPFR_RNDN); mpfr_set_d(r22025, a, MPFR_RNDN); mpfr_mul(r22026, r22018, r22025, MPFR_RNDN); mpfr_mul(r22027, r22016, r22026, MPFR_RNDN); mpfr_pow(r22028, r22027, r22023, MPFR_RNDN); mpfr_sub(r22029, r22024, r22028, MPFR_RNDN); mpfr_mul(r22030, r22016, r22018, MPFR_RNDN); mpfr_mul(r22031, r22030, r22025, MPFR_RNDN); mpfr_add(r22032, r22022, r22031, MPFR_RNDN); mpfr_mul(r22033, r22032, r22031, MPFR_RNDN); mpfr_mul(r22034, r22022, r22022, MPFR_RNDN); mpfr_add(r22035, r22033, r22034, MPFR_RNDN); mpfr_div(r22036, r22029, r22035, MPFR_RNDN); mpfr_sqrt(r22037, r22036, MPFR_RNDN); mpfr_sub(r22038, r22021, r22037, MPFR_RNDN); mpfr_div(r22039, r22019, r22038, MPFR_RNDN); mpfr_mul(r22040, r22015, r22039, MPFR_RNDN); return mpfr_get_d(r22040, MPFR_RNDN); } static mpfr_t r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052, r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r22041, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22042, "2", 10, MPFR_RNDN); mpfr_init(r22043); mpfr_init(r22044); mpfr_init(r22045); mpfr_init_set_str(r22046, "4", 10, MPFR_RNDN); mpfr_init(r22047); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init_set_str(r22051, "3", 10, MPFR_RNDN); mpfr_init(r22052); mpfr_init(r22053); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init(r22057); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); mpfr_init(r22068); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r22043, r22041, r22042, MPFR_RNDN); mpfr_set_d(r22044, c, MPFR_RNDN); mpfr_div(r22045, r22044, r22041, MPFR_RNDN); ; mpfr_mul(r22047, r22045, r22046, MPFR_RNDN); mpfr_set_d(r22048, b, MPFR_RNDN); mpfr_neg(r22049, r22048, MPFR_RNDN); mpfr_mul(r22050, r22048, r22048, MPFR_RNDN); ; mpfr_pow(r22052, r22050, r22051, MPFR_RNDN); mpfr_set_d(r22053, a, MPFR_RNDN); mpfr_mul(r22054, r22046, r22053, MPFR_RNDN); mpfr_mul(r22055, r22044, r22054, MPFR_RNDN); mpfr_pow(r22056, r22055, r22051, MPFR_RNDN); mpfr_sub(r22057, r22052, r22056, MPFR_RNDN); mpfr_mul(r22058, r22044, r22046, MPFR_RNDN); mpfr_mul(r22059, r22058, r22053, MPFR_RNDN); mpfr_add(r22060, r22050, r22059, MPFR_RNDN); mpfr_mul(r22061, r22060, r22059, MPFR_RNDN); mpfr_mul(r22062, r22050, r22050, MPFR_RNDN); mpfr_add(r22063, r22061, r22062, MPFR_RNDN); mpfr_div(r22064, r22057, r22063, MPFR_RNDN); mpfr_sqrt(r22065, r22064, MPFR_RNDN); mpfr_sub(r22066, r22049, r22065, MPFR_RNDN); mpfr_div(r22067, r22047, r22066, MPFR_RNDN); mpfr_mul(r22068, r22043, r22067, MPFR_RNDN); return mpfr_get_d(r22068, MPFR_RNDN); }