#include #include #include #include #include char *name = "Migdal et al, Equation (51)"; double f_if(float k, float n) { float r22536 = 1; float r22537 = k; float r22538 = sqrt(r22537); float r22539 = r22536 / r22538; float r22540 = 2; float r22541 = atan2(1.0, 0.0); float r22542 = r22540 * r22541; float r22543 = n; float r22544 = r22542 * r22543; float r22545 = r22536 - r22537; float r22546 = r22545 / r22540; float r22547 = pow(r22544, r22546); float r22548 = r22539 * r22547; return r22548; } double f_id(double k, double n) { double r22549 = 1; double r22550 = k; double r22551 = sqrt(r22550); double r22552 = r22549 / r22551; double r22553 = 2; double r22554 = atan2(1.0, 0.0); double r22555 = r22553 * r22554; double r22556 = n; double r22557 = r22555 * r22556; double r22558 = r22549 - r22550; double r22559 = r22558 / r22553; double r22560 = pow(r22557, r22559); double r22561 = r22552 * r22560; return r22561; } double f_of(float k, float n) { float r22562 = n; float r22563 = r22562 + r22562; float r22564 = sqrt(r22563); float r22565 = atan2(1.0, 0.0); float r22566 = r22564 * r22565; float r22567 = r22564 * r22566; float r22568 = 1; float r22569 = 2; float r22570 = r22568 / r22569; float r22571 = pow(r22567, r22570); float r22572 = r22563 * r22565; float r22573 = k; float r22574 = r22573 / r22569; float r22575 = pow(r22572, r22574); float r22576 = r22571 / r22575; float r22577 = sqrt(r22573); float r22578 = r22576 / r22577; return r22578; } double f_od(double k, double n) { double r22579 = n; double r22580 = r22579 + r22579; double r22581 = sqrt(r22580); double r22582 = atan2(1.0, 0.0); double r22583 = r22581 * r22582; double r22584 = r22581 * r22583; double r22585 = 1; double r22586 = 2; double r22587 = r22585 / r22586; double r22588 = pow(r22584, r22587); double r22589 = r22580 * r22582; double r22590 = k; double r22591 = r22590 / r22586; double r22592 = pow(r22589, r22591); double r22593 = r22588 / r22592; double r22594 = sqrt(r22590); double r22595 = r22593 / r22594; return r22595; } 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 r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r22596, "1", 10, MPFR_RNDN); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init_set_str(r22600, "2", 10, MPFR_RNDN); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); mpfr_init(r22606); mpfr_init(r22607); mpfr_init(r22608); } double f_im(double k, double n) { ; mpfr_set_d(r22597, k, MPFR_RNDN); mpfr_sqrt(r22598, r22597, MPFR_RNDN); mpfr_div(r22599, r22596, r22598, MPFR_RNDN); ; mpfr_const_pi(r22601, MPFR_RNDN); mpfr_mul(r22602, r22600, r22601, MPFR_RNDN); mpfr_set_d(r22603, n, MPFR_RNDN); mpfr_mul(r22604, r22602, r22603, MPFR_RNDN); mpfr_sub(r22605, r22596, r22597, MPFR_RNDN); mpfr_div(r22606, r22605, r22600, MPFR_RNDN); mpfr_pow(r22607, r22604, r22606, MPFR_RNDN); mpfr_mul(r22608, r22599, r22607, MPFR_RNDN); return mpfr_get_d(r22608, MPFR_RNDN); } static mpfr_t r22609, r22610, r22611, r22612, r22613, r22614, r22615, r22616, r22617, r22618, r22619, r22620, r22621, r22622, r22623, r22624, r22625; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22609); mpfr_init(r22610); mpfr_init(r22611); mpfr_init(r22612); mpfr_init(r22613); mpfr_init(r22614); mpfr_init_set_str(r22615, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22616, "2", 10, MPFR_RNDN); mpfr_init(r22617); mpfr_init(r22618); mpfr_init(r22619); mpfr_init(r22620); mpfr_init(r22621); mpfr_init(r22622); mpfr_init(r22623); mpfr_init(r22624); mpfr_init(r22625); } double f_fm(double k, double n) { mpfr_set_d(r22609, n, MPFR_RNDN); mpfr_add(r22610, r22609, r22609, MPFR_RNDN); mpfr_sqrt(r22611, r22610, MPFR_RNDN); mpfr_const_pi(r22612, MPFR_RNDN); mpfr_mul(r22613, r22611, r22612, MPFR_RNDN); mpfr_mul(r22614, r22611, r22613, MPFR_RNDN); ; ; mpfr_div(r22617, r22615, r22616, MPFR_RNDN); mpfr_pow(r22618, r22614, r22617, MPFR_RNDN); mpfr_mul(r22619, r22610, r22612, MPFR_RNDN); mpfr_set_d(r22620, k, MPFR_RNDN); mpfr_div(r22621, r22620, r22616, MPFR_RNDN); mpfr_pow(r22622, r22619, r22621, MPFR_RNDN); mpfr_div(r22623, r22618, r22622, MPFR_RNDN); mpfr_sqrt(r22624, r22620, MPFR_RNDN); mpfr_div(r22625, r22623, r22624, MPFR_RNDN); return mpfr_get_d(r22625, MPFR_RNDN); } static mpfr_t r22626, r22627, r22628, r22629, r22630, r22631, r22632, r22633, r22634, r22635, r22636, r22637, r22638, r22639, r22640, r22641, r22642; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22626); mpfr_init(r22627); mpfr_init(r22628); mpfr_init(r22629); mpfr_init(r22630); mpfr_init(r22631); mpfr_init_set_str(r22632, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22633, "2", 10, MPFR_RNDN); mpfr_init(r22634); mpfr_init(r22635); mpfr_init(r22636); mpfr_init(r22637); mpfr_init(r22638); mpfr_init(r22639); mpfr_init(r22640); mpfr_init(r22641); mpfr_init(r22642); } double f_dm(double k, double n) { mpfr_set_d(r22626, n, MPFR_RNDN); mpfr_add(r22627, r22626, r22626, MPFR_RNDN); mpfr_sqrt(r22628, r22627, MPFR_RNDN); mpfr_const_pi(r22629, MPFR_RNDN); mpfr_mul(r22630, r22628, r22629, MPFR_RNDN); mpfr_mul(r22631, r22628, r22630, MPFR_RNDN); ; ; mpfr_div(r22634, r22632, r22633, MPFR_RNDN); mpfr_pow(r22635, r22631, r22634, MPFR_RNDN); mpfr_mul(r22636, r22627, r22629, MPFR_RNDN); mpfr_set_d(r22637, k, MPFR_RNDN); mpfr_div(r22638, r22637, r22633, MPFR_RNDN); mpfr_pow(r22639, r22636, r22638, MPFR_RNDN); mpfr_div(r22640, r22635, r22639, MPFR_RNDN); mpfr_sqrt(r22641, r22637, MPFR_RNDN); mpfr_div(r22642, r22640, r22641, MPFR_RNDN); return mpfr_get_d(r22642, MPFR_RNDN); }