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