#include #include #include #include #include char *name = "Cubic critical"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r26208 = b; float r26209 = -r26208; float r26210 = r26208 * r26208; float r26211 = 3; float r26212 = a; float r26213 = r26211 * r26212; float r26214 = c; float r26215 = r26213 * r26214; float r26216 = r26210 - r26215; float r26217 = sqrt(r26216); float r26218 = r26209 + r26217; float r26219 = r26218 / r26213; return r26219; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26220 = b; double r26221 = -r26220; double r26222 = r26220 * r26220; double r26223 = 3; double r26224 = a; double r26225 = r26223 * r26224; double r26226 = c; double r26227 = r26225 * r26226; double r26228 = r26222 - r26227; double r26229 = sqrt(r26228); double r26230 = r26221 + r26229; double r26231 = r26230 / r26225; return r26231; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26232 = b; float r26233 = -r26232; float r26234 = -3.058392593590901e-134; bool r26235 = r26233 <= r26234; float r26236 = c; float r26237 = -r26236; float r26238 = cbrt(r26237); float r26239 = r26238 * r26238; float r26240 = 3; float r26241 = a; float r26242 = r26240 * r26241; float r26243 = r26232 * r26232; float r26244 = fma(r26237, r26242, r26243); float r26245 = sqrt(r26244); float r26246 = r26232 + r26245; float r26247 = sqrt(r26246); float r26248 = r26239 / r26247; float r26249 = r26245 + r26232; float r26250 = sqrt(r26249); float r26251 = r26238 / r26250; float r26252 = r26248 * r26251; float r26253 = fma(r26242, r26237, r26243); float r26254 = sqrt(r26253); float r26255 = r26254 - r26232; float r26256 = r26255 / r26242; float r26257 = r26235 ? r26252 : r26256; return r26257; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26258 = b; double r26259 = -r26258; double r26260 = -3.058392593590901e-134; bool r26261 = r26259 <= r26260; double r26262 = c; double r26263 = -r26262; double r26264 = cbrt(r26263); double r26265 = r26264 * r26264; double r26266 = 3; double r26267 = a; double r26268 = r26266 * r26267; double r26269 = r26258 * r26258; double r26270 = fma(r26263, r26268, r26269); double r26271 = sqrt(r26270); double r26272 = r26258 + r26271; double r26273 = sqrt(r26272); double r26274 = r26265 / r26273; double r26275 = r26271 + r26258; double r26276 = sqrt(r26275); double r26277 = r26264 / r26276; double r26278 = r26274 * r26277; double r26279 = fma(r26268, r26263, r26269); double r26280 = sqrt(r26279); double r26281 = r26280 - r26258; double r26282 = r26281 / r26268; double r26283 = r26261 ? r26278 : r26282; return r26283; } 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 r26284, r26285, r26286, r26287, r26288, r26289, r26290, r26291, r26292, r26293, r26294, r26295; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init(r26284); mpfr_init(r26285); mpfr_init(r26286); mpfr_init_set_str(r26287, "3", 10, MPFR_RNDN); mpfr_init(r26288); mpfr_init(r26289); mpfr_init(r26290); mpfr_init(r26291); mpfr_init(r26292); mpfr_init(r26293); mpfr_init(r26294); mpfr_init(r26295); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26284, b, MPFR_RNDN); mpfr_neg(r26285, r26284, MPFR_RNDN); mpfr_mul(r26286, r26284, r26284, MPFR_RNDN); ; mpfr_set_d(r26288, a, MPFR_RNDN); mpfr_mul(r26289, r26287, r26288, MPFR_RNDN); mpfr_set_d(r26290, c, MPFR_RNDN); mpfr_mul(r26291, r26289, r26290, MPFR_RNDN); mpfr_sub(r26292, r26286, r26291, MPFR_RNDN); mpfr_sqrt(r26293, r26292, MPFR_RNDN); mpfr_add(r26294, r26285, r26293, MPFR_RNDN); mpfr_div(r26295, r26294, r26289, MPFR_RNDN); return mpfr_get_d(r26295, MPFR_RNDN); } static mpfr_t r26296, r26297, r26298, r26299, r26300, r26301, r26302, r26303, r26304, r26305, r26306, r26307, r26308, r26309, r26310, r26311, r26312, r26313, r26314, r26315, r26316, r26317, r26318, r26319, r26320, r26321; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init(r26296); mpfr_init(r26297); mpfr_init_set_str(r26298, "-3.058392593590901e-134", 10, MPFR_RNDN); mpfr_init(r26299); mpfr_init(r26300); mpfr_init(r26301); mpfr_init(r26302); mpfr_init(r26303); mpfr_init_set_str(r26304, "3", 10, MPFR_RNDN); mpfr_init(r26305); mpfr_init(r26306); mpfr_init(r26307); mpfr_init(r26308); mpfr_init(r26309); mpfr_init(r26310); mpfr_init(r26311); mpfr_init(r26312); mpfr_init(r26313); mpfr_init(r26314); mpfr_init(r26315); mpfr_init(r26316); mpfr_init(r26317); mpfr_init(r26318); mpfr_init(r26319); mpfr_init(r26320); mpfr_init(r26321); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26296, b, MPFR_RNDN); mpfr_neg(r26297, r26296, MPFR_RNDN); ; mpfr_set_si(r26299, mpfr_cmp(r26297, r26298) <= 0, MPFR_RNDN); mpfr_set_d(r26300, c, MPFR_RNDN); mpfr_neg(r26301, r26300, MPFR_RNDN); mpfr_cbrt(r26302, r26301, MPFR_RNDN); mpfr_mul(r26303, r26302, r26302, MPFR_RNDN); ; mpfr_set_d(r26305, a, MPFR_RNDN); mpfr_mul(r26306, r26304, r26305, MPFR_RNDN); mpfr_mul(r26307, r26296, r26296, MPFR_RNDN); mpfr_fma(r26308, r26301, r26306, r26307, MPFR_RNDN); mpfr_sqrt(r26309, r26308, MPFR_RNDN); mpfr_add(r26310, r26296, r26309, MPFR_RNDN); mpfr_sqrt(r26311, r26310, MPFR_RNDN); mpfr_div(r26312, r26303, r26311, MPFR_RNDN); mpfr_add(r26313, r26309, r26296, MPFR_RNDN); mpfr_sqrt(r26314, r26313, MPFR_RNDN); mpfr_div(r26315, r26302, r26314, MPFR_RNDN); mpfr_mul(r26316, r26312, r26315, MPFR_RNDN); mpfr_fma(r26317, r26306, r26301, r26307, MPFR_RNDN); mpfr_sqrt(r26318, r26317, MPFR_RNDN); mpfr_sub(r26319, r26318, r26296, MPFR_RNDN); mpfr_div(r26320, r26319, r26306, MPFR_RNDN); if (mpfr_get_si(r26299, MPFR_RNDN)) { mpfr_set(r26321, r26316, MPFR_RNDN); } else { mpfr_set(r26321, r26320, MPFR_RNDN); }; return mpfr_get_d(r26321, MPFR_RNDN); } static mpfr_t r26322, r26323, r26324, r26325, r26326, r26327, r26328, r26329, r26330, r26331, r26332, r26333, r26334, r26335, r26336, r26337, r26338, r26339, r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init(r26322); mpfr_init(r26323); mpfr_init_set_str(r26324, "-3.058392593590901e-134", 10, MPFR_RNDN); mpfr_init(r26325); mpfr_init(r26326); mpfr_init(r26327); mpfr_init(r26328); mpfr_init(r26329); mpfr_init_set_str(r26330, "3", 10, MPFR_RNDN); mpfr_init(r26331); mpfr_init(r26332); mpfr_init(r26333); mpfr_init(r26334); mpfr_init(r26335); mpfr_init(r26336); mpfr_init(r26337); mpfr_init(r26338); mpfr_init(r26339); mpfr_init(r26340); mpfr_init(r26341); mpfr_init(r26342); mpfr_init(r26343); mpfr_init(r26344); mpfr_init(r26345); mpfr_init(r26346); mpfr_init(r26347); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26322, b, MPFR_RNDN); mpfr_neg(r26323, r26322, MPFR_RNDN); ; mpfr_set_si(r26325, mpfr_cmp(r26323, r26324) <= 0, MPFR_RNDN); mpfr_set_d(r26326, c, MPFR_RNDN); mpfr_neg(r26327, r26326, MPFR_RNDN); mpfr_cbrt(r26328, r26327, MPFR_RNDN); mpfr_mul(r26329, r26328, r26328, MPFR_RNDN); ; mpfr_set_d(r26331, a, MPFR_RNDN); mpfr_mul(r26332, r26330, r26331, MPFR_RNDN); mpfr_mul(r26333, r26322, r26322, MPFR_RNDN); mpfr_fma(r26334, r26327, r26332, r26333, MPFR_RNDN); mpfr_sqrt(r26335, r26334, MPFR_RNDN); mpfr_add(r26336, r26322, r26335, MPFR_RNDN); mpfr_sqrt(r26337, r26336, MPFR_RNDN); mpfr_div(r26338, r26329, r26337, MPFR_RNDN); mpfr_add(r26339, r26335, r26322, MPFR_RNDN); mpfr_sqrt(r26340, r26339, MPFR_RNDN); mpfr_div(r26341, r26328, r26340, MPFR_RNDN); mpfr_mul(r26342, r26338, r26341, MPFR_RNDN); mpfr_fma(r26343, r26332, r26327, r26333, MPFR_RNDN); mpfr_sqrt(r26344, r26343, MPFR_RNDN); mpfr_sub(r26345, r26344, r26322, MPFR_RNDN); mpfr_div(r26346, r26345, r26332, MPFR_RNDN); if (mpfr_get_si(r26325, MPFR_RNDN)) { mpfr_set(r26347, r26342, MPFR_RNDN); } else { mpfr_set(r26347, r26346, MPFR_RNDN); }; return mpfr_get_d(r26347, MPFR_RNDN); }