#include #include #include #include #include char *name = "Migdal et al, Equation (64)"; double f_if(float a1, float a2, float th) { float r25354 = th; float r25355 = cos(r25354); float r25356 = 2; float r25357 = sqrt(r25356); float r25358 = r25355 / r25357; float r25359 = a1; float r25360 = r25359 * r25359; float r25361 = r25358 * r25360; float r25362 = a2; float r25363 = r25362 * r25362; float r25364 = r25358 * r25363; float r25365 = r25361 + r25364; return r25365; } double f_id(double a1, double a2, double th) { double r25366 = th; double r25367 = cos(r25366); double r25368 = 2; double r25369 = sqrt(r25368); double r25370 = r25367 / r25369; double r25371 = a1; double r25372 = r25371 * r25371; double r25373 = r25370 * r25372; double r25374 = a2; double r25375 = r25374 * r25374; double r25376 = r25370 * r25375; double r25377 = r25373 + r25376; return r25377; } double f_of(float a1, float a2, float th) { float r25378 = 1; float r25379 = 2; float r25380 = sqrt(r25379); float r25381 = sqrt(r25380); float r25382 = sqrt(r25381); float r25383 = r25378 / r25382; float r25384 = r25383 / r25382; float r25385 = a1; float r25386 = r25385 * r25385; float r25387 = a2; float r25388 = r25387 * r25387; float r25389 = r25386 + r25388; float r25390 = r25389 / r25381; float r25391 = th; float r25392 = cos(r25391); float r25393 = r25390 * r25392; float r25394 = r25384 * r25393; return r25394; } double f_od(double a1, double a2, double th) { double r25395 = 1; double r25396 = 2; double r25397 = sqrt(r25396); double r25398 = sqrt(r25397); double r25399 = sqrt(r25398); double r25400 = r25395 / r25399; double r25401 = r25400 / r25399; double r25402 = a1; double r25403 = r25402 * r25402; double r25404 = a2; double r25405 = r25404 * r25404; double r25406 = r25403 + r25405; double r25407 = r25406 / r25398; double r25408 = th; double r25409 = cos(r25408); double r25410 = r25407 * r25409; double r25411 = r25401 * r25410; return r25411; } 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 r25412, r25413, r25414, r25415, r25416, r25417, r25418, r25419, r25420, r25421, r25422, r25423; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25412); mpfr_init(r25413); mpfr_init_set_str(r25414, "2", 10, MPFR_RNDN); mpfr_init(r25415); mpfr_init(r25416); mpfr_init(r25417); mpfr_init(r25418); mpfr_init(r25419); mpfr_init(r25420); mpfr_init(r25421); mpfr_init(r25422); mpfr_init(r25423); } double f_im(double a1, double a2, double th) { mpfr_set_d(r25412, th, MPFR_RNDN); mpfr_cos(r25413, r25412, MPFR_RNDN); ; mpfr_sqrt(r25415, r25414, MPFR_RNDN); mpfr_div(r25416, r25413, r25415, MPFR_RNDN); mpfr_set_d(r25417, a1, MPFR_RNDN); mpfr_mul(r25418, r25417, r25417, MPFR_RNDN); mpfr_mul(r25419, r25416, r25418, MPFR_RNDN); mpfr_set_d(r25420, a2, MPFR_RNDN); mpfr_mul(r25421, r25420, r25420, MPFR_RNDN); mpfr_mul(r25422, r25416, r25421, MPFR_RNDN); mpfr_add(r25423, r25419, r25422, MPFR_RNDN); return mpfr_get_d(r25423, MPFR_RNDN); } static mpfr_t r25424, r25425, r25426, r25427, r25428, r25429, r25430, r25431, r25432, r25433, r25434, r25435, r25436, r25437, r25438, r25439, r25440; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25424, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25425, "2", 10, MPFR_RNDN); mpfr_init(r25426); mpfr_init(r25427); mpfr_init(r25428); mpfr_init(r25429); mpfr_init(r25430); mpfr_init(r25431); mpfr_init(r25432); mpfr_init(r25433); mpfr_init(r25434); mpfr_init(r25435); mpfr_init(r25436); mpfr_init(r25437); mpfr_init(r25438); mpfr_init(r25439); mpfr_init(r25440); } double f_fm(double a1, double a2, double th) { ; ; mpfr_sqrt(r25426, r25425, MPFR_RNDN); mpfr_sqrt(r25427, r25426, MPFR_RNDN); mpfr_sqrt(r25428, r25427, MPFR_RNDN); mpfr_div(r25429, r25424, r25428, MPFR_RNDN); mpfr_div(r25430, r25429, r25428, MPFR_RNDN); mpfr_set_d(r25431, a1, MPFR_RNDN); mpfr_mul(r25432, r25431, r25431, MPFR_RNDN); mpfr_set_d(r25433, a2, MPFR_RNDN); mpfr_mul(r25434, r25433, r25433, MPFR_RNDN); mpfr_add(r25435, r25432, r25434, MPFR_RNDN); mpfr_div(r25436, r25435, r25427, MPFR_RNDN); mpfr_set_d(r25437, th, MPFR_RNDN); mpfr_cos(r25438, r25437, MPFR_RNDN); mpfr_mul(r25439, r25436, r25438, MPFR_RNDN); mpfr_mul(r25440, r25430, r25439, MPFR_RNDN); return mpfr_get_d(r25440, MPFR_RNDN); } static mpfr_t r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25441, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25442, "2", 10, MPFR_RNDN); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init(r25448); mpfr_init(r25449); mpfr_init(r25450); mpfr_init(r25451); mpfr_init(r25452); mpfr_init(r25453); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init(r25457); } double f_dm(double a1, double a2, double th) { ; ; mpfr_sqrt(r25443, r25442, MPFR_RNDN); mpfr_sqrt(r25444, r25443, MPFR_RNDN); mpfr_sqrt(r25445, r25444, MPFR_RNDN); mpfr_div(r25446, r25441, r25445, MPFR_RNDN); mpfr_div(r25447, r25446, r25445, MPFR_RNDN); mpfr_set_d(r25448, a1, MPFR_RNDN); mpfr_mul(r25449, r25448, r25448, MPFR_RNDN); mpfr_set_d(r25450, a2, MPFR_RNDN); mpfr_mul(r25451, r25450, r25450, MPFR_RNDN); mpfr_add(r25452, r25449, r25451, MPFR_RNDN); mpfr_div(r25453, r25452, r25444, MPFR_RNDN); mpfr_set_d(r25454, th, MPFR_RNDN); mpfr_cos(r25455, r25454, MPFR_RNDN); mpfr_mul(r25456, r25453, r25455, MPFR_RNDN); mpfr_mul(r25457, r25447, r25456, MPFR_RNDN); return mpfr_get_d(r25457, MPFR_RNDN); }