#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r25594 = N; float r25595 = 1; float r25596 = r25594 + r25595; float r25597 = log(r25596); float r25598 = log(r25594); float r25599 = r25597 - r25598; return r25599; } double f_id(double N) { double r25600 = N; double r25601 = 1; double r25602 = r25600 + r25601; double r25603 = log(r25602); double r25604 = log(r25600); double r25605 = r25603 - r25604; return r25605; } double f_of(float N) { float r25606 = N; float r25607 = 1; float r25608 = r25606 + r25607; float r25609 = log(r25608); float r25610 = log(r25606); float r25611 = r25609 - r25610; float r25612 = 0.0029377285253137393; bool r25613 = r25611 <= r25612; float r25614 = 1/3; float r25615 = r25606 * r25606; float r25616 = r25614 / r25615; float r25617 = 1/2; float r25618 = r25617 / r25606; float r25619 = r25607 - r25618; float r25620 = r25616 + r25619; float r25621 = r25606 / r25620; float r25622 = r25607 / r25621; float r25623 = r25608 / r25606; float r25624 = log(r25623); float r25625 = r25613 ? r25622 : r25624; return r25625; } double f_od(double N) { double r25626 = N; double r25627 = 1; double r25628 = r25626 + r25627; double r25629 = log(r25628); double r25630 = log(r25626); double r25631 = r25629 - r25630; double r25632 = 0.0029377285253137393; bool r25633 = r25631 <= r25632; double r25634 = 1/3; double r25635 = r25626 * r25626; double r25636 = r25634 / r25635; double r25637 = 1/2; double r25638 = r25637 / r25626; double r25639 = r25627 - r25638; double r25640 = r25636 + r25639; double r25641 = r25626 / r25640; double r25642 = r25627 / r25641; double r25643 = r25628 / r25626; double r25644 = log(r25643); double r25645 = r25633 ? r25642 : r25644; return r25645; } 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 r25646, r25647, r25648, r25649, r25650, r25651; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r25646); mpfr_init_set_str(r25647, "1", 10, MPFR_RNDN); mpfr_init(r25648); mpfr_init(r25649); mpfr_init(r25650); mpfr_init(r25651); } double f_im(double N) { mpfr_set_d(r25646, N, MPFR_RNDN); ; mpfr_add(r25648, r25646, r25647, MPFR_RNDN); mpfr_log(r25649, r25648, MPFR_RNDN); mpfr_log(r25650, r25646, MPFR_RNDN); mpfr_sub(r25651, r25649, r25650, MPFR_RNDN); return mpfr_get_d(r25651, MPFR_RNDN); } static mpfr_t r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659, r25660, r25661, r25662, r25663, r25664, r25665, r25666, r25667, r25668, r25669, r25670, r25671; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r25652); mpfr_init_set_str(r25653, "1", 10, MPFR_RNDN); mpfr_init(r25654); mpfr_init(r25655); mpfr_init(r25656); mpfr_init(r25657); mpfr_init_set_str(r25658, "0.0029377285253137393", 10, MPFR_RNDN); mpfr_init(r25659); mpfr_init_set_str(r25660, "1/3", 10, MPFR_RNDN); mpfr_init(r25661); mpfr_init(r25662); mpfr_init_set_str(r25663, "1/2", 10, MPFR_RNDN); mpfr_init(r25664); mpfr_init(r25665); mpfr_init(r25666); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); mpfr_init(r25670); mpfr_init(r25671); } double f_fm(double N) { mpfr_set_d(r25652, N, MPFR_RNDN); ; mpfr_add(r25654, r25652, r25653, MPFR_RNDN); mpfr_log(r25655, r25654, MPFR_RNDN); mpfr_log(r25656, r25652, MPFR_RNDN); mpfr_sub(r25657, r25655, r25656, MPFR_RNDN); ; mpfr_set_si(r25659, mpfr_cmp(r25657, r25658) <= 0, MPFR_RNDN); ; mpfr_mul(r25661, r25652, r25652, MPFR_RNDN); mpfr_div(r25662, r25660, r25661, MPFR_RNDN); ; mpfr_div(r25664, r25663, r25652, MPFR_RNDN); mpfr_sub(r25665, r25653, r25664, MPFR_RNDN); mpfr_add(r25666, r25662, r25665, MPFR_RNDN); mpfr_div(r25667, r25652, r25666, MPFR_RNDN); mpfr_div(r25668, r25653, r25667, MPFR_RNDN); mpfr_div(r25669, r25654, r25652, MPFR_RNDN); mpfr_log(r25670, r25669, MPFR_RNDN); if (mpfr_get_si(r25659, MPFR_RNDN)) { mpfr_set(r25671, r25668, MPFR_RNDN); } else { mpfr_set(r25671, r25670, MPFR_RNDN); }; return mpfr_get_d(r25671, MPFR_RNDN); } static mpfr_t r25672, r25673, r25674, r25675, r25676, r25677, r25678, r25679, r25680, r25681, r25682, r25683, r25684, r25685, r25686, r25687, r25688, r25689, r25690, r25691; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r25672); mpfr_init_set_str(r25673, "1", 10, MPFR_RNDN); mpfr_init(r25674); mpfr_init(r25675); mpfr_init(r25676); mpfr_init(r25677); mpfr_init_set_str(r25678, "0.0029377285253137393", 10, MPFR_RNDN); mpfr_init(r25679); mpfr_init_set_str(r25680, "1/3", 10, MPFR_RNDN); mpfr_init(r25681); mpfr_init(r25682); mpfr_init_set_str(r25683, "1/2", 10, MPFR_RNDN); mpfr_init(r25684); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init(r25689); mpfr_init(r25690); mpfr_init(r25691); } double f_dm(double N) { mpfr_set_d(r25672, N, MPFR_RNDN); ; mpfr_add(r25674, r25672, r25673, MPFR_RNDN); mpfr_log(r25675, r25674, MPFR_RNDN); mpfr_log(r25676, r25672, MPFR_RNDN); mpfr_sub(r25677, r25675, r25676, MPFR_RNDN); ; mpfr_set_si(r25679, mpfr_cmp(r25677, r25678) <= 0, MPFR_RNDN); ; mpfr_mul(r25681, r25672, r25672, MPFR_RNDN); mpfr_div(r25682, r25680, r25681, MPFR_RNDN); ; mpfr_div(r25684, r25683, r25672, MPFR_RNDN); mpfr_sub(r25685, r25673, r25684, MPFR_RNDN); mpfr_add(r25686, r25682, r25685, MPFR_RNDN); mpfr_div(r25687, r25672, r25686, MPFR_RNDN); mpfr_div(r25688, r25673, r25687, MPFR_RNDN); mpfr_div(r25689, r25674, r25672, MPFR_RNDN); mpfr_log(r25690, r25689, MPFR_RNDN); if (mpfr_get_si(r25679, MPFR_RNDN)) { mpfr_set(r25691, r25688, MPFR_RNDN); } else { mpfr_set(r25691, r25690, MPFR_RNDN); }; return mpfr_get_d(r25691, MPFR_RNDN); }