#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25690 = J; float r25691 = l; float r25692 = exp(r25691); float r25693 = -r25691; float r25694 = exp(r25693); float r25695 = r25692 - r25694; float r25696 = r25690 * r25695; float r25697 = K; float r25698 = 2; float r25699 = r25697 / r25698; float r25700 = cos(r25699); float r25701 = r25696 * r25700; float r25702 = U; float r25703 = r25701 + r25702; return r25703; } double f_id(double J, double l, double K, double U) { double r25704 = J; double r25705 = l; double r25706 = exp(r25705); double r25707 = -r25705; double r25708 = exp(r25707); double r25709 = r25706 - r25708; double r25710 = r25704 * r25709; double r25711 = K; double r25712 = 2; double r25713 = r25711 / r25712; double r25714 = cos(r25713); double r25715 = r25710 * r25714; double r25716 = U; double r25717 = r25715 + r25716; return r25717; } double f_of(float J, float l, float K, float U) { float r25718 = K; float r25719 = 2; float r25720 = r25718 / r25719; float r25721 = cos(r25720); float r25722 = J; float r25723 = r25721 * r25722; float r25724 = l; float r25725 = 1/3; float r25726 = r25724 * r25725; float r25727 = fma(r25726, r25724, r25719); float r25728 = 5; float r25729 = pow(r25724, r25728); float r25730 = 1/60; float r25731 = r25729 * r25730; float r25732 = fma(r25724, r25727, r25731); float r25733 = U; float r25734 = fma(r25723, r25732, r25733); return r25734; } double f_od(double J, double l, double K, double U) { double r25735 = K; double r25736 = 2; double r25737 = r25735 / r25736; double r25738 = cos(r25737); double r25739 = J; double r25740 = r25738 * r25739; double r25741 = l; double r25742 = 1/3; double r25743 = r25741 * r25742; double r25744 = fma(r25743, r25741, r25736); double r25745 = 5; double r25746 = pow(r25741, r25745); double r25747 = 1/60; double r25748 = r25746 * r25747; double r25749 = fma(r25741, r25744, r25748); double r25750 = U; double r25751 = fma(r25740, r25749, r25750); return r25751; } 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 r25752, r25753, r25754, r25755, r25756, r25757, r25758, r25759, r25760, r25761, r25762, r25763, r25764, r25765; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init(r25756); mpfr_init(r25757); mpfr_init(r25758); mpfr_init(r25759); mpfr_init_set_str(r25760, "2", 10, MPFR_RNDN); mpfr_init(r25761); mpfr_init(r25762); mpfr_init(r25763); mpfr_init(r25764); mpfr_init(r25765); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25752, J, MPFR_RNDN); mpfr_set_d(r25753, l, MPFR_RNDN); mpfr_exp(r25754, r25753, MPFR_RNDN); mpfr_neg(r25755, r25753, MPFR_RNDN); mpfr_exp(r25756, r25755, MPFR_RNDN); mpfr_sub(r25757, r25754, r25756, MPFR_RNDN); mpfr_mul(r25758, r25752, r25757, MPFR_RNDN); mpfr_set_d(r25759, K, MPFR_RNDN); ; mpfr_div(r25761, r25759, r25760, MPFR_RNDN); mpfr_cos(r25762, r25761, MPFR_RNDN); mpfr_mul(r25763, r25758, r25762, MPFR_RNDN); mpfr_set_d(r25764, U, MPFR_RNDN); mpfr_add(r25765, r25763, r25764, MPFR_RNDN); return mpfr_get_d(r25765, MPFR_RNDN); } static mpfr_t r25766, r25767, r25768, r25769, r25770, r25771, r25772, r25773, r25774, r25775, r25776, r25777, r25778, r25779, r25780, r25781, r25782; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25766); mpfr_init_set_str(r25767, "2", 10, MPFR_RNDN); mpfr_init(r25768); mpfr_init(r25769); mpfr_init(r25770); mpfr_init(r25771); mpfr_init(r25772); mpfr_init_set_str(r25773, "1/3", 10, MPFR_RNDN); mpfr_init(r25774); mpfr_init(r25775); mpfr_init_set_str(r25776, "5", 10, MPFR_RNDN); mpfr_init(r25777); mpfr_init_set_str(r25778, "1/60", 10, MPFR_RNDN); mpfr_init(r25779); mpfr_init(r25780); mpfr_init(r25781); mpfr_init(r25782); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25766, K, MPFR_RNDN); ; mpfr_div(r25768, r25766, r25767, MPFR_RNDN); mpfr_cos(r25769, r25768, MPFR_RNDN); mpfr_set_d(r25770, J, MPFR_RNDN); mpfr_mul(r25771, r25769, r25770, MPFR_RNDN); mpfr_set_d(r25772, l, MPFR_RNDN); ; mpfr_mul(r25774, r25772, r25773, MPFR_RNDN); mpfr_fma(r25775, r25774, r25772, r25767, MPFR_RNDN); ; mpfr_pow(r25777, r25772, r25776, MPFR_RNDN); ; mpfr_mul(r25779, r25777, r25778, MPFR_RNDN); mpfr_fma(r25780, r25772, r25775, r25779, MPFR_RNDN); mpfr_set_d(r25781, U, MPFR_RNDN); mpfr_fma(r25782, r25771, r25780, r25781, MPFR_RNDN); return mpfr_get_d(r25782, MPFR_RNDN); } static mpfr_t r25783, r25784, r25785, r25786, r25787, r25788, r25789, r25790, r25791, r25792, r25793, r25794, r25795, r25796, r25797, r25798, r25799; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25783); mpfr_init_set_str(r25784, "2", 10, MPFR_RNDN); mpfr_init(r25785); mpfr_init(r25786); mpfr_init(r25787); mpfr_init(r25788); mpfr_init(r25789); mpfr_init_set_str(r25790, "1/3", 10, MPFR_RNDN); mpfr_init(r25791); mpfr_init(r25792); mpfr_init_set_str(r25793, "5", 10, MPFR_RNDN); mpfr_init(r25794); mpfr_init_set_str(r25795, "1/60", 10, MPFR_RNDN); mpfr_init(r25796); mpfr_init(r25797); mpfr_init(r25798); mpfr_init(r25799); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r25783, K, MPFR_RNDN); ; mpfr_div(r25785, r25783, r25784, MPFR_RNDN); mpfr_cos(r25786, r25785, MPFR_RNDN); mpfr_set_d(r25787, J, MPFR_RNDN); mpfr_mul(r25788, r25786, r25787, MPFR_RNDN); mpfr_set_d(r25789, l, MPFR_RNDN); ; mpfr_mul(r25791, r25789, r25790, MPFR_RNDN); mpfr_fma(r25792, r25791, r25789, r25784, MPFR_RNDN); ; mpfr_pow(r25794, r25789, r25793, MPFR_RNDN); ; mpfr_mul(r25796, r25794, r25795, MPFR_RNDN); mpfr_fma(r25797, r25789, r25792, r25796, MPFR_RNDN); mpfr_set_d(r25798, U, MPFR_RNDN); mpfr_fma(r25799, r25788, r25797, r25798, MPFR_RNDN); return mpfr_get_d(r25799, MPFR_RNDN); }