#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r8685 = 2.0f; float r8686 = t; float r8687 = 3.0f; float r8688 = pow(r8686, r8687); float r8689 = l; float r8690 = r8689 * r8689; float r8691 = r8688 / r8690; float r8692 = k; float r8693 = sin(r8692); float r8694 = r8691 * r8693; float r8695 = tan(r8692); float r8696 = r8694 * r8695; float r8697 = 1.0f; float r8698 = r8692 / r8686; float r8699 = pow(r8698, r8685); float r8700 = r8697 + r8699; float r8701 = r8700 - r8697; float r8702 = r8696 * r8701; float r8703 = r8685 / r8702; return r8703; } double f_id(double t, double l, double k) { double r8704 = 2.0; double r8705 = t; double r8706 = 3.0; double r8707 = pow(r8705, r8706); double r8708 = l; double r8709 = r8708 * r8708; double r8710 = r8707 / r8709; double r8711 = k; double r8712 = sin(r8711); double r8713 = r8710 * r8712; double r8714 = tan(r8711); double r8715 = r8713 * r8714; double r8716 = 1.0; double r8717 = r8711 / r8705; double r8718 = pow(r8717, r8704); double r8719 = r8716 + r8718; double r8720 = r8719 - r8716; double r8721 = r8715 * r8720; double r8722 = r8704 / r8721; return r8722; } double f_of(float t, float l, float k) { float r8723 = l; float r8724 = t; float r8725 = r8723 / r8724; float r8726 = k; float r8727 = r8726 / r8724; float r8728 = fabs(r8727); float r8729 = r8725 / r8728; float r8730 = tan(r8726); float r8731 = r8725 / r8730; float r8732 = 2.0f; float r8733 = r8732 / r8724; float r8734 = sin(r8726); float r8735 = r8733 / r8734; float r8736 = r8735 / r8728; float r8737 = r8731 * r8736; float r8738 = r8729 * r8737; return r8738; } double f_od(double t, double l, double k) { double r8739 = l; double r8740 = t; double r8741 = r8739 / r8740; double r8742 = k; double r8743 = r8742 / r8740; double r8744 = fabs(r8743); double r8745 = r8741 / r8744; double r8746 = tan(r8742); double r8747 = r8741 / r8746; double r8748 = 2.0; double r8749 = r8748 / r8740; double r8750 = sin(r8742); double r8751 = r8749 / r8750; double r8752 = r8751 / r8744; double r8753 = r8747 * r8752; double r8754 = r8745 * r8753; return r8754; } 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 r8755, r8756, r8757, r8758, r8759, r8760, r8761, r8762, r8763, r8764, r8765, r8766, r8767, r8768, r8769, r8770, r8771, r8772, r8773; void setup_mpfr_f_im() { mpfr_set_default_prec(4176); mpfr_init_set_str(r8755, "2", 10, MPFR_RNDN); mpfr_init(r8756); mpfr_init_set_str(r8757, "3", 10, MPFR_RNDN); mpfr_init(r8758); mpfr_init(r8759); mpfr_init(r8760); mpfr_init(r8761); mpfr_init(r8762); mpfr_init(r8763); mpfr_init(r8764); mpfr_init(r8765); mpfr_init(r8766); mpfr_init_set_str(r8767, "1", 10, MPFR_RNDN); mpfr_init(r8768); mpfr_init(r8769); mpfr_init(r8770); mpfr_init(r8771); mpfr_init(r8772); mpfr_init(r8773); } double f_im(double t, double l, double k) { ; mpfr_set_d(r8756, t, MPFR_RNDN); ; mpfr_pow(r8758, r8756, r8757, MPFR_RNDN); mpfr_set_d(r8759, l, MPFR_RNDN); mpfr_mul(r8760, r8759, r8759, MPFR_RNDN); mpfr_div(r8761, r8758, r8760, MPFR_RNDN); mpfr_set_d(r8762, k, MPFR_RNDN); mpfr_sin(r8763, r8762, MPFR_RNDN); mpfr_mul(r8764, r8761, r8763, MPFR_RNDN); mpfr_tan(r8765, r8762, MPFR_RNDN); mpfr_mul(r8766, r8764, r8765, MPFR_RNDN); ; mpfr_div(r8768, r8762, r8756, MPFR_RNDN); mpfr_pow(r8769, r8768, r8755, MPFR_RNDN); mpfr_add(r8770, r8767, r8769, MPFR_RNDN); mpfr_sub(r8771, r8770, r8767, MPFR_RNDN); mpfr_mul(r8772, r8766, r8771, MPFR_RNDN); mpfr_div(r8773, r8755, r8772, MPFR_RNDN); return mpfr_get_d(r8773, MPFR_RNDN); } static mpfr_t r8774, r8775, r8776, r8777, r8778, r8779, r8780, r8781, r8782, r8783, r8784, r8785, r8786, r8787, r8788, r8789; void setup_mpfr_f_fm() { mpfr_set_default_prec(4176); mpfr_init(r8774); mpfr_init(r8775); mpfr_init(r8776); mpfr_init(r8777); mpfr_init(r8778); mpfr_init(r8779); mpfr_init(r8780); mpfr_init(r8781); mpfr_init(r8782); mpfr_init_set_str(r8783, "2", 10, MPFR_RNDN); mpfr_init(r8784); mpfr_init(r8785); mpfr_init(r8786); mpfr_init(r8787); mpfr_init(r8788); mpfr_init(r8789); } double f_fm(double t, double l, double k) { mpfr_set_d(r8774, l, MPFR_RNDN); mpfr_set_d(r8775, t, MPFR_RNDN); mpfr_div(r8776, r8774, r8775, MPFR_RNDN); mpfr_set_d(r8777, k, MPFR_RNDN); mpfr_div(r8778, r8777, r8775, MPFR_RNDN); mpfr_abs(r8779, r8778, MPFR_RNDN); mpfr_div(r8780, r8776, r8779, MPFR_RNDN); mpfr_tan(r8781, r8777, MPFR_RNDN); mpfr_div(r8782, r8776, r8781, MPFR_RNDN); ; mpfr_div(r8784, r8783, r8775, MPFR_RNDN); mpfr_sin(r8785, r8777, MPFR_RNDN); mpfr_div(r8786, r8784, r8785, MPFR_RNDN); mpfr_div(r8787, r8786, r8779, MPFR_RNDN); mpfr_mul(r8788, r8782, r8787, MPFR_RNDN); mpfr_mul(r8789, r8780, r8788, MPFR_RNDN); return mpfr_get_d(r8789, MPFR_RNDN); } static mpfr_t r8790, r8791, r8792, r8793, r8794, r8795, r8796, r8797, r8798, r8799, r8800, r8801, r8802, r8803, r8804, r8805; void setup_mpfr_f_dm() { mpfr_set_default_prec(4176); mpfr_init(r8790); mpfr_init(r8791); mpfr_init(r8792); mpfr_init(r8793); mpfr_init(r8794); mpfr_init(r8795); mpfr_init(r8796); mpfr_init(r8797); mpfr_init(r8798); mpfr_init_set_str(r8799, "2", 10, MPFR_RNDN); mpfr_init(r8800); mpfr_init(r8801); mpfr_init(r8802); mpfr_init(r8803); mpfr_init(r8804); mpfr_init(r8805); } double f_dm(double t, double l, double k) { mpfr_set_d(r8790, l, MPFR_RNDN); mpfr_set_d(r8791, t, MPFR_RNDN); mpfr_div(r8792, r8790, r8791, MPFR_RNDN); mpfr_set_d(r8793, k, MPFR_RNDN); mpfr_div(r8794, r8793, r8791, MPFR_RNDN); mpfr_abs(r8795, r8794, MPFR_RNDN); mpfr_div(r8796, r8792, r8795, MPFR_RNDN); mpfr_tan(r8797, r8793, MPFR_RNDN); mpfr_div(r8798, r8792, r8797, MPFR_RNDN); ; mpfr_div(r8800, r8799, r8791, MPFR_RNDN); mpfr_sin(r8801, r8793, MPFR_RNDN); mpfr_div(r8802, r8800, r8801, MPFR_RNDN); mpfr_div(r8803, r8802, r8795, MPFR_RNDN); mpfr_mul(r8804, r8798, r8803, MPFR_RNDN); mpfr_mul(r8805, r8796, r8804, MPFR_RNDN); return mpfr_get_d(r8805, MPFR_RNDN); }