#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r20538 = N; float r20539 = 1; float r20540 = r20538 + r20539; float r20541 = log(r20540); float r20542 = log(r20538); float r20543 = r20541 - r20542; return r20543; } double f_id(double N) { double r20544 = N; double r20545 = 1; double r20546 = r20544 + r20545; double r20547 = log(r20546); double r20548 = log(r20544); double r20549 = r20547 - r20548; return r20549; } double f_of(float N) { float r20550 = N; float r20551 = 5788.106078246906; bool r20552 = r20550 <= r20551; float r20553 = 1; float r20554 = r20550 + r20553; float r20555 = r20554 / r20550; float r20556 = sqrt(r20555); float r20557 = log(r20556); float r20558 = sqrt(r20556); float r20559 = log(r20558); float r20560 = r20559 + r20559; float r20561 = r20557 + r20560; float r20562 = r20553 / r20550; float r20563 = 1/3; float r20564 = r20563 / r20550; float r20565 = 1/2; float r20566 = r20564 - r20565; float r20567 = r20550 * r20550; float r20568 = r20566 / r20567; float r20569 = r20562 + r20568; float r20570 = r20552 ? r20561 : r20569; return r20570; } double f_od(double N) { double r20571 = N; double r20572 = 5788.106078246906; bool r20573 = r20571 <= r20572; double r20574 = 1; double r20575 = r20571 + r20574; double r20576 = r20575 / r20571; double r20577 = sqrt(r20576); double r20578 = log(r20577); double r20579 = sqrt(r20577); double r20580 = log(r20579); double r20581 = r20580 + r20580; double r20582 = r20578 + r20581; double r20583 = r20574 / r20571; double r20584 = 1/3; double r20585 = r20584 / r20571; double r20586 = 1/2; double r20587 = r20585 - r20586; double r20588 = r20571 * r20571; double r20589 = r20587 / r20588; double r20590 = r20583 + r20589; double r20591 = r20573 ? r20582 : r20590; return r20591; } 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 r20592, r20593, r20594, r20595, r20596, r20597; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r20592); mpfr_init_set_str(r20593, "1", 10, MPFR_RNDN); mpfr_init(r20594); mpfr_init(r20595); mpfr_init(r20596); mpfr_init(r20597); } double f_im(double N) { mpfr_set_d(r20592, N, MPFR_RNDN); ; mpfr_add(r20594, r20592, r20593, MPFR_RNDN); mpfr_log(r20595, r20594, MPFR_RNDN); mpfr_log(r20596, r20592, MPFR_RNDN); mpfr_sub(r20597, r20595, r20596, MPFR_RNDN); return mpfr_get_d(r20597, MPFR_RNDN); } static mpfr_t r20598, r20599, r20600, r20601, r20602, r20603, r20604, r20605, r20606, r20607, r20608, r20609, r20610, r20611, r20612, r20613, r20614, r20615, r20616, r20617, r20618; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r20598); mpfr_init_set_str(r20599, "5788.106078246906", 10, MPFR_RNDN); mpfr_init(r20600); mpfr_init_set_str(r20601, "1", 10, MPFR_RNDN); mpfr_init(r20602); mpfr_init(r20603); mpfr_init(r20604); mpfr_init(r20605); mpfr_init(r20606); mpfr_init(r20607); mpfr_init(r20608); mpfr_init(r20609); mpfr_init(r20610); mpfr_init_set_str(r20611, "1/3", 10, MPFR_RNDN); mpfr_init(r20612); mpfr_init_set_str(r20613, "1/2", 10, MPFR_RNDN); mpfr_init(r20614); mpfr_init(r20615); mpfr_init(r20616); mpfr_init(r20617); mpfr_init(r20618); } double f_fm(double N) { mpfr_set_d(r20598, N, MPFR_RNDN); ; mpfr_set_si(r20600, mpfr_cmp(r20598, r20599) <= 0, MPFR_RNDN); ; mpfr_add(r20602, r20598, r20601, MPFR_RNDN); mpfr_div(r20603, r20602, r20598, MPFR_RNDN); mpfr_sqrt(r20604, r20603, MPFR_RNDN); mpfr_log(r20605, r20604, MPFR_RNDN); mpfr_sqrt(r20606, r20604, MPFR_RNDN); mpfr_log(r20607, r20606, MPFR_RNDN); mpfr_add(r20608, r20607, r20607, MPFR_RNDN); mpfr_add(r20609, r20605, r20608, MPFR_RNDN); mpfr_div(r20610, r20601, r20598, MPFR_RNDN); ; mpfr_div(r20612, r20611, r20598, MPFR_RNDN); ; mpfr_sub(r20614, r20612, r20613, MPFR_RNDN); mpfr_mul(r20615, r20598, r20598, MPFR_RNDN); mpfr_div(r20616, r20614, r20615, MPFR_RNDN); mpfr_add(r20617, r20610, r20616, MPFR_RNDN); if (mpfr_get_si(r20600, MPFR_RNDN)) { mpfr_set(r20618, r20609, MPFR_RNDN); } else { mpfr_set(r20618, r20617, MPFR_RNDN); }; return mpfr_get_d(r20618, MPFR_RNDN); } static mpfr_t r20619, r20620, r20621, r20622, r20623, r20624, r20625, r20626, r20627, r20628, r20629, r20630, r20631, r20632, r20633, r20634, r20635, r20636, r20637, r20638, r20639; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r20619); mpfr_init_set_str(r20620, "5788.106078246906", 10, MPFR_RNDN); mpfr_init(r20621); mpfr_init_set_str(r20622, "1", 10, MPFR_RNDN); mpfr_init(r20623); mpfr_init(r20624); mpfr_init(r20625); mpfr_init(r20626); mpfr_init(r20627); mpfr_init(r20628); mpfr_init(r20629); mpfr_init(r20630); mpfr_init(r20631); mpfr_init_set_str(r20632, "1/3", 10, MPFR_RNDN); mpfr_init(r20633); mpfr_init_set_str(r20634, "1/2", 10, MPFR_RNDN); mpfr_init(r20635); mpfr_init(r20636); mpfr_init(r20637); mpfr_init(r20638); mpfr_init(r20639); } double f_dm(double N) { mpfr_set_d(r20619, N, MPFR_RNDN); ; mpfr_set_si(r20621, mpfr_cmp(r20619, r20620) <= 0, MPFR_RNDN); ; mpfr_add(r20623, r20619, r20622, MPFR_RNDN); mpfr_div(r20624, r20623, r20619, MPFR_RNDN); mpfr_sqrt(r20625, r20624, MPFR_RNDN); mpfr_log(r20626, r20625, MPFR_RNDN); mpfr_sqrt(r20627, r20625, MPFR_RNDN); mpfr_log(r20628, r20627, MPFR_RNDN); mpfr_add(r20629, r20628, r20628, MPFR_RNDN); mpfr_add(r20630, r20626, r20629, MPFR_RNDN); mpfr_div(r20631, r20622, r20619, MPFR_RNDN); ; mpfr_div(r20633, r20632, r20619, MPFR_RNDN); ; mpfr_sub(r20635, r20633, r20634, MPFR_RNDN); mpfr_mul(r20636, r20619, r20619, MPFR_RNDN); mpfr_div(r20637, r20635, r20636, MPFR_RNDN); mpfr_add(r20638, r20631, r20637, MPFR_RNDN); if (mpfr_get_si(r20621, MPFR_RNDN)) { mpfr_set(r20639, r20630, MPFR_RNDN); } else { mpfr_set(r20639, r20638, MPFR_RNDN); }; return mpfr_get_d(r20639, MPFR_RNDN); }