#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 r23470 = b; float r23471 = -r23470; float r23472 = r23470 * r23470; float r23473 = 3; float r23474 = a; float r23475 = r23473 * r23474; float r23476 = c; float r23477 = r23475 * r23476; float r23478 = r23472 - r23477; float r23479 = sqrt(r23478); float r23480 = r23471 + r23479; float r23481 = r23480 / r23475; return r23481; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23482 = b; double r23483 = -r23482; double r23484 = r23482 * r23482; double r23485 = 3; double r23486 = a; double r23487 = r23485 * r23486; double r23488 = c; double r23489 = r23487 * r23488; double r23490 = r23484 - r23489; double r23491 = sqrt(r23490); double r23492 = r23483 + r23491; double r23493 = r23492 / r23487; return r23493; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23494 = 1; float r23495 = c; float r23496 = r23494 / r23495; float r23497 = a; float r23498 = -r23497; float r23499 = r23497 / r23498; float r23500 = r23496 * r23499; float r23501 = b; float r23502 = r23501 * r23501; float r23503 = r23497 * r23495; float r23504 = 3; float r23505 = r23503 * r23504; float r23506 = r23502 - r23505; float r23507 = sqrt(r23506); float r23508 = r23507 + r23501; float r23509 = r23500 * r23508; float r23510 = r23494 / r23509; return r23510; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23511 = 1; double r23512 = c; double r23513 = r23511 / r23512; double r23514 = a; double r23515 = -r23514; double r23516 = r23514 / r23515; double r23517 = r23513 * r23516; double r23518 = b; double r23519 = r23518 * r23518; double r23520 = r23514 * r23512; double r23521 = 3; double r23522 = r23520 * r23521; double r23523 = r23519 - r23522; double r23524 = sqrt(r23523); double r23525 = r23524 + r23518; double r23526 = r23517 * r23525; double r23527 = r23511 / r23526; return r23527; } 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 r23528, r23529, r23530, r23531, r23532, r23533, r23534, r23535, r23536, r23537, r23538, r23539; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r23528); mpfr_init(r23529); mpfr_init(r23530); mpfr_init_set_str(r23531, "3", 10, MPFR_RNDN); mpfr_init(r23532); mpfr_init(r23533); mpfr_init(r23534); mpfr_init(r23535); mpfr_init(r23536); mpfr_init(r23537); mpfr_init(r23538); mpfr_init(r23539); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23528, b, MPFR_RNDN); mpfr_neg(r23529, r23528, MPFR_RNDN); mpfr_mul(r23530, r23528, r23528, MPFR_RNDN); ; mpfr_set_d(r23532, a, MPFR_RNDN); mpfr_mul(r23533, r23531, r23532, MPFR_RNDN); mpfr_set_d(r23534, c, MPFR_RNDN); mpfr_mul(r23535, r23533, r23534, MPFR_RNDN); mpfr_sub(r23536, r23530, r23535, MPFR_RNDN); mpfr_sqrt(r23537, r23536, MPFR_RNDN); mpfr_add(r23538, r23529, r23537, MPFR_RNDN); mpfr_div(r23539, r23538, r23533, MPFR_RNDN); return mpfr_get_d(r23539, MPFR_RNDN); } static mpfr_t r23540, r23541, r23542, r23543, r23544, r23545, r23546, r23547, r23548, r23549, r23550, r23551, r23552, r23553, r23554, r23555, r23556; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init_set_str(r23540, "1", 10, MPFR_RNDN); mpfr_init(r23541); mpfr_init(r23542); mpfr_init(r23543); mpfr_init(r23544); mpfr_init(r23545); mpfr_init(r23546); mpfr_init(r23547); mpfr_init(r23548); mpfr_init(r23549); mpfr_init_set_str(r23550, "3", 10, MPFR_RNDN); mpfr_init(r23551); mpfr_init(r23552); mpfr_init(r23553); mpfr_init(r23554); mpfr_init(r23555); mpfr_init(r23556); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r23541, c, MPFR_RNDN); mpfr_div(r23542, r23540, r23541, MPFR_RNDN); mpfr_set_d(r23543, a, MPFR_RNDN); mpfr_neg(r23544, r23543, MPFR_RNDN); mpfr_div(r23545, r23543, r23544, MPFR_RNDN); mpfr_mul(r23546, r23542, r23545, MPFR_RNDN); mpfr_set_d(r23547, b, MPFR_RNDN); mpfr_mul(r23548, r23547, r23547, MPFR_RNDN); mpfr_mul(r23549, r23543, r23541, MPFR_RNDN); ; mpfr_mul(r23551, r23549, r23550, MPFR_RNDN); mpfr_sub(r23552, r23548, r23551, MPFR_RNDN); mpfr_sqrt(r23553, r23552, MPFR_RNDN); mpfr_add(r23554, r23553, r23547, MPFR_RNDN); mpfr_mul(r23555, r23546, r23554, MPFR_RNDN); mpfr_div(r23556, r23540, r23555, MPFR_RNDN); return mpfr_get_d(r23556, MPFR_RNDN); } static mpfr_t r23557, r23558, r23559, r23560, r23561, r23562, r23563, r23564, r23565, r23566, r23567, r23568, r23569, r23570, r23571, r23572, r23573; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init_set_str(r23557, "1", 10, MPFR_RNDN); mpfr_init(r23558); mpfr_init(r23559); mpfr_init(r23560); mpfr_init(r23561); mpfr_init(r23562); mpfr_init(r23563); mpfr_init(r23564); mpfr_init(r23565); mpfr_init(r23566); mpfr_init_set_str(r23567, "3", 10, MPFR_RNDN); mpfr_init(r23568); mpfr_init(r23569); mpfr_init(r23570); mpfr_init(r23571); mpfr_init(r23572); mpfr_init(r23573); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { ; mpfr_set_d(r23558, c, MPFR_RNDN); mpfr_div(r23559, r23557, r23558, MPFR_RNDN); mpfr_set_d(r23560, a, MPFR_RNDN); mpfr_neg(r23561, r23560, MPFR_RNDN); mpfr_div(r23562, r23560, r23561, MPFR_RNDN); mpfr_mul(r23563, r23559, r23562, MPFR_RNDN); mpfr_set_d(r23564, b, MPFR_RNDN); mpfr_mul(r23565, r23564, r23564, MPFR_RNDN); mpfr_mul(r23566, r23560, r23558, MPFR_RNDN); ; mpfr_mul(r23568, r23566, r23567, MPFR_RNDN); mpfr_sub(r23569, r23565, r23568, MPFR_RNDN); mpfr_sqrt(r23570, r23569, MPFR_RNDN); mpfr_add(r23571, r23570, r23564, MPFR_RNDN); mpfr_mul(r23572, r23563, r23571, MPFR_RNDN); mpfr_div(r23573, r23557, r23572, MPFR_RNDN); return mpfr_get_d(r23573, MPFR_RNDN); }