#include #include #include #include #include char *name = "Maksimov and Kolovsky, Equation (32)"; double f_if(float K, float m, float n, float M, float l) { float r23483 = K; float r23484 = m; float r23485 = n; float r23486 = r23484 + r23485; float r23487 = r23483 * r23486; float r23488 = 2; float r23489 = r23487 / r23488; float r23490 = M; float r23491 = r23489 - r23490; float r23492 = cos(r23491); float r23493 = r23486 / r23488; float r23494 = r23493 - r23490; float r23495 = pow(r23494, r23488); float r23496 = -r23495; float r23497 = l; float r23498 = r23484 - r23485; float r23499 = fabs(r23498); float r23500 = r23497 - r23499; float r23501 = r23496 - r23500; float r23502 = exp(r23501); float r23503 = r23492 * r23502; return r23503; } double f_id(double K, double m, double n, double M, double l) { double r23504 = K; double r23505 = m; double r23506 = n; double r23507 = r23505 + r23506; double r23508 = r23504 * r23507; double r23509 = 2; double r23510 = r23508 / r23509; double r23511 = M; double r23512 = r23510 - r23511; double r23513 = cos(r23512); double r23514 = r23507 / r23509; double r23515 = r23514 - r23511; double r23516 = pow(r23515, r23509); double r23517 = -r23516; double r23518 = l; double r23519 = r23505 - r23506; double r23520 = fabs(r23519); double r23521 = r23518 - r23520; double r23522 = r23517 - r23521; double r23523 = exp(r23522); double r23524 = r23513 * r23523; return r23524; } double f_of(float __attribute__((unused)) K, float m, float n, float M, float l) { float r23525 = 1; float r23526 = m; float r23527 = n; float r23528 = r23526 + r23527; float r23529 = 2; float r23530 = r23528 / r23529; float r23531 = M; float r23532 = r23530 - r23531; float r23533 = pow(r23532, r23529); float r23534 = cbrt(r23533); float r23535 = cbrt(r23534); float r23536 = r23535 * r23535; float r23537 = r23534 * r23536; float r23538 = r23537 * r23535; float r23539 = r23536 * r23535; float r23540 = r23538 * r23539; float r23541 = -r23540; float r23542 = l; float r23543 = r23526 - r23527; float r23544 = fabs(r23543); float r23545 = r23542 - r23544; float r23546 = r23541 - r23545; float r23547 = exp(r23546); float r23548 = r23525 * r23547; return r23548; } double f_od(double __attribute__((unused)) K, double m, double n, double M, double l) { double r23549 = 1; double r23550 = m; double r23551 = n; double r23552 = r23550 + r23551; double r23553 = 2; double r23554 = r23552 / r23553; double r23555 = M; double r23556 = r23554 - r23555; double r23557 = pow(r23556, r23553); double r23558 = cbrt(r23557); double r23559 = cbrt(r23558); double r23560 = r23559 * r23559; double r23561 = r23558 * r23560; double r23562 = r23561 * r23559; double r23563 = r23560 * r23559; double r23564 = r23562 * r23563; double r23565 = -r23564; double r23566 = l; double r23567 = r23550 - r23551; double r23568 = fabs(r23567); double r23569 = r23566 - r23568; double r23570 = r23565 - r23569; double r23571 = exp(r23570); double r23572 = r23549 * r23571; return r23572; } 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 r23573, r23574, r23575, r23576, r23577, r23578, r23579, r23580, r23581, r23582, r23583, r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593; void setup_mpfr_f_im() { mpfr_set_default_prec(1616); mpfr_init(r23573); mpfr_init(r23574); mpfr_init(r23575); mpfr_init(r23576); mpfr_init(r23577); mpfr_init_set_str(r23578, "2", 10, MPFR_RNDN); mpfr_init(r23579); mpfr_init(r23580); mpfr_init(r23581); mpfr_init(r23582); mpfr_init(r23583); mpfr_init(r23584); mpfr_init(r23585); mpfr_init(r23586); mpfr_init(r23587); mpfr_init(r23588); mpfr_init(r23589); mpfr_init(r23590); mpfr_init(r23591); mpfr_init(r23592); mpfr_init(r23593); } double f_im(double K, double m, double n, double M, double l) { mpfr_set_d(r23573, K, MPFR_RNDN); mpfr_set_d(r23574, m, MPFR_RNDN); mpfr_set_d(r23575, n, MPFR_RNDN); mpfr_add(r23576, r23574, r23575, MPFR_RNDN); mpfr_mul(r23577, r23573, r23576, MPFR_RNDN); ; mpfr_div(r23579, r23577, r23578, MPFR_RNDN); mpfr_set_d(r23580, M, MPFR_RNDN); mpfr_sub(r23581, r23579, r23580, MPFR_RNDN); mpfr_cos(r23582, r23581, MPFR_RNDN); mpfr_div(r23583, r23576, r23578, MPFR_RNDN); mpfr_sub(r23584, r23583, r23580, MPFR_RNDN); mpfr_pow(r23585, r23584, r23578, MPFR_RNDN); mpfr_neg(r23586, r23585, MPFR_RNDN); mpfr_set_d(r23587, l, MPFR_RNDN); mpfr_sub(r23588, r23574, r23575, MPFR_RNDN); mpfr_abs(r23589, r23588, MPFR_RNDN); mpfr_sub(r23590, r23587, r23589, MPFR_RNDN); mpfr_sub(r23591, r23586, r23590, MPFR_RNDN); mpfr_exp(r23592, r23591, MPFR_RNDN); mpfr_mul(r23593, r23582, r23592, MPFR_RNDN); return mpfr_get_d(r23593, MPFR_RNDN); } static mpfr_t r23594, r23595, r23596, r23597, r23598, r23599, r23600, r23601, r23602, r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617; void setup_mpfr_f_fm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r23594, "1", 10, MPFR_RNDN); mpfr_init(r23595); mpfr_init(r23596); mpfr_init(r23597); mpfr_init_set_str(r23598, "2", 10, MPFR_RNDN); mpfr_init(r23599); mpfr_init(r23600); mpfr_init(r23601); mpfr_init(r23602); mpfr_init(r23603); mpfr_init(r23604); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init(r23611); mpfr_init(r23612); mpfr_init(r23613); mpfr_init(r23614); mpfr_init(r23615); mpfr_init(r23616); mpfr_init(r23617); } double f_fm(double __attribute__((unused)) K, double m, double n, double M, double l) { ; mpfr_set_d(r23595, m, MPFR_RNDN); mpfr_set_d(r23596, n, MPFR_RNDN); mpfr_add(r23597, r23595, r23596, MPFR_RNDN); ; mpfr_div(r23599, r23597, r23598, MPFR_RNDN); mpfr_set_d(r23600, M, MPFR_RNDN); mpfr_sub(r23601, r23599, r23600, MPFR_RNDN); mpfr_pow(r23602, r23601, r23598, MPFR_RNDN); mpfr_cbrt(r23603, r23602, MPFR_RNDN); mpfr_cbrt(r23604, r23603, MPFR_RNDN); mpfr_mul(r23605, r23604, r23604, MPFR_RNDN); mpfr_mul(r23606, r23603, r23605, MPFR_RNDN); mpfr_mul(r23607, r23606, r23604, MPFR_RNDN); mpfr_mul(r23608, r23605, r23604, MPFR_RNDN); mpfr_mul(r23609, r23607, r23608, MPFR_RNDN); mpfr_neg(r23610, r23609, MPFR_RNDN); mpfr_set_d(r23611, l, MPFR_RNDN); mpfr_sub(r23612, r23595, r23596, MPFR_RNDN); mpfr_abs(r23613, r23612, MPFR_RNDN); mpfr_sub(r23614, r23611, r23613, MPFR_RNDN); mpfr_sub(r23615, r23610, r23614, MPFR_RNDN); mpfr_exp(r23616, r23615, MPFR_RNDN); mpfr_mul(r23617, r23594, r23616, MPFR_RNDN); return mpfr_get_d(r23617, MPFR_RNDN); } static mpfr_t r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632, r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641; void setup_mpfr_f_dm() { mpfr_set_default_prec(1616); mpfr_init_set_str(r23618, "1", 10, MPFR_RNDN); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init_set_str(r23622, "2", 10, MPFR_RNDN); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); mpfr_init(r23628); mpfr_init(r23629); mpfr_init(r23630); mpfr_init(r23631); mpfr_init(r23632); mpfr_init(r23633); mpfr_init(r23634); mpfr_init(r23635); mpfr_init(r23636); mpfr_init(r23637); mpfr_init(r23638); mpfr_init(r23639); mpfr_init(r23640); mpfr_init(r23641); } double f_dm(double __attribute__((unused)) K, double m, double n, double M, double l) { ; mpfr_set_d(r23619, m, MPFR_RNDN); mpfr_set_d(r23620, n, MPFR_RNDN); mpfr_add(r23621, r23619, r23620, MPFR_RNDN); ; mpfr_div(r23623, r23621, r23622, MPFR_RNDN); mpfr_set_d(r23624, M, MPFR_RNDN); mpfr_sub(r23625, r23623, r23624, MPFR_RNDN); mpfr_pow(r23626, r23625, r23622, MPFR_RNDN); mpfr_cbrt(r23627, r23626, MPFR_RNDN); mpfr_cbrt(r23628, r23627, MPFR_RNDN); mpfr_mul(r23629, r23628, r23628, MPFR_RNDN); mpfr_mul(r23630, r23627, r23629, MPFR_RNDN); mpfr_mul(r23631, r23630, r23628, MPFR_RNDN); mpfr_mul(r23632, r23629, r23628, MPFR_RNDN); mpfr_mul(r23633, r23631, r23632, MPFR_RNDN); mpfr_neg(r23634, r23633, MPFR_RNDN); mpfr_set_d(r23635, l, MPFR_RNDN); mpfr_sub(r23636, r23619, r23620, MPFR_RNDN); mpfr_abs(r23637, r23636, MPFR_RNDN); mpfr_sub(r23638, r23635, r23637, MPFR_RNDN); mpfr_sub(r23639, r23634, r23638, MPFR_RNDN); mpfr_exp(r23640, r23639, MPFR_RNDN); mpfr_mul(r23641, r23618, r23640, MPFR_RNDN); return mpfr_get_d(r23641, MPFR_RNDN); }