#include #include #include #include #include char *name = "Migdal et al, Equation (64)"; double f_if(float a1, float a2, float th) { float r25398 = th; float r25399 = cos(r25398); float r25400 = 2; float r25401 = sqrt(r25400); float r25402 = r25399 / r25401; float r25403 = a1; float r25404 = r25403 * r25403; float r25405 = r25402 * r25404; float r25406 = a2; float r25407 = r25406 * r25406; float r25408 = r25402 * r25407; float r25409 = r25405 + r25408; return r25409; } double f_id(double a1, double a2, double th) { double r25410 = th; double r25411 = cos(r25410); double r25412 = 2; double r25413 = sqrt(r25412); double r25414 = r25411 / r25413; double r25415 = a1; double r25416 = r25415 * r25415; double r25417 = r25414 * r25416; double r25418 = a2; double r25419 = r25418 * r25418; double r25420 = r25414 * r25419; double r25421 = r25417 + r25420; return r25421; } double f_of(float a1, float a2, float th) { float r25422 = 1; float r25423 = 2; float r25424 = sqrt(r25423); float r25425 = sqrt(r25424); float r25426 = sqrt(r25425); float r25427 = r25422 / r25426; float r25428 = r25427 / r25426; float r25429 = a1; float r25430 = a2; float r25431 = r25430 * r25430; float r25432 = fma(r25429, r25429, r25431); float r25433 = r25432 / r25425; float r25434 = th; float r25435 = cos(r25434); float r25436 = r25433 * r25435; float r25437 = r25428 * r25436; return r25437; } double f_od(double a1, double a2, double th) { double r25438 = 1; double r25439 = 2; double r25440 = sqrt(r25439); double r25441 = sqrt(r25440); double r25442 = sqrt(r25441); double r25443 = r25438 / r25442; double r25444 = r25443 / r25442; double r25445 = a1; double r25446 = a2; double r25447 = r25446 * r25446; double r25448 = fma(r25445, r25445, r25447); double r25449 = r25448 / r25441; double r25450 = th; double r25451 = cos(r25450); double r25452 = r25449 * r25451; double r25453 = r25444 * r25452; return r25453; } 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 r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463, r25464, r25465; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25454); mpfr_init(r25455); mpfr_init_set_str(r25456, "2", 10, MPFR_RNDN); mpfr_init(r25457); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); mpfr_init(r25464); mpfr_init(r25465); } double f_im(double a1, double a2, double th) { mpfr_set_d(r25454, th, MPFR_RNDN); mpfr_cos(r25455, r25454, MPFR_RNDN); ; mpfr_sqrt(r25457, r25456, MPFR_RNDN); mpfr_div(r25458, r25455, r25457, MPFR_RNDN); mpfr_set_d(r25459, a1, MPFR_RNDN); mpfr_mul(r25460, r25459, r25459, MPFR_RNDN); mpfr_mul(r25461, r25458, r25460, MPFR_RNDN); mpfr_set_d(r25462, a2, MPFR_RNDN); mpfr_mul(r25463, r25462, r25462, MPFR_RNDN); mpfr_mul(r25464, r25458, r25463, MPFR_RNDN); mpfr_add(r25465, r25461, r25464, MPFR_RNDN); return mpfr_get_d(r25465, MPFR_RNDN); } static mpfr_t r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25466, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25467, "2", 10, MPFR_RNDN); mpfr_init(r25468); mpfr_init(r25469); mpfr_init(r25470); mpfr_init(r25471); mpfr_init(r25472); mpfr_init(r25473); mpfr_init(r25474); mpfr_init(r25475); mpfr_init(r25476); mpfr_init(r25477); mpfr_init(r25478); mpfr_init(r25479); mpfr_init(r25480); mpfr_init(r25481); } double f_fm(double a1, double a2, double th) { ; ; mpfr_sqrt(r25468, r25467, MPFR_RNDN); mpfr_sqrt(r25469, r25468, MPFR_RNDN); mpfr_sqrt(r25470, r25469, MPFR_RNDN); mpfr_div(r25471, r25466, r25470, MPFR_RNDN); mpfr_div(r25472, r25471, r25470, MPFR_RNDN); mpfr_set_d(r25473, a1, MPFR_RNDN); mpfr_set_d(r25474, a2, MPFR_RNDN); mpfr_mul(r25475, r25474, r25474, MPFR_RNDN); mpfr_fma(r25476, r25473, r25473, r25475, MPFR_RNDN); mpfr_div(r25477, r25476, r25469, MPFR_RNDN); mpfr_set_d(r25478, th, MPFR_RNDN); mpfr_cos(r25479, r25478, MPFR_RNDN); mpfr_mul(r25480, r25477, r25479, MPFR_RNDN); mpfr_mul(r25481, r25472, r25480, MPFR_RNDN); return mpfr_get_d(r25481, MPFR_RNDN); } static mpfr_t r25482, r25483, r25484, r25485, r25486, r25487, r25488, r25489, r25490, r25491, r25492, r25493, r25494, r25495, r25496, r25497; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r25482, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25483, "2", 10, MPFR_RNDN); mpfr_init(r25484); mpfr_init(r25485); mpfr_init(r25486); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init(r25490); mpfr_init(r25491); mpfr_init(r25492); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); mpfr_init(r25497); } double f_dm(double a1, double a2, double th) { ; ; mpfr_sqrt(r25484, r25483, MPFR_RNDN); mpfr_sqrt(r25485, r25484, MPFR_RNDN); mpfr_sqrt(r25486, r25485, MPFR_RNDN); mpfr_div(r25487, r25482, r25486, MPFR_RNDN); mpfr_div(r25488, r25487, r25486, MPFR_RNDN); mpfr_set_d(r25489, a1, MPFR_RNDN); mpfr_set_d(r25490, a2, MPFR_RNDN); mpfr_mul(r25491, r25490, r25490, MPFR_RNDN); mpfr_fma(r25492, r25489, r25489, r25491, MPFR_RNDN); mpfr_div(r25493, r25492, r25485, MPFR_RNDN); mpfr_set_d(r25494, th, MPFR_RNDN); mpfr_cos(r25495, r25494, MPFR_RNDN); mpfr_mul(r25496, r25493, r25495, MPFR_RNDN); mpfr_mul(r25497, r25488, r25496, MPFR_RNDN); return mpfr_get_d(r25497, MPFR_RNDN); }