#include #include #include #include #include char *name = "Quadratic roots, medium range"; double f_if(float a, float b, float c) { float r21050 = b; float r21051 = -r21050; float r21052 = r21050 * r21050; float r21053 = 4; float r21054 = a; float r21055 = r21053 * r21054; float r21056 = c; float r21057 = r21055 * r21056; float r21058 = r21052 - r21057; float r21059 = sqrt(r21058); float r21060 = r21051 + r21059; float r21061 = 2; float r21062 = r21061 * r21054; float r21063 = r21060 / r21062; return r21063; } double f_id(double a, double b, double c) { double r21064 = b; double r21065 = -r21064; double r21066 = r21064 * r21064; double r21067 = 4; double r21068 = a; double r21069 = r21067 * r21068; double r21070 = c; double r21071 = r21069 * r21070; double r21072 = r21066 - r21071; double r21073 = sqrt(r21072); double r21074 = r21065 + r21073; double r21075 = 2; double r21076 = r21075 * r21068; double r21077 = r21074 / r21076; return r21077; } double f_of(float a, float b, float c) { float r21078 = 4; float r21079 = 2; float r21080 = r21078 / r21079; float r21081 = c; float r21082 = b; float r21083 = -r21082; float r21084 = r21082 * r21082; float r21085 = a; float r21086 = r21081 * r21085; float r21087 = r21086 * r21078; float r21088 = r21084 - r21087; float r21089 = sqrt(r21088); float r21090 = r21083 - r21089; float r21091 = r21081 / r21090; float r21092 = r21080 * r21091; return r21092; } double f_od(double a, double b, double c) { double r21093 = 4; double r21094 = 2; double r21095 = r21093 / r21094; double r21096 = c; double r21097 = b; double r21098 = -r21097; double r21099 = r21097 * r21097; double r21100 = a; double r21101 = r21096 * r21100; double r21102 = r21101 * r21093; double r21103 = r21099 - r21102; double r21104 = sqrt(r21103); double r21105 = r21098 - r21104; double r21106 = r21096 / r21105; double r21107 = r21095 * r21106; return r21107; } 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 r21108, r21109, r21110, r21111, r21112, r21113, r21114, r21115, r21116, r21117, r21118, r21119, r21120, r21121; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r21108); mpfr_init(r21109); mpfr_init(r21110); mpfr_init_set_str(r21111, "4", 10, MPFR_RNDN); mpfr_init(r21112); mpfr_init(r21113); mpfr_init(r21114); mpfr_init(r21115); mpfr_init(r21116); mpfr_init(r21117); mpfr_init(r21118); mpfr_init_set_str(r21119, "2", 10, MPFR_RNDN); mpfr_init(r21120); mpfr_init(r21121); } double f_im(double a, double b, double c) { mpfr_set_d(r21108, b, MPFR_RNDN); mpfr_neg(r21109, r21108, MPFR_RNDN); mpfr_mul(r21110, r21108, r21108, MPFR_RNDN); ; mpfr_set_d(r21112, a, MPFR_RNDN); mpfr_mul(r21113, r21111, r21112, MPFR_RNDN); mpfr_set_d(r21114, c, MPFR_RNDN); mpfr_mul(r21115, r21113, r21114, MPFR_RNDN); mpfr_sub(r21116, r21110, r21115, MPFR_RNDN); mpfr_sqrt(r21117, r21116, MPFR_RNDN); mpfr_add(r21118, r21109, r21117, MPFR_RNDN); ; mpfr_mul(r21120, r21119, r21112, MPFR_RNDN); mpfr_div(r21121, r21118, r21120, MPFR_RNDN); return mpfr_get_d(r21121, MPFR_RNDN); } static mpfr_t r21122, r21123, r21124, r21125, r21126, r21127, r21128, r21129, r21130, r21131, r21132, r21133, r21134, r21135, r21136; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21122, "4", 10, MPFR_RNDN); mpfr_init_set_str(r21123, "2", 10, MPFR_RNDN); mpfr_init(r21124); mpfr_init(r21125); mpfr_init(r21126); mpfr_init(r21127); mpfr_init(r21128); mpfr_init(r21129); mpfr_init(r21130); mpfr_init(r21131); mpfr_init(r21132); mpfr_init(r21133); mpfr_init(r21134); mpfr_init(r21135); mpfr_init(r21136); } double f_fm(double a, double b, double c) { ; ; mpfr_div(r21124, r21122, r21123, MPFR_RNDN); mpfr_set_d(r21125, c, MPFR_RNDN); mpfr_set_d(r21126, b, MPFR_RNDN); mpfr_neg(r21127, r21126, MPFR_RNDN); mpfr_mul(r21128, r21126, r21126, MPFR_RNDN); mpfr_set_d(r21129, a, MPFR_RNDN); mpfr_mul(r21130, r21125, r21129, MPFR_RNDN); mpfr_mul(r21131, r21130, r21122, MPFR_RNDN); mpfr_sub(r21132, r21128, r21131, MPFR_RNDN); mpfr_sqrt(r21133, r21132, MPFR_RNDN); mpfr_sub(r21134, r21127, r21133, MPFR_RNDN); mpfr_div(r21135, r21125, r21134, MPFR_RNDN); mpfr_mul(r21136, r21124, r21135, MPFR_RNDN); return mpfr_get_d(r21136, MPFR_RNDN); } static mpfr_t r21137, r21138, r21139, r21140, r21141, r21142, r21143, r21144, r21145, r21146, r21147, r21148, r21149, r21150, r21151; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r21137, "4", 10, MPFR_RNDN); mpfr_init_set_str(r21138, "2", 10, MPFR_RNDN); mpfr_init(r21139); mpfr_init(r21140); mpfr_init(r21141); mpfr_init(r21142); mpfr_init(r21143); mpfr_init(r21144); mpfr_init(r21145); mpfr_init(r21146); mpfr_init(r21147); mpfr_init(r21148); mpfr_init(r21149); mpfr_init(r21150); mpfr_init(r21151); } double f_dm(double a, double b, double c) { ; ; mpfr_div(r21139, r21137, r21138, MPFR_RNDN); mpfr_set_d(r21140, c, MPFR_RNDN); mpfr_set_d(r21141, b, MPFR_RNDN); mpfr_neg(r21142, r21141, MPFR_RNDN); mpfr_mul(r21143, r21141, r21141, MPFR_RNDN); mpfr_set_d(r21144, a, MPFR_RNDN); mpfr_mul(r21145, r21140, r21144, MPFR_RNDN); mpfr_mul(r21146, r21145, r21137, MPFR_RNDN); mpfr_sub(r21147, r21143, r21146, MPFR_RNDN); mpfr_sqrt(r21148, r21147, MPFR_RNDN); mpfr_sub(r21149, r21142, r21148, MPFR_RNDN); mpfr_div(r21150, r21140, r21149, MPFR_RNDN); mpfr_mul(r21151, r21139, r21150, MPFR_RNDN); return mpfr_get_d(r21151, MPFR_RNDN); }