#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r22365 = 1; float r22366 = 2; float r22367 = r22365 / r22366; float r22368 = l; float r22369 = r22366 * r22368; float r22370 = Om; float r22371 = r22369 / r22370; float r22372 = pow(r22371, r22366); float r22373 = kx; float r22374 = sin(r22373); float r22375 = pow(r22374, r22366); float r22376 = ky; float r22377 = sin(r22376); float r22378 = pow(r22377, r22366); float r22379 = r22375 + r22378; float r22380 = r22372 * r22379; float r22381 = r22365 + r22380; float r22382 = sqrt(r22381); float r22383 = r22365 / r22382; float r22384 = r22365 + r22383; float r22385 = r22367 * r22384; float r22386 = sqrt(r22385); return r22386; } double f_id(double l, double Om, double kx, double ky) { double r22387 = 1; double r22388 = 2; double r22389 = r22387 / r22388; double r22390 = l; double r22391 = r22388 * r22390; double r22392 = Om; double r22393 = r22391 / r22392; double r22394 = pow(r22393, r22388); double r22395 = kx; double r22396 = sin(r22395); double r22397 = pow(r22396, r22388); double r22398 = ky; double r22399 = sin(r22398); double r22400 = pow(r22399, r22388); double r22401 = r22397 + r22400; double r22402 = r22394 * r22401; double r22403 = r22387 + r22402; double r22404 = sqrt(r22403); double r22405 = r22387 / r22404; double r22406 = r22387 + r22405; double r22407 = r22389 * r22406; double r22408 = sqrt(r22407); return r22408; } double f_of(float l, float Om, float kx, float ky) { float r22409 = 1; float r22410 = 2; float r22411 = r22409 / r22410; float r22412 = l; float r22413 = Om; float r22414 = r22413 / r22410; float r22415 = r22412 / r22414; float r22416 = kx; float r22417 = sin(r22416); float r22418 = ky; float r22419 = sin(r22418); float r22420 = hypot(r22417, r22419); float r22421 = r22415 * r22420; float r22422 = hypot(r22409, r22421); float r22423 = expm1(r22422); float r22424 = log1p(r22423); float r22425 = r22409 / r22424; float r22426 = r22409 + r22425; float r22427 = r22411 * r22426; float r22428 = sqrt(r22427); return r22428; } double f_od(double l, double Om, double kx, double ky) { double r22429 = 1; double r22430 = 2; double r22431 = r22429 / r22430; double r22432 = l; double r22433 = Om; double r22434 = r22433 / r22430; double r22435 = r22432 / r22434; double r22436 = kx; double r22437 = sin(r22436); double r22438 = ky; double r22439 = sin(r22438); double r22440 = hypot(r22437, r22439); double r22441 = r22435 * r22440; double r22442 = hypot(r22429, r22441); double r22443 = expm1(r22442); double r22444 = log1p(r22443); double r22445 = r22429 / r22444; double r22446 = r22429 + r22445; double r22447 = r22431 * r22446; double r22448 = sqrt(r22447); return r22448; } 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 r22449, r22450, r22451, r22452, r22453, r22454, r22455, r22456, r22457, r22458, r22459, r22460, r22461, r22462, r22463, r22464, r22465, r22466, r22467, r22468, r22469, r22470; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22449, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22450, "2", 10, MPFR_RNDN); mpfr_init(r22451); mpfr_init(r22452); mpfr_init(r22453); mpfr_init(r22454); mpfr_init(r22455); mpfr_init(r22456); mpfr_init(r22457); mpfr_init(r22458); mpfr_init(r22459); mpfr_init(r22460); mpfr_init(r22461); mpfr_init(r22462); mpfr_init(r22463); mpfr_init(r22464); mpfr_init(r22465); mpfr_init(r22466); mpfr_init(r22467); mpfr_init(r22468); mpfr_init(r22469); mpfr_init(r22470); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22451, r22449, r22450, MPFR_RNDN); mpfr_set_d(r22452, l, MPFR_RNDN); mpfr_mul(r22453, r22450, r22452, MPFR_RNDN); mpfr_set_d(r22454, Om, MPFR_RNDN); mpfr_div(r22455, r22453, r22454, MPFR_RNDN); mpfr_pow(r22456, r22455, r22450, MPFR_RNDN); mpfr_set_d(r22457, kx, MPFR_RNDN); mpfr_sin(r22458, r22457, MPFR_RNDN); mpfr_pow(r22459, r22458, r22450, MPFR_RNDN); mpfr_set_d(r22460, ky, MPFR_RNDN); mpfr_sin(r22461, r22460, MPFR_RNDN); mpfr_pow(r22462, r22461, r22450, MPFR_RNDN); mpfr_add(r22463, r22459, r22462, MPFR_RNDN); mpfr_mul(r22464, r22456, r22463, MPFR_RNDN); mpfr_add(r22465, r22449, r22464, MPFR_RNDN); mpfr_sqrt(r22466, r22465, MPFR_RNDN); mpfr_div(r22467, r22449, r22466, MPFR_RNDN); mpfr_add(r22468, r22449, r22467, MPFR_RNDN); mpfr_mul(r22469, r22451, r22468, MPFR_RNDN); mpfr_sqrt(r22470, r22469, MPFR_RNDN); return mpfr_get_d(r22470, MPFR_RNDN); } static mpfr_t r22471, r22472, r22473, r22474, r22475, r22476, r22477, r22478, r22479, r22480, r22481, r22482, r22483, r22484, r22485, r22486, r22487, r22488, r22489, r22490; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22471, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22472, "2", 10, MPFR_RNDN); mpfr_init(r22473); mpfr_init(r22474); mpfr_init(r22475); mpfr_init(r22476); mpfr_init(r22477); mpfr_init(r22478); mpfr_init(r22479); mpfr_init(r22480); mpfr_init(r22481); mpfr_init(r22482); mpfr_init(r22483); mpfr_init(r22484); mpfr_init(r22485); mpfr_init(r22486); mpfr_init(r22487); mpfr_init(r22488); mpfr_init(r22489); mpfr_init(r22490); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22473, r22471, r22472, MPFR_RNDN); mpfr_set_d(r22474, l, MPFR_RNDN); mpfr_set_d(r22475, Om, MPFR_RNDN); mpfr_div(r22476, r22475, r22472, MPFR_RNDN); mpfr_div(r22477, r22474, r22476, MPFR_RNDN); mpfr_set_d(r22478, kx, MPFR_RNDN); mpfr_sin(r22479, r22478, MPFR_RNDN); mpfr_set_d(r22480, ky, MPFR_RNDN); mpfr_sin(r22481, r22480, MPFR_RNDN); mpfr_hypot(r22482, r22479, r22481, MPFR_RNDN); mpfr_mul(r22483, r22477, r22482, MPFR_RNDN); mpfr_hypot(r22484, r22471, r22483, MPFR_RNDN); mpfr_expm1(r22485, r22484, MPFR_RNDN); mpfr_log1p(r22486, r22485, MPFR_RNDN); mpfr_div(r22487, r22471, r22486, MPFR_RNDN); mpfr_add(r22488, r22471, r22487, MPFR_RNDN); mpfr_mul(r22489, r22473, r22488, MPFR_RNDN); mpfr_sqrt(r22490, r22489, MPFR_RNDN); return mpfr_get_d(r22490, MPFR_RNDN); } static mpfr_t r22491, r22492, r22493, r22494, r22495, r22496, r22497, r22498, r22499, r22500, r22501, r22502, r22503, r22504, r22505, r22506, r22507, r22508, r22509, r22510; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22491, "1", 10, MPFR_RNDN); mpfr_init_set_str(r22492, "2", 10, MPFR_RNDN); mpfr_init(r22493); mpfr_init(r22494); mpfr_init(r22495); mpfr_init(r22496); mpfr_init(r22497); mpfr_init(r22498); mpfr_init(r22499); mpfr_init(r22500); mpfr_init(r22501); mpfr_init(r22502); mpfr_init(r22503); mpfr_init(r22504); mpfr_init(r22505); mpfr_init(r22506); mpfr_init(r22507); mpfr_init(r22508); mpfr_init(r22509); mpfr_init(r22510); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r22493, r22491, r22492, MPFR_RNDN); mpfr_set_d(r22494, l, MPFR_RNDN); mpfr_set_d(r22495, Om, MPFR_RNDN); mpfr_div(r22496, r22495, r22492, MPFR_RNDN); mpfr_div(r22497, r22494, r22496, MPFR_RNDN); mpfr_set_d(r22498, kx, MPFR_RNDN); mpfr_sin(r22499, r22498, MPFR_RNDN); mpfr_set_d(r22500, ky, MPFR_RNDN); mpfr_sin(r22501, r22500, MPFR_RNDN); mpfr_hypot(r22502, r22499, r22501, MPFR_RNDN); mpfr_mul(r22503, r22497, r22502, MPFR_RNDN); mpfr_hypot(r22504, r22491, r22503, MPFR_RNDN); mpfr_expm1(r22505, r22504, MPFR_RNDN); mpfr_log1p(r22506, r22505, MPFR_RNDN); mpfr_div(r22507, r22491, r22506, MPFR_RNDN); mpfr_add(r22508, r22491, r22507, MPFR_RNDN); mpfr_mul(r22509, r22493, r22508, MPFR_RNDN); mpfr_sqrt(r22510, r22509, MPFR_RNDN); return mpfr_get_d(r22510, MPFR_RNDN); }