#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r23504 = x; float r23505 = 1; float r23506 = B; float r23507 = tan(r23506); float r23508 = r23505 / r23507; float r23509 = r23504 * r23508; float r23510 = -r23509; float r23511 = F; float r23512 = sin(r23506); float r23513 = r23511 / r23512; float r23514 = r23511 * r23511; float r23515 = 2; float r23516 = r23514 + r23515; float r23517 = r23515 * r23504; float r23518 = r23516 + r23517; float r23519 = r23505 / r23515; float r23520 = -r23519; float r23521 = pow(r23518, r23520); float r23522 = r23513 * r23521; float r23523 = r23510 + r23522; return r23523; } double f_id(double F, double B, double x) { double r23524 = x; double r23525 = 1; double r23526 = B; double r23527 = tan(r23526); double r23528 = r23525 / r23527; double r23529 = r23524 * r23528; double r23530 = -r23529; double r23531 = F; double r23532 = sin(r23526); double r23533 = r23531 / r23532; double r23534 = r23531 * r23531; double r23535 = 2; double r23536 = r23534 + r23535; double r23537 = r23535 * r23524; double r23538 = r23536 + r23537; double r23539 = r23525 / r23535; double r23540 = -r23539; double r23541 = pow(r23538, r23540); double r23542 = r23533 * r23541; double r23543 = r23530 + r23542; return r23543; } double f_of(float F, float B, float x) { float r23544 = x; float r23545 = -r23544; float r23546 = B; float r23547 = tan(r23546); float r23548 = r23545 / r23547; float r23549 = r23544 + r23544; float r23550 = 2; float r23551 = F; float r23552 = r23551 * r23551; float r23553 = r23550 + r23552; float r23554 = r23549 + r23553; float r23555 = 1; float r23556 = r23555 / r23550; float r23557 = -r23556; float r23558 = pow(r23554, r23557); float r23559 = r23558 * r23551; float r23560 = sin(r23546); float r23561 = r23555 / r23560; float r23562 = r23559 * r23561; float r23563 = r23548 + r23562; return r23563; } double f_od(double F, double B, double x) { double r23564 = x; double r23565 = -r23564; double r23566 = B; double r23567 = tan(r23566); double r23568 = r23565 / r23567; double r23569 = r23564 + r23564; double r23570 = 2; double r23571 = F; double r23572 = r23571 * r23571; double r23573 = r23570 + r23572; double r23574 = r23569 + r23573; double r23575 = 1; double r23576 = r23575 / r23570; double r23577 = -r23576; double r23578 = pow(r23574, r23577); double r23579 = r23578 * r23571; double r23580 = sin(r23566); double r23581 = r23575 / r23580; double r23582 = r23579 * r23581; double r23583 = r23568 + r23582; return r23583; } 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 r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593, r23594, r23595, r23596, r23597, r23598, r23599, r23600, r23601, r23602, r23603; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23584); mpfr_init_set_str(r23585, "1", 10, MPFR_RNDN); mpfr_init(r23586); mpfr_init(r23587); mpfr_init(r23588); mpfr_init(r23589); mpfr_init(r23590); mpfr_init(r23591); mpfr_init(r23592); mpfr_init(r23593); mpfr_init(r23594); mpfr_init_set_str(r23595, "2", 10, MPFR_RNDN); mpfr_init(r23596); mpfr_init(r23597); mpfr_init(r23598); mpfr_init(r23599); mpfr_init(r23600); mpfr_init(r23601); mpfr_init(r23602); mpfr_init(r23603); } double f_im(double F, double B, double x) { mpfr_set_d(r23584, x, MPFR_RNDN); ; mpfr_set_d(r23586, B, MPFR_RNDN); mpfr_tan(r23587, r23586, MPFR_RNDN); mpfr_div(r23588, r23585, r23587, MPFR_RNDN); mpfr_mul(r23589, r23584, r23588, MPFR_RNDN); mpfr_neg(r23590, r23589, MPFR_RNDN); mpfr_set_d(r23591, F, MPFR_RNDN); mpfr_sin(r23592, r23586, MPFR_RNDN); mpfr_div(r23593, r23591, r23592, MPFR_RNDN); mpfr_mul(r23594, r23591, r23591, MPFR_RNDN); ; mpfr_add(r23596, r23594, r23595, MPFR_RNDN); mpfr_mul(r23597, r23595, r23584, MPFR_RNDN); mpfr_add(r23598, r23596, r23597, MPFR_RNDN); mpfr_div(r23599, r23585, r23595, MPFR_RNDN); mpfr_neg(r23600, r23599, MPFR_RNDN); mpfr_pow(r23601, r23598, r23600, MPFR_RNDN); mpfr_mul(r23602, r23593, r23601, MPFR_RNDN); mpfr_add(r23603, r23590, r23602, MPFR_RNDN); return mpfr_get_d(r23603, MPFR_RNDN); } static mpfr_t r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23604); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init_set_str(r23610, "2", 10, MPFR_RNDN); mpfr_init(r23611); mpfr_init(r23612); mpfr_init(r23613); mpfr_init(r23614); mpfr_init_set_str(r23615, "1", 10, MPFR_RNDN); mpfr_init(r23616); mpfr_init(r23617); mpfr_init(r23618); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); } double f_fm(double F, double B, double x) { mpfr_set_d(r23604, x, MPFR_RNDN); mpfr_neg(r23605, r23604, MPFR_RNDN); mpfr_set_d(r23606, B, MPFR_RNDN); mpfr_tan(r23607, r23606, MPFR_RNDN); mpfr_div(r23608, r23605, r23607, MPFR_RNDN); mpfr_add(r23609, r23604, r23604, MPFR_RNDN); ; mpfr_set_d(r23611, F, MPFR_RNDN); mpfr_mul(r23612, r23611, r23611, MPFR_RNDN); mpfr_add(r23613, r23610, r23612, MPFR_RNDN); mpfr_add(r23614, r23609, r23613, MPFR_RNDN); ; mpfr_div(r23616, r23615, r23610, MPFR_RNDN); mpfr_neg(r23617, r23616, MPFR_RNDN); mpfr_pow(r23618, r23614, r23617, MPFR_RNDN); mpfr_mul(r23619, r23618, r23611, MPFR_RNDN); mpfr_sin(r23620, r23606, MPFR_RNDN); mpfr_div(r23621, r23615, r23620, MPFR_RNDN); mpfr_mul(r23622, r23619, r23621, MPFR_RNDN); mpfr_add(r23623, r23608, r23622, MPFR_RNDN); return mpfr_get_d(r23623, MPFR_RNDN); } static mpfr_t r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23624); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); mpfr_init(r23628); mpfr_init(r23629); mpfr_init_set_str(r23630, "2", 10, MPFR_RNDN); mpfr_init(r23631); mpfr_init(r23632); mpfr_init(r23633); mpfr_init(r23634); mpfr_init_set_str(r23635, "1", 10, MPFR_RNDN); mpfr_init(r23636); mpfr_init(r23637); mpfr_init(r23638); mpfr_init(r23639); mpfr_init(r23640); mpfr_init(r23641); mpfr_init(r23642); mpfr_init(r23643); } double f_dm(double F, double B, double x) { mpfr_set_d(r23624, x, MPFR_RNDN); mpfr_neg(r23625, r23624, MPFR_RNDN); mpfr_set_d(r23626, B, MPFR_RNDN); mpfr_tan(r23627, r23626, MPFR_RNDN); mpfr_div(r23628, r23625, r23627, MPFR_RNDN); mpfr_add(r23629, r23624, r23624, MPFR_RNDN); ; mpfr_set_d(r23631, F, MPFR_RNDN); mpfr_mul(r23632, r23631, r23631, MPFR_RNDN); mpfr_add(r23633, r23630, r23632, MPFR_RNDN); mpfr_add(r23634, r23629, r23633, MPFR_RNDN); ; mpfr_div(r23636, r23635, r23630, MPFR_RNDN); mpfr_neg(r23637, r23636, MPFR_RNDN); mpfr_pow(r23638, r23634, r23637, MPFR_RNDN); mpfr_mul(r23639, r23638, r23631, MPFR_RNDN); mpfr_sin(r23640, r23626, MPFR_RNDN); mpfr_div(r23641, r23635, r23640, MPFR_RNDN); mpfr_mul(r23642, r23639, r23641, MPFR_RNDN); mpfr_add(r23643, r23628, r23642, MPFR_RNDN); return mpfr_get_d(r23643, MPFR_RNDN); }