#include #include #include #include #include char *name = "Cubic critical, narrow range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r23381 = b; float r23382 = -r23381; float r23383 = r23381 * r23381; float r23384 = 3; float r23385 = a; float r23386 = r23384 * r23385; float r23387 = c; float r23388 = r23386 * r23387; float r23389 = r23383 - r23388; float r23390 = sqrt(r23389); float r23391 = r23382 + r23390; float r23392 = r23391 / r23386; return r23392; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23393 = b; double r23394 = -r23393; double r23395 = r23393 * r23393; double r23396 = 3; double r23397 = a; double r23398 = r23396 * r23397; double r23399 = c; double r23400 = r23398 * r23399; double r23401 = r23395 - r23400; double r23402 = sqrt(r23401); double r23403 = r23394 + r23402; double r23404 = r23403 / r23398; return r23404; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23405 = a; float r23406 = 3; float r23407 = r23405 / r23406; float r23408 = c; float r23409 = r23407 * r23408; float r23410 = -r23406; float r23411 = r23406 * r23405; float r23412 = -r23408; float r23413 = b; float r23414 = r23413 * r23413; float r23415 = fma(r23411, r23412, r23414); float r23416 = sqrt(r23415); float r23417 = r23416 + r23413; float r23418 = r23410 / r23417; float r23419 = r23418 / r23405; float r23420 = r23409 * r23419; return r23420; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23421 = a; double r23422 = 3; double r23423 = r23421 / r23422; double r23424 = c; double r23425 = r23423 * r23424; double r23426 = -r23422; double r23427 = r23422 * r23421; double r23428 = -r23424; double r23429 = b; double r23430 = r23429 * r23429; double r23431 = fma(r23427, r23428, r23430); double r23432 = sqrt(r23431); double r23433 = r23432 + r23429; double r23434 = r23426 / r23433; double r23435 = r23434 / r23421; double r23436 = r23425 * r23435; return r23436; } 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 r23437, r23438, r23439, r23440, r23441, r23442, r23443, r23444, r23445, r23446, r23447, r23448; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r23437); mpfr_init(r23438); mpfr_init(r23439); mpfr_init_set_str(r23440, "3", 10, MPFR_RNDN); mpfr_init(r23441); mpfr_init(r23442); mpfr_init(r23443); mpfr_init(r23444); mpfr_init(r23445); mpfr_init(r23446); mpfr_init(r23447); mpfr_init(r23448); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23437, b, MPFR_RNDN); mpfr_neg(r23438, r23437, MPFR_RNDN); mpfr_mul(r23439, r23437, r23437, MPFR_RNDN); ; mpfr_set_d(r23441, a, MPFR_RNDN); mpfr_mul(r23442, r23440, r23441, MPFR_RNDN); mpfr_set_d(r23443, c, MPFR_RNDN); mpfr_mul(r23444, r23442, r23443, MPFR_RNDN); mpfr_sub(r23445, r23439, r23444, MPFR_RNDN); mpfr_sqrt(r23446, r23445, MPFR_RNDN); mpfr_add(r23447, r23438, r23446, MPFR_RNDN); mpfr_div(r23448, r23447, r23442, MPFR_RNDN); return mpfr_get_d(r23448, MPFR_RNDN); } static mpfr_t r23449, r23450, r23451, r23452, r23453, r23454, r23455, r23456, r23457, r23458, r23459, r23460, r23461, r23462, r23463, r23464; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23449); mpfr_init_set_str(r23450, "3", 10, MPFR_RNDN); mpfr_init(r23451); mpfr_init(r23452); mpfr_init(r23453); mpfr_init(r23454); mpfr_init(r23455); mpfr_init(r23456); mpfr_init(r23457); mpfr_init(r23458); mpfr_init(r23459); mpfr_init(r23460); mpfr_init(r23461); mpfr_init(r23462); mpfr_init(r23463); mpfr_init(r23464); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23449, a, MPFR_RNDN); ; mpfr_div(r23451, r23449, r23450, MPFR_RNDN); mpfr_set_d(r23452, c, MPFR_RNDN); mpfr_mul(r23453, r23451, r23452, MPFR_RNDN); mpfr_neg(r23454, r23450, MPFR_RNDN); mpfr_mul(r23455, r23450, r23449, MPFR_RNDN); mpfr_neg(r23456, r23452, MPFR_RNDN); mpfr_set_d(r23457, b, MPFR_RNDN); mpfr_mul(r23458, r23457, r23457, MPFR_RNDN); mpfr_fma(r23459, r23455, r23456, r23458, MPFR_RNDN); mpfr_sqrt(r23460, r23459, MPFR_RNDN); mpfr_add(r23461, r23460, r23457, MPFR_RNDN); mpfr_div(r23462, r23454, r23461, MPFR_RNDN); mpfr_div(r23463, r23462, r23449, MPFR_RNDN); mpfr_mul(r23464, r23453, r23463, MPFR_RNDN); return mpfr_get_d(r23464, MPFR_RNDN); } static mpfr_t r23465, r23466, r23467, r23468, r23469, r23470, r23471, r23472, r23473, r23474, r23475, r23476, r23477, r23478, r23479, r23480; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23465); mpfr_init_set_str(r23466, "3", 10, MPFR_RNDN); mpfr_init(r23467); mpfr_init(r23468); mpfr_init(r23469); mpfr_init(r23470); mpfr_init(r23471); mpfr_init(r23472); mpfr_init(r23473); mpfr_init(r23474); mpfr_init(r23475); mpfr_init(r23476); mpfr_init(r23477); mpfr_init(r23478); mpfr_init(r23479); mpfr_init(r23480); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23465, a, MPFR_RNDN); ; mpfr_div(r23467, r23465, r23466, MPFR_RNDN); mpfr_set_d(r23468, c, MPFR_RNDN); mpfr_mul(r23469, r23467, r23468, MPFR_RNDN); mpfr_neg(r23470, r23466, MPFR_RNDN); mpfr_mul(r23471, r23466, r23465, MPFR_RNDN); mpfr_neg(r23472, r23468, MPFR_RNDN); mpfr_set_d(r23473, b, MPFR_RNDN); mpfr_mul(r23474, r23473, r23473, MPFR_RNDN); mpfr_fma(r23475, r23471, r23472, r23474, MPFR_RNDN); mpfr_sqrt(r23476, r23475, MPFR_RNDN); mpfr_add(r23477, r23476, r23473, MPFR_RNDN); mpfr_div(r23478, r23470, r23477, MPFR_RNDN); mpfr_div(r23479, r23478, r23465, MPFR_RNDN); mpfr_mul(r23480, r23469, r23479, MPFR_RNDN); return mpfr_get_d(r23480, MPFR_RNDN); }