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