#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (23)"; double f_if(float F, float B, float x) { float r23648 = x; float r23649 = 1; float r23650 = B; float r23651 = tan(r23650); float r23652 = r23649 / r23651; float r23653 = r23648 * r23652; float r23654 = -r23653; float r23655 = F; float r23656 = sin(r23650); float r23657 = r23655 / r23656; float r23658 = r23655 * r23655; float r23659 = 2; float r23660 = r23658 + r23659; float r23661 = r23659 * r23648; float r23662 = r23660 + r23661; float r23663 = r23649 / r23659; float r23664 = -r23663; float r23665 = pow(r23662, r23664); float r23666 = r23657 * r23665; float r23667 = r23654 + r23666; return r23667; } double f_id(double F, double B, double x) { double r23668 = x; double r23669 = 1; double r23670 = B; double r23671 = tan(r23670); double r23672 = r23669 / r23671; double r23673 = r23668 * r23672; double r23674 = -r23673; double r23675 = F; double r23676 = sin(r23670); double r23677 = r23675 / r23676; double r23678 = r23675 * r23675; double r23679 = 2; double r23680 = r23678 + r23679; double r23681 = r23679 * r23668; double r23682 = r23680 + r23681; double r23683 = r23669 / r23679; double r23684 = -r23683; double r23685 = pow(r23682, r23684); double r23686 = r23677 * r23685; double r23687 = r23674 + r23686; return r23687; } double f_of(float F, float B, float x) { float r23688 = x; float r23689 = -r23688; float r23690 = B; float r23691 = tan(r23690); float r23692 = r23689 / r23691; float r23693 = r23688 + r23688; float r23694 = 2; float r23695 = F; float r23696 = r23695 * r23695; float r23697 = r23694 + r23696; float r23698 = r23693 + r23697; float r23699 = 1; float r23700 = r23699 / r23694; float r23701 = -r23700; float r23702 = pow(r23698, r23701); float r23703 = r23702 * r23695; float r23704 = sin(r23690); float r23705 = r23699 / r23704; float r23706 = r23703 * r23705; float r23707 = r23692 + r23706; return r23707; } double f_od(double F, double B, double x) { double r23708 = x; double r23709 = -r23708; double r23710 = B; double r23711 = tan(r23710); double r23712 = r23709 / r23711; double r23713 = r23708 + r23708; double r23714 = 2; double r23715 = F; double r23716 = r23715 * r23715; double r23717 = r23714 + r23716; double r23718 = r23713 + r23717; double r23719 = 1; double r23720 = r23719 / r23714; double r23721 = -r23720; double r23722 = pow(r23718, r23721); double r23723 = r23722 * r23715; double r23724 = sin(r23710); double r23725 = r23719 / r23724; double r23726 = r23723 * r23725; double r23727 = r23712 + r23726; return r23727; } 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 r23728, r23729, r23730, r23731, r23732, r23733, r23734, r23735, r23736, r23737, r23738, r23739, r23740, r23741, r23742, r23743, r23744, r23745, r23746, r23747; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23728); mpfr_init_set_str(r23729, "1", 10, MPFR_RNDN); mpfr_init(r23730); mpfr_init(r23731); mpfr_init(r23732); mpfr_init(r23733); mpfr_init(r23734); mpfr_init(r23735); mpfr_init(r23736); mpfr_init(r23737); mpfr_init(r23738); mpfr_init_set_str(r23739, "2", 10, MPFR_RNDN); mpfr_init(r23740); mpfr_init(r23741); mpfr_init(r23742); mpfr_init(r23743); mpfr_init(r23744); mpfr_init(r23745); mpfr_init(r23746); mpfr_init(r23747); } double f_im(double F, double B, double x) { mpfr_set_d(r23728, x, MPFR_RNDN); ; mpfr_set_d(r23730, B, MPFR_RNDN); mpfr_tan(r23731, r23730, MPFR_RNDN); mpfr_div(r23732, r23729, r23731, MPFR_RNDN); mpfr_mul(r23733, r23728, r23732, MPFR_RNDN); mpfr_neg(r23734, r23733, MPFR_RNDN); mpfr_set_d(r23735, F, MPFR_RNDN); mpfr_sin(r23736, r23730, MPFR_RNDN); mpfr_div(r23737, r23735, r23736, MPFR_RNDN); mpfr_mul(r23738, r23735, r23735, MPFR_RNDN); ; mpfr_add(r23740, r23738, r23739, MPFR_RNDN); mpfr_mul(r23741, r23739, r23728, MPFR_RNDN); mpfr_add(r23742, r23740, r23741, MPFR_RNDN); mpfr_div(r23743, r23729, r23739, MPFR_RNDN); mpfr_neg(r23744, r23743, MPFR_RNDN); mpfr_pow(r23745, r23742, r23744, MPFR_RNDN); mpfr_mul(r23746, r23737, r23745, MPFR_RNDN); mpfr_add(r23747, r23734, r23746, MPFR_RNDN); return mpfr_get_d(r23747, MPFR_RNDN); } static mpfr_t r23748, r23749, r23750, r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23748); mpfr_init(r23749); mpfr_init(r23750); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init_set_str(r23754, "2", 10, MPFR_RNDN); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init_set_str(r23759, "1", 10, MPFR_RNDN); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); } double f_fm(double F, double B, double x) { mpfr_set_d(r23748, x, MPFR_RNDN); mpfr_neg(r23749, r23748, MPFR_RNDN); mpfr_set_d(r23750, B, MPFR_RNDN); mpfr_tan(r23751, r23750, MPFR_RNDN); mpfr_div(r23752, r23749, r23751, MPFR_RNDN); mpfr_add(r23753, r23748, r23748, MPFR_RNDN); ; mpfr_set_d(r23755, F, MPFR_RNDN); mpfr_mul(r23756, r23755, r23755, MPFR_RNDN); mpfr_add(r23757, r23754, r23756, MPFR_RNDN); mpfr_add(r23758, r23753, r23757, MPFR_RNDN); ; mpfr_div(r23760, r23759, r23754, MPFR_RNDN); mpfr_neg(r23761, r23760, MPFR_RNDN); mpfr_pow(r23762, r23758, r23761, MPFR_RNDN); mpfr_mul(r23763, r23762, r23755, MPFR_RNDN); mpfr_sin(r23764, r23750, MPFR_RNDN); mpfr_div(r23765, r23759, r23764, MPFR_RNDN); mpfr_mul(r23766, r23763, r23765, MPFR_RNDN); mpfr_add(r23767, r23752, r23766, MPFR_RNDN); return mpfr_get_d(r23767, MPFR_RNDN); } static mpfr_t r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778, r23779, r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23768); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init_set_str(r23774, "2", 10, MPFR_RNDN); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); mpfr_init_set_str(r23779, "1", 10, MPFR_RNDN); mpfr_init(r23780); mpfr_init(r23781); mpfr_init(r23782); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); } double f_dm(double F, double B, double x) { mpfr_set_d(r23768, x, MPFR_RNDN); mpfr_neg(r23769, r23768, MPFR_RNDN); mpfr_set_d(r23770, B, MPFR_RNDN); mpfr_tan(r23771, r23770, MPFR_RNDN); mpfr_div(r23772, r23769, r23771, MPFR_RNDN); mpfr_add(r23773, r23768, r23768, MPFR_RNDN); ; mpfr_set_d(r23775, F, MPFR_RNDN); mpfr_mul(r23776, r23775, r23775, MPFR_RNDN); mpfr_add(r23777, r23774, r23776, MPFR_RNDN); mpfr_add(r23778, r23773, r23777, MPFR_RNDN); ; mpfr_div(r23780, r23779, r23774, MPFR_RNDN); mpfr_neg(r23781, r23780, MPFR_RNDN); mpfr_pow(r23782, r23778, r23781, MPFR_RNDN); mpfr_mul(r23783, r23782, r23775, MPFR_RNDN); mpfr_sin(r23784, r23770, MPFR_RNDN); mpfr_div(r23785, r23779, r23784, MPFR_RNDN); mpfr_mul(r23786, r23783, r23785, MPFR_RNDN); mpfr_add(r23787, r23772, r23786, MPFR_RNDN); return mpfr_get_d(r23787, MPFR_RNDN); }