#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (6)"; double f_if(float F, float l) { float r19531 = atan2(1.0, 0.0); float r19532 = l; float r19533 = r19531 * r19532; float r19534 = 1; float r19535 = F; float r19536 = r19535 * r19535; float r19537 = r19534 / r19536; float r19538 = tan(r19533); float r19539 = r19537 * r19538; float r19540 = r19533 - r19539; return r19540; } double f_id(double F, double l) { double r19541 = atan2(1.0, 0.0); double r19542 = l; double r19543 = r19541 * r19542; double r19544 = 1; double r19545 = F; double r19546 = r19545 * r19545; double r19547 = r19544 / r19546; double r19548 = tan(r19543); double r19549 = r19547 * r19548; double r19550 = r19543 - r19549; return r19550; } double f_of(float F, float l) { float r19551 = atan2(1.0, 0.0); float r19552 = l; float r19553 = r19551 * r19552; float r19554 = 1; float r19555 = F; float r19556 = r19554 / r19555; float r19557 = r19555 / r19551; float r19558 = r19557 / r19552; float r19559 = r19552 * r19555; float r19560 = 1/3; float r19561 = r19551 * r19560; float r19562 = 3; float r19563 = pow(r19561, r19562); float r19564 = cbrt(r19563); float r19565 = r19559 * r19564; float r19566 = r19558 - r19565; float r19567 = r19556 / r19566; float r19568 = r19553 - r19567; return r19568; } double f_od(double F, double l) { double r19569 = atan2(1.0, 0.0); double r19570 = l; double r19571 = r19569 * r19570; double r19572 = 1; double r19573 = F; double r19574 = r19572 / r19573; double r19575 = r19573 / r19569; double r19576 = r19575 / r19570; double r19577 = r19570 * r19573; double r19578 = 1/3; double r19579 = r19569 * r19578; double r19580 = 3; double r19581 = pow(r19579, r19580); double r19582 = cbrt(r19581); double r19583 = r19577 * r19582; double r19584 = r19576 - r19583; double r19585 = r19574 / r19584; double r19586 = r19571 - 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; void setup_mpfr_f_im() { mpfr_set_default_prec(3472); mpfr_init(r19587); mpfr_init(r19588); mpfr_init(r19589); mpfr_init_set_str(r19590, "1", 10, MPFR_RNDN); mpfr_init(r19591); mpfr_init(r19592); mpfr_init(r19593); mpfr_init(r19594); mpfr_init(r19595); mpfr_init(r19596); } double f_im(double F, double l) { mpfr_const_pi(r19587, MPFR_RNDN); mpfr_set_d(r19588, l, MPFR_RNDN); mpfr_mul(r19589, r19587, r19588, MPFR_RNDN); ; mpfr_set_d(r19591, F, MPFR_RNDN); mpfr_mul(r19592, r19591, r19591, MPFR_RNDN); mpfr_div(r19593, r19590, r19592, MPFR_RNDN); mpfr_tan(r19594, r19589, MPFR_RNDN); mpfr_mul(r19595, r19593, r19594, MPFR_RNDN); mpfr_sub(r19596, r19589, r19595, MPFR_RNDN); return mpfr_get_d(r19596, MPFR_RNDN); } static mpfr_t r19597, r19598, r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607, r19608, r19609, r19610, r19611, r19612, r19613, r19614; void setup_mpfr_f_fm() { mpfr_set_default_prec(3472); mpfr_init(r19597); mpfr_init(r19598); mpfr_init(r19599); mpfr_init_set_str(r19600, "1", 10, MPFR_RNDN); mpfr_init(r19601); mpfr_init(r19602); mpfr_init(r19603); mpfr_init(r19604); mpfr_init(r19605); mpfr_init_set_str(r19606, "1/3", 10, MPFR_RNDN); mpfr_init(r19607); mpfr_init_set_str(r19608, "3", 10, MPFR_RNDN); mpfr_init(r19609); mpfr_init(r19610); mpfr_init(r19611); mpfr_init(r19612); mpfr_init(r19613); mpfr_init(r19614); } double f_fm(double F, double l) { mpfr_const_pi(r19597, MPFR_RNDN); mpfr_set_d(r19598, l, MPFR_RNDN); mpfr_mul(r19599, r19597, r19598, MPFR_RNDN); ; mpfr_set_d(r19601, F, MPFR_RNDN); mpfr_div(r19602, r19600, r19601, MPFR_RNDN); mpfr_div(r19603, r19601, r19597, MPFR_RNDN); mpfr_div(r19604, r19603, r19598, MPFR_RNDN); mpfr_mul(r19605, r19598, r19601, MPFR_RNDN); ; mpfr_mul(r19607, r19597, r19606, MPFR_RNDN); ; mpfr_pow(r19609, r19607, r19608, MPFR_RNDN); mpfr_cbrt(r19610, r19609, MPFR_RNDN); mpfr_mul(r19611, r19605, r19610, MPFR_RNDN); mpfr_sub(r19612, r19604, r19611, MPFR_RNDN); mpfr_div(r19613, r19602, r19612, MPFR_RNDN); mpfr_sub(r19614, r19599, r19613, MPFR_RNDN); return mpfr_get_d(r19614, MPFR_RNDN); } static mpfr_t r19615, r19616, r19617, r19618, r19619, r19620, r19621, r19622, r19623, r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632; void setup_mpfr_f_dm() { mpfr_set_default_prec(3472); mpfr_init(r19615); mpfr_init(r19616); mpfr_init(r19617); mpfr_init_set_str(r19618, "1", 10, MPFR_RNDN); mpfr_init(r19619); mpfr_init(r19620); mpfr_init(r19621); mpfr_init(r19622); mpfr_init(r19623); mpfr_init_set_str(r19624, "1/3", 10, MPFR_RNDN); mpfr_init(r19625); mpfr_init_set_str(r19626, "3", 10, MPFR_RNDN); mpfr_init(r19627); mpfr_init(r19628); mpfr_init(r19629); mpfr_init(r19630); mpfr_init(r19631); mpfr_init(r19632); } double f_dm(double F, double l) { mpfr_const_pi(r19615, MPFR_RNDN); mpfr_set_d(r19616, l, MPFR_RNDN); mpfr_mul(r19617, r19615, r19616, MPFR_RNDN); ; mpfr_set_d(r19619, F, MPFR_RNDN); mpfr_div(r19620, r19618, r19619, MPFR_RNDN); mpfr_div(r19621, r19619, r19615, MPFR_RNDN); mpfr_div(r19622, r19621, r19616, MPFR_RNDN); mpfr_mul(r19623, r19616, r19619, MPFR_RNDN); ; mpfr_mul(r19625, r19615, r19624, MPFR_RNDN); ; mpfr_pow(r19627, r19625, r19626, MPFR_RNDN); mpfr_cbrt(r19628, r19627, MPFR_RNDN); mpfr_mul(r19629, r19623, r19628, MPFR_RNDN); mpfr_sub(r19630, r19622, r19629, MPFR_RNDN); mpfr_div(r19631, r19620, r19630, MPFR_RNDN); mpfr_sub(r19632, r19617, r19631, MPFR_RNDN); return mpfr_get_d(r19632, MPFR_RNDN); }