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