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