#include #include #include #include #include char *name = "Falkner and Boettcher, Equation (22+)"; double f_if(float v) { float r23523 = 4; float r23524 = 3; float r23525 = atan2(1.0, 0.0); float r23526 = r23524 * r23525; float r23527 = 1; float r23528 = v; float r23529 = r23528 * r23528; float r23530 = r23527 - r23529; float r23531 = r23526 * r23530; float r23532 = 2; float r23533 = 6; float r23534 = r23533 * r23529; float r23535 = r23532 - r23534; float r23536 = sqrt(r23535); float r23537 = r23531 * r23536; float r23538 = r23523 / r23537; return r23538; } double f_id(double v) { double r23539 = 4; double r23540 = 3; double r23541 = atan2(1.0, 0.0); double r23542 = r23540 * r23541; double r23543 = 1; double r23544 = v; double r23545 = r23544 * r23544; double r23546 = r23543 - r23545; double r23547 = r23542 * r23546; double r23548 = 2; double r23549 = 6; double r23550 = r23549 * r23545; double r23551 = r23548 - r23550; double r23552 = sqrt(r23551); double r23553 = r23547 * r23552; double r23554 = r23539 / r23553; return r23554; } double f_of(float v) { float r23555 = 4; float r23556 = atan2(1.0, 0.0); float r23557 = 3; float r23558 = r23556 * r23557; float r23559 = r23555 / r23558; float r23560 = 1; float r23561 = v; float r23562 = r23561 * r23561; float r23563 = r23560 - r23562; float r23564 = r23559 / r23563; float r23565 = 2; float r23566 = 6; float r23567 = r23566 * r23561; float r23568 = r23561 * r23567; float r23569 = r23565 - r23568; float r23570 = r23564 / r23569; float r23571 = r23570 * r23564; float r23572 = sqrt(r23569); float r23573 = r23572 / r23564; float r23574 = r23571 / r23573; float r23575 = cbrt(r23574); return r23575; } double f_od(double v) { double r23576 = 4; double r23577 = atan2(1.0, 0.0); double r23578 = 3; double r23579 = r23577 * r23578; double r23580 = r23576 / r23579; double r23581 = 1; double r23582 = v; double r23583 = r23582 * r23582; double r23584 = r23581 - r23583; double r23585 = r23580 / r23584; double r23586 = 2; double r23587 = 6; double r23588 = r23587 * r23582; double r23589 = r23582 * r23588; double r23590 = r23586 - r23589; double r23591 = r23585 / r23590; double r23592 = r23591 * r23585; double r23593 = sqrt(r23590); double r23594 = r23593 / r23585; double r23595 = r23592 / r23594; double r23596 = cbrt(r23595); return r23596; } 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 r23597, r23598, r23599, r23600, r23601, r23602, r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23597, "4", 10, MPFR_RNDN); mpfr_init_set_str(r23598, "3", 10, MPFR_RNDN); mpfr_init(r23599); mpfr_init(r23600); mpfr_init_set_str(r23601, "1", 10, MPFR_RNDN); mpfr_init(r23602); mpfr_init(r23603); mpfr_init(r23604); mpfr_init(r23605); mpfr_init_set_str(r23606, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23607, "6", 10, MPFR_RNDN); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init(r23611); mpfr_init(r23612); } double f_im(double v) { ; ; mpfr_const_pi(r23599, MPFR_RNDN); mpfr_mul(r23600, r23598, r23599, MPFR_RNDN); ; mpfr_set_d(r23602, v, MPFR_RNDN); mpfr_mul(r23603, r23602, r23602, MPFR_RNDN); mpfr_sub(r23604, r23601, r23603, MPFR_RNDN); mpfr_mul(r23605, r23600, r23604, MPFR_RNDN); ; ; mpfr_mul(r23608, r23607, r23603, MPFR_RNDN); mpfr_sub(r23609, r23606, r23608, MPFR_RNDN); mpfr_sqrt(r23610, r23609, MPFR_RNDN); mpfr_mul(r23611, r23605, r23610, MPFR_RNDN); mpfr_div(r23612, r23597, r23611, MPFR_RNDN); return mpfr_get_d(r23612, MPFR_RNDN); } static mpfr_t r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23613, "4", 10, MPFR_RNDN); mpfr_init(r23614); mpfr_init_set_str(r23615, "3", 10, MPFR_RNDN); mpfr_init(r23616); mpfr_init(r23617); mpfr_init_set_str(r23618, "1", 10, MPFR_RNDN); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init_set_str(r23623, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23624, "6", 10, MPFR_RNDN); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); mpfr_init(r23628); mpfr_init(r23629); mpfr_init(r23630); mpfr_init(r23631); mpfr_init(r23632); mpfr_init(r23633); } double f_fm(double v) { ; mpfr_const_pi(r23614, MPFR_RNDN); ; mpfr_mul(r23616, r23614, r23615, MPFR_RNDN); mpfr_div(r23617, r23613, r23616, MPFR_RNDN); ; mpfr_set_d(r23619, v, MPFR_RNDN); mpfr_mul(r23620, r23619, r23619, MPFR_RNDN); mpfr_sub(r23621, r23618, r23620, MPFR_RNDN); mpfr_div(r23622, r23617, r23621, MPFR_RNDN); ; ; mpfr_mul(r23625, r23624, r23619, MPFR_RNDN); mpfr_mul(r23626, r23619, r23625, MPFR_RNDN); mpfr_sub(r23627, r23623, r23626, MPFR_RNDN); mpfr_div(r23628, r23622, r23627, MPFR_RNDN); mpfr_mul(r23629, r23628, r23622, MPFR_RNDN); mpfr_sqrt(r23630, r23627, MPFR_RNDN); mpfr_div(r23631, r23630, r23622, MPFR_RNDN); mpfr_div(r23632, r23629, r23631, MPFR_RNDN); mpfr_cbrt(r23633, r23632, MPFR_RNDN); return mpfr_get_d(r23633, MPFR_RNDN); } static mpfr_t r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647, r23648, r23649, r23650, r23651, r23652, r23653, r23654; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23634, "4", 10, MPFR_RNDN); mpfr_init(r23635); mpfr_init_set_str(r23636, "3", 10, MPFR_RNDN); mpfr_init(r23637); mpfr_init(r23638); mpfr_init_set_str(r23639, "1", 10, MPFR_RNDN); mpfr_init(r23640); mpfr_init(r23641); mpfr_init(r23642); mpfr_init(r23643); mpfr_init_set_str(r23644, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23645, "6", 10, MPFR_RNDN); mpfr_init(r23646); mpfr_init(r23647); mpfr_init(r23648); mpfr_init(r23649); mpfr_init(r23650); mpfr_init(r23651); mpfr_init(r23652); mpfr_init(r23653); mpfr_init(r23654); } double f_dm(double v) { ; mpfr_const_pi(r23635, MPFR_RNDN); ; mpfr_mul(r23637, r23635, r23636, MPFR_RNDN); mpfr_div(r23638, r23634, r23637, MPFR_RNDN); ; mpfr_set_d(r23640, v, MPFR_RNDN); mpfr_mul(r23641, r23640, r23640, MPFR_RNDN); mpfr_sub(r23642, r23639, r23641, MPFR_RNDN); mpfr_div(r23643, r23638, r23642, MPFR_RNDN); ; ; mpfr_mul(r23646, r23645, r23640, MPFR_RNDN); mpfr_mul(r23647, r23640, r23646, MPFR_RNDN); mpfr_sub(r23648, r23644, r23647, MPFR_RNDN); mpfr_div(r23649, r23643, r23648, MPFR_RNDN); mpfr_mul(r23650, r23649, r23643, MPFR_RNDN); mpfr_sqrt(r23651, r23648, MPFR_RNDN); mpfr_div(r23652, r23651, r23643, MPFR_RNDN); mpfr_div(r23653, r23650, r23652, MPFR_RNDN); mpfr_cbrt(r23654, r23653, MPFR_RNDN); return mpfr_get_d(r23654, MPFR_RNDN); }