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