#include #include #include #include #include char *name = "Migdal et al, Equation (51)"; double f_if(float k, float n) { float r9439 = 1; float r9440 = k; float r9441 = sqrt(r9440); float r9442 = r9439 / r9441; float r9443 = 2; float r9444 = atan2(1.0, 0.0); float r9445 = r9443 * r9444; float r9446 = n; float r9447 = r9445 * r9446; float r9448 = r9439 - r9440; float r9449 = r9448 / r9443; float r9450 = pow(r9447, r9449); float r9451 = r9442 * r9450; return r9451; } double f_id(double k, double n) { double r9452 = 1; double r9453 = k; double r9454 = sqrt(r9453); double r9455 = r9452 / r9454; double r9456 = 2; double r9457 = atan2(1.0, 0.0); double r9458 = r9456 * r9457; double r9459 = n; double r9460 = r9458 * r9459; double r9461 = r9452 - r9453; double r9462 = r9461 / r9456; double r9463 = pow(r9460, r9462); double r9464 = r9455 * r9463; return r9464; } double f_of(float k, float n) { float r9465 = n; float r9466 = atan2(1.0, 0.0); float r9467 = r9465 * r9466; float r9468 = sqrt(r9467); float r9469 = r9468 * r9468; float r9470 = 2; float r9471 = r9469 * r9470; float r9472 = 1; float r9473 = r9472 / r9470; float r9474 = pow(r9471, r9473); float r9475 = r9467 * r9470; float r9476 = k; float r9477 = r9476 / r9470; float r9478 = pow(r9475, r9477); float r9479 = r9474 / r9478; float r9480 = sqrt(r9476); float r9481 = r9479 / r9480; return r9481; } double f_od(double k, double n) { double r9482 = n; double r9483 = atan2(1.0, 0.0); double r9484 = r9482 * r9483; double r9485 = sqrt(r9484); double r9486 = r9485 * r9485; double r9487 = 2; double r9488 = r9486 * r9487; double r9489 = 1; double r9490 = r9489 / r9487; double r9491 = pow(r9488, r9490); double r9492 = r9484 * r9487; double r9493 = k; double r9494 = r9493 / r9487; double r9495 = pow(r9492, r9494); double r9496 = r9491 / r9495; double r9497 = sqrt(r9493); double r9498 = r9496 / r9497; return r9498; } 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 r9499, r9500, r9501, r9502, r9503, r9504, r9505, r9506, r9507, r9508, r9509, r9510, r9511; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init_set_str(r9499, "1", 10, MPFR_RNDN); mpfr_init(r9500); mpfr_init(r9501); mpfr_init(r9502); mpfr_init_set_str(r9503, "2", 10, MPFR_RNDN); mpfr_init(r9504); mpfr_init(r9505); mpfr_init(r9506); mpfr_init(r9507); mpfr_init(r9508); mpfr_init(r9509); mpfr_init(r9510); mpfr_init(r9511); } double f_im(double k, double n) { ; mpfr_set_d(r9500, k, MPFR_RNDN); mpfr_sqrt(r9501, r9500, MPFR_RNDN); mpfr_div(r9502, r9499, r9501, MPFR_RNDN); ; mpfr_const_pi(r9504, MPFR_RNDN); mpfr_mul(r9505, r9503, r9504, MPFR_RNDN); mpfr_set_d(r9506, n, MPFR_RNDN); mpfr_mul(r9507, r9505, r9506, MPFR_RNDN); mpfr_sub(r9508, r9499, r9500, MPFR_RNDN); mpfr_div(r9509, r9508, r9503, MPFR_RNDN); mpfr_pow(r9510, r9507, r9509, MPFR_RNDN); mpfr_mul(r9511, r9502, r9510, MPFR_RNDN); return mpfr_get_d(r9511, MPFR_RNDN); } static mpfr_t r9512, r9513, r9514, r9515, r9516, r9517, r9518, r9519, r9520, r9521, r9522, r9523, r9524, r9525, r9526, r9527, r9528; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r9512); mpfr_init(r9513); mpfr_init(r9514); mpfr_init(r9515); mpfr_init(r9516); mpfr_init_set_str(r9517, "2", 10, MPFR_RNDN); mpfr_init(r9518); mpfr_init_set_str(r9519, "1", 10, MPFR_RNDN); mpfr_init(r9520); mpfr_init(r9521); mpfr_init(r9522); mpfr_init(r9523); mpfr_init(r9524); mpfr_init(r9525); mpfr_init(r9526); mpfr_init(r9527); mpfr_init(r9528); } double f_fm(double k, double n) { mpfr_set_d(r9512, n, MPFR_RNDN); mpfr_const_pi(r9513, MPFR_RNDN); mpfr_mul(r9514, r9512, r9513, MPFR_RNDN); mpfr_sqrt(r9515, r9514, MPFR_RNDN); mpfr_mul(r9516, r9515, r9515, MPFR_RNDN); ; mpfr_mul(r9518, r9516, r9517, MPFR_RNDN); ; mpfr_div(r9520, r9519, r9517, MPFR_RNDN); mpfr_pow(r9521, r9518, r9520, MPFR_RNDN); mpfr_mul(r9522, r9514, r9517, MPFR_RNDN); mpfr_set_d(r9523, k, MPFR_RNDN); mpfr_div(r9524, r9523, r9517, MPFR_RNDN); mpfr_pow(r9525, r9522, r9524, MPFR_RNDN); mpfr_div(r9526, r9521, r9525, MPFR_RNDN); mpfr_sqrt(r9527, r9523, MPFR_RNDN); mpfr_div(r9528, r9526, r9527, MPFR_RNDN); return mpfr_get_d(r9528, MPFR_RNDN); } static mpfr_t r9529, r9530, r9531, r9532, r9533, r9534, r9535, r9536, r9537, r9538, r9539, r9540, r9541, r9542, r9543, r9544, r9545; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r9529); mpfr_init(r9530); mpfr_init(r9531); mpfr_init(r9532); mpfr_init(r9533); mpfr_init_set_str(r9534, "2", 10, MPFR_RNDN); mpfr_init(r9535); mpfr_init_set_str(r9536, "1", 10, MPFR_RNDN); mpfr_init(r9537); mpfr_init(r9538); mpfr_init(r9539); mpfr_init(r9540); mpfr_init(r9541); mpfr_init(r9542); mpfr_init(r9543); mpfr_init(r9544); mpfr_init(r9545); } double f_dm(double k, double n) { mpfr_set_d(r9529, n, MPFR_RNDN); mpfr_const_pi(r9530, MPFR_RNDN); mpfr_mul(r9531, r9529, r9530, MPFR_RNDN); mpfr_sqrt(r9532, r9531, MPFR_RNDN); mpfr_mul(r9533, r9532, r9532, MPFR_RNDN); ; mpfr_mul(r9535, r9533, r9534, MPFR_RNDN); ; mpfr_div(r9537, r9536, r9534, MPFR_RNDN); mpfr_pow(r9538, r9535, r9537, MPFR_RNDN); mpfr_mul(r9539, r9531, r9534, MPFR_RNDN); mpfr_set_d(r9540, k, MPFR_RNDN); mpfr_div(r9541, r9540, r9534, MPFR_RNDN); mpfr_pow(r9542, r9539, r9541, MPFR_RNDN); mpfr_div(r9543, r9538, r9542, MPFR_RNDN); mpfr_sqrt(r9544, r9540, MPFR_RNDN); mpfr_div(r9545, r9543, r9544, MPFR_RNDN); return mpfr_get_d(r9545, MPFR_RNDN); }