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