#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r20339 = 1; float r20340 = 2; float r20341 = r20339 / r20340; float r20342 = l; float r20343 = r20340 * r20342; float r20344 = Om; float r20345 = r20343 / r20344; float r20346 = pow(r20345, r20340); float r20347 = kx; float r20348 = sin(r20347); float r20349 = pow(r20348, r20340); float r20350 = ky; float r20351 = sin(r20350); float r20352 = pow(r20351, r20340); float r20353 = r20349 + r20352; float r20354 = r20346 * r20353; float r20355 = r20339 + r20354; float r20356 = sqrt(r20355); float r20357 = r20339 / r20356; float r20358 = r20339 + r20357; float r20359 = r20341 * r20358; float r20360 = sqrt(r20359); return r20360; } double f_id(double l, double Om, double kx, double ky) { double r20361 = 1; double r20362 = 2; double r20363 = r20361 / r20362; double r20364 = l; double r20365 = r20362 * r20364; double r20366 = Om; double r20367 = r20365 / r20366; double r20368 = pow(r20367, r20362); double r20369 = kx; double r20370 = sin(r20369); double r20371 = pow(r20370, r20362); double r20372 = ky; double r20373 = sin(r20372); double r20374 = pow(r20373, r20362); double r20375 = r20371 + r20374; double r20376 = r20368 * r20375; double r20377 = r20361 + r20376; double r20378 = sqrt(r20377); double r20379 = r20361 / r20378; double r20380 = r20361 + r20379; double r20381 = r20363 * r20380; double r20382 = sqrt(r20381); return r20382; } double f_of(float l, float Om, float kx, float ky) { float r20383 = 1; float r20384 = 2; float r20385 = r20383 / r20384; float r20386 = l; float r20387 = r20384 * r20386; float r20388 = Om; float r20389 = r20387 / r20388; float r20390 = kx; float r20391 = sin(r20390); float r20392 = pow(r20391, r20384); float r20393 = ky; float r20394 = sin(r20393); float r20395 = pow(r20394, r20384); float r20396 = r20392 + r20395; float r20397 = r20389 * r20396; float r20398 = r20389 * r20397; float r20399 = r20383 + r20398; float r20400 = sqrt(r20399); float r20401 = r20383 / r20400; float r20402 = r20383 + r20401; float r20403 = r20385 * r20402; float r20404 = sqrt(r20403); return r20404; } double f_od(double l, double Om, double kx, double ky) { double r20405 = 1; double r20406 = 2; double r20407 = r20405 / r20406; double r20408 = l; double r20409 = r20406 * r20408; double r20410 = Om; double r20411 = r20409 / r20410; double r20412 = kx; double r20413 = sin(r20412); double r20414 = pow(r20413, r20406); double r20415 = ky; double r20416 = sin(r20415); double r20417 = pow(r20416, r20406); double r20418 = r20414 + r20417; double r20419 = r20411 * r20418; double r20420 = r20411 * r20419; double r20421 = r20405 + r20420; double r20422 = sqrt(r20421); double r20423 = r20405 / r20422; double r20424 = r20405 + r20423; double r20425 = r20407 * r20424; double r20426 = sqrt(r20425); return r20426; } 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 r20427, r20428, r20429, r20430, r20431, r20432, r20433, r20434, r20435, r20436, r20437, r20438, r20439, r20440, r20441, r20442, r20443, r20444, r20445, r20446, r20447, r20448; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20427, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20428, "2", 10, MPFR_RNDN); mpfr_init(r20429); mpfr_init(r20430); mpfr_init(r20431); mpfr_init(r20432); mpfr_init(r20433); mpfr_init(r20434); mpfr_init(r20435); mpfr_init(r20436); mpfr_init(r20437); mpfr_init(r20438); mpfr_init(r20439); mpfr_init(r20440); mpfr_init(r20441); mpfr_init(r20442); mpfr_init(r20443); mpfr_init(r20444); mpfr_init(r20445); mpfr_init(r20446); mpfr_init(r20447); mpfr_init(r20448); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r20429, r20427, r20428, MPFR_RNDN); mpfr_set_d(r20430, l, MPFR_RNDN); mpfr_mul(r20431, r20428, r20430, MPFR_RNDN); mpfr_set_d(r20432, Om, MPFR_RNDN); mpfr_div(r20433, r20431, r20432, MPFR_RNDN); mpfr_pow(r20434, r20433, r20428, MPFR_RNDN); mpfr_set_d(r20435, kx, MPFR_RNDN); mpfr_sin(r20436, r20435, MPFR_RNDN); mpfr_pow(r20437, r20436, r20428, MPFR_RNDN); mpfr_set_d(r20438, ky, MPFR_RNDN); mpfr_sin(r20439, r20438, MPFR_RNDN); mpfr_pow(r20440, r20439, r20428, MPFR_RNDN); mpfr_add(r20441, r20437, r20440, MPFR_RNDN); mpfr_mul(r20442, r20434, r20441, MPFR_RNDN); mpfr_add(r20443, r20427, r20442, MPFR_RNDN); mpfr_sqrt(r20444, r20443, MPFR_RNDN); mpfr_div(r20445, r20427, r20444, MPFR_RNDN); mpfr_add(r20446, r20427, r20445, MPFR_RNDN); mpfr_mul(r20447, r20429, r20446, MPFR_RNDN); mpfr_sqrt(r20448, r20447, MPFR_RNDN); return mpfr_get_d(r20448, MPFR_RNDN); } static mpfr_t r20449, r20450, r20451, r20452, r20453, r20454, r20455, r20456, r20457, r20458, r20459, r20460, r20461, r20462, r20463, r20464, r20465, r20466, r20467, r20468, r20469, r20470; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20449, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20450, "2", 10, MPFR_RNDN); mpfr_init(r20451); mpfr_init(r20452); mpfr_init(r20453); mpfr_init(r20454); mpfr_init(r20455); mpfr_init(r20456); mpfr_init(r20457); mpfr_init(r20458); mpfr_init(r20459); mpfr_init(r20460); mpfr_init(r20461); mpfr_init(r20462); mpfr_init(r20463); mpfr_init(r20464); mpfr_init(r20465); mpfr_init(r20466); mpfr_init(r20467); mpfr_init(r20468); mpfr_init(r20469); mpfr_init(r20470); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r20451, r20449, r20450, MPFR_RNDN); mpfr_set_d(r20452, l, MPFR_RNDN); mpfr_mul(r20453, r20450, r20452, MPFR_RNDN); mpfr_set_d(r20454, Om, MPFR_RNDN); mpfr_div(r20455, r20453, r20454, MPFR_RNDN); mpfr_set_d(r20456, kx, MPFR_RNDN); mpfr_sin(r20457, r20456, MPFR_RNDN); mpfr_pow(r20458, r20457, r20450, MPFR_RNDN); mpfr_set_d(r20459, ky, MPFR_RNDN); mpfr_sin(r20460, r20459, MPFR_RNDN); mpfr_pow(r20461, r20460, r20450, MPFR_RNDN); mpfr_add(r20462, r20458, r20461, MPFR_RNDN); mpfr_mul(r20463, r20455, r20462, MPFR_RNDN); mpfr_mul(r20464, r20455, r20463, MPFR_RNDN); mpfr_add(r20465, r20449, r20464, MPFR_RNDN); mpfr_sqrt(r20466, r20465, MPFR_RNDN); mpfr_div(r20467, r20449, r20466, MPFR_RNDN); mpfr_add(r20468, r20449, r20467, MPFR_RNDN); mpfr_mul(r20469, r20451, r20468, MPFR_RNDN); mpfr_sqrt(r20470, r20469, MPFR_RNDN); return mpfr_get_d(r20470, MPFR_RNDN); } static mpfr_t r20471, r20472, r20473, r20474, r20475, r20476, r20477, r20478, r20479, r20480, r20481, r20482, r20483, r20484, r20485, r20486, r20487, r20488, r20489, r20490, r20491, r20492; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20471, "1", 10, MPFR_RNDN); mpfr_init_set_str(r20472, "2", 10, MPFR_RNDN); mpfr_init(r20473); mpfr_init(r20474); mpfr_init(r20475); mpfr_init(r20476); mpfr_init(r20477); mpfr_init(r20478); mpfr_init(r20479); mpfr_init(r20480); mpfr_init(r20481); mpfr_init(r20482); mpfr_init(r20483); mpfr_init(r20484); mpfr_init(r20485); mpfr_init(r20486); mpfr_init(r20487); mpfr_init(r20488); mpfr_init(r20489); mpfr_init(r20490); mpfr_init(r20491); mpfr_init(r20492); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r20473, r20471, r20472, MPFR_RNDN); mpfr_set_d(r20474, l, MPFR_RNDN); mpfr_mul(r20475, r20472, r20474, MPFR_RNDN); mpfr_set_d(r20476, Om, MPFR_RNDN); mpfr_div(r20477, r20475, r20476, MPFR_RNDN); mpfr_set_d(r20478, kx, MPFR_RNDN); mpfr_sin(r20479, r20478, MPFR_RNDN); mpfr_pow(r20480, r20479, r20472, MPFR_RNDN); mpfr_set_d(r20481, ky, MPFR_RNDN); mpfr_sin(r20482, r20481, MPFR_RNDN); mpfr_pow(r20483, r20482, r20472, MPFR_RNDN); mpfr_add(r20484, r20480, r20483, MPFR_RNDN); mpfr_mul(r20485, r20477, r20484, MPFR_RNDN); mpfr_mul(r20486, r20477, r20485, MPFR_RNDN); mpfr_add(r20487, r20471, r20486, MPFR_RNDN); mpfr_sqrt(r20488, r20487, MPFR_RNDN); mpfr_div(r20489, r20471, r20488, MPFR_RNDN); mpfr_add(r20490, r20471, r20489, MPFR_RNDN); mpfr_mul(r20491, r20473, r20490, MPFR_RNDN); mpfr_sqrt(r20492, r20491, MPFR_RNDN); return mpfr_get_d(r20492, MPFR_RNDN); }