#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (3)"; double f_if(float J, float K, float U) { float r9576 = -2.0f; float r9577 = J; float r9578 = r9576 * r9577; float r9579 = K; float r9580 = 2.0f; float r9581 = r9579 / r9580; float r9582 = cos(r9581); float r9583 = r9578 * r9582; float r9584 = 1.0f; float r9585 = U; float r9586 = r9580 * r9577; float r9587 = r9586 * r9582; float r9588 = r9585 / r9587; float r9589 = pow(r9588, r9580); float r9590 = r9584 + r9589; float r9591 = sqrt(r9590); float r9592 = r9583 * r9591; return r9592; } double f_id(double J, double K, double U) { double r9593 = -2.0; double r9594 = J; double r9595 = r9593 * r9594; double r9596 = K; double r9597 = 2.0; double r9598 = r9596 / r9597; double r9599 = cos(r9598); double r9600 = r9595 * r9599; double r9601 = 1.0; double r9602 = U; double r9603 = r9597 * r9594; double r9604 = r9603 * r9599; double r9605 = r9602 / r9604; double r9606 = pow(r9605, r9597); double r9607 = r9601 + r9606; double r9608 = sqrt(r9607); double r9609 = r9600 * r9608; return r9609; } double f_of(float J, float K, float U) { float r9610 = U; float r9611 = J; float r9612 = 2.0f; float r9613 = r9611 * r9612; float r9614 = K; float r9615 = r9614 / r9612; float r9616 = cos(r9615); float r9617 = r9613 * r9616; float r9618 = r9610 / r9617; float r9619 = pow(r9618, r9612); float r9620 = 1.0f; float r9621 = r9619 + r9620; float r9622 = sqrt(r9621); float r9623 = -2.0f; float r9624 = r9623 * r9611; float r9625 = r9616 * r9624; float r9626 = r9622 * r9625; return r9626; } double f_od(double J, double K, double U) { double r9627 = U; double r9628 = J; double r9629 = 2.0; double r9630 = r9628 * r9629; double r9631 = K; double r9632 = r9631 / r9629; double r9633 = cos(r9632); double r9634 = r9630 * r9633; double r9635 = r9627 / r9634; double r9636 = pow(r9635, r9629); double r9637 = 1.0; double r9638 = r9636 + r9637; double r9639 = sqrt(r9638); double r9640 = -2.0; double r9641 = r9640 * r9628; double r9642 = r9633 * r9641; double r9643 = r9639 * r9642; return r9643; } 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 r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658, r9659, r9660; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r9644, "-2", 10, MPFR_RNDN); mpfr_init(r9645); mpfr_init(r9646); mpfr_init(r9647); mpfr_init_set_str(r9648, "2", 10, MPFR_RNDN); mpfr_init(r9649); mpfr_init(r9650); mpfr_init(r9651); mpfr_init_set_str(r9652, "1", 10, MPFR_RNDN); mpfr_init(r9653); mpfr_init(r9654); mpfr_init(r9655); mpfr_init(r9656); mpfr_init(r9657); mpfr_init(r9658); mpfr_init(r9659); mpfr_init(r9660); } double f_im(double J, double K, double U) { ; mpfr_set_d(r9645, J, MPFR_RNDN); mpfr_mul(r9646, r9644, r9645, MPFR_RNDN); mpfr_set_d(r9647, K, MPFR_RNDN); ; mpfr_div(r9649, r9647, r9648, MPFR_RNDN); mpfr_cos(r9650, r9649, MPFR_RNDN); mpfr_mul(r9651, r9646, r9650, MPFR_RNDN); ; mpfr_set_d(r9653, U, MPFR_RNDN); mpfr_mul(r9654, r9648, r9645, MPFR_RNDN); mpfr_mul(r9655, r9654, r9650, MPFR_RNDN); mpfr_div(r9656, r9653, r9655, MPFR_RNDN); mpfr_pow(r9657, r9656, r9648, MPFR_RNDN); mpfr_add(r9658, r9652, r9657, MPFR_RNDN); mpfr_sqrt(r9659, r9658, MPFR_RNDN); mpfr_mul(r9660, r9651, r9659, MPFR_RNDN); return mpfr_get_d(r9660, MPFR_RNDN); } static mpfr_t r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r9661); mpfr_init(r9662); mpfr_init_set_str(r9663, "2", 10, MPFR_RNDN); mpfr_init(r9664); mpfr_init(r9665); mpfr_init(r9666); mpfr_init(r9667); mpfr_init(r9668); mpfr_init(r9669); mpfr_init(r9670); mpfr_init_set_str(r9671, "1", 10, MPFR_RNDN); mpfr_init(r9672); mpfr_init(r9673); mpfr_init_set_str(r9674, "-2", 10, MPFR_RNDN); mpfr_init(r9675); mpfr_init(r9676); mpfr_init(r9677); } double f_fm(double J, double K, double U) { mpfr_set_d(r9661, U, MPFR_RNDN); mpfr_set_d(r9662, J, MPFR_RNDN); ; mpfr_mul(r9664, r9662, r9663, MPFR_RNDN); mpfr_set_d(r9665, K, MPFR_RNDN); mpfr_div(r9666, r9665, r9663, MPFR_RNDN); mpfr_cos(r9667, r9666, MPFR_RNDN); mpfr_mul(r9668, r9664, r9667, MPFR_RNDN); mpfr_div(r9669, r9661, r9668, MPFR_RNDN); mpfr_pow(r9670, r9669, r9663, MPFR_RNDN); ; mpfr_add(r9672, r9670, r9671, MPFR_RNDN); mpfr_sqrt(r9673, r9672, MPFR_RNDN); ; mpfr_mul(r9675, r9674, r9662, MPFR_RNDN); mpfr_mul(r9676, r9667, r9675, MPFR_RNDN); mpfr_mul(r9677, r9673, r9676, MPFR_RNDN); return mpfr_get_d(r9677, MPFR_RNDN); } static mpfr_t r9678, r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689, r9690, r9691, r9692, r9693, r9694; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r9678); mpfr_init(r9679); mpfr_init_set_str(r9680, "2", 10, MPFR_RNDN); mpfr_init(r9681); mpfr_init(r9682); mpfr_init(r9683); mpfr_init(r9684); mpfr_init(r9685); mpfr_init(r9686); mpfr_init(r9687); mpfr_init_set_str(r9688, "1", 10, MPFR_RNDN); mpfr_init(r9689); mpfr_init(r9690); mpfr_init_set_str(r9691, "-2", 10, MPFR_RNDN); mpfr_init(r9692); mpfr_init(r9693); mpfr_init(r9694); } double f_dm(double J, double K, double U) { mpfr_set_d(r9678, U, MPFR_RNDN); mpfr_set_d(r9679, J, MPFR_RNDN); ; mpfr_mul(r9681, r9679, r9680, MPFR_RNDN); mpfr_set_d(r9682, K, MPFR_RNDN); mpfr_div(r9683, r9682, r9680, MPFR_RNDN); mpfr_cos(r9684, r9683, MPFR_RNDN); mpfr_mul(r9685, r9681, r9684, MPFR_RNDN); mpfr_div(r9686, r9678, r9685, MPFR_RNDN); mpfr_pow(r9687, r9686, r9680, MPFR_RNDN); ; mpfr_add(r9689, r9687, r9688, MPFR_RNDN); mpfr_sqrt(r9690, r9689, MPFR_RNDN); ; mpfr_mul(r9692, r9691, r9679, MPFR_RNDN); mpfr_mul(r9693, r9684, r9692, MPFR_RNDN); mpfr_mul(r9694, r9690, r9693, MPFR_RNDN); return mpfr_get_d(r9694, MPFR_RNDN); }