#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r9528 = N; float r9529 = 1; float r9530 = r9528 + r9529; float r9531 = log(r9530); float r9532 = log(r9528); float r9533 = r9531 - r9532; return r9533; } double f_id(double N) { double r9534 = N; double r9535 = 1; double r9536 = r9534 + r9535; double r9537 = log(r9536); double r9538 = log(r9534); double r9539 = r9537 - r9538; return r9539; } double f_of(float N) { float r9540 = N; float r9541 = 1; float r9542 = r9540 + r9541; float r9543 = log(r9542); float r9544 = log(r9540); float r9545 = r9543 - r9544; float r9546 = 4.798456013914258e-05; bool r9547 = r9545 <= r9546; float r9548 = 1/3; float r9549 = r9548 / r9540; float r9550 = 1/2; float r9551 = r9549 - r9550; float r9552 = r9540 * r9540; float r9553 = r9541 / r9552; float r9554 = r9541 / r9540; float r9555 = fma(r9551, r9553, r9554); float r9556 = log1p(r9540); float r9557 = sqrt(r9556); float r9558 = -r9544; float r9559 = fma(r9557, r9557, r9558); float r9560 = r9547 ? r9555 : r9559; return r9560; } double f_od(double N) { double r9561 = N; double r9562 = 1; double r9563 = r9561 + r9562; double r9564 = log(r9563); double r9565 = log(r9561); double r9566 = r9564 - r9565; double r9567 = 4.798456013914258e-05; bool r9568 = r9566 <= r9567; double r9569 = 1/3; double r9570 = r9569 / r9561; double r9571 = 1/2; double r9572 = r9570 - r9571; double r9573 = r9561 * r9561; double r9574 = r9562 / r9573; double r9575 = r9562 / r9561; double r9576 = fma(r9572, r9574, r9575); double r9577 = log1p(r9561); double r9578 = sqrt(r9577); double r9579 = -r9565; double r9580 = fma(r9578, r9578, r9579); double r9581 = r9568 ? r9576 : r9580; return r9581; } 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 r9582, r9583, r9584, r9585, r9586, r9587; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r9582); mpfr_init_set_str(r9583, "1", 10, MPFR_RNDN); mpfr_init(r9584); mpfr_init(r9585); mpfr_init(r9586); mpfr_init(r9587); } double f_im(double N) { mpfr_set_d(r9582, N, MPFR_RNDN); ; mpfr_add(r9584, r9582, r9583, MPFR_RNDN); mpfr_log(r9585, r9584, MPFR_RNDN); mpfr_log(r9586, r9582, MPFR_RNDN); mpfr_sub(r9587, r9585, r9586, MPFR_RNDN); return mpfr_get_d(r9587, MPFR_RNDN); } static mpfr_t r9588, r9589, r9590, r9591, r9592, r9593, r9594, r9595, r9596, r9597, r9598, r9599, r9600, r9601, r9602, r9603, r9604, r9605, r9606, r9607, r9608; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r9588); mpfr_init_set_str(r9589, "1", 10, MPFR_RNDN); mpfr_init(r9590); mpfr_init(r9591); mpfr_init(r9592); mpfr_init(r9593); mpfr_init_set_str(r9594, "4.798456013914258e-05", 10, MPFR_RNDN); mpfr_init(r9595); mpfr_init_set_str(r9596, "1/3", 10, MPFR_RNDN); mpfr_init(r9597); mpfr_init_set_str(r9598, "1/2", 10, MPFR_RNDN); mpfr_init(r9599); mpfr_init(r9600); mpfr_init(r9601); mpfr_init(r9602); mpfr_init(r9603); mpfr_init(r9604); mpfr_init(r9605); mpfr_init(r9606); mpfr_init(r9607); mpfr_init(r9608); } double f_fm(double N) { mpfr_set_d(r9588, N, MPFR_RNDN); ; mpfr_add(r9590, r9588, r9589, MPFR_RNDN); mpfr_log(r9591, r9590, MPFR_RNDN); mpfr_log(r9592, r9588, MPFR_RNDN); mpfr_sub(r9593, r9591, r9592, MPFR_RNDN); ; mpfr_set_si(r9595, mpfr_cmp(r9593, r9594) <= 0, MPFR_RNDN); ; mpfr_div(r9597, r9596, r9588, MPFR_RNDN); ; mpfr_sub(r9599, r9597, r9598, MPFR_RNDN); mpfr_mul(r9600, r9588, r9588, MPFR_RNDN); mpfr_div(r9601, r9589, r9600, MPFR_RNDN); mpfr_div(r9602, r9589, r9588, MPFR_RNDN); mpfr_fma(r9603, r9599, r9601, r9602, MPFR_RNDN); mpfr_log1p(r9604, r9588, MPFR_RNDN); mpfr_sqrt(r9605, r9604, MPFR_RNDN); mpfr_neg(r9606, r9592, MPFR_RNDN); mpfr_fma(r9607, r9605, r9605, r9606, MPFR_RNDN); if (mpfr_get_si(r9595, MPFR_RNDN)) { mpfr_set(r9608, r9603, MPFR_RNDN); } else { mpfr_set(r9608, r9607, MPFR_RNDN); }; return mpfr_get_d(r9608, MPFR_RNDN); } static mpfr_t r9609, r9610, r9611, r9612, r9613, r9614, r9615, r9616, r9617, r9618, r9619, r9620, r9621, r9622, r9623, r9624, r9625, r9626, r9627, r9628, r9629; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r9609); mpfr_init_set_str(r9610, "1", 10, MPFR_RNDN); mpfr_init(r9611); mpfr_init(r9612); mpfr_init(r9613); mpfr_init(r9614); mpfr_init_set_str(r9615, "4.798456013914258e-05", 10, MPFR_RNDN); mpfr_init(r9616); mpfr_init_set_str(r9617, "1/3", 10, MPFR_RNDN); mpfr_init(r9618); mpfr_init_set_str(r9619, "1/2", 10, MPFR_RNDN); mpfr_init(r9620); mpfr_init(r9621); mpfr_init(r9622); mpfr_init(r9623); mpfr_init(r9624); mpfr_init(r9625); mpfr_init(r9626); mpfr_init(r9627); mpfr_init(r9628); mpfr_init(r9629); } double f_dm(double N) { mpfr_set_d(r9609, N, MPFR_RNDN); ; mpfr_add(r9611, r9609, r9610, MPFR_RNDN); mpfr_log(r9612, r9611, MPFR_RNDN); mpfr_log(r9613, r9609, MPFR_RNDN); mpfr_sub(r9614, r9612, r9613, MPFR_RNDN); ; mpfr_set_si(r9616, mpfr_cmp(r9614, r9615) <= 0, MPFR_RNDN); ; mpfr_div(r9618, r9617, r9609, MPFR_RNDN); ; mpfr_sub(r9620, r9618, r9619, MPFR_RNDN); mpfr_mul(r9621, r9609, r9609, MPFR_RNDN); mpfr_div(r9622, r9610, r9621, MPFR_RNDN); mpfr_div(r9623, r9610, r9609, MPFR_RNDN); mpfr_fma(r9624, r9620, r9622, r9623, MPFR_RNDN); mpfr_log1p(r9625, r9609, MPFR_RNDN); mpfr_sqrt(r9626, r9625, MPFR_RNDN); mpfr_neg(r9627, r9613, MPFR_RNDN); mpfr_fma(r9628, r9626, r9626, r9627, MPFR_RNDN); if (mpfr_get_si(r9616, MPFR_RNDN)) { mpfr_set(r9629, r9624, MPFR_RNDN); } else { mpfr_set(r9629, r9628, MPFR_RNDN); }; return mpfr_get_d(r9629, MPFR_RNDN); }