#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r22490 = 2; float r22491 = t; float r22492 = 3; float r22493 = pow(r22491, r22492); float r22494 = l; float r22495 = r22494 * r22494; float r22496 = r22493 / r22495; float r22497 = k; float r22498 = sin(r22497); float r22499 = r22496 * r22498; float r22500 = tan(r22497); float r22501 = r22499 * r22500; float r22502 = 1; float r22503 = r22497 / r22491; float r22504 = pow(r22503, r22490); float r22505 = r22502 + r22504; float r22506 = r22505 - r22502; float r22507 = r22501 * r22506; float r22508 = r22490 / r22507; return r22508; } double f_id(double t, double l, double k) { double r22509 = 2; double r22510 = t; double r22511 = 3; double r22512 = pow(r22510, r22511); double r22513 = l; double r22514 = r22513 * r22513; double r22515 = r22512 / r22514; double r22516 = k; double r22517 = sin(r22516); double r22518 = r22515 * r22517; double r22519 = tan(r22516); double r22520 = r22518 * r22519; double r22521 = 1; double r22522 = r22516 / r22510; double r22523 = pow(r22522, r22509); double r22524 = r22521 + r22523; double r22525 = r22524 - r22521; double r22526 = r22520 * r22525; double r22527 = r22509 / r22526; return r22527; } double f_of(float t, float l, float k) { float r22528 = l; float r22529 = -r22528; float r22530 = k; float r22531 = r22529 / r22530; float r22532 = 2; float r22533 = sin(r22530); float r22534 = r22532 / r22533; float r22535 = r22531 * r22534; float r22536 = cos(r22530); float r22537 = t; float r22538 = r22536 / r22537; float r22539 = r22530 / r22529; float r22540 = r22538 / r22539; float r22541 = r22540 / r22533; float r22542 = r22535 * r22541; return r22542; } double f_od(double t, double l, double k) { double r22543 = l; double r22544 = -r22543; double r22545 = k; double r22546 = r22544 / r22545; double r22547 = 2; double r22548 = sin(r22545); double r22549 = r22547 / r22548; double r22550 = r22546 * r22549; double r22551 = cos(r22545); double r22552 = t; double r22553 = r22551 / r22552; double r22554 = r22545 / r22544; double r22555 = r22553 / r22554; double r22556 = r22555 / r22548; double r22557 = r22550 * r22556; return r22557; } 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 r22558, r22559, r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576; void setup_mpfr_f_im() { mpfr_set_default_prec(4432); mpfr_init_set_str(r22558, "2", 10, MPFR_RNDN); mpfr_init(r22559); mpfr_init_set_str(r22560, "3", 10, MPFR_RNDN); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init_set_str(r22570, "1", 10, MPFR_RNDN); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); } double f_im(double t, double l, double k) { ; mpfr_set_d(r22559, t, MPFR_RNDN); ; mpfr_pow(r22561, r22559, r22560, MPFR_RNDN); mpfr_set_d(r22562, l, MPFR_RNDN); mpfr_mul(r22563, r22562, r22562, MPFR_RNDN); mpfr_div(r22564, r22561, r22563, MPFR_RNDN); mpfr_set_d(r22565, k, MPFR_RNDN); mpfr_sin(r22566, r22565, MPFR_RNDN); mpfr_mul(r22567, r22564, r22566, MPFR_RNDN); mpfr_tan(r22568, r22565, MPFR_RNDN); mpfr_mul(r22569, r22567, r22568, MPFR_RNDN); ; mpfr_div(r22571, r22565, r22559, MPFR_RNDN); mpfr_pow(r22572, r22571, r22558, MPFR_RNDN); mpfr_add(r22573, r22570, r22572, MPFR_RNDN); mpfr_sub(r22574, r22573, r22570, MPFR_RNDN); mpfr_mul(r22575, r22569, r22574, MPFR_RNDN); mpfr_div(r22576, r22558, r22575, MPFR_RNDN); return mpfr_get_d(r22576, MPFR_RNDN); } static mpfr_t r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591; void setup_mpfr_f_fm() { mpfr_set_default_prec(4432); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init_set_str(r22581, "2", 10, MPFR_RNDN); mpfr_init(r22582); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); } double f_fm(double t, double l, double k) { mpfr_set_d(r22577, l, MPFR_RNDN); mpfr_neg(r22578, r22577, MPFR_RNDN); mpfr_set_d(r22579, k, MPFR_RNDN); mpfr_div(r22580, r22578, r22579, MPFR_RNDN); ; mpfr_sin(r22582, r22579, MPFR_RNDN); mpfr_div(r22583, r22581, r22582, MPFR_RNDN); mpfr_mul(r22584, r22580, r22583, MPFR_RNDN); mpfr_cos(r22585, r22579, MPFR_RNDN); mpfr_set_d(r22586, t, MPFR_RNDN); mpfr_div(r22587, r22585, r22586, MPFR_RNDN); mpfr_div(r22588, r22579, r22578, MPFR_RNDN); mpfr_div(r22589, r22587, r22588, MPFR_RNDN); mpfr_div(r22590, r22589, r22582, MPFR_RNDN); mpfr_mul(r22591, r22584, r22590, MPFR_RNDN); return mpfr_get_d(r22591, MPFR_RNDN); } static mpfr_t r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606; void setup_mpfr_f_dm() { mpfr_set_default_prec(4432); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init(r22595); mpfr_init_set_str(r22596, "2", 10, MPFR_RNDN); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); mpfr_init(r22606); } double f_dm(double t, double l, double k) { mpfr_set_d(r22592, l, MPFR_RNDN); mpfr_neg(r22593, r22592, MPFR_RNDN); mpfr_set_d(r22594, k, MPFR_RNDN); mpfr_div(r22595, r22593, r22594, MPFR_RNDN); ; mpfr_sin(r22597, r22594, MPFR_RNDN); mpfr_div(r22598, r22596, r22597, MPFR_RNDN); mpfr_mul(r22599, r22595, r22598, MPFR_RNDN); mpfr_cos(r22600, r22594, MPFR_RNDN); mpfr_set_d(r22601, t, MPFR_RNDN); mpfr_div(r22602, r22600, r22601, MPFR_RNDN); mpfr_div(r22603, r22594, r22593, MPFR_RNDN); mpfr_div(r22604, r22602, r22603, MPFR_RNDN); mpfr_div(r22605, r22604, r22597, MPFR_RNDN); mpfr_mul(r22606, r22599, r22605, MPFR_RNDN); return mpfr_get_d(r22606, MPFR_RNDN); }