#include #include #include #include #include char *name = "Migdal et al, Equation (51)"; double f_if(float k, float n) { float r23339 = 1; float r23340 = k; float r23341 = sqrt(r23340); float r23342 = r23339 / r23341; float r23343 = 2; float r23344 = atan2(1.0, 0.0); float r23345 = r23343 * r23344; float r23346 = n; float r23347 = r23345 * r23346; float r23348 = r23339 - r23340; float r23349 = r23348 / r23343; float r23350 = pow(r23347, r23349); float r23351 = r23342 * r23350; return r23351; } double f_id(double k, double n) { double r23352 = 1; double r23353 = k; double r23354 = sqrt(r23353); double r23355 = r23352 / r23354; double r23356 = 2; double r23357 = atan2(1.0, 0.0); double r23358 = r23356 * r23357; double r23359 = n; double r23360 = r23358 * r23359; double r23361 = r23352 - r23353; double r23362 = r23361 / r23356; double r23363 = pow(r23360, r23362); double r23364 = r23355 * r23363; return r23364; } double f_of(float k, float n) { float r23365 = k; float r23366 = 1/2; float r23367 = -r23366; float r23368 = pow(r23365, r23367); float r23369 = 2; float r23370 = atan2(1.0, 0.0); float r23371 = r23369 * r23370; float r23372 = n; float r23373 = r23371 * r23372; float r23374 = 1; float r23375 = r23374 / r23369; float r23376 = pow(r23373, r23375); float r23377 = r23365 / r23369; float r23378 = pow(r23373, r23377); float r23379 = r23376 / r23378; float r23380 = r23368 * r23379; return r23380; } double f_od(double k, double n) { double r23381 = k; double r23382 = 1/2; double r23383 = -r23382; double r23384 = pow(r23381, r23383); double r23385 = 2; double r23386 = atan2(1.0, 0.0); double r23387 = r23385 * r23386; double r23388 = n; double r23389 = r23387 * r23388; double r23390 = 1; double r23391 = r23390 / r23385; double r23392 = pow(r23389, r23391); double r23393 = r23381 / r23385; double r23394 = pow(r23389, r23393); double r23395 = r23392 / r23394; double r23396 = r23384 * r23395; return r23396; } 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 r23397, r23398, r23399, r23400, r23401, r23402, r23403, r23404, r23405, r23406, r23407, r23408, r23409; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r23397, "1", 10, MPFR_RNDN); mpfr_init(r23398); mpfr_init(r23399); mpfr_init(r23400); mpfr_init_set_str(r23401, "2", 10, MPFR_RNDN); mpfr_init(r23402); mpfr_init(r23403); mpfr_init(r23404); mpfr_init(r23405); mpfr_init(r23406); mpfr_init(r23407); mpfr_init(r23408); mpfr_init(r23409); } double f_im(double k, double n) { ; mpfr_set_d(r23398, k, MPFR_RNDN); mpfr_sqrt(r23399, r23398, MPFR_RNDN); mpfr_div(r23400, r23397, r23399, MPFR_RNDN); ; mpfr_const_pi(r23402, MPFR_RNDN); mpfr_mul(r23403, r23401, r23402, MPFR_RNDN); mpfr_set_d(r23404, n, MPFR_RNDN); mpfr_mul(r23405, r23403, r23404, MPFR_RNDN); mpfr_sub(r23406, r23397, r23398, MPFR_RNDN); mpfr_div(r23407, r23406, r23401, MPFR_RNDN); mpfr_pow(r23408, r23405, r23407, MPFR_RNDN); mpfr_mul(r23409, r23400, r23408, MPFR_RNDN); return mpfr_get_d(r23409, MPFR_RNDN); } static mpfr_t r23410, r23411, r23412, r23413, r23414, r23415, r23416, r23417, r23418, r23419, r23420, r23421, r23422, r23423, r23424, r23425; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r23410); mpfr_init_set_str(r23411, "1/2", 10, MPFR_RNDN); mpfr_init(r23412); mpfr_init(r23413); mpfr_init_set_str(r23414, "2", 10, MPFR_RNDN); mpfr_init(r23415); mpfr_init(r23416); mpfr_init(r23417); mpfr_init(r23418); mpfr_init_set_str(r23419, "1", 10, MPFR_RNDN); mpfr_init(r23420); mpfr_init(r23421); mpfr_init(r23422); mpfr_init(r23423); mpfr_init(r23424); mpfr_init(r23425); } double f_fm(double k, double n) { mpfr_set_d(r23410, k, MPFR_RNDN); ; mpfr_neg(r23412, r23411, MPFR_RNDN); mpfr_pow(r23413, r23410, r23412, MPFR_RNDN); ; mpfr_const_pi(r23415, MPFR_RNDN); mpfr_mul(r23416, r23414, r23415, MPFR_RNDN); mpfr_set_d(r23417, n, MPFR_RNDN); mpfr_mul(r23418, r23416, r23417, MPFR_RNDN); ; mpfr_div(r23420, r23419, r23414, MPFR_RNDN); mpfr_pow(r23421, r23418, r23420, MPFR_RNDN); mpfr_div(r23422, r23410, r23414, MPFR_RNDN); mpfr_pow(r23423, r23418, r23422, MPFR_RNDN); mpfr_div(r23424, r23421, r23423, MPFR_RNDN); mpfr_mul(r23425, r23413, r23424, MPFR_RNDN); return mpfr_get_d(r23425, MPFR_RNDN); } static mpfr_t r23426, r23427, r23428, r23429, r23430, r23431, r23432, r23433, r23434, r23435, r23436, r23437, r23438, r23439, r23440, r23441; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r23426); mpfr_init_set_str(r23427, "1/2", 10, MPFR_RNDN); mpfr_init(r23428); mpfr_init(r23429); mpfr_init_set_str(r23430, "2", 10, MPFR_RNDN); mpfr_init(r23431); mpfr_init(r23432); mpfr_init(r23433); mpfr_init(r23434); mpfr_init_set_str(r23435, "1", 10, MPFR_RNDN); mpfr_init(r23436); mpfr_init(r23437); mpfr_init(r23438); mpfr_init(r23439); mpfr_init(r23440); mpfr_init(r23441); } double f_dm(double k, double n) { mpfr_set_d(r23426, k, MPFR_RNDN); ; mpfr_neg(r23428, r23427, MPFR_RNDN); mpfr_pow(r23429, r23426, r23428, MPFR_RNDN); ; mpfr_const_pi(r23431, MPFR_RNDN); mpfr_mul(r23432, r23430, r23431, MPFR_RNDN); mpfr_set_d(r23433, n, MPFR_RNDN); mpfr_mul(r23434, r23432, r23433, MPFR_RNDN); ; mpfr_div(r23436, r23435, r23430, MPFR_RNDN); mpfr_pow(r23437, r23434, r23436, MPFR_RNDN); mpfr_div(r23438, r23426, r23430, MPFR_RNDN); mpfr_pow(r23439, r23434, r23438, MPFR_RNDN); mpfr_div(r23440, r23437, r23439, MPFR_RNDN); mpfr_mul(r23441, r23429, r23440, MPFR_RNDN); return mpfr_get_d(r23441, MPFR_RNDN); }