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