#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r20934 = b; float r20935 = -r20934; float r20936 = r20934 * r20934; float r20937 = 4; float r20938 = a; float r20939 = r20937 * r20938; float r20940 = c; float r20941 = r20939 * r20940; float r20942 = r20936 - r20941; float r20943 = sqrt(r20942); float r20944 = r20935 + r20943; float r20945 = 2; float r20946 = r20945 * r20938; float r20947 = r20944 / r20946; return r20947; } double f_id(double a, double b, double c) { double r20948 = b; double r20949 = -r20948; double r20950 = r20948 * r20948; double r20951 = 4; double r20952 = a; double r20953 = r20951 * r20952; double r20954 = c; double r20955 = r20953 * r20954; double r20956 = r20950 - r20955; double r20957 = sqrt(r20956); double r20958 = r20949 + r20957; double r20959 = 2; double r20960 = r20959 * r20952; double r20961 = r20958 / r20960; return r20961; } double f_of(float a, float b, float c) { float r20962 = c; float r20963 = 2; float r20964 = r20962 / r20963; float r20965 = 4; float r20966 = b; float r20967 = -r20966; float r20968 = r20966 * r20966; float r20969 = a; float r20970 = r20962 * r20965; float r20971 = r20969 * r20970; float r20972 = r20968 - r20971; float r20973 = sqrt(r20972); float r20974 = r20967 - r20973; float r20975 = r20965 / r20974; float r20976 = r20964 * r20975; return r20976; } double f_od(double a, double b, double c) { double r20977 = c; double r20978 = 2; double r20979 = r20977 / r20978; double r20980 = 4; double r20981 = b; double r20982 = -r20981; double r20983 = r20981 * r20981; double r20984 = a; double r20985 = r20977 * r20980; double r20986 = r20984 * r20985; double r20987 = r20983 - r20986; double r20988 = sqrt(r20987); double r20989 = r20982 - r20988; double r20990 = r20980 / r20989; double r20991 = r20979 * r20990; return r20991; } 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 r20992, r20993, r20994, r20995, r20996, r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r20992); mpfr_init(r20993); mpfr_init(r20994); mpfr_init_set_str(r20995, "4", 10, MPFR_RNDN); mpfr_init(r20996); mpfr_init(r20997); mpfr_init(r20998); mpfr_init(r20999); mpfr_init(r21000); mpfr_init(r21001); mpfr_init(r21002); mpfr_init_set_str(r21003, "2", 10, MPFR_RNDN); mpfr_init(r21004); mpfr_init(r21005); } double f_im(double a, double b, double c) { mpfr_set_d(r20992, b, MPFR_RNDN); mpfr_neg(r20993, r20992, MPFR_RNDN); mpfr_mul(r20994, r20992, r20992, MPFR_RNDN); ; mpfr_set_d(r20996, a, MPFR_RNDN); mpfr_mul(r20997, r20995, r20996, MPFR_RNDN); mpfr_set_d(r20998, c, MPFR_RNDN); mpfr_mul(r20999, r20997, r20998, MPFR_RNDN); mpfr_sub(r21000, r20994, r20999, MPFR_RNDN); mpfr_sqrt(r21001, r21000, MPFR_RNDN); mpfr_add(r21002, r20993, r21001, MPFR_RNDN); ; mpfr_mul(r21004, r21003, r20996, MPFR_RNDN); mpfr_div(r21005, r21002, r21004, MPFR_RNDN); return mpfr_get_d(r21005, MPFR_RNDN); } static mpfr_t r21006, r21007, r21008, r21009, r21010, r21011, r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r21006); mpfr_init_set_str(r21007, "2", 10, MPFR_RNDN); mpfr_init(r21008); mpfr_init_set_str(r21009, "4", 10, MPFR_RNDN); mpfr_init(r21010); mpfr_init(r21011); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init(r21015); mpfr_init(r21016); mpfr_init(r21017); mpfr_init(r21018); mpfr_init(r21019); mpfr_init(r21020); } double f_fm(double a, double b, double c) { mpfr_set_d(r21006, c, MPFR_RNDN); ; mpfr_div(r21008, r21006, r21007, MPFR_RNDN); ; mpfr_set_d(r21010, b, MPFR_RNDN); mpfr_neg(r21011, r21010, MPFR_RNDN); mpfr_mul(r21012, r21010, r21010, MPFR_RNDN); mpfr_set_d(r21013, a, MPFR_RNDN); mpfr_mul(r21014, r21006, r21009, MPFR_RNDN); mpfr_mul(r21015, r21013, r21014, MPFR_RNDN); mpfr_sub(r21016, r21012, r21015, MPFR_RNDN); mpfr_sqrt(r21017, r21016, MPFR_RNDN); mpfr_sub(r21018, r21011, r21017, MPFR_RNDN); mpfr_div(r21019, r21009, r21018, MPFR_RNDN); mpfr_mul(r21020, r21008, r21019, MPFR_RNDN); return mpfr_get_d(r21020, MPFR_RNDN); } static mpfr_t r21021, r21022, r21023, r21024, r21025, r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r21021); mpfr_init_set_str(r21022, "2", 10, MPFR_RNDN); mpfr_init(r21023); mpfr_init_set_str(r21024, "4", 10, MPFR_RNDN); mpfr_init(r21025); mpfr_init(r21026); mpfr_init(r21027); mpfr_init(r21028); mpfr_init(r21029); mpfr_init(r21030); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); mpfr_init(r21035); } double f_dm(double a, double b, double c) { mpfr_set_d(r21021, c, MPFR_RNDN); ; mpfr_div(r21023, r21021, r21022, MPFR_RNDN); ; mpfr_set_d(r21025, b, MPFR_RNDN); mpfr_neg(r21026, r21025, MPFR_RNDN); mpfr_mul(r21027, r21025, r21025, MPFR_RNDN); mpfr_set_d(r21028, a, MPFR_RNDN); mpfr_mul(r21029, r21021, r21024, MPFR_RNDN); mpfr_mul(r21030, r21028, r21029, MPFR_RNDN); mpfr_sub(r21031, r21027, r21030, MPFR_RNDN); mpfr_sqrt(r21032, r21031, MPFR_RNDN); mpfr_sub(r21033, r21026, r21032, MPFR_RNDN); mpfr_div(r21034, r21024, r21033, MPFR_RNDN); mpfr_mul(r21035, r21023, r21034, MPFR_RNDN); return mpfr_get_d(r21035, MPFR_RNDN); }