#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (3)"; double f_if(float J, float K, float U) { float r22455 = -2; float r22456 = J; float r22457 = r22455 * r22456; float r22458 = K; float r22459 = 2; float r22460 = r22458 / r22459; float r22461 = cos(r22460); float r22462 = r22457 * r22461; float r22463 = 1; float r22464 = U; float r22465 = r22459 * r22456; float r22466 = r22465 * r22461; float r22467 = r22464 / r22466; float r22468 = pow(r22467, r22459); float r22469 = r22463 + r22468; float r22470 = sqrt(r22469); float r22471 = r22462 * r22470; return r22471; } double f_id(double J, double K, double U) { double r22472 = -2; double r22473 = J; double r22474 = r22472 * r22473; double r22475 = K; double r22476 = 2; double r22477 = r22475 / r22476; double r22478 = cos(r22477); double r22479 = r22474 * r22478; double r22480 = 1; double r22481 = U; double r22482 = r22476 * r22473; double r22483 = r22482 * r22478; double r22484 = r22481 / r22483; double r22485 = pow(r22484, r22476); double r22486 = r22480 + r22485; double r22487 = sqrt(r22486); double r22488 = r22479 * r22487; return r22488; } double f_of(float J, float K, float U) { float r22489 = J; float r22490 = -2; float r22491 = r22489 * r22490; float r22492 = K; float r22493 = 2; float r22494 = r22492 / r22493; float r22495 = cos(r22494); float r22496 = 1; float r22497 = U; float r22498 = r22489 + r22489; float r22499 = r22497 / r22498; float r22500 = r22499 / r22495; float r22501 = hypot(r22496, r22500); float r22502 = r22495 * r22501; float r22503 = r22491 * r22502; return r22503; } double f_od(double J, double K, double U) { double r22504 = J; double r22505 = -2; double r22506 = r22504 * r22505; double r22507 = K; double r22508 = 2; double r22509 = r22507 / r22508; double r22510 = cos(r22509); double r22511 = 1; double r22512 = U; double r22513 = r22504 + r22504; double r22514 = r22512 / r22513; double r22515 = r22514 / r22510; double r22516 = hypot(r22511, r22515); double r22517 = r22510 * r22516; double r22518 = r22506 * r22517; return r22518; } 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 r22519, r22520, r22521, r22522, r22523, r22524, r22525, r22526, r22527, r22528, r22529, r22530, r22531, r22532, r22533, r22534, r22535; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22519, "-2", 10, MPFR_RNDN); mpfr_init(r22520); mpfr_init(r22521); mpfr_init(r22522); mpfr_init_set_str(r22523, "2", 10, MPFR_RNDN); mpfr_init(r22524); mpfr_init(r22525); mpfr_init(r22526); mpfr_init_set_str(r22527, "1", 10, MPFR_RNDN); mpfr_init(r22528); mpfr_init(r22529); mpfr_init(r22530); mpfr_init(r22531); mpfr_init(r22532); mpfr_init(r22533); mpfr_init(r22534); mpfr_init(r22535); } double f_im(double J, double K, double U) { ; mpfr_set_d(r22520, J, MPFR_RNDN); mpfr_mul(r22521, r22519, r22520, MPFR_RNDN); mpfr_set_d(r22522, K, MPFR_RNDN); ; mpfr_div(r22524, r22522, r22523, MPFR_RNDN); mpfr_cos(r22525, r22524, MPFR_RNDN); mpfr_mul(r22526, r22521, r22525, MPFR_RNDN); ; mpfr_set_d(r22528, U, MPFR_RNDN); mpfr_mul(r22529, r22523, r22520, MPFR_RNDN); mpfr_mul(r22530, r22529, r22525, MPFR_RNDN); mpfr_div(r22531, r22528, r22530, MPFR_RNDN); mpfr_pow(r22532, r22531, r22523, MPFR_RNDN); mpfr_add(r22533, r22527, r22532, MPFR_RNDN); mpfr_sqrt(r22534, r22533, MPFR_RNDN); mpfr_mul(r22535, r22526, r22534, MPFR_RNDN); return mpfr_get_d(r22535, MPFR_RNDN); } static mpfr_t r22536, r22537, r22538, r22539, r22540, r22541, r22542, r22543, r22544, r22545, r22546, r22547, r22548, r22549, r22550; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22536); mpfr_init_set_str(r22537, "-2", 10, MPFR_RNDN); mpfr_init(r22538); mpfr_init(r22539); mpfr_init_set_str(r22540, "2", 10, MPFR_RNDN); mpfr_init(r22541); mpfr_init(r22542); mpfr_init_set_str(r22543, "1", 10, MPFR_RNDN); mpfr_init(r22544); mpfr_init(r22545); mpfr_init(r22546); mpfr_init(r22547); mpfr_init(r22548); mpfr_init(r22549); mpfr_init(r22550); } double f_fm(double J, double K, double U) { mpfr_set_d(r22536, J, MPFR_RNDN); ; mpfr_mul(r22538, r22536, r22537, MPFR_RNDN); mpfr_set_d(r22539, K, MPFR_RNDN); ; mpfr_div(r22541, r22539, r22540, MPFR_RNDN); mpfr_cos(r22542, r22541, MPFR_RNDN); ; mpfr_set_d(r22544, U, MPFR_RNDN); mpfr_add(r22545, r22536, r22536, MPFR_RNDN); mpfr_div(r22546, r22544, r22545, MPFR_RNDN); mpfr_div(r22547, r22546, r22542, MPFR_RNDN); mpfr_hypot(r22548, r22543, r22547, MPFR_RNDN); mpfr_mul(r22549, r22542, r22548, MPFR_RNDN); mpfr_mul(r22550, r22538, r22549, MPFR_RNDN); return mpfr_get_d(r22550, MPFR_RNDN); } static mpfr_t r22551, r22552, r22553, r22554, r22555, r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22551); mpfr_init_set_str(r22552, "-2", 10, MPFR_RNDN); mpfr_init(r22553); mpfr_init(r22554); mpfr_init_set_str(r22555, "2", 10, MPFR_RNDN); mpfr_init(r22556); mpfr_init(r22557); mpfr_init_set_str(r22558, "1", 10, MPFR_RNDN); mpfr_init(r22559); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); } double f_dm(double J, double K, double U) { mpfr_set_d(r22551, J, MPFR_RNDN); ; mpfr_mul(r22553, r22551, r22552, MPFR_RNDN); mpfr_set_d(r22554, K, MPFR_RNDN); ; mpfr_div(r22556, r22554, r22555, MPFR_RNDN); mpfr_cos(r22557, r22556, MPFR_RNDN); ; mpfr_set_d(r22559, U, MPFR_RNDN); mpfr_add(r22560, r22551, r22551, MPFR_RNDN); mpfr_div(r22561, r22559, r22560, MPFR_RNDN); mpfr_div(r22562, r22561, r22557, MPFR_RNDN); mpfr_hypot(r22563, r22558, r22562, MPFR_RNDN); mpfr_mul(r22564, r22557, r22563, MPFR_RNDN); mpfr_mul(r22565, r22553, r22564, MPFR_RNDN); return mpfr_get_d(r22565, MPFR_RNDN); }