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