#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (3)"; double f_if(float J, float K, float U) { float r22525 = -2; float r22526 = J; float r22527 = r22525 * r22526; float r22528 = K; float r22529 = 2; float r22530 = r22528 / r22529; float r22531 = cos(r22530); float r22532 = r22527 * r22531; float r22533 = 1; float r22534 = U; float r22535 = r22529 * r22526; float r22536 = r22535 * r22531; float r22537 = r22534 / r22536; float r22538 = pow(r22537, r22529); float r22539 = r22533 + r22538; float r22540 = sqrt(r22539); float r22541 = r22532 * r22540; return r22541; } double f_id(double J, double K, double U) { double r22542 = -2; double r22543 = J; double r22544 = r22542 * r22543; double r22545 = K; double r22546 = 2; double r22547 = r22545 / r22546; double r22548 = cos(r22547); double r22549 = r22544 * r22548; double r22550 = 1; double r22551 = U; double r22552 = r22546 * r22543; double r22553 = r22552 * r22548; double r22554 = r22551 / r22553; double r22555 = pow(r22554, r22546); double r22556 = r22550 + r22555; double r22557 = sqrt(r22556); double r22558 = r22549 * r22557; return r22558; } double f_of(float J, float K, float U) { float r22559 = J; float r22560 = -2; float r22561 = r22559 * r22560; float r22562 = K; float r22563 = 2; float r22564 = r22562 / r22563; float r22565 = cos(r22564); float r22566 = 1; float r22567 = U; float r22568 = r22559 + r22559; float r22569 = r22567 / r22568; float r22570 = r22569 / r22565; float r22571 = hypot(r22566, r22570); float r22572 = r22565 * r22571; float r22573 = r22561 * r22572; return r22573; } double f_od(double J, double K, double U) { double r22574 = J; double r22575 = -2; double r22576 = r22574 * r22575; double r22577 = K; double r22578 = 2; double r22579 = r22577 / r22578; double r22580 = cos(r22579); double r22581 = 1; double r22582 = U; double r22583 = r22574 + r22574; double r22584 = r22582 / r22583; double r22585 = r22584 / r22580; double r22586 = hypot(r22581, r22585); double r22587 = r22580 * r22586; double r22588 = r22576 * r22587; return r22588; } 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 r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22589, "-2", 10, MPFR_RNDN); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init_set_str(r22593, "2", 10, MPFR_RNDN); mpfr_init(r22594); mpfr_init(r22595); mpfr_init(r22596); mpfr_init_set_str(r22597, "1", 10, MPFR_RNDN); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); } double f_im(double J, double K, double U) { ; mpfr_set_d(r22590, J, MPFR_RNDN); mpfr_mul(r22591, r22589, r22590, MPFR_RNDN); mpfr_set_d(r22592, K, MPFR_RNDN); ; mpfr_div(r22594, r22592, r22593, MPFR_RNDN); mpfr_cos(r22595, r22594, MPFR_RNDN); mpfr_mul(r22596, r22591, r22595, MPFR_RNDN); ; mpfr_set_d(r22598, U, MPFR_RNDN); mpfr_mul(r22599, r22593, r22590, MPFR_RNDN); mpfr_mul(r22600, r22599, r22595, MPFR_RNDN); mpfr_div(r22601, r22598, r22600, MPFR_RNDN); mpfr_pow(r22602, r22601, r22593, MPFR_RNDN); mpfr_add(r22603, r22597, r22602, MPFR_RNDN); mpfr_sqrt(r22604, r22603, MPFR_RNDN); mpfr_mul(r22605, r22596, r22604, MPFR_RNDN); return mpfr_get_d(r22605, MPFR_RNDN); } static mpfr_t r22606, r22607, r22608, r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619, r22620; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22606); mpfr_init_set_str(r22607, "-2", 10, MPFR_RNDN); mpfr_init(r22608); mpfr_init(r22609); mpfr_init_set_str(r22610, "2", 10, MPFR_RNDN); mpfr_init(r22611); mpfr_init(r22612); mpfr_init_set_str(r22613, "1", 10, MPFR_RNDN); mpfr_init(r22614); mpfr_init(r22615); mpfr_init(r22616); mpfr_init(r22617); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); } double f_fm(double J, double K, double U) { mpfr_set_d(r22606, J, MPFR_RNDN); ; mpfr_mul(r22608, r22606, r22607, MPFR_RNDN); mpfr_set_d(r22609, K, MPFR_RNDN); ; mpfr_div(r22611, r22609, r22610, MPFR_RNDN); mpfr_cos(r22612, r22611, MPFR_RNDN); ; mpfr_set_d(r22614, U, MPFR_RNDN); mpfr_add(r22615, r22606, r22606, MPFR_RNDN); mpfr_div(r22616, r22614, r22615, MPFR_RNDN); mpfr_div(r22617, r22616, r22612, MPFR_RNDN); mpfr_hypot(r22618, r22613, r22617, MPFR_RNDN); mpfr_mul(r22619, r22612, r22618, MPFR_RNDN); mpfr_mul(r22620, r22608, r22619, MPFR_RNDN); return mpfr_get_d(r22620, MPFR_RNDN); } static mpfr_t r22621, r22622, r22623, r22624, r22625, r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22621); mpfr_init_set_str(r22622, "-2", 10, MPFR_RNDN); mpfr_init(r22623); mpfr_init(r22624); mpfr_init_set_str(r22625, "2", 10, MPFR_RNDN); mpfr_init(r22626); mpfr_init(r22627); mpfr_init_set_str(r22628, "1", 10, MPFR_RNDN); mpfr_init(r22629); mpfr_init(r22630); mpfr_init(r22631); mpfr_init(r22632); mpfr_init(r22633); mpfr_init(r22634); mpfr_init(r22635); } double f_dm(double J, double K, double U) { mpfr_set_d(r22621, J, MPFR_RNDN); ; mpfr_mul(r22623, r22621, r22622, MPFR_RNDN); mpfr_set_d(r22624, K, MPFR_RNDN); ; mpfr_div(r22626, r22624, r22625, MPFR_RNDN); mpfr_cos(r22627, r22626, MPFR_RNDN); ; mpfr_set_d(r22629, U, MPFR_RNDN); mpfr_add(r22630, r22621, r22621, MPFR_RNDN); mpfr_div(r22631, r22629, r22630, MPFR_RNDN); mpfr_div(r22632, r22631, r22627, MPFR_RNDN); mpfr_hypot(r22633, r22628, r22632, MPFR_RNDN); mpfr_mul(r22634, r22627, r22633, MPFR_RNDN); mpfr_mul(r22635, r22623, r22634, MPFR_RNDN); return mpfr_get_d(r22635, MPFR_RNDN); }