#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r21139 = N; float r21140 = 1; float r21141 = r21139 + r21140; float r21142 = log(r21141); float r21143 = log(r21139); float r21144 = r21142 - r21143; return r21144; } double f_id(double N) { double r21145 = N; double r21146 = 1; double r21147 = r21145 + r21146; double r21148 = log(r21147); double r21149 = log(r21145); double r21150 = r21148 - r21149; return r21150; } double f_of(float N) { float r21151 = N; float r21152 = 1; float r21153 = r21151 + r21152; float r21154 = log(r21153); float r21155 = log(r21151); float r21156 = r21154 - r21155; float r21157 = 2.8879688583671737e-05; bool r21158 = r21156 <= r21157; float r21159 = r21152 / r21151; float r21160 = 1/3; float r21161 = r21160 / r21151; float r21162 = 1/2; float r21163 = r21161 - r21162; float r21164 = r21151 * r21151; float r21165 = r21152 / r21164; float r21166 = r21163 * r21165; float r21167 = r21159 + r21166; float r21168 = r21153 / r21151; float r21169 = log(r21168); float r21170 = r21158 ? r21167 : r21169; return r21170; } double f_od(double N) { double r21171 = N; double r21172 = 1; double r21173 = r21171 + r21172; double r21174 = log(r21173); double r21175 = log(r21171); double r21176 = r21174 - r21175; double r21177 = 2.8879688583671737e-05; bool r21178 = r21176 <= r21177; double r21179 = r21172 / r21171; double r21180 = 1/3; double r21181 = r21180 / r21171; double r21182 = 1/2; double r21183 = r21181 - r21182; double r21184 = r21171 * r21171; double r21185 = r21172 / r21184; double r21186 = r21183 * r21185; double r21187 = r21179 + r21186; double r21188 = r21173 / r21171; double r21189 = log(r21188); double r21190 = r21178 ? r21187 : r21189; return r21190; } 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 r21191, r21192, r21193, r21194, r21195, r21196; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r21191); mpfr_init_set_str(r21192, "1", 10, MPFR_RNDN); mpfr_init(r21193); mpfr_init(r21194); mpfr_init(r21195); mpfr_init(r21196); } double f_im(double N) { mpfr_set_d(r21191, N, MPFR_RNDN); ; mpfr_add(r21193, r21191, r21192, MPFR_RNDN); mpfr_log(r21194, r21193, MPFR_RNDN); mpfr_log(r21195, r21191, MPFR_RNDN); mpfr_sub(r21196, r21194, r21195, MPFR_RNDN); return mpfr_get_d(r21196, MPFR_RNDN); } static mpfr_t r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r21197); mpfr_init_set_str(r21198, "1", 10, MPFR_RNDN); mpfr_init(r21199); mpfr_init(r21200); mpfr_init(r21201); mpfr_init(r21202); mpfr_init_set_str(r21203, "2.8879688583671737e-05", 10, MPFR_RNDN); mpfr_init(r21204); mpfr_init(r21205); mpfr_init_set_str(r21206, "1/3", 10, MPFR_RNDN); mpfr_init(r21207); mpfr_init_set_str(r21208, "1/2", 10, MPFR_RNDN); mpfr_init(r21209); mpfr_init(r21210); mpfr_init(r21211); mpfr_init(r21212); mpfr_init(r21213); mpfr_init(r21214); mpfr_init(r21215); mpfr_init(r21216); } double f_fm(double N) { mpfr_set_d(r21197, N, MPFR_RNDN); ; mpfr_add(r21199, r21197, r21198, MPFR_RNDN); mpfr_log(r21200, r21199, MPFR_RNDN); mpfr_log(r21201, r21197, MPFR_RNDN); mpfr_sub(r21202, r21200, r21201, MPFR_RNDN); ; mpfr_set_si(r21204, mpfr_cmp(r21202, r21203) <= 0, MPFR_RNDN); mpfr_div(r21205, r21198, r21197, MPFR_RNDN); ; mpfr_div(r21207, r21206, r21197, MPFR_RNDN); ; mpfr_sub(r21209, r21207, r21208, MPFR_RNDN); mpfr_mul(r21210, r21197, r21197, MPFR_RNDN); mpfr_div(r21211, r21198, r21210, MPFR_RNDN); mpfr_mul(r21212, r21209, r21211, MPFR_RNDN); mpfr_add(r21213, r21205, r21212, MPFR_RNDN); mpfr_div(r21214, r21199, r21197, MPFR_RNDN); mpfr_log(r21215, r21214, MPFR_RNDN); if (mpfr_get_si(r21204, MPFR_RNDN)) { mpfr_set(r21216, r21213, MPFR_RNDN); } else { mpfr_set(r21216, r21215, MPFR_RNDN); }; return mpfr_get_d(r21216, MPFR_RNDN); } static mpfr_t r21217, r21218, r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229, r21230, r21231, r21232, r21233, r21234, r21235, r21236; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r21217); mpfr_init_set_str(r21218, "1", 10, MPFR_RNDN); mpfr_init(r21219); mpfr_init(r21220); mpfr_init(r21221); mpfr_init(r21222); mpfr_init_set_str(r21223, "2.8879688583671737e-05", 10, MPFR_RNDN); mpfr_init(r21224); mpfr_init(r21225); mpfr_init_set_str(r21226, "1/3", 10, MPFR_RNDN); mpfr_init(r21227); mpfr_init_set_str(r21228, "1/2", 10, MPFR_RNDN); mpfr_init(r21229); mpfr_init(r21230); mpfr_init(r21231); mpfr_init(r21232); mpfr_init(r21233); mpfr_init(r21234); mpfr_init(r21235); mpfr_init(r21236); } double f_dm(double N) { mpfr_set_d(r21217, N, MPFR_RNDN); ; mpfr_add(r21219, r21217, r21218, MPFR_RNDN); mpfr_log(r21220, r21219, MPFR_RNDN); mpfr_log(r21221, r21217, MPFR_RNDN); mpfr_sub(r21222, r21220, r21221, MPFR_RNDN); ; mpfr_set_si(r21224, mpfr_cmp(r21222, r21223) <= 0, MPFR_RNDN); mpfr_div(r21225, r21218, r21217, MPFR_RNDN); ; mpfr_div(r21227, r21226, r21217, MPFR_RNDN); ; mpfr_sub(r21229, r21227, r21228, MPFR_RNDN); mpfr_mul(r21230, r21217, r21217, MPFR_RNDN); mpfr_div(r21231, r21218, r21230, MPFR_RNDN); mpfr_mul(r21232, r21229, r21231, MPFR_RNDN); mpfr_add(r21233, r21225, r21232, MPFR_RNDN); mpfr_div(r21234, r21219, r21217, MPFR_RNDN); mpfr_log(r21235, r21234, MPFR_RNDN); if (mpfr_get_si(r21224, MPFR_RNDN)) { mpfr_set(r21236, r21233, MPFR_RNDN); } else { mpfr_set(r21236, r21235, MPFR_RNDN); }; return mpfr_get_d(r21236, MPFR_RNDN); }