#include #include #include #include #include char *name = "Toniolo and Linder, Equation (10-)"; double f_if(float t, float l, float k) { float r19523 = 2.0f; float r19524 = t; float r19525 = 3.0f; float r19526 = pow(r19524, r19525); float r19527 = l; float r19528 = r19527 * r19527; float r19529 = r19526 / r19528; float r19530 = k; float r19531 = sin(r19530); float r19532 = r19529 * r19531; float r19533 = tan(r19530); float r19534 = r19532 * r19533; float r19535 = 1.0f; float r19536 = r19530 / r19524; float r19537 = r19536 * r19536; float r19538 = r19535 + r19537; float r19539 = r19538 - r19535; float r19540 = r19534 * r19539; float r19541 = r19523 / r19540; return r19541; } double f_id(double t, double l, double k) { double r19542 = 2.0; double r19543 = t; double r19544 = 3.0; double r19545 = pow(r19543, r19544); double r19546 = l; double r19547 = r19546 * r19546; double r19548 = r19545 / r19547; double r19549 = k; double r19550 = sin(r19549); double r19551 = r19548 * r19550; double r19552 = tan(r19549); double r19553 = r19551 * r19552; double r19554 = 1.0; double r19555 = r19549 / r19543; double r19556 = r19555 * r19555; double r19557 = r19554 + r19556; double r19558 = r19557 - r19554; double r19559 = r19553 * r19558; double r19560 = r19542 / r19559; return r19560; } double f_of(float t, float l, float k) { float r19561 = 2.0f; float r19562 = k; float r19563 = tan(r19562); float r19564 = cbrt(r19563); float r19565 = sin(r19562); float r19566 = cbrt(r19565); float r19567 = r19564 * r19566; float r19568 = t; float r19569 = r19562 / r19568; float r19570 = cbrt(r19569); float r19571 = l; float r19572 = cbrt(r19571); float r19573 = r19572 * r19572; float r19574 = r19568 / r19573; float r19575 = r19570 * r19574; float r19576 = r19570 * r19575; float r19577 = r19567 * r19576; float r19578 = r19577 * (r19577 * r19577); float r19579 = r19561 / r19578; return r19579; } double f_od(double t, double l, double k) { double r19580 = 2.0; double r19581 = k; double r19582 = tan(r19581); double r19583 = cbrt(r19582); double r19584 = sin(r19581); double r19585 = cbrt(r19584); double r19586 = r19583 * r19585; double r19587 = t; double r19588 = r19581 / r19587; double r19589 = cbrt(r19588); double r19590 = l; double r19591 = cbrt(r19590); double r19592 = r19591 * r19591; double r19593 = r19587 / r19592; double r19594 = r19589 * r19593; double r19595 = r19589 * r19594; double r19596 = r19586 * r19595; double r19597 = r19596 * (r19596 * r19596); double r19598 = r19580 / r19597; return r19598; } 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 r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607, r19608, r19609, r19610, r19611, r19612, r19613, r19614, r19615, r19616, r19617; void setup_mpfr_f_im() { mpfr_set_default_prec(3216); mpfr_init_set_str(r19599, "2", 10, MPFR_RNDN); mpfr_init(r19600); mpfr_init_set_str(r19601, "3", 10, MPFR_RNDN); mpfr_init(r19602); mpfr_init(r19603); mpfr_init(r19604); mpfr_init(r19605); mpfr_init(r19606); mpfr_init(r19607); mpfr_init(r19608); mpfr_init(r19609); mpfr_init(r19610); mpfr_init_set_str(r19611, "1", 10, MPFR_RNDN); mpfr_init(r19612); mpfr_init(r19613); mpfr_init(r19614); mpfr_init(r19615); mpfr_init(r19616); mpfr_init(r19617); } double f_im(double t, double l, double k) { ; mpfr_set_d(r19600, t, MPFR_RNDN); ; mpfr_pow(r19602, r19600, r19601, MPFR_RNDN); mpfr_set_d(r19603, l, MPFR_RNDN); mpfr_sqr(r19604, r19603, MPFR_RNDN); mpfr_div(r19605, r19602, r19604, MPFR_RNDN); mpfr_set_d(r19606, k, MPFR_RNDN); mpfr_sin(r19607, r19606, MPFR_RNDN); mpfr_mul(r19608, r19605, r19607, MPFR_RNDN); mpfr_tan(r19609, r19606, MPFR_RNDN); mpfr_mul(r19610, r19608, r19609, MPFR_RNDN); ; mpfr_div(r19612, r19606, r19600, MPFR_RNDN); mpfr_sqr(r19613, r19612, MPFR_RNDN); mpfr_add(r19614, r19611, r19613, MPFR_RNDN); mpfr_sub(r19615, r19614, r19611, MPFR_RNDN); mpfr_mul(r19616, r19610, r19615, MPFR_RNDN); mpfr_div(r19617, r19599, r19616, MPFR_RNDN); return mpfr_get_d(r19617, MPFR_RNDN); } static mpfr_t r19618, r19619, r19620, r19621, r19622, r19623, r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636; void setup_mpfr_f_fm() { mpfr_set_default_prec(3216); mpfr_init_set_str(r19618, "2", 10, MPFR_RNDN); mpfr_init(r19619); mpfr_init(r19620); mpfr_init(r19621); mpfr_init(r19622); mpfr_init(r19623); mpfr_init(r19624); mpfr_init(r19625); mpfr_init(r19626); mpfr_init(r19627); mpfr_init(r19628); mpfr_init(r19629); mpfr_init(r19630); mpfr_init(r19631); mpfr_init(r19632); mpfr_init(r19633); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); } double f_fm(double t, double l, double k) { ; mpfr_set_d(r19619, k, MPFR_RNDN); mpfr_tan(r19620, r19619, MPFR_RNDN); mpfr_cbrt(r19621, r19620, MPFR_RNDN); mpfr_sin(r19622, r19619, MPFR_RNDN); mpfr_cbrt(r19623, r19622, MPFR_RNDN); mpfr_mul(r19624, r19621, r19623, MPFR_RNDN); mpfr_set_d(r19625, t, MPFR_RNDN); mpfr_div(r19626, r19619, r19625, MPFR_RNDN); mpfr_cbrt(r19627, r19626, MPFR_RNDN); mpfr_set_d(r19628, l, MPFR_RNDN); mpfr_cbrt(r19629, r19628, MPFR_RNDN); mpfr_mul(r19630, r19629, r19629, MPFR_RNDN); mpfr_div(r19631, r19625, r19630, MPFR_RNDN); mpfr_mul(r19632, r19627, r19631, MPFR_RNDN); mpfr_mul(r19633, r19627, r19632, MPFR_RNDN); mpfr_mul(r19634, r19624, r19633, MPFR_RNDN); mpfr_mul(r19635, r19634, r19634, MPFR_RNDN); mpfr_mul(r19635, r19635, r19634, MPFR_RNDN); mpfr_div(r19636, r19618, r19635, MPFR_RNDN); return mpfr_get_d(r19636, MPFR_RNDN); } static mpfr_t r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650, r19651, r19652, r19653, r19654, r19655; void setup_mpfr_f_dm() { mpfr_set_default_prec(3216); mpfr_init_set_str(r19637, "2", 10, MPFR_RNDN); mpfr_init(r19638); mpfr_init(r19639); mpfr_init(r19640); mpfr_init(r19641); mpfr_init(r19642); mpfr_init(r19643); mpfr_init(r19644); mpfr_init(r19645); mpfr_init(r19646); mpfr_init(r19647); mpfr_init(r19648); mpfr_init(r19649); mpfr_init(r19650); mpfr_init(r19651); mpfr_init(r19652); mpfr_init(r19653); mpfr_init(r19654); mpfr_init(r19655); } double f_dm(double t, double l, double k) { ; mpfr_set_d(r19638, k, MPFR_RNDN); mpfr_tan(r19639, r19638, MPFR_RNDN); mpfr_cbrt(r19640, r19639, MPFR_RNDN); mpfr_sin(r19641, r19638, MPFR_RNDN); mpfr_cbrt(r19642, r19641, MPFR_RNDN); mpfr_mul(r19643, r19640, r19642, MPFR_RNDN); mpfr_set_d(r19644, t, MPFR_RNDN); mpfr_div(r19645, r19638, r19644, MPFR_RNDN); mpfr_cbrt(r19646, r19645, MPFR_RNDN); mpfr_set_d(r19647, l, MPFR_RNDN); mpfr_cbrt(r19648, r19647, MPFR_RNDN); mpfr_mul(r19649, r19648, r19648, MPFR_RNDN); mpfr_div(r19650, r19644, r19649, MPFR_RNDN); mpfr_mul(r19651, r19646, r19650, MPFR_RNDN); mpfr_mul(r19652, r19646, r19651, MPFR_RNDN); mpfr_mul(r19653, r19643, r19652, MPFR_RNDN); mpfr_mul(r19654, r19653, r19653, MPFR_RNDN); mpfr_mul(r19654, r19654, r19653, MPFR_RNDN); mpfr_div(r19655, r19637, r19654, MPFR_RNDN); return mpfr_get_d(r19655, MPFR_RNDN); }