#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 = kx; float r21484 = sin(r21483); float r21485 = ky; float r21486 = sin(r21485); float r21487 = hypot(r21484, r21486); float r21488 = l; float r21489 = r21488 + r21488; float r21490 = Om; float r21491 = r21489 / r21490; float r21492 = r21487 * r21491; float r21493 = hypot(r21480, r21492); float r21494 = expm1(r21493); float r21495 = log1p(r21494); float r21496 = r21480 / r21495; float r21497 = r21480 + r21496; float r21498 = r21482 * r21497; float r21499 = sqrt(r21498); return r21499; } double f_od(double l, double Om, double kx, double ky) { double r21500 = 1; double r21501 = 2; double r21502 = r21500 / r21501; double r21503 = kx; double r21504 = sin(r21503); double r21505 = ky; double r21506 = sin(r21505); double r21507 = hypot(r21504, r21506); double r21508 = l; double r21509 = r21508 + r21508; double r21510 = Om; double r21511 = r21509 / r21510; double r21512 = r21507 * r21511; double r21513 = hypot(r21500, r21512); double r21514 = expm1(r21513); double r21515 = log1p(r21514); double r21516 = r21500 / r21515; double r21517 = r21500 + r21516; double r21518 = r21502 * r21517; double r21519 = sqrt(r21518); return r21519; } 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 r21520, r21521, r21522, r21523, r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535, r21536, r21537, r21538, r21539, r21540, r21541; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21520, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21521, "2", 10, MPFR_RNDN); mpfr_init(r21522); mpfr_init(r21523); mpfr_init(r21524); mpfr_init(r21525); 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); } double f_im(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21522, r21520, r21521, MPFR_RNDN); mpfr_set_d(r21523, l, MPFR_RNDN); mpfr_mul(r21524, r21521, r21523, MPFR_RNDN); mpfr_set_d(r21525, Om, MPFR_RNDN); mpfr_div(r21526, r21524, r21525, MPFR_RNDN); mpfr_pow(r21527, r21526, r21521, MPFR_RNDN); mpfr_set_d(r21528, kx, MPFR_RNDN); mpfr_sin(r21529, r21528, MPFR_RNDN); mpfr_pow(r21530, r21529, r21521, MPFR_RNDN); mpfr_set_d(r21531, ky, MPFR_RNDN); mpfr_sin(r21532, r21531, MPFR_RNDN); mpfr_pow(r21533, r21532, r21521, MPFR_RNDN); mpfr_add(r21534, r21530, r21533, MPFR_RNDN); mpfr_mul(r21535, r21527, r21534, MPFR_RNDN); mpfr_add(r21536, r21520, r21535, MPFR_RNDN); mpfr_sqrt(r21537, r21536, MPFR_RNDN); mpfr_div(r21538, r21520, r21537, MPFR_RNDN); mpfr_add(r21539, r21520, r21538, MPFR_RNDN); mpfr_mul(r21540, r21522, r21539, MPFR_RNDN); mpfr_sqrt(r21541, r21540, MPFR_RNDN); return mpfr_get_d(r21541, MPFR_RNDN); } static mpfr_t r21542, r21543, r21544, r21545, r21546, r21547, r21548, r21549, r21550, r21551, r21552, r21553, r21554, r21555, r21556, r21557, r21558, r21559, r21560, r21561; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21542, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21543, "2", 10, MPFR_RNDN); mpfr_init(r21544); mpfr_init(r21545); mpfr_init(r21546); mpfr_init(r21547); 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); } double f_fm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21544, r21542, r21543, MPFR_RNDN); mpfr_set_d(r21545, kx, MPFR_RNDN); mpfr_sin(r21546, r21545, MPFR_RNDN); mpfr_set_d(r21547, ky, MPFR_RNDN); mpfr_sin(r21548, r21547, MPFR_RNDN); mpfr_hypot(r21549, r21546, r21548, MPFR_RNDN); mpfr_set_d(r21550, l, MPFR_RNDN); mpfr_add(r21551, r21550, r21550, MPFR_RNDN); mpfr_set_d(r21552, Om, MPFR_RNDN); mpfr_div(r21553, r21551, r21552, MPFR_RNDN); mpfr_mul(r21554, r21549, r21553, MPFR_RNDN); mpfr_hypot(r21555, r21542, r21554, MPFR_RNDN); mpfr_expm1(r21556, r21555, MPFR_RNDN); mpfr_log1p(r21557, r21556, MPFR_RNDN); mpfr_div(r21558, r21542, r21557, MPFR_RNDN); mpfr_add(r21559, r21542, r21558, MPFR_RNDN); mpfr_mul(r21560, r21544, r21559, MPFR_RNDN); mpfr_sqrt(r21561, r21560, MPFR_RNDN); return mpfr_get_d(r21561, MPFR_RNDN); } static mpfr_t r21562, r21563, r21564, r21565, r21566, r21567, r21568, r21569, r21570, r21571, r21572, r21573, r21574, r21575, r21576, r21577, r21578, r21579, r21580, r21581; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21562, "1", 10, MPFR_RNDN); mpfr_init_set_str(r21563, "2", 10, MPFR_RNDN); mpfr_init(r21564); mpfr_init(r21565); mpfr_init(r21566); mpfr_init(r21567); mpfr_init(r21568); mpfr_init(r21569); 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); } double f_dm(double l, double Om, double kx, double ky) { ; ; mpfr_div(r21564, r21562, r21563, MPFR_RNDN); mpfr_set_d(r21565, kx, MPFR_RNDN); mpfr_sin(r21566, r21565, MPFR_RNDN); mpfr_set_d(r21567, ky, MPFR_RNDN); mpfr_sin(r21568, r21567, MPFR_RNDN); mpfr_hypot(r21569, r21566, r21568, MPFR_RNDN); mpfr_set_d(r21570, l, MPFR_RNDN); mpfr_add(r21571, r21570, r21570, MPFR_RNDN); mpfr_set_d(r21572, Om, MPFR_RNDN); mpfr_div(r21573, r21571, r21572, MPFR_RNDN); mpfr_mul(r21574, r21569, r21573, MPFR_RNDN); mpfr_hypot(r21575, r21562, r21574, MPFR_RNDN); mpfr_expm1(r21576, r21575, MPFR_RNDN); mpfr_log1p(r21577, r21576, MPFR_RNDN); mpfr_div(r21578, r21562, r21577, MPFR_RNDN); mpfr_add(r21579, r21562, r21578, MPFR_RNDN); mpfr_mul(r21580, r21564, r21579, MPFR_RNDN); mpfr_sqrt(r21581, r21580, MPFR_RNDN); return mpfr_get_d(r21581, MPFR_RNDN); }