#include #include #include #include #include char *name = "Toniolo and Linder, Equation (3a)"; double f_if(float l, float Om, float kx, float ky) { float r21436 = 1; float r21437 = 2; float r21438 = r21436 / r21437; float r21439 = l; float r21440 = r21437 * r21439; float r21441 = Om; float r21442 = r21440 / r21441; float r21443 = pow(r21442, r21437); float r21444 = kx; float r21445 = sin(r21444); float r21446 = pow(r21445, r21437); float r21447 = ky; float r21448 = sin(r21447); float r21449 = pow(r21448, r21437); float r21450 = r21446 + r21449; float r21451 = r21443 * r21450; float r21452 = r21436 + r21451; float r21453 = sqrt(r21452); float r21454 = r21436 / r21453; float r21455 = r21436 + r21454; float r21456 = r21438 * r21455; float r21457 = sqrt(r21456); return r21457; } double f_id(double l, double Om, double kx, double ky) { double r21458 = 1; double r21459 = 2; double r21460 = r21458 / r21459; double r21461 = l; double r21462 = r21459 * r21461; double r21463 = Om; double r21464 = r21462 / r21463; double r21465 = pow(r21464, r21459); double r21466 = kx; double r21467 = sin(r21466); double r21468 = pow(r21467, r21459); double r21469 = ky; double r21470 = sin(r21469); double r21471 = pow(r21470, r21459); double r21472 = r21468 + r21471; double r21473 = r21465 * r21472; double r21474 = r21458 + r21473; double r21475 = sqrt(r21474); double r21476 = r21458 / r21475; double r21477 = r21458 + r21476; double r21478 = r21460 * r21477; double r21479 = sqrt(r21478); return r21479; } double f_of(float l, float Om, float kx, float ky) { float r21480 = 1; float r21481 = 2; float r21482 = r21480 / r21481; float r21483 = l; float r21484 = r21481 * r21483; float r21485 = Om; float r21486 = r21484 / r21485; float r21487 = pow(r21486, r21481); float r21488 = kx; float r21489 = sin(r21488); float r21490 = pow(r21489, r21481); float r21491 = ky; float r21492 = sin(r21491); float r21493 = pow(r21492, r21481); float r21494 = r21490 + r21493; float r21495 = r21487 * r21494; float r21496 = r21480 + r21495; float r21497 = sqrt(r21496); float r21498 = r21480 / r21497; float r21499 = r21480 + r21498; float r21500 = r21482 * r21499; float r21501 = sqrt(r21500); return r21501; } double f_od(double l, double Om, double kx, double ky) { double r21502 = 1; double r21503 = 2; double r21504 = r21502 / r21503; double r21505 = l; double r21506 = r21503 * r21505; double r21507 = Om; double r21508 = r21506 / r21507; double r21509 = pow(r21508, r21503); double r21510 = kx; double r21511 = sin(r21510); double r21512 = pow(r21511, r21503); double r21513 = ky; double r21514 = sin(r21513); double r21515 = pow(r21514, r21503); double r21516 = r21512 + r21515; double r21517 = r21509 * r21516; double r21518 = r21502 + r21517; double r21519 = sqrt(r21518); double r21520 = r21502 / r21519; double r21521 = r21502 + r21520; double r21522 = r21504 * r21521; double r21523 = sqrt(r21522); return r21523; } 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 r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535, r21536, r21537, r21538, r21539, r21540, r21541, r21542, r21543, r21544, r21545; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21524, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21525, "2", 10, MPFR_RNDN); mpfr_init(r21526); mpfr_init(r21527); mpfr_init(r21528); mpfr_init(r21529); mpfr_init(r21530); mpfr_init(r21531); mpfr_init(r21532); mpfr_init(r21533); mpfr_init(r21534); mpfr_init(r21535); mpfr_init(r21536); mpfr_init(r21537); mpfr_init(r21538); mpfr_init(r21539); mpfr_init(r21540); mpfr_init(r21541); mpfr_init(r21542); mpfr_init(r21543); mpfr_init(r21544); mpfr_init(r21545); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21526, r21524, r21525, MPFR_RNDN); mpfr_set_d(r21527, l, MPFR_RNDN); mpfr_mul(r21528, r21525, r21527, MPFR_RNDN); mpfr_set_d(r21529, Om, MPFR_RNDN); mpfr_div(r21530, r21528, r21529, MPFR_RNDN); mpfr_pow(r21531, r21530, r21525, MPFR_RNDN); mpfr_set_d(r21532, kx, MPFR_RNDN); mpfr_sin(r21533, r21532, MPFR_RNDN); mpfr_pow(r21534, r21533, r21525, MPFR_RNDN); mpfr_set_d(r21535, ky, MPFR_RNDN); mpfr_sin(r21536, r21535, MPFR_RNDN); mpfr_pow(r21537, r21536, r21525, MPFR_RNDN); mpfr_add(r21538, r21534, r21537, MPFR_RNDN); mpfr_mul(r21539, r21531, r21538, MPFR_RNDN); mpfr_add(r21540, r21524, r21539, MPFR_RNDN); mpfr_sqrt(r21541, r21540, MPFR_RNDN); mpfr_div(r21542, r21524, r21541, MPFR_RNDN); mpfr_add(r21543, r21524, r21542, MPFR_RNDN); mpfr_mul(r21544, r21526, r21543, MPFR_RNDN); mpfr_sqrt(r21545, r21544, MPFR_RNDN); return mpfr_get_d(r21545, MPFR_RNDN); } static mpfr_t r21546, r21547, r21548, r21549, r21550, r21551, r21552, r21553, r21554, r21555, r21556, r21557, r21558, r21559, r21560, r21561, r21562, r21563, r21564, r21565, r21566, r21567; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21546, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21547, "2", 10, MPFR_RNDN); mpfr_init(r21548); mpfr_init(r21549); mpfr_init(r21550); mpfr_init(r21551); mpfr_init(r21552); mpfr_init(r21553); mpfr_init(r21554); mpfr_init(r21555); mpfr_init(r21556); mpfr_init(r21557); mpfr_init(r21558); mpfr_init(r21559); mpfr_init(r21560); mpfr_init(r21561); mpfr_init(r21562); mpfr_init(r21563); mpfr_init(r21564); mpfr_init(r21565); mpfr_init(r21566); mpfr_init(r21567); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21548, r21546, r21547, MPFR_RNDN); mpfr_set_d(r21549, l, MPFR_RNDN); mpfr_mul(r21550, r21547, r21549, MPFR_RNDN); mpfr_set_d(r21551, Om, MPFR_RNDN); mpfr_div(r21552, r21550, r21551, MPFR_RNDN); mpfr_pow(r21553, r21552, r21547, MPFR_RNDN); mpfr_set_d(r21554, kx, MPFR_RNDN); mpfr_sin(r21555, r21554, MPFR_RNDN); mpfr_pow(r21556, r21555, r21547, MPFR_RNDN); mpfr_set_d(r21557, ky, MPFR_RNDN); mpfr_sin(r21558, r21557, MPFR_RNDN); mpfr_pow(r21559, r21558, r21547, MPFR_RNDN); mpfr_add(r21560, r21556, r21559, MPFR_RNDN); mpfr_mul(r21561, r21553, r21560, MPFR_RNDN); mpfr_add(r21562, r21546, r21561, MPFR_RNDN); mpfr_sqrt(r21563, r21562, MPFR_RNDN); mpfr_div(r21564, r21546, r21563, MPFR_RNDN); mpfr_add(r21565, r21546, r21564, MPFR_RNDN); mpfr_mul(r21566, r21548, r21565, MPFR_RNDN); mpfr_sqrt(r21567, r21566, MPFR_RNDN); return mpfr_get_d(r21567, MPFR_RNDN); } static mpfr_t r21568, r21569, r21570, r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579, r21580, r21581, r21582, r21583, r21584, r21585, r21586, r21587, r21588, r21589; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21568, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21569, "2", 10, MPFR_RNDN); mpfr_init(r21570); mpfr_init(r21571); mpfr_init(r21572); mpfr_init(r21573); mpfr_init(r21574); mpfr_init(r21575); mpfr_init(r21576); mpfr_init(r21577); mpfr_init(r21578); mpfr_init(r21579); mpfr_init(r21580); mpfr_init(r21581); mpfr_init(r21582); mpfr_init(r21583); mpfr_init(r21584); mpfr_init(r21585); mpfr_init(r21586); mpfr_init(r21587); mpfr_init(r21588); mpfr_init(r21589); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21570, r21568, r21569, MPFR_RNDN); mpfr_set_d(r21571, l, MPFR_RNDN); mpfr_mul(r21572, r21569, r21571, MPFR_RNDN); mpfr_set_d(r21573, Om, MPFR_RNDN); mpfr_div(r21574, r21572, r21573, MPFR_RNDN); mpfr_pow(r21575, r21574, r21569, MPFR_RNDN); mpfr_set_d(r21576, kx, MPFR_RNDN); mpfr_sin(r21577, r21576, MPFR_RNDN); mpfr_pow(r21578, r21577, r21569, MPFR_RNDN); mpfr_set_d(r21579, ky, MPFR_RNDN); mpfr_sin(r21580, r21579, MPFR_RNDN); mpfr_pow(r21581, r21580, r21569, MPFR_RNDN); mpfr_add(r21582, r21578, r21581, MPFR_RNDN); mpfr_mul(r21583, r21575, r21582, MPFR_RNDN); mpfr_add(r21584, r21568, r21583, MPFR_RNDN); mpfr_sqrt(r21585, r21584, MPFR_RNDN); mpfr_div(r21586, r21568, r21585, MPFR_RNDN); mpfr_add(r21587, r21568, r21586, MPFR_RNDN); mpfr_mul(r21588, r21570, r21587, MPFR_RNDN); mpfr_sqrt(r21589, r21588, MPFR_RNDN); return mpfr_get_d(r21589, MPFR_RNDN); }