#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (6)"; double f_if(float F, float l) { float r20554 = atan2(1.0, 0.0); float r20555 = l; float r20556 = r20554 * r20555; float r20557 = 1; float r20558 = F; float r20559 = r20558 * r20558; float r20560 = r20557 / r20559; float r20561 = tan(r20556); float r20562 = r20560 * r20561; float r20563 = r20556 - r20562; return r20563; } double f_id(double F, double l) { double r20564 = atan2(1.0, 0.0); double r20565 = l; double r20566 = r20564 * r20565; double r20567 = 1; double r20568 = F; double r20569 = r20568 * r20568; double r20570 = r20567 / r20569; double r20571 = tan(r20566); double r20572 = r20570 * r20571; double r20573 = r20566 - r20572; return r20573; } double f_of(float F, float l) { float r20574 = atan2(1.0, 0.0); float r20575 = l; float r20576 = r20574 * r20575; float r20577 = 1; float r20578 = F; float r20579 = r20577 / r20578; float r20580 = r20578 / r20574; float r20581 = r20580 / r20575; float r20582 = r20575 * r20578; float r20583 = 1/3; float r20584 = r20583 * r20574; float r20585 = r20582 * r20584; float r20586 = cbrt(r20585); float r20587 = r20586 * r20586; float r20588 = r20587 * r20586; float r20589 = r20581 - r20588; float r20590 = r20579 / r20589; float r20591 = r20576 - r20590; return r20591; } double f_od(double F, double l) { double r20592 = atan2(1.0, 0.0); double r20593 = l; double r20594 = r20592 * r20593; double r20595 = 1; double r20596 = F; double r20597 = r20595 / r20596; double r20598 = r20596 / r20592; double r20599 = r20598 / r20593; double r20600 = r20593 * r20596; double r20601 = 1/3; double r20602 = r20601 * r20592; double r20603 = r20600 * r20602; double r20604 = cbrt(r20603); double r20605 = r20604 * r20604; double r20606 = r20605 * r20604; double r20607 = r20599 - r20606; double r20608 = r20597 / r20607; double r20609 = r20594 - r20608; return r20609; } 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 r20610, r20611, r20612, r20613, r20614, r20615, r20616, r20617, r20618, r20619; void setup_mpfr_f_im() { mpfr_set_default_prec(3728); mpfr_init(r20610); mpfr_init(r20611); mpfr_init(r20612); mpfr_init_set_str(r20613, "1", 10, MPFR_RNDN); mpfr_init(r20614); mpfr_init(r20615); mpfr_init(r20616); mpfr_init(r20617); mpfr_init(r20618); mpfr_init(r20619); } double f_im(double F, double l) { mpfr_const_pi(r20610, MPFR_RNDN); mpfr_set_d(r20611, l, MPFR_RNDN); mpfr_mul(r20612, r20610, r20611, MPFR_RNDN); ; mpfr_set_d(r20614, F, MPFR_RNDN); mpfr_mul(r20615, r20614, r20614, MPFR_RNDN); mpfr_div(r20616, r20613, r20615, MPFR_RNDN); mpfr_tan(r20617, r20612, MPFR_RNDN); mpfr_mul(r20618, r20616, r20617, MPFR_RNDN); mpfr_sub(r20619, r20612, r20618, MPFR_RNDN); return mpfr_get_d(r20619, MPFR_RNDN); } static mpfr_t r20620, r20621, r20622, r20623, r20624, r20625, r20626, r20627, r20628, r20629, r20630, r20631, r20632, r20633, r20634, r20635, r20636, r20637; void setup_mpfr_f_fm() { mpfr_set_default_prec(3728); mpfr_init(r20620); mpfr_init(r20621); mpfr_init(r20622); mpfr_init_set_str(r20623, "1", 10, MPFR_RNDN); mpfr_init(r20624); mpfr_init(r20625); mpfr_init(r20626); mpfr_init(r20627); mpfr_init(r20628); mpfr_init_set_str(r20629, "1/3", 10, MPFR_RNDN); mpfr_init(r20630); mpfr_init(r20631); mpfr_init(r20632); mpfr_init(r20633); mpfr_init(r20634); mpfr_init(r20635); mpfr_init(r20636); mpfr_init(r20637); } double f_fm(double F, double l) { mpfr_const_pi(r20620, MPFR_RNDN); mpfr_set_d(r20621, l, MPFR_RNDN); mpfr_mul(r20622, r20620, r20621, MPFR_RNDN); ; mpfr_set_d(r20624, F, MPFR_RNDN); mpfr_div(r20625, r20623, r20624, MPFR_RNDN); mpfr_div(r20626, r20624, r20620, MPFR_RNDN); mpfr_div(r20627, r20626, r20621, MPFR_RNDN); mpfr_mul(r20628, r20621, r20624, MPFR_RNDN); ; mpfr_mul(r20630, r20629, r20620, MPFR_RNDN); mpfr_mul(r20631, r20628, r20630, MPFR_RNDN); mpfr_cbrt(r20632, r20631, MPFR_RNDN); mpfr_mul(r20633, r20632, r20632, MPFR_RNDN); mpfr_mul(r20634, r20633, r20632, MPFR_RNDN); mpfr_sub(r20635, r20627, r20634, MPFR_RNDN); mpfr_div(r20636, r20625, r20635, MPFR_RNDN); mpfr_sub(r20637, r20622, r20636, MPFR_RNDN); return mpfr_get_d(r20637, MPFR_RNDN); } static mpfr_t r20638, r20639, r20640, r20641, r20642, r20643, r20644, r20645, r20646, r20647, r20648, r20649, r20650, r20651, r20652, r20653, r20654, r20655; void setup_mpfr_f_dm() { mpfr_set_default_prec(3728); mpfr_init(r20638); mpfr_init(r20639); mpfr_init(r20640); mpfr_init_set_str(r20641, "1", 10, MPFR_RNDN); mpfr_init(r20642); mpfr_init(r20643); mpfr_init(r20644); mpfr_init(r20645); mpfr_init(r20646); mpfr_init_set_str(r20647, "1/3", 10, MPFR_RNDN); mpfr_init(r20648); mpfr_init(r20649); mpfr_init(r20650); mpfr_init(r20651); mpfr_init(r20652); mpfr_init(r20653); mpfr_init(r20654); mpfr_init(r20655); } double f_dm(double F, double l) { mpfr_const_pi(r20638, MPFR_RNDN); mpfr_set_d(r20639, l, MPFR_RNDN); mpfr_mul(r20640, r20638, r20639, MPFR_RNDN); ; mpfr_set_d(r20642, F, MPFR_RNDN); mpfr_div(r20643, r20641, r20642, MPFR_RNDN); mpfr_div(r20644, r20642, r20638, MPFR_RNDN); mpfr_div(r20645, r20644, r20639, MPFR_RNDN); mpfr_mul(r20646, r20639, r20642, MPFR_RNDN); ; mpfr_mul(r20648, r20647, r20638, MPFR_RNDN); mpfr_mul(r20649, r20646, r20648, MPFR_RNDN); mpfr_cbrt(r20650, r20649, MPFR_RNDN); mpfr_mul(r20651, r20650, r20650, MPFR_RNDN); mpfr_mul(r20652, r20651, r20650, MPFR_RNDN); mpfr_sub(r20653, r20645, r20652, MPFR_RNDN); mpfr_div(r20654, r20643, r20653, MPFR_RNDN); mpfr_sub(r20655, r20640, r20654, MPFR_RNDN); return mpfr_get_d(r20655, MPFR_RNDN); }