#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r20909 = b; float r20910 = -r20909; float r20911 = r20909 * r20909; float r20912 = 4; float r20913 = a; float r20914 = r20912 * r20913; float r20915 = c; float r20916 = r20914 * r20915; float r20917 = r20911 - r20916; float r20918 = sqrt(r20917); float r20919 = r20910 + r20918; float r20920 = 2; float r20921 = r20920 * r20913; float r20922 = r20919 / r20921; return r20922; } double f_id(double a, double b, double c) { double r20923 = b; double r20924 = -r20923; double r20925 = r20923 * r20923; double r20926 = 4; double r20927 = a; double r20928 = r20926 * r20927; double r20929 = c; double r20930 = r20928 * r20929; double r20931 = r20925 - r20930; double r20932 = sqrt(r20931); double r20933 = r20924 + r20932; double r20934 = 2; double r20935 = r20934 * r20927; double r20936 = r20933 / r20935; return r20936; } double f_of(float a, float b, float c) { float r20937 = 4; float r20938 = 2; float r20939 = r20937 / r20938; float r20940 = c; float r20941 = b; float r20942 = -r20941; float r20943 = r20941 * r20941; float r20944 = a; float r20945 = r20940 * r20944; float r20946 = r20945 * r20937; float r20947 = r20943 - r20946; float r20948 = sqrt(r20947); float r20949 = r20942 - r20948; float r20950 = r20940 / r20949; float r20951 = r20939 * r20950; return r20951; } double f_od(double a, double b, double c) { double r20952 = 4; double r20953 = 2; double r20954 = r20952 / r20953; double r20955 = c; double r20956 = b; double r20957 = -r20956; double r20958 = r20956 * r20956; double r20959 = a; double r20960 = r20955 * r20959; double r20961 = r20960 * r20952; double r20962 = r20958 - r20961; double r20963 = sqrt(r20962); double r20964 = r20957 - r20963; double r20965 = r20955 / r20964; double r20966 = r20954 * r20965; return r20966; } 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 r20967, r20968, r20969, r20970, r20971, r20972, r20973, r20974, r20975, r20976, r20977, r20978, r20979, r20980; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r20967); mpfr_init(r20968); mpfr_init(r20969); mpfr_init_set_str(r20970, "4", 10, MPFR_RNDN); mpfr_init(r20971); mpfr_init(r20972); mpfr_init(r20973); mpfr_init(r20974); mpfr_init(r20975); mpfr_init(r20976); mpfr_init(r20977); mpfr_init_set_str(r20978, "2", 10, MPFR_RNDN); mpfr_init(r20979); mpfr_init(r20980); } double f_im(double a, double b, double c) { mpfr_set_d(r20967, b, MPFR_RNDN); mpfr_neg(r20968, r20967, MPFR_RNDN); mpfr_mul(r20969, r20967, r20967, MPFR_RNDN); ; mpfr_set_d(r20971, a, MPFR_RNDN); mpfr_mul(r20972, r20970, r20971, MPFR_RNDN); mpfr_set_d(r20973, c, MPFR_RNDN); mpfr_mul(r20974, r20972, r20973, MPFR_RNDN); mpfr_sub(r20975, r20969, r20974, MPFR_RNDN); mpfr_sqrt(r20976, r20975, MPFR_RNDN); mpfr_add(r20977, r20968, r20976, MPFR_RNDN); ; mpfr_mul(r20979, r20978, r20971, MPFR_RNDN); mpfr_div(r20980, r20977, r20979, MPFR_RNDN); return mpfr_get_d(r20980, MPFR_RNDN); } static mpfr_t r20981, r20982, r20983, r20984, r20985, r20986, r20987, r20988, r20989, r20990, r20991, r20992, r20993, r20994, r20995; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r20981, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20982, "2", 10, MPFR_RNDN); mpfr_init(r20983); mpfr_init(r20984); mpfr_init(r20985); mpfr_init(r20986); mpfr_init(r20987); mpfr_init(r20988); mpfr_init(r20989); mpfr_init(r20990); mpfr_init(r20991); mpfr_init(r20992); mpfr_init(r20993); mpfr_init(r20994); mpfr_init(r20995); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r20983, r20981, r20982, MPFR_RNDN); mpfr_set_d(r20984, c, MPFR_RNDN); mpfr_set_d(r20985, b, MPFR_RNDN); mpfr_neg(r20986, r20985, MPFR_RNDN); mpfr_mul(r20987, r20985, r20985, MPFR_RNDN); mpfr_set_d(r20988, a, MPFR_RNDN); mpfr_mul(r20989, r20984, r20988, MPFR_RNDN); mpfr_mul(r20990, r20989, r20981, MPFR_RNDN); mpfr_sub(r20991, r20987, r20990, MPFR_RNDN); mpfr_sqrt(r20992, r20991, MPFR_RNDN); mpfr_sub(r20993, r20986, r20992, MPFR_RNDN); mpfr_div(r20994, r20984, r20993, MPFR_RNDN); mpfr_mul(r20995, r20983, r20994, MPFR_RNDN); return mpfr_get_d(r20995, MPFR_RNDN); } static mpfr_t r20996, r20997, r20998, r20999, r21000, r21001, r21002, r21003, r21004, r21005, r21006, r21007, r21008, r21009, r21010; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r20996, "4", 10, MPFR_RNDN); mpfr_init_set_str(r20997, "2", 10, MPFR_RNDN); mpfr_init(r20998); mpfr_init(r20999); mpfr_init(r21000); mpfr_init(r21001); mpfr_init(r21002); mpfr_init(r21003); mpfr_init(r21004); mpfr_init(r21005); mpfr_init(r21006); mpfr_init(r21007); mpfr_init(r21008); mpfr_init(r21009); mpfr_init(r21010); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r20998, r20996, r20997, MPFR_RNDN); mpfr_set_d(r20999, c, MPFR_RNDN); mpfr_set_d(r21000, b, MPFR_RNDN); mpfr_neg(r21001, r21000, MPFR_RNDN); mpfr_mul(r21002, r21000, r21000, MPFR_RNDN); mpfr_set_d(r21003, a, MPFR_RNDN); mpfr_mul(r21004, r20999, r21003, MPFR_RNDN); mpfr_mul(r21005, r21004, r20996, MPFR_RNDN); mpfr_sub(r21006, r21002, r21005, MPFR_RNDN); mpfr_sqrt(r21007, r21006, MPFR_RNDN); mpfr_sub(r21008, r21001, r21007, MPFR_RNDN); mpfr_div(r21009, r20999, r21008, MPFR_RNDN); mpfr_mul(r21010, r20998, r21009, MPFR_RNDN); return mpfr_get_d(r21010, MPFR_RNDN); }