#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (3)"; double f_if(float J, float K, float U) { float r22493 = -2; float r22494 = J; float r22495 = r22493 * r22494; float r22496 = K; float r22497 = 2; float r22498 = r22496 / r22497; float r22499 = cos(r22498); float r22500 = r22495 * r22499; float r22501 = 1; float r22502 = U; float r22503 = r22497 * r22494; float r22504 = r22503 * r22499; float r22505 = r22502 / r22504; float r22506 = pow(r22505, r22497); float r22507 = r22501 + r22506; float r22508 = sqrt(r22507); float r22509 = r22500 * r22508; return r22509; } double f_id(double J, double K, double U) { double r22510 = -2; double r22511 = J; double r22512 = r22510 * r22511; double r22513 = K; double r22514 = 2; double r22515 = r22513 / r22514; double r22516 = cos(r22515); double r22517 = r22512 * r22516; double r22518 = 1; double r22519 = U; double r22520 = r22514 * r22511; double r22521 = r22520 * r22516; double r22522 = r22519 / r22521; double r22523 = pow(r22522, r22514); double r22524 = r22518 + r22523; double r22525 = sqrt(r22524); double r22526 = r22517 * r22525; return r22526; } double f_of(float J, float K, float U) { float r22527 = J; float r22528 = -2; float r22529 = r22527 * r22528; float r22530 = K; float r22531 = 2; float r22532 = r22530 / r22531; float r22533 = cos(r22532); float r22534 = 1; float r22535 = U; float r22536 = r22527 + r22527; float r22537 = r22535 / r22536; float r22538 = r22537 / r22533; float r22539 = hypot(r22534, r22538); float r22540 = r22533 * r22539; float r22541 = r22529 * r22540; return r22541; } double f_od(double J, double K, double U) { double r22542 = J; double r22543 = -2; double r22544 = r22542 * r22543; double r22545 = K; double r22546 = 2; double r22547 = r22545 / r22546; double r22548 = cos(r22547); double r22549 = 1; double r22550 = U; double r22551 = r22542 + r22542; double r22552 = r22550 / r22551; double r22553 = r22552 / r22548; double r22554 = hypot(r22549, r22553); double r22555 = r22548 * r22554; double r22556 = r22544 * r22555; return r22556; } 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 r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22557, "-2", 10, MPFR_RNDN); mpfr_init(r22558); mpfr_init(r22559); mpfr_init(r22560); mpfr_init_set_str(r22561, "2", 10, MPFR_RNDN); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init_set_str(r22565, "1", 10, MPFR_RNDN); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); } double f_im(double J, double K, double U) { ; mpfr_set_d(r22558, J, MPFR_RNDN); mpfr_mul(r22559, r22557, r22558, MPFR_RNDN); mpfr_set_d(r22560, K, MPFR_RNDN); ; mpfr_div(r22562, r22560, r22561, MPFR_RNDN); mpfr_cos(r22563, r22562, MPFR_RNDN); mpfr_mul(r22564, r22559, r22563, MPFR_RNDN); ; mpfr_set_d(r22566, U, MPFR_RNDN); mpfr_mul(r22567, r22561, r22558, MPFR_RNDN); mpfr_mul(r22568, r22567, r22563, MPFR_RNDN); mpfr_div(r22569, r22566, r22568, MPFR_RNDN); mpfr_pow(r22570, r22569, r22561, MPFR_RNDN); mpfr_add(r22571, r22565, r22570, MPFR_RNDN); mpfr_sqrt(r22572, r22571, MPFR_RNDN); mpfr_mul(r22573, r22564, r22572, MPFR_RNDN); return mpfr_get_d(r22573, MPFR_RNDN); } static mpfr_t r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22574); mpfr_init_set_str(r22575, "-2", 10, MPFR_RNDN); mpfr_init(r22576); mpfr_init(r22577); mpfr_init_set_str(r22578, "2", 10, MPFR_RNDN); mpfr_init(r22579); mpfr_init(r22580); mpfr_init_set_str(r22581, "1", 10, MPFR_RNDN); mpfr_init(r22582); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); } double f_fm(double J, double K, double U) { mpfr_set_d(r22574, J, MPFR_RNDN); ; mpfr_mul(r22576, r22574, r22575, MPFR_RNDN); mpfr_set_d(r22577, K, MPFR_RNDN); ; mpfr_div(r22579, r22577, r22578, MPFR_RNDN); mpfr_cos(r22580, r22579, MPFR_RNDN); ; mpfr_set_d(r22582, U, MPFR_RNDN); mpfr_add(r22583, r22574, r22574, MPFR_RNDN); mpfr_div(r22584, r22582, r22583, MPFR_RNDN); mpfr_div(r22585, r22584, r22580, MPFR_RNDN); mpfr_hypot(r22586, r22581, r22585, MPFR_RNDN); mpfr_mul(r22587, r22580, r22586, MPFR_RNDN); mpfr_mul(r22588, r22576, r22587, MPFR_RNDN); return mpfr_get_d(r22588, MPFR_RNDN); } static mpfr_t r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22589); mpfr_init_set_str(r22590, "-2", 10, MPFR_RNDN); mpfr_init(r22591); mpfr_init(r22592); mpfr_init_set_str(r22593, "2", 10, MPFR_RNDN); mpfr_init(r22594); mpfr_init(r22595); mpfr_init_set_str(r22596, "1", 10, MPFR_RNDN); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); } double f_dm(double J, double K, double U) { mpfr_set_d(r22589, 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_set_d(r22597, U, MPFR_RNDN); mpfr_add(r22598, r22589, r22589, MPFR_RNDN); mpfr_div(r22599, r22597, r22598, MPFR_RNDN); mpfr_div(r22600, r22599, r22595, MPFR_RNDN); mpfr_hypot(r22601, r22596, r22600, MPFR_RNDN); mpfr_mul(r22602, r22595, r22601, MPFR_RNDN); mpfr_mul(r22603, r22591, r22602, MPFR_RNDN); return mpfr_get_d(r22603, MPFR_RNDN); }