#include #include #include #include #include char *name = "Cubic critical, wide range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r26056 = b; float r26057 = -r26056; float r26058 = r26056 * r26056; float r26059 = 3; float r26060 = a; float r26061 = r26059 * r26060; float r26062 = c; float r26063 = r26061 * r26062; float r26064 = r26058 - r26063; float r26065 = sqrt(r26064); float r26066 = r26057 + r26065; float r26067 = r26066 / r26061; return r26067; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26068 = b; double r26069 = -r26068; double r26070 = r26068 * r26068; double r26071 = 3; double r26072 = a; double r26073 = r26071 * r26072; double r26074 = c; double r26075 = r26073 * r26074; double r26076 = r26070 - r26075; double r26077 = sqrt(r26076); double r26078 = r26069 + r26077; double r26079 = r26078 / r26073; return r26079; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26080 = a; float r26081 = 3; float r26082 = r26080 / r26081; float r26083 = c; float r26084 = r26082 * r26083; float r26085 = -r26081; float r26086 = r26081 * r26080; float r26087 = -r26083; float r26088 = b; float r26089 = r26088 * r26088; float r26090 = fma(r26086, r26087, r26089); float r26091 = sqrt(r26090); float r26092 = r26091 + r26088; float r26093 = r26085 / r26092; float r26094 = r26093 / r26080; float r26095 = r26084 * r26094; return r26095; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26096 = a; double r26097 = 3; double r26098 = r26096 / r26097; double r26099 = c; double r26100 = r26098 * r26099; double r26101 = -r26097; double r26102 = r26097 * r26096; double r26103 = -r26099; double r26104 = b; double r26105 = r26104 * r26104; double r26106 = fma(r26102, r26103, r26105); double r26107 = sqrt(r26106); double r26108 = r26107 + r26104; double r26109 = r26101 / r26108; double r26110 = r26109 / r26096; double r26111 = r26100 * r26110; return r26111; } 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 r26112, r26113, r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r26112); mpfr_init(r26113); mpfr_init(r26114); mpfr_init_set_str(r26115, "3", 10, MPFR_RNDN); mpfr_init(r26116); mpfr_init(r26117); mpfr_init(r26118); mpfr_init(r26119); mpfr_init(r26120); mpfr_init(r26121); mpfr_init(r26122); mpfr_init(r26123); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26112, b, MPFR_RNDN); mpfr_neg(r26113, r26112, MPFR_RNDN); mpfr_mul(r26114, r26112, r26112, MPFR_RNDN); ; mpfr_set_d(r26116, a, MPFR_RNDN); mpfr_mul(r26117, r26115, r26116, MPFR_RNDN); mpfr_set_d(r26118, c, MPFR_RNDN); mpfr_mul(r26119, r26117, r26118, MPFR_RNDN); mpfr_sub(r26120, r26114, r26119, MPFR_RNDN); mpfr_sqrt(r26121, r26120, MPFR_RNDN); mpfr_add(r26122, r26113, r26121, MPFR_RNDN); mpfr_div(r26123, r26122, r26117, MPFR_RNDN); return mpfr_get_d(r26123, MPFR_RNDN); } static mpfr_t r26124, r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r26124); mpfr_init_set_str(r26125, "3", 10, MPFR_RNDN); mpfr_init(r26126); mpfr_init(r26127); mpfr_init(r26128); mpfr_init(r26129); mpfr_init(r26130); mpfr_init(r26131); mpfr_init(r26132); mpfr_init(r26133); mpfr_init(r26134); mpfr_init(r26135); mpfr_init(r26136); mpfr_init(r26137); mpfr_init(r26138); mpfr_init(r26139); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26124, a, MPFR_RNDN); ; mpfr_div(r26126, r26124, r26125, MPFR_RNDN); mpfr_set_d(r26127, c, MPFR_RNDN); mpfr_mul(r26128, r26126, r26127, MPFR_RNDN); mpfr_neg(r26129, r26125, MPFR_RNDN); mpfr_mul(r26130, r26125, r26124, MPFR_RNDN); mpfr_neg(r26131, r26127, MPFR_RNDN); mpfr_set_d(r26132, b, MPFR_RNDN); mpfr_mul(r26133, r26132, r26132, MPFR_RNDN); mpfr_fma(r26134, r26130, r26131, r26133, MPFR_RNDN); mpfr_sqrt(r26135, r26134, MPFR_RNDN); mpfr_add(r26136, r26135, r26132, MPFR_RNDN); mpfr_div(r26137, r26129, r26136, MPFR_RNDN); mpfr_div(r26138, r26137, r26124, MPFR_RNDN); mpfr_mul(r26139, r26128, r26138, MPFR_RNDN); return mpfr_get_d(r26139, MPFR_RNDN); } static mpfr_t r26140, r26141, r26142, r26143, r26144, r26145, r26146, r26147, r26148, r26149, r26150, r26151, r26152, r26153, r26154, r26155; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r26140); mpfr_init_set_str(r26141, "3", 10, MPFR_RNDN); mpfr_init(r26142); mpfr_init(r26143); mpfr_init(r26144); mpfr_init(r26145); mpfr_init(r26146); mpfr_init(r26147); mpfr_init(r26148); mpfr_init(r26149); mpfr_init(r26150); mpfr_init(r26151); mpfr_init(r26152); mpfr_init(r26153); mpfr_init(r26154); mpfr_init(r26155); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26140, a, MPFR_RNDN); ; mpfr_div(r26142, r26140, r26141, MPFR_RNDN); mpfr_set_d(r26143, c, MPFR_RNDN); mpfr_mul(r26144, r26142, r26143, MPFR_RNDN); mpfr_neg(r26145, r26141, MPFR_RNDN); mpfr_mul(r26146, r26141, r26140, MPFR_RNDN); mpfr_neg(r26147, r26143, MPFR_RNDN); mpfr_set_d(r26148, b, MPFR_RNDN); mpfr_mul(r26149, r26148, r26148, MPFR_RNDN); mpfr_fma(r26150, r26146, r26147, r26149, MPFR_RNDN); mpfr_sqrt(r26151, r26150, MPFR_RNDN); mpfr_add(r26152, r26151, r26148, MPFR_RNDN); mpfr_div(r26153, r26145, r26152, MPFR_RNDN); mpfr_div(r26154, r26153, r26140, MPFR_RNDN); mpfr_mul(r26155, r26144, r26154, MPFR_RNDN); return mpfr_get_d(r26155, MPFR_RNDN); }