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