#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r20532 = 1; float r20533 = atan2(1.0, 0.0); float r20534 = 4; float r20535 = r20533 / r20534; float r20536 = r20532 / r20535; float r20537 = f; float r20538 = r20535 * r20537; float r20539 = exp(r20538); float r20540 = -r20538; float r20541 = exp(r20540); float r20542 = r20539 + r20541; float r20543 = r20539 - r20541; float r20544 = r20542 / r20543; float r20545 = log(r20544); float r20546 = r20536 * r20545; float r20547 = -r20546; return r20547; } double f_id(double f) { double r20548 = 1; double r20549 = atan2(1.0, 0.0); double r20550 = 4; double r20551 = r20549 / r20550; double r20552 = r20548 / r20551; double r20553 = f; double r20554 = r20551 * r20553; double r20555 = exp(r20554); double r20556 = -r20554; double r20557 = exp(r20556); double r20558 = r20555 + r20557; double r20559 = r20555 - r20557; double r20560 = r20558 / r20559; double r20561 = log(r20560); double r20562 = r20552 * r20561; double r20563 = -r20562; return r20563; } double f_of(float f) { float r20564 = atan2(1.0, 0.0); float r20565 = f; float r20566 = 1/12; float r20567 = r20565 * r20566; float r20568 = r20567 * r20565; float r20569 = 7/5760; float r20570 = 4; float r20571 = pow(r20565, r20570); float r20572 = r20569 * r20571; float r20573 = r20564 * r20564; float r20574 = r20572 * r20573; float r20575 = r20568 - r20574; float r20576 = r20570 / r20564; float r20577 = log(r20576); float r20578 = r20576 * r20577; float r20579 = log(r20565); float r20580 = r20576 * r20579; float r20581 = r20578 - r20580; float r20582 = fma(r20564, r20575, r20581); float r20583 = -r20582; return r20583; } double f_od(double f) { double r20584 = atan2(1.0, 0.0); double r20585 = f; double r20586 = 1/12; double r20587 = r20585 * r20586; double r20588 = r20587 * r20585; double r20589 = 7/5760; double r20590 = 4; double r20591 = pow(r20585, r20590); double r20592 = r20589 * r20591; double r20593 = r20584 * r20584; double r20594 = r20592 * r20593; double r20595 = r20588 - r20594; double r20596 = r20590 / r20584; double r20597 = log(r20596); double r20598 = r20596 * r20597; double r20599 = log(r20585); double r20600 = r20596 * r20599; double r20601 = r20598 - r20600; double r20602 = fma(r20584, r20595, r20601); double r20603 = -r20602; return r20603; } 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 r20604, r20605, r20606, r20607, r20608, r20609, r20610, r20611, r20612, r20613, r20614, r20615, r20616, r20617, r20618, r20619; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20604, "1", 10, MPFR_RNDN); mpfr_init(r20605); mpfr_init_set_str(r20606, "4", 10, MPFR_RNDN); mpfr_init(r20607); mpfr_init(r20608); mpfr_init(r20609); mpfr_init(r20610); mpfr_init(r20611); mpfr_init(r20612); mpfr_init(r20613); mpfr_init(r20614); mpfr_init(r20615); mpfr_init(r20616); mpfr_init(r20617); mpfr_init(r20618); mpfr_init(r20619); } double f_im(double f) { ; mpfr_const_pi(r20605, MPFR_RNDN); ; mpfr_div(r20607, r20605, r20606, MPFR_RNDN); mpfr_div(r20608, r20604, r20607, MPFR_RNDN); mpfr_set_d(r20609, f, MPFR_RNDN); mpfr_mul(r20610, r20607, r20609, MPFR_RNDN); mpfr_exp(r20611, r20610, MPFR_RNDN); mpfr_neg(r20612, r20610, MPFR_RNDN); mpfr_exp(r20613, r20612, MPFR_RNDN); mpfr_add(r20614, r20611, r20613, MPFR_RNDN); mpfr_sub(r20615, r20611, r20613, MPFR_RNDN); mpfr_div(r20616, r20614, r20615, MPFR_RNDN); mpfr_log(r20617, r20616, MPFR_RNDN); mpfr_mul(r20618, r20608, r20617, MPFR_RNDN); mpfr_neg(r20619, 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, r20638, r20639; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20620); mpfr_init(r20621); mpfr_init_set_str(r20622, "1/12", 10, MPFR_RNDN); mpfr_init(r20623); mpfr_init(r20624); mpfr_init_set_str(r20625, "7/5760", 10, MPFR_RNDN); mpfr_init_set_str(r20626, "4", 10, MPFR_RNDN); mpfr_init(r20627); mpfr_init(r20628); mpfr_init(r20629); mpfr_init(r20630); mpfr_init(r20631); mpfr_init(r20632); mpfr_init(r20633); mpfr_init(r20634); mpfr_init(r20635); mpfr_init(r20636); mpfr_init(r20637); mpfr_init(r20638); mpfr_init(r20639); } double f_fm(double f) { mpfr_const_pi(r20620, MPFR_RNDN); mpfr_set_d(r20621, f, MPFR_RNDN); ; mpfr_mul(r20623, r20621, r20622, MPFR_RNDN); mpfr_mul(r20624, r20623, r20621, MPFR_RNDN); ; ; mpfr_pow(r20627, r20621, r20626, MPFR_RNDN); mpfr_mul(r20628, r20625, r20627, MPFR_RNDN); mpfr_mul(r20629, r20620, r20620, MPFR_RNDN); mpfr_mul(r20630, r20628, r20629, MPFR_RNDN); mpfr_sub(r20631, r20624, r20630, MPFR_RNDN); mpfr_div(r20632, r20626, r20620, MPFR_RNDN); mpfr_log(r20633, r20632, MPFR_RNDN); mpfr_mul(r20634, r20632, r20633, MPFR_RNDN); mpfr_log(r20635, r20621, MPFR_RNDN); mpfr_mul(r20636, r20632, r20635, MPFR_RNDN); mpfr_sub(r20637, r20634, r20636, MPFR_RNDN); mpfr_fma(r20638, r20620, r20631, r20637, MPFR_RNDN); mpfr_neg(r20639, r20638, MPFR_RNDN); return mpfr_get_d(r20639, MPFR_RNDN); } static mpfr_t r20640, r20641, r20642, r20643, r20644, r20645, r20646, r20647, r20648, r20649, r20650, r20651, r20652, r20653, r20654, r20655, r20656, r20657, r20658, r20659; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20640); mpfr_init(r20641); mpfr_init_set_str(r20642, "1/12", 10, MPFR_RNDN); mpfr_init(r20643); mpfr_init(r20644); mpfr_init_set_str(r20645, "7/5760", 10, MPFR_RNDN); mpfr_init_set_str(r20646, "4", 10, MPFR_RNDN); mpfr_init(r20647); mpfr_init(r20648); mpfr_init(r20649); mpfr_init(r20650); mpfr_init(r20651); mpfr_init(r20652); mpfr_init(r20653); mpfr_init(r20654); mpfr_init(r20655); mpfr_init(r20656); mpfr_init(r20657); mpfr_init(r20658); mpfr_init(r20659); } double f_dm(double f) { mpfr_const_pi(r20640, MPFR_RNDN); mpfr_set_d(r20641, f, MPFR_RNDN); ; mpfr_mul(r20643, r20641, r20642, MPFR_RNDN); mpfr_mul(r20644, r20643, r20641, MPFR_RNDN); ; ; mpfr_pow(r20647, r20641, r20646, MPFR_RNDN); mpfr_mul(r20648, r20645, r20647, MPFR_RNDN); mpfr_mul(r20649, r20640, r20640, MPFR_RNDN); mpfr_mul(r20650, r20648, r20649, MPFR_RNDN); mpfr_sub(r20651, r20644, r20650, MPFR_RNDN); mpfr_div(r20652, r20646, r20640, MPFR_RNDN); mpfr_log(r20653, r20652, MPFR_RNDN); mpfr_mul(r20654, r20652, r20653, MPFR_RNDN); mpfr_log(r20655, r20641, MPFR_RNDN); mpfr_mul(r20656, r20652, r20655, MPFR_RNDN); mpfr_sub(r20657, r20654, r20656, MPFR_RNDN); mpfr_fma(r20658, r20640, r20651, r20657, MPFR_RNDN); mpfr_neg(r20659, r20658, MPFR_RNDN); return mpfr_get_d(r20659, MPFR_RNDN); }