#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 r26232 = b; float r26233 = -r26232; float r26234 = r26232 * r26232; float r26235 = 3; float r26236 = a; float r26237 = r26235 * r26236; float r26238 = c; float r26239 = r26237 * r26238; float r26240 = r26234 - r26239; float r26241 = sqrt(r26240); float r26242 = r26233 + r26241; float r26243 = r26242 / r26237; return r26243; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26244 = b; double r26245 = -r26244; double r26246 = r26244 * r26244; double r26247 = 3; double r26248 = a; double r26249 = r26247 * r26248; double r26250 = c; double r26251 = r26249 * r26250; double r26252 = r26246 - r26251; double r26253 = sqrt(r26252); double r26254 = r26245 + r26253; double r26255 = r26254 / r26249; return r26255; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26256 = c; float r26257 = -r26256; float r26258 = a; float r26259 = r26257 * r26258; float r26260 = 1; float r26261 = 3; float r26262 = r26261 * r26258; float r26263 = b; float r26264 = r26263 * r26263; float r26265 = fma(r26262, r26257, r26264); float r26266 = sqrt(r26265); float r26267 = r26266 + r26263; float r26268 = r26260 / r26267; float r26269 = r26268 / r26258; float r26270 = r26259 * r26269; return r26270; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26271 = c; double r26272 = -r26271; double r26273 = a; double r26274 = r26272 * r26273; double r26275 = 1; double r26276 = 3; double r26277 = r26276 * r26273; double r26278 = b; double r26279 = r26278 * r26278; double r26280 = fma(r26277, r26272, r26279); double r26281 = sqrt(r26280); double r26282 = r26281 + r26278; double r26283 = r26275 / r26282; double r26284 = r26283 / r26273; double r26285 = r26274 * r26284; return r26285; } 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 r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295, r26296, r26297; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26286); mpfr_init(r26287); mpfr_init(r26288); mpfr_init_set_str(r26289, "3", 10, MPFR_RNDN); mpfr_init(r26290); mpfr_init(r26291); mpfr_init(r26292); mpfr_init(r26293); mpfr_init(r26294); mpfr_init(r26295); mpfr_init(r26296); mpfr_init(r26297); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26286, b, MPFR_RNDN); mpfr_neg(r26287, r26286, MPFR_RNDN); mpfr_mul(r26288, r26286, r26286, MPFR_RNDN); ; mpfr_set_d(r26290, a, MPFR_RNDN); mpfr_mul(r26291, r26289, r26290, MPFR_RNDN); mpfr_set_d(r26292, c, MPFR_RNDN); mpfr_mul(r26293, r26291, r26292, MPFR_RNDN); mpfr_sub(r26294, r26288, r26293, MPFR_RNDN); mpfr_sqrt(r26295, r26294, MPFR_RNDN); mpfr_add(r26296, r26287, r26295, MPFR_RNDN); mpfr_div(r26297, r26296, r26291, MPFR_RNDN); return mpfr_get_d(r26297, MPFR_RNDN); } static mpfr_t r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26298); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init_set_str(r26302, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26303, "3", 10, MPFR_RNDN); mpfr_init(r26304); mpfr_init(r26305); mpfr_init(r26306); mpfr_init(r26307); mpfr_init(r26308); mpfr_init(r26309); mpfr_init(r26310); mpfr_init(r26311); mpfr_init(r26312); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26298, c, MPFR_RNDN); mpfr_neg(r26299, r26298, MPFR_RNDN); mpfr_set_d(r26300, a, MPFR_RNDN); mpfr_mul(r26301, r26299, r26300, MPFR_RNDN); ; ; mpfr_mul(r26304, r26303, r26300, MPFR_RNDN); mpfr_set_d(r26305, b, MPFR_RNDN); mpfr_mul(r26306, r26305, r26305, MPFR_RNDN); mpfr_fma(r26307, r26304, r26299, r26306, MPFR_RNDN); mpfr_sqrt(r26308, r26307, MPFR_RNDN); mpfr_add(r26309, r26308, r26305, MPFR_RNDN); mpfr_div(r26310, r26302, r26309, MPFR_RNDN); mpfr_div(r26311, r26310, r26300, MPFR_RNDN); mpfr_mul(r26312, r26301, r26311, MPFR_RNDN); return mpfr_get_d(r26312, MPFR_RNDN); } static mpfr_t r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321, r26322, r26323, r26324, r26325, r26326, r26327; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init(r26316); mpfr_init_set_str(r26317, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26318, "3", 10, MPFR_RNDN); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); mpfr_init(r26322); mpfr_init(r26323); mpfr_init(r26324); mpfr_init(r26325); mpfr_init(r26326); mpfr_init(r26327); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26313, c, MPFR_RNDN); mpfr_neg(r26314, r26313, MPFR_RNDN); mpfr_set_d(r26315, a, MPFR_RNDN); mpfr_mul(r26316, r26314, r26315, MPFR_RNDN); ; ; mpfr_mul(r26319, r26318, r26315, MPFR_RNDN); mpfr_set_d(r26320, b, MPFR_RNDN); mpfr_mul(r26321, r26320, r26320, MPFR_RNDN); mpfr_fma(r26322, r26319, r26314, r26321, MPFR_RNDN); mpfr_sqrt(r26323, r26322, MPFR_RNDN); mpfr_add(r26324, r26323, r26320, MPFR_RNDN); mpfr_div(r26325, r26317, r26324, MPFR_RNDN); mpfr_div(r26326, r26325, r26315, MPFR_RNDN); mpfr_mul(r26327, r26316, r26326, MPFR_RNDN); return mpfr_get_d(r26327, MPFR_RNDN); }