#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21350 = 1; float r21351 = 2; float r21352 = r21350 / r21351; float r21353 = l; float r21354 = r21351 * r21353; float r21355 = Om; float r21356 = r21354 / r21355; float r21357 = pow(r21356, r21351); float r21358 = kx; float r21359 = sin(r21358); float r21360 = pow(r21359, r21351); float r21361 = ky; float r21362 = sin(r21361); float r21363 = pow(r21362, r21351); float r21364 = r21360 + r21363; float r21365 = r21357 * r21364; float r21366 = r21350 + r21365; float r21367 = sqrt(r21366); float r21368 = r21350 / r21367; float r21369 = r21350 + r21368; float r21370 = r21352 * r21369; float r21371 = sqrt(r21370); return r21371; } double f_id(double l, double Om, double kx, double ky) { double r21372 = 1; double r21373 = 2; double r21374 = r21372 / r21373; double r21375 = l; double r21376 = r21373 * r21375; double r21377 = Om; double r21378 = r21376 / r21377; double r21379 = pow(r21378, r21373); double r21380 = kx; double r21381 = sin(r21380); double r21382 = pow(r21381, r21373); double r21383 = ky; double r21384 = sin(r21383); double r21385 = pow(r21384, r21373); double r21386 = r21382 + r21385; double r21387 = r21379 * r21386; double r21388 = r21372 + r21387; double r21389 = sqrt(r21388); double r21390 = r21372 / r21389; double r21391 = r21372 + r21390; double r21392 = r21374 * r21391; double r21393 = sqrt(r21392); return r21393; } double f_of(float l, float Om, float kx, float ky) { float r21394 = 1; float r21395 = 2; float r21396 = r21394 / r21395; float r21397 = l; float r21398 = r21395 * r21397; float r21399 = Om; float r21400 = r21398 / r21399; float r21401 = pow(r21400, r21395); float r21402 = kx; float r21403 = sin(r21402); float r21404 = pow(r21403, r21395); float r21405 = ky; float r21406 = sin(r21405); float r21407 = pow(r21406, r21395); float r21408 = r21404 + r21407; float r21409 = r21401 * r21408; float r21410 = r21394 + r21409; float r21411 = sqrt(r21410); float r21412 = r21394 / r21411; float r21413 = r21394 + r21412; float r21414 = r21396 * r21413; float r21415 = sqrt(r21414); return r21415; } double f_od(double l, double Om, double kx, double ky) { double r21416 = 1; double r21417 = 2; double r21418 = r21416 / r21417; double r21419 = l; double r21420 = r21417 * r21419; double r21421 = Om; double r21422 = r21420 / r21421; double r21423 = pow(r21422, r21417); double r21424 = kx; double r21425 = sin(r21424); double r21426 = pow(r21425, r21417); double r21427 = ky; double r21428 = sin(r21427); double r21429 = pow(r21428, r21417); double r21430 = r21426 + r21429; double r21431 = r21423 * r21430; double r21432 = r21416 + r21431; double r21433 = sqrt(r21432); double r21434 = r21416 / r21433; double r21435 = r21416 + r21434; double r21436 = r21418 * r21435; double r21437 = sqrt(r21436); return r21437; } 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 r21438, r21439, r21440, r21441, r21442, r21443, r21444, r21445, r21446, r21447, r21448, r21449, r21450, r21451, r21452, r21453, r21454, r21455, r21456, r21457, r21458, r21459; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21438, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21439, "2", 10, MPFR_RNDN); mpfr_init(r21440); mpfr_init(r21441); mpfr_init(r21442); mpfr_init(r21443); mpfr_init(r21444); mpfr_init(r21445); mpfr_init(r21446); mpfr_init(r21447); mpfr_init(r21448); mpfr_init(r21449); mpfr_init(r21450); mpfr_init(r21451); mpfr_init(r21452); mpfr_init(r21453); mpfr_init(r21454); mpfr_init(r21455); mpfr_init(r21456); mpfr_init(r21457); mpfr_init(r21458); mpfr_init(r21459); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21440, r21438, r21439, MPFR_RNDN); mpfr_set_d(r21441, l, MPFR_RNDN); mpfr_mul(r21442, r21439, r21441, MPFR_RNDN); mpfr_set_d(r21443, Om, MPFR_RNDN); mpfr_div(r21444, r21442, r21443, MPFR_RNDN); mpfr_pow(r21445, r21444, r21439, MPFR_RNDN); mpfr_set_d(r21446, kx, MPFR_RNDN); mpfr_sin(r21447, r21446, MPFR_RNDN); mpfr_pow(r21448, r21447, r21439, MPFR_RNDN); mpfr_set_d(r21449, ky, MPFR_RNDN); mpfr_sin(r21450, r21449, MPFR_RNDN); mpfr_pow(r21451, r21450, r21439, MPFR_RNDN); mpfr_add(r21452, r21448, r21451, MPFR_RNDN); mpfr_mul(r21453, r21445, r21452, MPFR_RNDN); mpfr_add(r21454, r21438, r21453, MPFR_RNDN); mpfr_sqrt(r21455, r21454, MPFR_RNDN); mpfr_div(r21456, r21438, r21455, MPFR_RNDN); mpfr_add(r21457, r21438, r21456, MPFR_RNDN); mpfr_mul(r21458, r21440, r21457, MPFR_RNDN); mpfr_sqrt(r21459, r21458, MPFR_RNDN); return mpfr_get_d(r21459, MPFR_RNDN); } static mpfr_t r21460, r21461, r21462, r21463, r21464, r21465, r21466, r21467, r21468, r21469, r21470, r21471, r21472, r21473, r21474, r21475, r21476, r21477, r21478, r21479, r21480, r21481; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21460, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21461, "2", 10, MPFR_RNDN); mpfr_init(r21462); mpfr_init(r21463); mpfr_init(r21464); mpfr_init(r21465); mpfr_init(r21466); mpfr_init(r21467); mpfr_init(r21468); mpfr_init(r21469); mpfr_init(r21470); mpfr_init(r21471); mpfr_init(r21472); mpfr_init(r21473); mpfr_init(r21474); mpfr_init(r21475); mpfr_init(r21476); mpfr_init(r21477); mpfr_init(r21478); mpfr_init(r21479); mpfr_init(r21480); mpfr_init(r21481); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21462, r21460, r21461, MPFR_RNDN); mpfr_set_d(r21463, l, MPFR_RNDN); mpfr_mul(r21464, r21461, r21463, MPFR_RNDN); mpfr_set_d(r21465, Om, MPFR_RNDN); mpfr_div(r21466, r21464, r21465, MPFR_RNDN); mpfr_pow(r21467, r21466, r21461, MPFR_RNDN); mpfr_set_d(r21468, kx, MPFR_RNDN); mpfr_sin(r21469, r21468, MPFR_RNDN); mpfr_pow(r21470, r21469, r21461, MPFR_RNDN); mpfr_set_d(r21471, ky, MPFR_RNDN); mpfr_sin(r21472, r21471, MPFR_RNDN); mpfr_pow(r21473, r21472, r21461, MPFR_RNDN); mpfr_add(r21474, r21470, r21473, MPFR_RNDN); mpfr_mul(r21475, r21467, r21474, MPFR_RNDN); mpfr_add(r21476, r21460, r21475, MPFR_RNDN); mpfr_sqrt(r21477, r21476, MPFR_RNDN); mpfr_div(r21478, r21460, r21477, MPFR_RNDN); mpfr_add(r21479, r21460, r21478, MPFR_RNDN); mpfr_mul(r21480, r21462, r21479, MPFR_RNDN); mpfr_sqrt(r21481, r21480, MPFR_RNDN); return mpfr_get_d(r21481, MPFR_RNDN); } static mpfr_t r21482, r21483, r21484, r21485, r21486, r21487, r21488, r21489, r21490, r21491, r21492, r21493, r21494, r21495, r21496, r21497, r21498, r21499, r21500, r21501, r21502, r21503; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21482, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21483, "2", 10, MPFR_RNDN); mpfr_init(r21484); mpfr_init(r21485); mpfr_init(r21486); mpfr_init(r21487); mpfr_init(r21488); mpfr_init(r21489); mpfr_init(r21490); mpfr_init(r21491); mpfr_init(r21492); mpfr_init(r21493); mpfr_init(r21494); mpfr_init(r21495); mpfr_init(r21496); mpfr_init(r21497); mpfr_init(r21498); mpfr_init(r21499); mpfr_init(r21500); mpfr_init(r21501); mpfr_init(r21502); mpfr_init(r21503); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21484, r21482, r21483, MPFR_RNDN); mpfr_set_d(r21485, l, MPFR_RNDN); mpfr_mul(r21486, r21483, r21485, MPFR_RNDN); mpfr_set_d(r21487, Om, MPFR_RNDN); mpfr_div(r21488, r21486, r21487, MPFR_RNDN); mpfr_pow(r21489, r21488, r21483, MPFR_RNDN); mpfr_set_d(r21490, kx, MPFR_RNDN); mpfr_sin(r21491, r21490, MPFR_RNDN); mpfr_pow(r21492, r21491, r21483, MPFR_RNDN); mpfr_set_d(r21493, ky, MPFR_RNDN); mpfr_sin(r21494, r21493, MPFR_RNDN); mpfr_pow(r21495, r21494, r21483, MPFR_RNDN); mpfr_add(r21496, r21492, r21495, MPFR_RNDN); mpfr_mul(r21497, r21489, r21496, MPFR_RNDN); mpfr_add(r21498, r21482, r21497, MPFR_RNDN); mpfr_sqrt(r21499, r21498, MPFR_RNDN); mpfr_div(r21500, r21482, r21499, MPFR_RNDN); mpfr_add(r21501, r21482, r21500, MPFR_RNDN); mpfr_mul(r21502, r21484, r21501, MPFR_RNDN); mpfr_sqrt(r21503, r21502, MPFR_RNDN); return mpfr_get_d(r21503, MPFR_RNDN); }