#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r25614 = J; float r25615 = l; float r25616 = exp(r25615); float r25617 = -r25615; float r25618 = exp(r25617); float r25619 = r25616 - r25618; float r25620 = r25614 * r25619; float r25621 = K; float r25622 = 2; float r25623 = r25621 / r25622; float r25624 = cos(r25623); float r25625 = r25620 * r25624; float r25626 = U; float r25627 = r25625 + r25626; return r25627; } double f_id(double J, double l, double K, double U) { double r25628 = J; double r25629 = l; double r25630 = exp(r25629); double r25631 = -r25629; double r25632 = exp(r25631); double r25633 = r25630 - r25632; double r25634 = r25628 * r25633; double r25635 = K; double r25636 = 2; double r25637 = r25635 / r25636; double r25638 = cos(r25637); double r25639 = r25634 * r25638; double r25640 = U; double r25641 = r25639 + r25640; return r25641; } double f_of(float J, float l, float K, float U) { float r25642 = K; float r25643 = 2; float r25644 = r25642 / r25643; float r25645 = cos(r25644); float r25646 = J; float r25647 = r25645 * r25646; float r25648 = l; float r25649 = 1/3; float r25650 = r25648 * r25649; float r25651 = fma(r25650, r25648, r25643); float r25652 = 5; float r25653 = pow(r25648, r25652); float r25654 = 1/60; float r25655 = r25653 * r25654; float r25656 = fma(r25648, r25651, r25655); float r25657 = U; float r25658 = fma(r25647, r25656, r25657); return r25658; } double f_od(double J, double l, double K, double U) { double r25659 = K; double r25660 = 2; double r25661 = r25659 / r25660; double r25662 = cos(r25661); double r25663 = J; double r25664 = r25662 * r25663; double r25665 = l; double r25666 = 1/3; double r25667 = r25665 * r25666; double r25668 = fma(r25667, r25665, r25660); double r25669 = 5; double r25670 = pow(r25665, r25669); double r25671 = 1/60; double r25672 = r25670 * r25671; double r25673 = fma(r25665, r25668, r25672); double r25674 = U; double r25675 = fma(r25664, r25673, r25674); return r25675; } 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 r25676, r25677, r25678, r25679, r25680, r25681, r25682, r25683, r25684, r25685, r25686, r25687, r25688, r25689; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25676); mpfr_init(r25677); mpfr_init(r25678); mpfr_init(r25679); mpfr_init(r25680); mpfr_init(r25681); mpfr_init(r25682); mpfr_init(r25683); mpfr_init_set_str(r25684, "2", 10, MPFR_RNDN); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init(r25689); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r25676, J, MPFR_RNDN); mpfr_set_d(r25677, l, MPFR_RNDN); mpfr_exp(r25678, r25677, MPFR_RNDN); mpfr_neg(r25679, r25677, MPFR_RNDN); mpfr_exp(r25680, r25679, MPFR_RNDN); mpfr_sub(r25681, r25678, r25680, MPFR_RNDN); mpfr_mul(r25682, r25676, r25681, MPFR_RNDN); mpfr_set_d(r25683, K, MPFR_RNDN); ; mpfr_div(r25685, r25683, r25684, MPFR_RNDN); mpfr_cos(r25686, r25685, MPFR_RNDN); mpfr_mul(r25687, r25682, r25686, MPFR_RNDN); mpfr_set_d(r25688, U, MPFR_RNDN); mpfr_add(r25689, r25687, r25688, MPFR_RNDN); return mpfr_get_d(r25689, MPFR_RNDN); } static mpfr_t r25690, r25691, r25692, r25693, r25694, r25695, r25696, r25697, r25698, r25699, r25700, r25701, r25702, r25703, r25704, r25705, r25706; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25690); mpfr_init_set_str(r25691, "2", 10, MPFR_RNDN); mpfr_init(r25692); mpfr_init(r25693); mpfr_init(r25694); mpfr_init(r25695); mpfr_init(r25696); mpfr_init_set_str(r25697, "1/3", 10, MPFR_RNDN); mpfr_init(r25698); mpfr_init(r25699); mpfr_init_set_str(r25700, "5", 10, MPFR_RNDN); mpfr_init(r25701); mpfr_init_set_str(r25702, "1/60", 10, MPFR_RNDN); mpfr_init(r25703); mpfr_init(r25704); mpfr_init(r25705); mpfr_init(r25706); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r25690, K, MPFR_RNDN); ; mpfr_div(r25692, r25690, r25691, MPFR_RNDN); mpfr_cos(r25693, r25692, MPFR_RNDN); mpfr_set_d(r25694, J, MPFR_RNDN); mpfr_mul(r25695, r25693, r25694, MPFR_RNDN); mpfr_set_d(r25696, l, MPFR_RNDN); ; mpfr_mul(r25698, r25696, r25697, MPFR_RNDN); mpfr_fma(r25699, r25698, r25696, r25691, MPFR_RNDN); ; mpfr_pow(r25701, r25696, r25700, MPFR_RNDN); ; mpfr_mul(r25703, r25701, r25702, MPFR_RNDN); mpfr_fma(r25704, r25696, r25699, r25703, MPFR_RNDN); mpfr_set_d(r25705, U, MPFR_RNDN); mpfr_fma(r25706, r25695, r25704, r25705, MPFR_RNDN); return mpfr_get_d(r25706, MPFR_RNDN); } static mpfr_t r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25707); mpfr_init_set_str(r25708, "2", 10, MPFR_RNDN); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init(r25712); mpfr_init(r25713); mpfr_init_set_str(r25714, "1/3", 10, MPFR_RNDN); mpfr_init(r25715); mpfr_init(r25716); mpfr_init_set_str(r25717, "5", 10, MPFR_RNDN); mpfr_init(r25718); mpfr_init_set_str(r25719, "1/60", 10, MPFR_RNDN); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r25707, K, MPFR_RNDN); ; mpfr_div(r25709, r25707, r25708, MPFR_RNDN); mpfr_cos(r25710, r25709, MPFR_RNDN); mpfr_set_d(r25711, J, MPFR_RNDN); mpfr_mul(r25712, r25710, r25711, MPFR_RNDN); mpfr_set_d(r25713, l, MPFR_RNDN); ; mpfr_mul(r25715, r25713, r25714, MPFR_RNDN); mpfr_fma(r25716, r25715, r25713, r25708, MPFR_RNDN); ; mpfr_pow(r25718, r25713, r25717, MPFR_RNDN); ; mpfr_mul(r25720, r25718, r25719, MPFR_RNDN); mpfr_fma(r25721, r25713, r25716, r25720, MPFR_RNDN); mpfr_set_d(r25722, U, MPFR_RNDN); mpfr_fma(r25723, r25712, r25721, r25722, MPFR_RNDN); return mpfr_get_d(r25723, MPFR_RNDN); }