#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 r26286 = b; float r26287 = -r26286; float r26288 = r26286 * r26286; float r26289 = 3; float r26290 = a; float r26291 = r26289 * r26290; float r26292 = c; float r26293 = r26291 * r26292; float r26294 = r26288 - r26293; float r26295 = sqrt(r26294); float r26296 = r26287 + r26295; float r26297 = r26296 / r26291; return r26297; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26298 = b; double r26299 = -r26298; double r26300 = r26298 * r26298; double r26301 = 3; double r26302 = a; double r26303 = r26301 * r26302; double r26304 = c; double r26305 = r26303 * r26304; double r26306 = r26300 - r26305; double r26307 = sqrt(r26306); double r26308 = r26299 + r26307; double r26309 = r26308 / r26303; return r26309; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26310 = c; float r26311 = -r26310; float r26312 = a; float r26313 = r26311 * r26312; float r26314 = 1; float r26315 = 3; float r26316 = r26315 * r26312; float r26317 = b; float r26318 = r26317 * r26317; float r26319 = fma(r26316, r26311, r26318); float r26320 = sqrt(r26319); float r26321 = r26320 + r26317; float r26322 = r26314 / r26321; float r26323 = r26322 / r26312; float r26324 = r26313 * r26323; return r26324; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26325 = c; double r26326 = -r26325; double r26327 = a; double r26328 = r26326 * r26327; double r26329 = 1; double r26330 = 3; double r26331 = r26330 * r26327; double r26332 = b; double r26333 = r26332 * r26332; double r26334 = fma(r26331, r26326, r26333); double r26335 = sqrt(r26334); double r26336 = r26335 + r26332; double r26337 = r26329 / r26336; double r26338 = r26337 / r26327; double r26339 = r26328 * r26338; return r26339; } 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 r26340, r26341, r26342, r26343, r26344, r26345, r26346, r26347, r26348, r26349, r26350, r26351; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26340); mpfr_init(r26341); mpfr_init(r26342); mpfr_init_set_str(r26343, "3", 10, MPFR_RNDN); mpfr_init(r26344); mpfr_init(r26345); mpfr_init(r26346); mpfr_init(r26347); mpfr_init(r26348); mpfr_init(r26349); mpfr_init(r26350); mpfr_init(r26351); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26340, b, MPFR_RNDN); mpfr_neg(r26341, r26340, MPFR_RNDN); mpfr_mul(r26342, r26340, r26340, MPFR_RNDN); ; mpfr_set_d(r26344, a, MPFR_RNDN); mpfr_mul(r26345, r26343, r26344, MPFR_RNDN); mpfr_set_d(r26346, c, MPFR_RNDN); mpfr_mul(r26347, r26345, r26346, MPFR_RNDN); mpfr_sub(r26348, r26342, r26347, MPFR_RNDN); mpfr_sqrt(r26349, r26348, MPFR_RNDN); mpfr_add(r26350, r26341, r26349, MPFR_RNDN); mpfr_div(r26351, r26350, r26345, MPFR_RNDN); return mpfr_get_d(r26351, MPFR_RNDN); } static mpfr_t r26352, r26353, r26354, r26355, r26356, r26357, r26358, r26359, r26360, r26361, r26362, r26363, r26364, r26365, r26366; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26352); mpfr_init(r26353); mpfr_init(r26354); mpfr_init(r26355); mpfr_init_set_str(r26356, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26357, "3", 10, MPFR_RNDN); mpfr_init(r26358); mpfr_init(r26359); mpfr_init(r26360); mpfr_init(r26361); mpfr_init(r26362); mpfr_init(r26363); mpfr_init(r26364); mpfr_init(r26365); mpfr_init(r26366); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26352, c, MPFR_RNDN); mpfr_neg(r26353, r26352, MPFR_RNDN); mpfr_set_d(r26354, a, MPFR_RNDN); mpfr_mul(r26355, r26353, r26354, MPFR_RNDN); ; ; mpfr_mul(r26358, r26357, r26354, MPFR_RNDN); mpfr_set_d(r26359, b, MPFR_RNDN); mpfr_mul(r26360, r26359, r26359, MPFR_RNDN); mpfr_fma(r26361, r26358, r26353, r26360, MPFR_RNDN); mpfr_sqrt(r26362, r26361, MPFR_RNDN); mpfr_add(r26363, r26362, r26359, MPFR_RNDN); mpfr_div(r26364, r26356, r26363, MPFR_RNDN); mpfr_div(r26365, r26364, r26354, MPFR_RNDN); mpfr_mul(r26366, r26355, r26365, MPFR_RNDN); return mpfr_get_d(r26366, MPFR_RNDN); } static mpfr_t r26367, r26368, r26369, r26370, r26371, r26372, r26373, r26374, r26375, r26376, r26377, r26378, r26379, r26380, r26381; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26367); mpfr_init(r26368); mpfr_init(r26369); mpfr_init(r26370); mpfr_init_set_str(r26371, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26372, "3", 10, MPFR_RNDN); mpfr_init(r26373); mpfr_init(r26374); mpfr_init(r26375); mpfr_init(r26376); mpfr_init(r26377); mpfr_init(r26378); mpfr_init(r26379); mpfr_init(r26380); mpfr_init(r26381); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26367, c, MPFR_RNDN); mpfr_neg(r26368, r26367, MPFR_RNDN); mpfr_set_d(r26369, a, MPFR_RNDN); mpfr_mul(r26370, r26368, r26369, MPFR_RNDN); ; ; mpfr_mul(r26373, r26372, r26369, MPFR_RNDN); mpfr_set_d(r26374, b, MPFR_RNDN); mpfr_mul(r26375, r26374, r26374, MPFR_RNDN); mpfr_fma(r26376, r26373, r26368, r26375, MPFR_RNDN); mpfr_sqrt(r26377, r26376, MPFR_RNDN); mpfr_add(r26378, r26377, r26374, MPFR_RNDN); mpfr_div(r26379, r26371, r26378, MPFR_RNDN); mpfr_div(r26380, r26379, r26369, MPFR_RNDN); mpfr_mul(r26381, r26370, r26380, MPFR_RNDN); return mpfr_get_d(r26381, MPFR_RNDN); }