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