#include #include #include #include #include char *name = "Quadratic roots, narrow range"; double f_if(float a, float b, float c) { float r20950 = b; float r20951 = -r20950; float r20952 = r20950 * r20950; float r20953 = 4; float r20954 = a; float r20955 = r20953 * r20954; float r20956 = c; float r20957 = r20955 * r20956; float r20958 = r20952 - r20957; float r20959 = sqrt(r20958); float r20960 = r20951 + r20959; float r20961 = 2; float r20962 = r20961 * r20954; float r20963 = r20960 / r20962; return r20963; } double f_id(double a, double b, double c) { double r20964 = b; double r20965 = -r20964; double r20966 = r20964 * r20964; double r20967 = 4; double r20968 = a; double r20969 = r20967 * r20968; double r20970 = c; double r20971 = r20969 * r20970; double r20972 = r20966 - r20971; double r20973 = sqrt(r20972); double r20974 = r20965 + r20973; double r20975 = 2; double r20976 = r20975 * r20968; double r20977 = r20974 / r20976; return r20977; } double f_of(float a, float b, float c) { float r20978 = 1; float r20979 = 2; float r20980 = r20978 / r20979; float r20981 = c; float r20982 = 4; float r20983 = r20981 * r20982; float r20984 = b; float r20985 = -r20984; float r20986 = r20984 * r20984; float r20987 = a; float r20988 = r20982 * r20987; float r20989 = r20981 * r20988; float r20990 = r20986 - r20989; float r20991 = sqrt(r20990); float r20992 = r20985 - r20991; float r20993 = r20983 / r20992; float r20994 = r20980 * r20993; return r20994; } double f_od(double a, double b, double c) { double r20995 = 1; double r20996 = 2; double r20997 = r20995 / r20996; double r20998 = c; double r20999 = 4; double r21000 = r20998 * r20999; double r21001 = b; double r21002 = -r21001; double r21003 = r21001 * r21001; double r21004 = a; double r21005 = r20999 * r21004; double r21006 = r20998 * r21005; double r21007 = r21003 - r21006; double r21008 = sqrt(r21007); double r21009 = r21002 - r21008; double r21010 = r21000 / r21009; double r21011 = r20997 * r21010; return r21011; } 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 r21012, r21013, r21014, r21015, r21016, r21017, r21018, r21019, r21020, r21021, r21022, r21023, r21024, r21025; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21012); mpfr_init(r21013); mpfr_init(r21014); mpfr_init_set_str(r21015, "4", 10, MPFR_RNDN); mpfr_init(r21016); mpfr_init(r21017); mpfr_init(r21018); mpfr_init(r21019); mpfr_init(r21020); mpfr_init(r21021); mpfr_init(r21022); mpfr_init_set_str(r21023, "2", 10, MPFR_RNDN); mpfr_init(r21024); mpfr_init(r21025); } double f_im(double a, double b, double c) { mpfr_set_d(r21012, b, MPFR_RNDN); mpfr_neg(r21013, r21012, MPFR_RNDN); mpfr_mul(r21014, r21012, r21012, MPFR_RNDN); ; mpfr_set_d(r21016, a, MPFR_RNDN); mpfr_mul(r21017, r21015, r21016, MPFR_RNDN); mpfr_set_d(r21018, c, MPFR_RNDN); mpfr_mul(r21019, r21017, r21018, MPFR_RNDN); mpfr_sub(r21020, r21014, r21019, MPFR_RNDN); mpfr_sqrt(r21021, r21020, MPFR_RNDN); mpfr_add(r21022, r21013, r21021, MPFR_RNDN); ; mpfr_mul(r21024, r21023, r21016, MPFR_RNDN); mpfr_div(r21025, r21022, r21024, MPFR_RNDN); return mpfr_get_d(r21025, MPFR_RNDN); } static mpfr_t r21026, r21027, r21028, r21029, r21030, r21031, r21032, r21033, r21034, r21035, r21036, r21037, r21038, r21039, r21040, r21041, r21042; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21026, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21027, "2", 10, MPFR_RNDN); mpfr_init(r21028); mpfr_init(r21029); mpfr_init_set_str(r21030, "4", 10, MPFR_RNDN); mpfr_init(r21031); mpfr_init(r21032); mpfr_init(r21033); mpfr_init(r21034); mpfr_init(r21035); mpfr_init(r21036); mpfr_init(r21037); mpfr_init(r21038); mpfr_init(r21039); mpfr_init(r21040); mpfr_init(r21041); mpfr_init(r21042); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r21028, r21026, r21027, MPFR_RNDN); mpfr_set_d(r21029, c, MPFR_RNDN); ; mpfr_mul(r21031, r21029, r21030, MPFR_RNDN); mpfr_set_d(r21032, b, MPFR_RNDN); mpfr_neg(r21033, r21032, MPFR_RNDN); mpfr_mul(r21034, r21032, r21032, MPFR_RNDN); mpfr_set_d(r21035, a, MPFR_RNDN); mpfr_mul(r21036, r21030, r21035, MPFR_RNDN); mpfr_mul(r21037, r21029, r21036, MPFR_RNDN); mpfr_sub(r21038, r21034, r21037, MPFR_RNDN); mpfr_sqrt(r21039, r21038, MPFR_RNDN); mpfr_sub(r21040, r21033, r21039, MPFR_RNDN); mpfr_div(r21041, r21031, r21040, MPFR_RNDN); mpfr_mul(r21042, r21028, r21041, MPFR_RNDN); return mpfr_get_d(r21042, MPFR_RNDN); } static mpfr_t r21043, r21044, r21045, r21046, r21047, r21048, r21049, r21050, r21051, r21052, r21053, r21054, r21055, r21056, r21057, r21058, r21059; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21043, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21044, "2", 10, MPFR_RNDN); mpfr_init(r21045); mpfr_init(r21046); mpfr_init_set_str(r21047, "4", 10, MPFR_RNDN); mpfr_init(r21048); mpfr_init(r21049); mpfr_init(r21050); mpfr_init(r21051); mpfr_init(r21052); mpfr_init(r21053); mpfr_init(r21054); mpfr_init(r21055); mpfr_init(r21056); mpfr_init(r21057); mpfr_init(r21058); mpfr_init(r21059); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r21045, r21043, r21044, MPFR_RNDN); mpfr_set_d(r21046, c, MPFR_RNDN); ; mpfr_mul(r21048, r21046, r21047, MPFR_RNDN); mpfr_set_d(r21049, b, MPFR_RNDN); mpfr_neg(r21050, r21049, MPFR_RNDN); mpfr_mul(r21051, r21049, r21049, MPFR_RNDN); mpfr_set_d(r21052, a, MPFR_RNDN); mpfr_mul(r21053, r21047, r21052, MPFR_RNDN); mpfr_mul(r21054, r21046, r21053, MPFR_RNDN); mpfr_sub(r21055, r21051, r21054, MPFR_RNDN); mpfr_sqrt(r21056, r21055, MPFR_RNDN); mpfr_sub(r21057, r21050, r21056, MPFR_RNDN); mpfr_div(r21058, r21048, r21057, MPFR_RNDN); mpfr_mul(r21059, r21045, r21058, MPFR_RNDN); return mpfr_get_d(r21059, MPFR_RNDN); }