#include #include #include #include #include char *name = "Cubic critical, medium range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r10063 = b; float r10064 = -r10063; float r10065 = r10063 * r10063; float r10066 = 3; float r10067 = a; float r10068 = r10066 * r10067; float r10069 = c; float r10070 = r10068 * r10069; float r10071 = r10065 - r10070; float r10072 = sqrt(r10071); float r10073 = r10064 + r10072; float r10074 = r10073 / r10068; return r10074; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r10075 = b; double r10076 = -r10075; double r10077 = r10075 * r10075; double r10078 = 3; double r10079 = a; double r10080 = r10078 * r10079; double r10081 = c; double r10082 = r10080 * r10081; double r10083 = r10077 - r10082; double r10084 = sqrt(r10083); double r10085 = r10076 + r10084; double r10086 = r10085 / r10080; return r10086; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r10087 = 3; float r10088 = r10087 / r10087; float r10089 = c; float r10090 = b; float r10091 = -r10090; float r10092 = -r10089; float r10093 = a; float r10094 = r10087 * r10093; float r10095 = r10090 * r10090; float r10096 = fma(r10092, r10094, r10095); float r10097 = pow(r10096, r10087); float r10098 = cbrt(r10097); float r10099 = sqrt(r10098); float r10100 = r10091 - r10099; float r10101 = r10089 / r10100; float r10102 = r10088 * r10101; return r10102; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r10103 = 3; double r10104 = r10103 / r10103; double r10105 = c; double r10106 = b; double r10107 = -r10106; double r10108 = -r10105; double r10109 = a; double r10110 = r10103 * r10109; double r10111 = r10106 * r10106; double r10112 = fma(r10108, r10110, r10111); double r10113 = pow(r10112, r10103); double r10114 = cbrt(r10113); double r10115 = sqrt(r10114); double r10116 = r10107 - r10115; double r10117 = r10105 / r10116; double r10118 = r10104 * r10117; return r10118; } 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 r10119, r10120, r10121, r10122, r10123, r10124, r10125, r10126, r10127, r10128, r10129, r10130; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10119); mpfr_init(r10120); mpfr_init(r10121); mpfr_init_set_str(r10122, "3", 10, MPFR_RNDN); mpfr_init(r10123); mpfr_init(r10124); mpfr_init(r10125); mpfr_init(r10126); mpfr_init(r10127); mpfr_init(r10128); mpfr_init(r10129); mpfr_init(r10130); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r10119, b, MPFR_RNDN); mpfr_neg(r10120, r10119, MPFR_RNDN); mpfr_mul(r10121, r10119, r10119, MPFR_RNDN); ; mpfr_set_d(r10123, a, MPFR_RNDN); mpfr_mul(r10124, r10122, r10123, MPFR_RNDN); mpfr_set_d(r10125, c, MPFR_RNDN); mpfr_mul(r10126, r10124, r10125, MPFR_RNDN); mpfr_sub(r10127, r10121, r10126, MPFR_RNDN); mpfr_sqrt(r10128, r10127, MPFR_RNDN); mpfr_add(r10129, r10120, r10128, MPFR_RNDN); mpfr_div(r10130, r10129, r10124, MPFR_RNDN); return mpfr_get_d(r10130, MPFR_RNDN); } static mpfr_t r10131, r10132, r10133, r10134, r10135, r10136, r10137, r10138, r10139, r10140, r10141, r10142, r10143, r10144, r10145, r10146; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10131, "3", 10, MPFR_RNDN); mpfr_init(r10132); mpfr_init(r10133); mpfr_init(r10134); mpfr_init(r10135); mpfr_init(r10136); mpfr_init(r10137); mpfr_init(r10138); mpfr_init(r10139); mpfr_init(r10140); mpfr_init(r10141); mpfr_init(r10142); mpfr_init(r10143); mpfr_init(r10144); mpfr_init(r10145); mpfr_init(r10146); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_div(r10132, r10131, r10131, MPFR_RNDN); mpfr_set_d(r10133, c, MPFR_RNDN); mpfr_set_d(r10134, b, MPFR_RNDN); mpfr_neg(r10135, r10134, MPFR_RNDN); mpfr_neg(r10136, r10133, MPFR_RNDN); mpfr_set_d(r10137, a, MPFR_RNDN); mpfr_mul(r10138, r10131, r10137, MPFR_RNDN); mpfr_mul(r10139, r10134, r10134, MPFR_RNDN); mpfr_fma(r10140, r10136, r10138, r10139, MPFR_RNDN); mpfr_pow(r10141, r10140, r10131, MPFR_RNDN); mpfr_cbrt(r10142, r10141, MPFR_RNDN); mpfr_sqrt(r10143, r10142, MPFR_RNDN); mpfr_sub(r10144, r10135, r10143, MPFR_RNDN); mpfr_div(r10145, r10133, r10144, MPFR_RNDN); mpfr_mul(r10146, r10132, r10145, MPFR_RNDN); return mpfr_get_d(r10146, MPFR_RNDN); } static mpfr_t r10147, r10148, r10149, r10150, r10151, r10152, r10153, r10154, r10155, r10156, r10157, r10158, r10159, r10160, r10161, r10162; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10147, "3", 10, MPFR_RNDN); mpfr_init(r10148); mpfr_init(r10149); mpfr_init(r10150); mpfr_init(r10151); mpfr_init(r10152); mpfr_init(r10153); mpfr_init(r10154); mpfr_init(r10155); mpfr_init(r10156); mpfr_init(r10157); mpfr_init(r10158); mpfr_init(r10159); mpfr_init(r10160); mpfr_init(r10161); mpfr_init(r10162); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_div(r10148, r10147, r10147, MPFR_RNDN); mpfr_set_d(r10149, c, MPFR_RNDN); mpfr_set_d(r10150, b, MPFR_RNDN); mpfr_neg(r10151, r10150, MPFR_RNDN); mpfr_neg(r10152, r10149, MPFR_RNDN); mpfr_set_d(r10153, a, MPFR_RNDN); mpfr_mul(r10154, r10147, r10153, MPFR_RNDN); mpfr_mul(r10155, r10150, r10150, MPFR_RNDN); mpfr_fma(r10156, r10152, r10154, r10155, MPFR_RNDN); mpfr_pow(r10157, r10156, r10147, MPFR_RNDN); mpfr_cbrt(r10158, r10157, MPFR_RNDN); mpfr_sqrt(r10159, r10158, MPFR_RNDN); mpfr_sub(r10160, r10151, r10159, MPFR_RNDN); mpfr_div(r10161, r10149, r10160, MPFR_RNDN); mpfr_mul(r10162, r10148, r10161, MPFR_RNDN); return mpfr_get_d(r10162, MPFR_RNDN); }