#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (4)"; double f_if(float J, float l, float K, float U) { float r23599 = J; float r23600 = l; float r23601 = exp(r23600); float r23602 = -r23600; float r23603 = exp(r23602); float r23604 = r23601 - r23603; float r23605 = r23599 * r23604; float r23606 = K; float r23607 = 2; float r23608 = r23606 / r23607; float r23609 = cos(r23608); float r23610 = r23605 * r23609; float r23611 = U; float r23612 = r23610 + r23611; return r23612; } double f_id(double J, double l, double K, double U) { double r23613 = J; double r23614 = l; double r23615 = exp(r23614); double r23616 = -r23614; double r23617 = exp(r23616); double r23618 = r23615 - r23617; double r23619 = r23613 * r23618; double r23620 = K; double r23621 = 2; double r23622 = r23620 / r23621; double r23623 = cos(r23622); double r23624 = r23619 * r23623; double r23625 = U; double r23626 = r23624 + r23625; return r23626; } double f_of(float J, float l, float K, float U) { float r23627 = K; float r23628 = 2; float r23629 = r23627 / r23628; float r23630 = cos(r23629); float r23631 = J; float r23632 = r23630 * r23631; float r23633 = l; float r23634 = 1/3; float r23635 = r23633 * r23634; float r23636 = fma(r23635, r23633, r23628); float r23637 = 5; float r23638 = pow(r23633, r23637); float r23639 = 1/60; float r23640 = r23638 * r23639; float r23641 = fma(r23633, r23636, r23640); float r23642 = U; float r23643 = fma(r23632, r23641, r23642); return r23643; } double f_od(double J, double l, double K, double U) { double r23644 = K; double r23645 = 2; double r23646 = r23644 / r23645; double r23647 = cos(r23646); double r23648 = J; double r23649 = r23647 * r23648; double r23650 = l; double r23651 = 1/3; double r23652 = r23650 * r23651; double r23653 = fma(r23652, r23650, r23645); double r23654 = 5; double r23655 = pow(r23650, r23654); double r23656 = 1/60; double r23657 = r23655 * r23656; double r23658 = fma(r23650, r23653, r23657); double r23659 = U; double r23660 = fma(r23649, r23658, r23659); return r23660; } 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 r23661, r23662, r23663, r23664, r23665, r23666, r23667, r23668, r23669, r23670, r23671, r23672, r23673, r23674; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r23661); mpfr_init(r23662); mpfr_init(r23663); mpfr_init(r23664); mpfr_init(r23665); mpfr_init(r23666); mpfr_init(r23667); mpfr_init(r23668); mpfr_init_set_str(r23669, "2", 10, MPFR_RNDN); mpfr_init(r23670); mpfr_init(r23671); mpfr_init(r23672); mpfr_init(r23673); mpfr_init(r23674); } double f_im(double J, double l, double K, double U) { mpfr_set_d(r23661, J, MPFR_RNDN); mpfr_set_d(r23662, l, MPFR_RNDN); mpfr_exp(r23663, r23662, MPFR_RNDN); mpfr_neg(r23664, r23662, MPFR_RNDN); mpfr_exp(r23665, r23664, MPFR_RNDN); mpfr_sub(r23666, r23663, r23665, MPFR_RNDN); mpfr_mul(r23667, r23661, r23666, MPFR_RNDN); mpfr_set_d(r23668, K, MPFR_RNDN); ; mpfr_div(r23670, r23668, r23669, MPFR_RNDN); mpfr_cos(r23671, r23670, MPFR_RNDN); mpfr_mul(r23672, r23667, r23671, MPFR_RNDN); mpfr_set_d(r23673, U, MPFR_RNDN); mpfr_add(r23674, r23672, r23673, MPFR_RNDN); return mpfr_get_d(r23674, MPFR_RNDN); } static mpfr_t r23675, r23676, r23677, r23678, r23679, r23680, r23681, r23682, r23683, r23684, r23685, r23686, r23687, r23688, r23689, r23690, r23691; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23675); mpfr_init_set_str(r23676, "2", 10, MPFR_RNDN); mpfr_init(r23677); mpfr_init(r23678); mpfr_init(r23679); mpfr_init(r23680); mpfr_init(r23681); mpfr_init_set_str(r23682, "1/3", 10, MPFR_RNDN); mpfr_init(r23683); mpfr_init(r23684); mpfr_init_set_str(r23685, "5", 10, MPFR_RNDN); mpfr_init(r23686); mpfr_init_set_str(r23687, "1/60", 10, MPFR_RNDN); mpfr_init(r23688); mpfr_init(r23689); mpfr_init(r23690); mpfr_init(r23691); } double f_fm(double J, double l, double K, double U) { mpfr_set_d(r23675, K, MPFR_RNDN); ; mpfr_div(r23677, r23675, r23676, MPFR_RNDN); mpfr_cos(r23678, r23677, MPFR_RNDN); mpfr_set_d(r23679, J, MPFR_RNDN); mpfr_mul(r23680, r23678, r23679, MPFR_RNDN); mpfr_set_d(r23681, l, MPFR_RNDN); ; mpfr_mul(r23683, r23681, r23682, MPFR_RNDN); mpfr_fma(r23684, r23683, r23681, r23676, MPFR_RNDN); ; mpfr_pow(r23686, r23681, r23685, MPFR_RNDN); ; mpfr_mul(r23688, r23686, r23687, MPFR_RNDN); mpfr_fma(r23689, r23681, r23684, r23688, MPFR_RNDN); mpfr_set_d(r23690, U, MPFR_RNDN); mpfr_fma(r23691, r23680, r23689, r23690, MPFR_RNDN); return mpfr_get_d(r23691, MPFR_RNDN); } static mpfr_t r23692, r23693, r23694, r23695, r23696, r23697, r23698, r23699, r23700, r23701, r23702, r23703, r23704, r23705, r23706, r23707, r23708; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23692); mpfr_init_set_str(r23693, "2", 10, MPFR_RNDN); mpfr_init(r23694); mpfr_init(r23695); mpfr_init(r23696); mpfr_init(r23697); mpfr_init(r23698); mpfr_init_set_str(r23699, "1/3", 10, MPFR_RNDN); mpfr_init(r23700); mpfr_init(r23701); mpfr_init_set_str(r23702, "5", 10, MPFR_RNDN); mpfr_init(r23703); mpfr_init_set_str(r23704, "1/60", 10, MPFR_RNDN); mpfr_init(r23705); mpfr_init(r23706); mpfr_init(r23707); mpfr_init(r23708); } double f_dm(double J, double l, double K, double U) { mpfr_set_d(r23692, K, MPFR_RNDN); ; mpfr_div(r23694, r23692, r23693, MPFR_RNDN); mpfr_cos(r23695, r23694, MPFR_RNDN); mpfr_set_d(r23696, J, MPFR_RNDN); mpfr_mul(r23697, r23695, r23696, MPFR_RNDN); mpfr_set_d(r23698, l, MPFR_RNDN); ; mpfr_mul(r23700, r23698, r23699, MPFR_RNDN); mpfr_fma(r23701, r23700, r23698, r23693, MPFR_RNDN); ; mpfr_pow(r23703, r23698, r23702, MPFR_RNDN); ; mpfr_mul(r23705, r23703, r23704, MPFR_RNDN); mpfr_fma(r23706, r23698, r23701, r23705, MPFR_RNDN); mpfr_set_d(r23707, U, MPFR_RNDN); mpfr_fma(r23708, r23697, r23706, r23707, MPFR_RNDN); return mpfr_get_d(r23708, MPFR_RNDN); }