#include #include #include #include #include char *name = "Cubic critical, wide range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r23409 = b; float r23410 = -r23409; float r23411 = r23409 * r23409; float r23412 = 3; float r23413 = a; float r23414 = r23412 * r23413; float r23415 = c; float r23416 = r23414 * r23415; float r23417 = r23411 - r23416; float r23418 = sqrt(r23417); float r23419 = r23410 + r23418; float r23420 = r23419 / r23414; return r23420; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23421 = b; double r23422 = -r23421; double r23423 = r23421 * r23421; double r23424 = 3; double r23425 = a; double r23426 = r23424 * r23425; double r23427 = c; double r23428 = r23426 * r23427; double r23429 = r23423 - r23428; double r23430 = sqrt(r23429); double r23431 = r23422 + r23430; double r23432 = r23431 / r23426; return r23432; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23433 = c; float r23434 = a; float r23435 = 3; float r23436 = r23434 * r23435; float r23437 = r23433 * r23436; float r23438 = b; float r23439 = -r23438; float r23440 = r23438 * r23438; float r23441 = r23435 * r23434; float r23442 = r23441 * r23433; float r23443 = r23440 - r23442; float r23444 = sqrt(r23443); float r23445 = r23439 - r23444; float r23446 = r23437 / r23445; float r23447 = r23446 / r23441; return r23447; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23448 = c; double r23449 = a; double r23450 = 3; double r23451 = r23449 * r23450; double r23452 = r23448 * r23451; double r23453 = b; double r23454 = -r23453; double r23455 = r23453 * r23453; double r23456 = r23450 * r23449; double r23457 = r23456 * r23448; double r23458 = r23455 - r23457; double r23459 = sqrt(r23458); double r23460 = r23454 - r23459; double r23461 = r23452 / r23460; double r23462 = r23461 / r23456; return r23462; } 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 r23463, r23464, r23465, r23466, r23467, r23468, r23469, r23470, r23471, r23472, r23473, r23474; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r23463); mpfr_init(r23464); 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); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23463, b, MPFR_RNDN); mpfr_neg(r23464, r23463, MPFR_RNDN); mpfr_mul(r23465, r23463, r23463, MPFR_RNDN); ; mpfr_set_d(r23467, a, MPFR_RNDN); mpfr_mul(r23468, r23466, r23467, MPFR_RNDN); mpfr_set_d(r23469, c, MPFR_RNDN); mpfr_mul(r23470, r23468, r23469, MPFR_RNDN); mpfr_sub(r23471, r23465, r23470, MPFR_RNDN); mpfr_sqrt(r23472, r23471, MPFR_RNDN); mpfr_add(r23473, r23464, r23472, MPFR_RNDN); mpfr_div(r23474, r23473, r23468, MPFR_RNDN); return mpfr_get_d(r23474, MPFR_RNDN); } static mpfr_t r23475, r23476, r23477, r23478, r23479, r23480, r23481, r23482, r23483, r23484, r23485, r23486, r23487, r23488, r23489; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r23475); mpfr_init(r23476); mpfr_init_set_str(r23477, "3", 10, MPFR_RNDN); mpfr_init(r23478); 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); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23475, c, MPFR_RNDN); mpfr_set_d(r23476, a, MPFR_RNDN); ; mpfr_mul(r23478, r23476, r23477, MPFR_RNDN); mpfr_mul(r23479, r23475, 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, r23477, r23476, MPFR_RNDN); mpfr_mul(r23484, r23483, r23475, 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_div(r23489, r23488, r23483, MPFR_RNDN); return mpfr_get_d(r23489, MPFR_RNDN); } static mpfr_t r23490, r23491, r23492, r23493, r23494, r23495, r23496, r23497, r23498, r23499, r23500, r23501, r23502, r23503, r23504; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r23490); mpfr_init(r23491); mpfr_init_set_str(r23492, "3", 10, MPFR_RNDN); mpfr_init(r23493); mpfr_init(r23494); mpfr_init(r23495); mpfr_init(r23496); mpfr_init(r23497); mpfr_init(r23498); mpfr_init(r23499); mpfr_init(r23500); mpfr_init(r23501); mpfr_init(r23502); mpfr_init(r23503); mpfr_init(r23504); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23490, c, MPFR_RNDN); mpfr_set_d(r23491, a, MPFR_RNDN); ; mpfr_mul(r23493, r23491, r23492, MPFR_RNDN); mpfr_mul(r23494, r23490, r23493, MPFR_RNDN); mpfr_set_d(r23495, b, MPFR_RNDN); mpfr_neg(r23496, r23495, MPFR_RNDN); mpfr_mul(r23497, r23495, r23495, MPFR_RNDN); mpfr_mul(r23498, r23492, r23491, MPFR_RNDN); mpfr_mul(r23499, r23498, r23490, MPFR_RNDN); mpfr_sub(r23500, r23497, r23499, MPFR_RNDN); mpfr_sqrt(r23501, r23500, MPFR_RNDN); mpfr_sub(r23502, r23496, r23501, MPFR_RNDN); mpfr_div(r23503, r23494, r23502, MPFR_RNDN); mpfr_div(r23504, r23503, r23498, MPFR_RNDN); return mpfr_get_d(r23504, MPFR_RNDN); }