#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r21007 = b; float r21008 = -r21007; float r21009 = r21007 * r21007; float r21010 = 4; float r21011 = a; float r21012 = r21010 * r21011; float r21013 = c; float r21014 = r21012 * r21013; float r21015 = r21009 - r21014; float r21016 = sqrt(r21015); float r21017 = r21008 + r21016; float r21018 = 2; float r21019 = r21018 * r21011; float r21020 = r21017 / r21019; return r21020; } double f_id(double a, double b, double c) { double r21021 = b; double r21022 = -r21021; double r21023 = r21021 * r21021; double r21024 = 4; double r21025 = a; double r21026 = r21024 * r21025; double r21027 = c; double r21028 = r21026 * r21027; double r21029 = r21023 - r21028; double r21030 = sqrt(r21029); double r21031 = r21022 + r21030; double r21032 = 2; double r21033 = r21032 * r21025; double r21034 = r21031 / r21033; return r21034; } double f_of(float a, float b, float c) { float r21035 = 1; float r21036 = 2; float r21037 = r21035 / r21036; float r21038 = c; float r21039 = 4; float r21040 = r21038 * r21039; float r21041 = b; float r21042 = -r21041; float r21043 = r21041 * r21041; float r21044 = a; float r21045 = r21039 * r21044; float r21046 = r21038 * r21045; float r21047 = r21043 - r21046; float r21048 = sqrt(r21047); float r21049 = r21042 - r21048; float r21050 = r21040 / r21049; float r21051 = r21037 * r21050; return r21051; } double f_od(double a, double b, double c) { double r21052 = 1; double r21053 = 2; double r21054 = r21052 / r21053; double r21055 = c; double r21056 = 4; double r21057 = r21055 * r21056; double r21058 = b; double r21059 = -r21058; double r21060 = r21058 * r21058; double r21061 = a; double r21062 = r21056 * r21061; double r21063 = r21055 * r21062; double r21064 = r21060 - r21063; double r21065 = sqrt(r21064); double r21066 = r21059 - r21065; double r21067 = r21057 / r21066; double r21068 = r21054 * r21067; return r21068; } 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 r21069, r21070, r21071, r21072, r21073, r21074, r21075, r21076, r21077, r21078, r21079, r21080, r21081, r21082; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21069); mpfr_init(r21070); mpfr_init(r21071); mpfr_init_set_str(r21072, "4", 10, MPFR_RNDN); mpfr_init(r21073); mpfr_init(r21074); mpfr_init(r21075); mpfr_init(r21076); mpfr_init(r21077); mpfr_init(r21078); mpfr_init(r21079); mpfr_init_set_str(r21080, "2", 10, MPFR_RNDN); mpfr_init(r21081); mpfr_init(r21082); } double f_im(double a, double b, double c) { mpfr_set_d(r21069, b, MPFR_RNDN); mpfr_neg(r21070, r21069, MPFR_RNDN); mpfr_mul(r21071, r21069, r21069, MPFR_RNDN); ; mpfr_set_d(r21073, a, MPFR_RNDN); mpfr_mul(r21074, r21072, r21073, MPFR_RNDN); mpfr_set_d(r21075, c, MPFR_RNDN); mpfr_mul(r21076, r21074, r21075, MPFR_RNDN); mpfr_sub(r21077, r21071, r21076, MPFR_RNDN); mpfr_sqrt(r21078, r21077, MPFR_RNDN); mpfr_add(r21079, r21070, r21078, MPFR_RNDN); ; mpfr_mul(r21081, r21080, r21073, MPFR_RNDN); mpfr_div(r21082, r21079, r21081, MPFR_RNDN); return mpfr_get_d(r21082, MPFR_RNDN); } static mpfr_t r21083, r21084, r21085, r21086, r21087, r21088, r21089, r21090, r21091, r21092, r21093, r21094, r21095, r21096, r21097, r21098, r21099; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21083, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21084, "2", 10, MPFR_RNDN); mpfr_init(r21085); mpfr_init(r21086); mpfr_init_set_str(r21087, "4", 10, MPFR_RNDN); mpfr_init(r21088); mpfr_init(r21089); mpfr_init(r21090); mpfr_init(r21091); mpfr_init(r21092); mpfr_init(r21093); mpfr_init(r21094); mpfr_init(r21095); mpfr_init(r21096); mpfr_init(r21097); mpfr_init(r21098); mpfr_init(r21099); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r21085, r21083, r21084, MPFR_RNDN); mpfr_set_d(r21086, c, MPFR_RNDN); ; mpfr_mul(r21088, r21086, r21087, MPFR_RNDN); mpfr_set_d(r21089, b, MPFR_RNDN); mpfr_neg(r21090, r21089, MPFR_RNDN); mpfr_mul(r21091, r21089, r21089, MPFR_RNDN); mpfr_set_d(r21092, a, MPFR_RNDN); mpfr_mul(r21093, r21087, r21092, MPFR_RNDN); mpfr_mul(r21094, r21086, r21093, MPFR_RNDN); mpfr_sub(r21095, r21091, r21094, MPFR_RNDN); mpfr_sqrt(r21096, r21095, MPFR_RNDN); mpfr_sub(r21097, r21090, r21096, MPFR_RNDN); mpfr_div(r21098, r21088, r21097, MPFR_RNDN); mpfr_mul(r21099, r21085, r21098, MPFR_RNDN); return mpfr_get_d(r21099, MPFR_RNDN); } static mpfr_t r21100, r21101, r21102, r21103, r21104, r21105, r21106, r21107, r21108, r21109, r21110, r21111, r21112, r21113, r21114, r21115, r21116; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21100, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21101, "2", 10, MPFR_RNDN); mpfr_init(r21102); mpfr_init(r21103); mpfr_init_set_str(r21104, "4", 10, MPFR_RNDN); mpfr_init(r21105); mpfr_init(r21106); mpfr_init(r21107); mpfr_init(r21108); mpfr_init(r21109); mpfr_init(r21110); mpfr_init(r21111); mpfr_init(r21112); mpfr_init(r21113); mpfr_init(r21114); mpfr_init(r21115); mpfr_init(r21116); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r21102, r21100, r21101, MPFR_RNDN); mpfr_set_d(r21103, c, MPFR_RNDN); ; mpfr_mul(r21105, r21103, r21104, MPFR_RNDN); mpfr_set_d(r21106, b, MPFR_RNDN); mpfr_neg(r21107, r21106, MPFR_RNDN); mpfr_mul(r21108, r21106, r21106, MPFR_RNDN); mpfr_set_d(r21109, a, MPFR_RNDN); mpfr_mul(r21110, r21104, r21109, MPFR_RNDN); mpfr_mul(r21111, r21103, r21110, MPFR_RNDN); mpfr_sub(r21112, r21108, r21111, MPFR_RNDN); mpfr_sqrt(r21113, r21112, MPFR_RNDN); mpfr_sub(r21114, r21107, r21113, MPFR_RNDN); mpfr_div(r21115, r21105, r21114, MPFR_RNDN); mpfr_mul(r21116, r21102, r21115, MPFR_RNDN); return mpfr_get_d(r21116, MPFR_RNDN); }