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