#include #include #include #include #include char *name = "Cubic critical"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r26279 = b; float r26280 = -r26279; float r26281 = r26279 * r26279; float r26282 = 3; float r26283 = a; float r26284 = r26282 * r26283; float r26285 = c; float r26286 = r26284 * r26285; float r26287 = r26281 - r26286; float r26288 = sqrt(r26287); float r26289 = r26280 + r26288; float r26290 = r26289 / r26284; return r26290; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26291 = b; double r26292 = -r26291; double r26293 = r26291 * r26291; double r26294 = 3; double r26295 = a; double r26296 = r26294 * r26295; double r26297 = c; double r26298 = r26296 * r26297; double r26299 = r26293 - r26298; double r26300 = sqrt(r26299); double r26301 = r26292 + r26300; double r26302 = r26301 / r26296; return r26302; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26303 = 1; float r26304 = a; float r26305 = 3; float r26306 = r26304 * r26305; float r26307 = r26303 / r26306; float r26308 = c; float r26309 = -r26308; float r26310 = b; float r26311 = r26310 * r26310; float r26312 = fma(r26306, r26309, r26311); float r26313 = sqrt(r26312); float r26314 = r26313 - r26310; float r26315 = r26307 * r26314; float r26316 = -inf.0; bool r26317 = r26315 <= r26316; float r26318 = r26310 + r26310; float r26319 = r26309 / r26318; float r26320 = -1.5982339424674708e-255; bool r26321 = r26315 <= r26320; float r26322 = r26314 / r26305; float r26323 = r26322 / r26304; float r26324 = 2.855099414002002e-254; bool r26325 = r26315 <= r26324; float r26326 = 3.4191119758542845e+293; bool r26327 = r26315 <= r26326; float r26328 = r26327 ? r26315 : r26319; float r26329 = r26325 ? r26319 : r26328; float r26330 = r26321 ? r26323 : r26329; float r26331 = r26317 ? r26319 : r26330; return r26331; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26332 = 1; double r26333 = a; double r26334 = 3; double r26335 = r26333 * r26334; double r26336 = r26332 / r26335; double r26337 = c; double r26338 = -r26337; double r26339 = b; double r26340 = r26339 * r26339; double r26341 = fma(r26335, r26338, r26340); double r26342 = sqrt(r26341); double r26343 = r26342 - r26339; double r26344 = r26336 * r26343; double r26345 = -inf.0; bool r26346 = r26344 <= r26345; double r26347 = r26339 + r26339; double r26348 = r26338 / r26347; double r26349 = -1.5982339424674708e-255; bool r26350 = r26344 <= r26349; double r26351 = r26343 / r26334; double r26352 = r26351 / r26333; double r26353 = 2.855099414002002e-254; bool r26354 = r26344 <= r26353; double r26355 = 3.4191119758542845e+293; bool r26356 = r26344 <= r26355; double r26357 = r26356 ? r26344 : r26348; double r26358 = r26354 ? r26348 : r26357; double r26359 = r26350 ? r26352 : r26358; double r26360 = r26346 ? r26348 : r26359; return r26360; } 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 r26361, r26362, r26363, r26364, r26365, r26366, r26367, r26368, r26369, r26370, r26371, r26372; void setup_mpfr_f_im() { mpfr_set_default_prec(3408); mpfr_init(r26361); mpfr_init(r26362); mpfr_init(r26363); mpfr_init_set_str(r26364, "3", 10, MPFR_RNDN); mpfr_init(r26365); mpfr_init(r26366); mpfr_init(r26367); mpfr_init(r26368); mpfr_init(r26369); mpfr_init(r26370); mpfr_init(r26371); mpfr_init(r26372); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26361, b, MPFR_RNDN); mpfr_neg(r26362, r26361, MPFR_RNDN); mpfr_mul(r26363, r26361, r26361, MPFR_RNDN); ; mpfr_set_d(r26365, a, MPFR_RNDN); mpfr_mul(r26366, r26364, r26365, MPFR_RNDN); mpfr_set_d(r26367, c, MPFR_RNDN); mpfr_mul(r26368, r26366, r26367, MPFR_RNDN); mpfr_sub(r26369, r26363, r26368, MPFR_RNDN); mpfr_sqrt(r26370, r26369, MPFR_RNDN); mpfr_add(r26371, r26362, r26370, MPFR_RNDN); mpfr_div(r26372, r26371, r26366, MPFR_RNDN); return mpfr_get_d(r26372, MPFR_RNDN); } static mpfr_t r26373, r26374, r26375, r26376, r26377, r26378, r26379, r26380, r26381, r26382, r26383, r26384, r26385, r26386, r26387, r26388, r26389, r26390, r26391, r26392, r26393, r26394, r26395, r26396, r26397, r26398, r26399, r26400, r26401; void setup_mpfr_f_fm() { mpfr_set_default_prec(3408); mpfr_init_set_str(r26373, "1", 10, MPFR_RNDN); mpfr_init(r26374); mpfr_init_set_str(r26375, "3", 10, MPFR_RNDN); mpfr_init(r26376); mpfr_init(r26377); mpfr_init(r26378); mpfr_init(r26379); mpfr_init(r26380); mpfr_init(r26381); mpfr_init(r26382); mpfr_init(r26383); mpfr_init(r26384); mpfr_init(r26385); mpfr_init_set_str(r26386, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26387); mpfr_init(r26388); mpfr_init(r26389); mpfr_init_set_str(r26390, "-1.5982339424674708e-255", 10, MPFR_RNDN); mpfr_init(r26391); mpfr_init(r26392); mpfr_init(r26393); mpfr_init_set_str(r26394, "2.855099414002002e-254", 10, MPFR_RNDN); mpfr_init(r26395); mpfr_init_set_str(r26396, "3.4191119758542845e+293", 10, MPFR_RNDN); mpfr_init(r26397); mpfr_init(r26398); mpfr_init(r26399); mpfr_init(r26400); mpfr_init(r26401); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r26374, a, MPFR_RNDN); ; mpfr_mul(r26376, r26374, r26375, MPFR_RNDN); mpfr_div(r26377, r26373, r26376, MPFR_RNDN); mpfr_set_d(r26378, c, MPFR_RNDN); mpfr_neg(r26379, r26378, MPFR_RNDN); mpfr_set_d(r26380, b, MPFR_RNDN); mpfr_mul(r26381, r26380, r26380, MPFR_RNDN); mpfr_fma(r26382, r26376, r26379, r26381, MPFR_RNDN); mpfr_sqrt(r26383, r26382, MPFR_RNDN); mpfr_sub(r26384, r26383, r26380, MPFR_RNDN); mpfr_mul(r26385, r26377, r26384, MPFR_RNDN); ; mpfr_set_si(r26387, mpfr_cmp(r26385, r26386) <= 0, MPFR_RNDN); mpfr_add(r26388, r26380, r26380, MPFR_RNDN); mpfr_div(r26389, r26379, r26388, MPFR_RNDN); ; mpfr_set_si(r26391, mpfr_cmp(r26385, r26390) <= 0, MPFR_RNDN); mpfr_div(r26392, r26384, r26375, MPFR_RNDN); mpfr_div(r26393, r26392, r26374, MPFR_RNDN); ; mpfr_set_si(r26395, mpfr_cmp(r26385, r26394) <= 0, MPFR_RNDN); ; mpfr_set_si(r26397, mpfr_cmp(r26385, r26396) <= 0, MPFR_RNDN); if (mpfr_get_si(r26397, MPFR_RNDN)) { mpfr_set(r26398, r26385, MPFR_RNDN); } else { mpfr_set(r26398, r26389, MPFR_RNDN); }; if (mpfr_get_si(r26395, MPFR_RNDN)) { mpfr_set(r26399, r26389, MPFR_RNDN); } else { mpfr_set(r26399, r26398, MPFR_RNDN); }; if (mpfr_get_si(r26391, MPFR_RNDN)) { mpfr_set(r26400, r26393, MPFR_RNDN); } else { mpfr_set(r26400, r26399, MPFR_RNDN); }; if (mpfr_get_si(r26387, MPFR_RNDN)) { mpfr_set(r26401, r26389, MPFR_RNDN); } else { mpfr_set(r26401, r26400, MPFR_RNDN); }; return mpfr_get_d(r26401, MPFR_RNDN); } static mpfr_t r26402, r26403, r26404, r26405, r26406, r26407, r26408, r26409, r26410, r26411, r26412, r26413, r26414, r26415, r26416, r26417, r26418, r26419, r26420, r26421, r26422, r26423, r26424, r26425, r26426, r26427, r26428, r26429, r26430; void setup_mpfr_f_dm() { mpfr_set_default_prec(3408); mpfr_init_set_str(r26402, "1", 10, MPFR_RNDN); mpfr_init(r26403); mpfr_init_set_str(r26404, "3", 10, MPFR_RNDN); mpfr_init(r26405); mpfr_init(r26406); mpfr_init(r26407); mpfr_init(r26408); mpfr_init(r26409); mpfr_init(r26410); mpfr_init(r26411); mpfr_init(r26412); mpfr_init(r26413); mpfr_init(r26414); mpfr_init_set_str(r26415, "-inf.0", 10, MPFR_RNDN); mpfr_init(r26416); mpfr_init(r26417); mpfr_init(r26418); mpfr_init_set_str(r26419, "-1.5982339424674708e-255", 10, MPFR_RNDN); mpfr_init(r26420); mpfr_init(r26421); mpfr_init(r26422); mpfr_init_set_str(r26423, "2.855099414002002e-254", 10, MPFR_RNDN); mpfr_init(r26424); mpfr_init_set_str(r26425, "3.4191119758542845e+293", 10, MPFR_RNDN); mpfr_init(r26426); mpfr_init(r26427); mpfr_init(r26428); mpfr_init(r26429); mpfr_init(r26430); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r26403, a, MPFR_RNDN); ; mpfr_mul(r26405, r26403, r26404, MPFR_RNDN); mpfr_div(r26406, r26402, r26405, MPFR_RNDN); mpfr_set_d(r26407, c, MPFR_RNDN); mpfr_neg(r26408, r26407, MPFR_RNDN); mpfr_set_d(r26409, b, MPFR_RNDN); mpfr_mul(r26410, r26409, r26409, MPFR_RNDN); mpfr_fma(r26411, r26405, r26408, r26410, MPFR_RNDN); mpfr_sqrt(r26412, r26411, MPFR_RNDN); mpfr_sub(r26413, r26412, r26409, MPFR_RNDN); mpfr_mul(r26414, r26406, r26413, MPFR_RNDN); ; mpfr_set_si(r26416, mpfr_cmp(r26414, r26415) <= 0, MPFR_RNDN); mpfr_add(r26417, r26409, r26409, MPFR_RNDN); mpfr_div(r26418, r26408, r26417, MPFR_RNDN); ; mpfr_set_si(r26420, mpfr_cmp(r26414, r26419) <= 0, MPFR_RNDN); mpfr_div(r26421, r26413, r26404, MPFR_RNDN); mpfr_div(r26422, r26421, r26403, MPFR_RNDN); ; mpfr_set_si(r26424, mpfr_cmp(r26414, r26423) <= 0, MPFR_RNDN); ; mpfr_set_si(r26426, mpfr_cmp(r26414, r26425) <= 0, MPFR_RNDN); if (mpfr_get_si(r26426, MPFR_RNDN)) { mpfr_set(r26427, r26414, MPFR_RNDN); } else { mpfr_set(r26427, r26418, MPFR_RNDN); }; if (mpfr_get_si(r26424, MPFR_RNDN)) { mpfr_set(r26428, r26418, MPFR_RNDN); } else { mpfr_set(r26428, r26427, MPFR_RNDN); }; if (mpfr_get_si(r26420, MPFR_RNDN)) { mpfr_set(r26429, r26422, MPFR_RNDN); } else { mpfr_set(r26429, r26428, MPFR_RNDN); }; if (mpfr_get_si(r26416, MPFR_RNDN)) { mpfr_set(r26430, r26418, MPFR_RNDN); } else { mpfr_set(r26430, r26429, MPFR_RNDN); }; return mpfr_get_d(r26430, MPFR_RNDN); }