#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25657 = J; float r25658 = l; float r25659 = exp(r25658); float r25660 = -r25658; float r25661 = exp(r25660); float r25662 = r25659 - r25661; float r25663 = r25657 * r25662; float r25664 = K; float r25665 = 2; float r25666 = r25664 / r25665; float r25667 = cos(r25666); float r25668 = r25663 * r25667; float r25669 = U; float r25670 = r25668 + r25669; return r25670; } double f_id(double J, double l, double K, double U) { double r25671 = J; double r25672 = l; double r25673 = exp(r25672); double r25674 = -r25672; double r25675 = exp(r25674); double r25676 = r25673 - r25675; double r25677 = r25671 * r25676; double r25678 = K; double r25679 = 2; double r25680 = r25678 / r25679; double r25681 = cos(r25680); double r25682 = r25677 * r25681; double r25683 = U; double r25684 = r25682 + r25683; return r25684; } double f_of(float J, float l, float K, float U) { float r25685 = K; float r25686 = 2; float r25687 = r25685 / r25686; float r25688 = cos(r25687); float r25689 = J; float r25690 = r25688 * r25689; float r25691 = l; float r25692 = 1/3; float r25693 = r25691 * r25692; float r25694 = fma(r25693, r25691, r25686); float r25695 = 5; float r25696 = pow(r25691, r25695); float r25697 = 1/60; float r25698 = r25696 * r25697; float r25699 = fma(r25691, r25694, r25698); float r25700 = U; float r25701 = fma(r25690, r25699, r25700); return r25701; } double f_od(double J, double l, double K, double U) { double r25702 = K; double r25703 = 2; double r25704 = r25702 / r25703; double r25705 = cos(r25704); double r25706 = J; double r25707 = r25705 * r25706; double r25708 = l; double r25709 = 1/3; double r25710 = r25708 * r25709; double r25711 = fma(r25710, r25708, r25703); double r25712 = 5; double r25713 = pow(r25708, r25712); double r25714 = 1/60; double r25715 = r25713 * r25714; double r25716 = fma(r25708, r25711, r25715); double r25717 = U; double r25718 = fma(r25707, r25716, r25717); return r25718; } 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 r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729, r25730, r25731, r25732; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init_set_str(r25727, "2", 10, MPFR_RNDN); mpfr_init(r25728); mpfr_init(r25729); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25719, J, MPFR_RNDN); mpfr_set_d(r25720, l, MPFR_RNDN); mpfr_exp(r25721, r25720, MPFR_RNDN); mpfr_neg(r25722, r25720, MPFR_RNDN); mpfr_exp(r25723, r25722, MPFR_RNDN); mpfr_sub(r25724, r25721, r25723, MPFR_RNDN); mpfr_mul(r25725, r25719, r25724, MPFR_RNDN); mpfr_set_d(r25726, K, MPFR_RNDN); ; mpfr_div(r25728, r25726, r25727, MPFR_RNDN); mpfr_cos(r25729, r25728, MPFR_RNDN); mpfr_mul(r25730, r25725, r25729, MPFR_RNDN); mpfr_set_d(r25731, U, MPFR_RNDN); mpfr_add(r25732, r25730, r25731, MPFR_RNDN); return mpfr_get_d(r25732, MPFR_RNDN); } static mpfr_t r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25733); mpfr_init_set_str(r25734, "2", 10, MPFR_RNDN); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); mpfr_init(r25738); mpfr_init(r25739); mpfr_init_set_str(r25740, "1/3", 10, MPFR_RNDN); mpfr_init(r25741); mpfr_init(r25742); mpfr_init_set_str(r25743, "5", 10, MPFR_RNDN); mpfr_init(r25744); mpfr_init_set_str(r25745, "1/60", 10, MPFR_RNDN); mpfr_init(r25746); mpfr_init(r25747); mpfr_init(r25748); mpfr_init(r25749); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25733, K, MPFR_RNDN); ; mpfr_div(r25735, r25733, r25734, MPFR_RNDN); mpfr_cos(r25736, r25735, MPFR_RNDN); mpfr_set_d(r25737, J, MPFR_RNDN); mpfr_mul(r25738, r25736, r25737, MPFR_RNDN); mpfr_set_d(r25739, l, MPFR_RNDN); ; mpfr_mul(r25741, r25739, r25740, MPFR_RNDN); mpfr_fma(r25742, r25741, r25739, r25734, MPFR_RNDN); ; mpfr_pow(r25744, r25739, r25743, MPFR_RNDN); ; mpfr_mul(r25746, r25744, r25745, MPFR_RNDN); mpfr_fma(r25747, r25739, r25742, r25746, MPFR_RNDN); mpfr_set_d(r25748, U, MPFR_RNDN); mpfr_fma(r25749, r25738, r25747, r25748, MPFR_RNDN); return mpfr_get_d(r25749, MPFR_RNDN); } static mpfr_t r25750, r25751, r25752, r25753, r25754, r25755, r25756, r25757, r25758, r25759, r25760, r25761, r25762, r25763, r25764, r25765, r25766; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25750); mpfr_init_set_str(r25751, "2", 10, MPFR_RNDN); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init(r25756); mpfr_init_set_str(r25757, "1/3", 10, MPFR_RNDN); mpfr_init(r25758); mpfr_init(r25759); mpfr_init_set_str(r25760, "5", 10, MPFR_RNDN); mpfr_init(r25761); mpfr_init_set_str(r25762, "1/60", 10, MPFR_RNDN); mpfr_init(r25763); mpfr_init(r25764); mpfr_init(r25765); mpfr_init(r25766); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r25750, K, MPFR_RNDN); ; mpfr_div(r25752, r25750, r25751, MPFR_RNDN); mpfr_cos(r25753, r25752, MPFR_RNDN); mpfr_set_d(r25754, J, MPFR_RNDN); mpfr_mul(r25755, r25753, r25754, MPFR_RNDN); mpfr_set_d(r25756, l, MPFR_RNDN); ; mpfr_mul(r25758, r25756, r25757, MPFR_RNDN); mpfr_fma(r25759, r25758, r25756, r25751, MPFR_RNDN); ; mpfr_pow(r25761, r25756, r25760, MPFR_RNDN); ; mpfr_mul(r25763, r25761, r25762, MPFR_RNDN); mpfr_fma(r25764, r25756, r25759, r25763, MPFR_RNDN); mpfr_set_d(r25765, U, MPFR_RNDN); mpfr_fma(r25766, r25755, r25764, r25765, MPFR_RNDN); return mpfr_get_d(r25766, MPFR_RNDN); }