#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25589 = J; float r25590 = l; float r25591 = exp(r25590); float r25592 = -r25590; float r25593 = exp(r25592); float r25594 = r25591 - r25593; float r25595 = r25589 * r25594; float r25596 = K; float r25597 = 2; float r25598 = r25596 / r25597; float r25599 = cos(r25598); float r25600 = r25595 * r25599; float r25601 = U; float r25602 = r25600 + r25601; return r25602; } double f_id(double J, double l, double K, double U) { double r25603 = J; double r25604 = l; double r25605 = exp(r25604); double r25606 = -r25604; double r25607 = exp(r25606); double r25608 = r25605 - r25607; double r25609 = r25603 * r25608; double r25610 = K; double r25611 = 2; double r25612 = r25610 / r25611; double r25613 = cos(r25612); double r25614 = r25609 * r25613; double r25615 = U; double r25616 = r25614 + r25615; return r25616; } double f_of(float J, float l, float K, float U) { float r25617 = K; float r25618 = 2; float r25619 = r25617 / r25618; float r25620 = cos(r25619); float r25621 = J; float r25622 = r25620 * r25621; float r25623 = l; float r25624 = 1/3; float r25625 = r25623 * r25624; float r25626 = fma(r25625, r25623, r25618); float r25627 = 5; float r25628 = pow(r25623, r25627); float r25629 = 1/60; float r25630 = r25628 * r25629; float r25631 = fma(r25623, r25626, r25630); float r25632 = U; float r25633 = fma(r25622, r25631, r25632); return r25633; } double f_od(double J, double l, double K, double U) { double r25634 = K; double r25635 = 2; double r25636 = r25634 / r25635; double r25637 = cos(r25636); double r25638 = J; double r25639 = r25637 * r25638; double r25640 = l; double r25641 = 1/3; double r25642 = r25640 * r25641; double r25643 = fma(r25642, r25640, r25635); double r25644 = 5; double r25645 = pow(r25640, r25644); double r25646 = 1/60; double r25647 = r25645 * r25646; double r25648 = fma(r25640, r25643, r25647); double r25649 = U; double r25650 = fma(r25639, r25648, r25649); return r25650; } 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 r25651, r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659, r25660, r25661, r25662, r25663, r25664; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25651); mpfr_init(r25652); mpfr_init(r25653); mpfr_init(r25654); mpfr_init(r25655); mpfr_init(r25656); mpfr_init(r25657); mpfr_init(r25658); mpfr_init_set_str(r25659, "2", 10, MPFR_RNDN); mpfr_init(r25660); mpfr_init(r25661); mpfr_init(r25662); mpfr_init(r25663); mpfr_init(r25664); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25651, J, MPFR_RNDN); mpfr_set_d(r25652, l, MPFR_RNDN); mpfr_exp(r25653, r25652, MPFR_RNDN); mpfr_neg(r25654, r25652, MPFR_RNDN); mpfr_exp(r25655, r25654, MPFR_RNDN); mpfr_sub(r25656, r25653, r25655, MPFR_RNDN); mpfr_mul(r25657, r25651, r25656, MPFR_RNDN); mpfr_set_d(r25658, K, MPFR_RNDN); ; mpfr_div(r25660, r25658, r25659, MPFR_RNDN); mpfr_cos(r25661, r25660, MPFR_RNDN); mpfr_mul(r25662, r25657, r25661, MPFR_RNDN); mpfr_set_d(r25663, U, MPFR_RNDN); mpfr_add(r25664, r25662, r25663, MPFR_RNDN); return mpfr_get_d(r25664, MPFR_RNDN); } static mpfr_t r25665, r25666, r25667, r25668, r25669, r25670, r25671, r25672, r25673, r25674, r25675, r25676, r25677, r25678, r25679, r25680, r25681; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25665); mpfr_init_set_str(r25666, "2", 10, MPFR_RNDN); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); mpfr_init(r25670); mpfr_init(r25671); mpfr_init_set_str(r25672, "1/3", 10, MPFR_RNDN); mpfr_init(r25673); mpfr_init(r25674); mpfr_init_set_str(r25675, "5", 10, MPFR_RNDN); mpfr_init(r25676); mpfr_init_set_str(r25677, "1/60", 10, MPFR_RNDN); mpfr_init(r25678); mpfr_init(r25679); mpfr_init(r25680); mpfr_init(r25681); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25665, K, MPFR_RNDN); ; mpfr_div(r25667, r25665, r25666, MPFR_RNDN); mpfr_cos(r25668, r25667, MPFR_RNDN); mpfr_set_d(r25669, J, MPFR_RNDN); mpfr_mul(r25670, r25668, r25669, MPFR_RNDN); mpfr_set_d(r25671, l, MPFR_RNDN); ; mpfr_mul(r25673, r25671, r25672, MPFR_RNDN); mpfr_fma(r25674, r25673, r25671, r25666, MPFR_RNDN); ; mpfr_pow(r25676, r25671, r25675, MPFR_RNDN); ; mpfr_mul(r25678, r25676, r25677, MPFR_RNDN); mpfr_fma(r25679, r25671, r25674, r25678, MPFR_RNDN); mpfr_set_d(r25680, U, MPFR_RNDN); mpfr_fma(r25681, r25670, r25679, r25680, MPFR_RNDN); return mpfr_get_d(r25681, MPFR_RNDN); } static mpfr_t r25682, r25683, r25684, r25685, r25686, r25687, r25688, r25689, r25690, r25691, r25692, r25693, r25694, r25695, r25696, r25697, r25698; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25682); mpfr_init_set_str(r25683, "2", 10, MPFR_RNDN); mpfr_init(r25684); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init_set_str(r25689, "1/3", 10, MPFR_RNDN); mpfr_init(r25690); mpfr_init(r25691); mpfr_init_set_str(r25692, "5", 10, MPFR_RNDN); mpfr_init(r25693); mpfr_init_set_str(r25694, "1/60", 10, MPFR_RNDN); mpfr_init(r25695); mpfr_init(r25696); mpfr_init(r25697); mpfr_init(r25698); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r25682, K, MPFR_RNDN); ; mpfr_div(r25684, r25682, r25683, MPFR_RNDN); mpfr_cos(r25685, r25684, MPFR_RNDN); mpfr_set_d(r25686, J, MPFR_RNDN); mpfr_mul(r25687, r25685, r25686, MPFR_RNDN); mpfr_set_d(r25688, l, MPFR_RNDN); ; mpfr_mul(r25690, r25688, r25689, MPFR_RNDN); mpfr_fma(r25691, r25690, r25688, r25683, MPFR_RNDN); ; mpfr_pow(r25693, r25688, r25692, MPFR_RNDN); ; mpfr_mul(r25695, r25693, r25694, MPFR_RNDN); mpfr_fma(r25696, r25688, r25691, r25695, MPFR_RNDN); mpfr_set_d(r25697, U, MPFR_RNDN); mpfr_fma(r25698, r25687, r25696, r25697, MPFR_RNDN); return mpfr_get_d(r25698, MPFR_RNDN); }