#include #include #include #include #include char *name = "2nthrt (problem 3.4.6)"; double f_if(float x, float n) { float r21070 = x; float r21071 = 1; float r21072 = r21070 + r21071; float r21073 = n; float r21074 = r21071 / r21073; float r21075 = pow(r21072, r21074); float r21076 = pow(r21070, r21074); float r21077 = r21075 - r21076; return r21077; } double f_id(double x, double n) { double r21078 = x; double r21079 = 1; double r21080 = r21078 + r21079; double r21081 = n; double r21082 = r21079 / r21081; double r21083 = pow(r21080, r21082); double r21084 = pow(r21078, r21082); double r21085 = r21083 - r21084; return r21085; } double f_of(float x, float n) { float r21086 = 1; float r21087 = n; float r21088 = r21086 / r21087; float r21089 = -5.363094132516128e-12; bool r21090 = r21088 <= r21089; float r21091 = x; float r21092 = r21091 + r21086; float r21093 = pow(r21092, r21088); float r21094 = pow(r21091, r21088); float r21095 = r21093 - r21094; float r21096 = cbrt(r21095); float r21097 = r21096 * r21096; float r21098 = sqrt(r21092); float r21099 = pow(r21098, r21088); float r21100 = sqrt(r21094); float r21101 = r21099 + r21100; float r21102 = sqrt(r21098); float r21103 = pow(r21102, r21088); float r21104 = sqrt(r21100); float r21105 = r21103 + r21104; float r21106 = r21086 + r21091; float r21107 = sqrt(r21106); float r21108 = sqrt(r21107); float r21109 = pow(r21108, r21088); float r21110 = r21109 - r21104; float r21111 = exp(r21110); float r21112 = log(r21111); float r21113 = r21105 * r21112; float r21114 = r21101 * r21113; float r21115 = cbrt(r21114); float r21116 = r21097 * r21115; float r21117 = -1.4555907476059826e-116; bool r21118 = r21088 <= r21117; float r21119 = pow(r21106, r21088); float r21120 = r21119 - r21086; float r21121 = 1/2; float r21122 = r21121 / r21087; float r21123 = r21122 / r21087; float r21124 = log(r21091); float r21125 = r21124 * r21124; float r21126 = r21124 / r21087; float r21127 = fma(r21123, r21125, r21126); float r21128 = r21120 - r21127; float r21129 = -1.5275620307579182e-258; bool r21130 = r21088 <= r21129; float r21131 = r21091 * r21087; float r21132 = r21086 / r21131; float r21133 = r21091 * r21091; float r21134 = r21122 / r21133; float r21135 = r21132 - r21134; float r21136 = -r21124; float r21137 = r21087 * r21087; float r21138 = r21137 * r21091; float r21139 = r21136 / r21138; float r21140 = r21135 - r21139; float r21141 = 3.711468536197641e-292; bool r21142 = r21088 <= r21141; float r21143 = 1.5724217142242113e-266; bool r21144 = r21088 <= r21143; float r21145 = 6.762441816806913e-258; bool r21146 = r21088 <= r21145; float r21147 = 5.85591473108959e-110; bool r21148 = r21088 <= r21147; float r21149 = 3.594664189067294e-15; bool r21150 = r21088 <= r21149; float r21151 = log1p(r21091); float r21152 = r21151 / r21087; float r21153 = exp(r21152); float r21154 = r21153 - r21094; float r21155 = r21150 ? r21128 : r21154; float r21156 = r21148 ? r21140 : r21155; float r21157 = r21146 ? r21128 : r21156; float r21158 = r21144 ? r21140 : r21157; float r21159 = r21142 ? r21128 : r21158; float r21160 = r21130 ? r21140 : r21159; float r21161 = r21118 ? r21128 : r21160; float r21162 = r21090 ? r21116 : r21161; return r21162; } double f_od(double x, double n) { double r21163 = 1; double r21164 = n; double r21165 = r21163 / r21164; double r21166 = -5.363094132516128e-12; bool r21167 = r21165 <= r21166; double r21168 = x; double r21169 = r21168 + r21163; double r21170 = pow(r21169, r21165); double r21171 = pow(r21168, r21165); double r21172 = r21170 - r21171; double r21173 = cbrt(r21172); double r21174 = r21173 * r21173; double r21175 = sqrt(r21169); double r21176 = pow(r21175, r21165); double r21177 = sqrt(r21171); double r21178 = r21176 + r21177; double r21179 = sqrt(r21175); double r21180 = pow(r21179, r21165); double r21181 = sqrt(r21177); double r21182 = r21180 + r21181; double r21183 = r21163 + r21168; double r21184 = sqrt(r21183); double r21185 = sqrt(r21184); double r21186 = pow(r21185, r21165); double r21187 = r21186 - r21181; double r21188 = exp(r21187); double r21189 = log(r21188); double r21190 = r21182 * r21189; double r21191 = r21178 * r21190; double r21192 = cbrt(r21191); double r21193 = r21174 * r21192; double r21194 = -1.4555907476059826e-116; bool r21195 = r21165 <= r21194; double r21196 = pow(r21183, r21165); double r21197 = r21196 - r21163; double r21198 = 1/2; double r21199 = r21198 / r21164; double r21200 = r21199 / r21164; double r21201 = log(r21168); double r21202 = r21201 * r21201; double r21203 = r21201 / r21164; double r21204 = fma(r21200, r21202, r21203); double r21205 = r21197 - r21204; double r21206 = -1.5275620307579182e-258; bool r21207 = r21165 <= r21206; double r21208 = r21168 * r21164; double r21209 = r21163 / r21208; double r21210 = r21168 * r21168; double r21211 = r21199 / r21210; double r21212 = r21209 - r21211; double r21213 = -r21201; double r21214 = r21164 * r21164; double r21215 = r21214 * r21168; double r21216 = r21213 / r21215; double r21217 = r21212 - r21216; double r21218 = 3.711468536197641e-292; bool r21219 = r21165 <= r21218; double r21220 = 1.5724217142242113e-266; bool r21221 = r21165 <= r21220; double r21222 = 6.762441816806913e-258; bool r21223 = r21165 <= r21222; double r21224 = 5.85591473108959e-110; bool r21225 = r21165 <= r21224; double r21226 = 3.594664189067294e-15; bool r21227 = r21165 <= r21226; double r21228 = log1p(r21168); double r21229 = r21228 / r21164; double r21230 = exp(r21229); double r21231 = r21230 - r21171; double r21232 = r21227 ? r21205 : r21231; double r21233 = r21225 ? r21217 : r21232; double r21234 = r21223 ? r21205 : r21233; double r21235 = r21221 ? r21217 : r21234; double r21236 = r21219 ? r21205 : r21235; double r21237 = r21207 ? r21217 : r21236; double r21238 = r21195 ? r21205 : r21237; double r21239 = r21167 ? r21193 : r21238; return r21239; } 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 r21240, r21241, r21242, r21243, r21244, r21245, r21246, r21247; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r21240); mpfr_init_set_str(r21241, "1", 10, MPFR_RNDN); mpfr_init(r21242); mpfr_init(r21243); mpfr_init(r21244); mpfr_init(r21245); mpfr_init(r21246); mpfr_init(r21247); } double f_im(double x, double n) { mpfr_set_d(r21240, x, MPFR_RNDN); ; mpfr_add(r21242, r21240, r21241, MPFR_RNDN); mpfr_set_d(r21243, n, MPFR_RNDN); mpfr_div(r21244, r21241, r21243, MPFR_RNDN); mpfr_pow(r21245, r21242, r21244, MPFR_RNDN); mpfr_pow(r21246, r21240, r21244, MPFR_RNDN); mpfr_sub(r21247, r21245, r21246, MPFR_RNDN); return mpfr_get_d(r21247, MPFR_RNDN); } static mpfr_t r21248, r21249, r21250, r21251, r21252, r21253, r21254, r21255, r21256, r21257, r21258, r21259, r21260, r21261, r21262, r21263, r21264, r21265, r21266, r21267, r21268, r21269, r21270, r21271, r21272, r21273, r21274, r21275, r21276, r21277, r21278, r21279, r21280, r21281, r21282, r21283, r21284, r21285, r21286, r21287, r21288, r21289, r21290, r21291, r21292, r21293, r21294, r21295, r21296, r21297, r21298, r21299, r21300, r21301, r21302, r21303, r21304, r21305, r21306, r21307, r21308, r21309, r21310, r21311, r21312, r21313, r21314, r21315, r21316, r21317, r21318, r21319, r21320, r21321, r21322, r21323, r21324; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21248, "1", 10, MPFR_RNDN); mpfr_init(r21249); mpfr_init(r21250); mpfr_init_set_str(r21251, "-5.363094132516128e-12", 10, MPFR_RNDN); mpfr_init(r21252); mpfr_init(r21253); mpfr_init(r21254); mpfr_init(r21255); mpfr_init(r21256); mpfr_init(r21257); mpfr_init(r21258); mpfr_init(r21259); mpfr_init(r21260); mpfr_init(r21261); mpfr_init(r21262); mpfr_init(r21263); mpfr_init(r21264); mpfr_init(r21265); mpfr_init(r21266); mpfr_init(r21267); mpfr_init(r21268); mpfr_init(r21269); mpfr_init(r21270); mpfr_init(r21271); mpfr_init(r21272); mpfr_init(r21273); mpfr_init(r21274); mpfr_init(r21275); mpfr_init(r21276); mpfr_init(r21277); mpfr_init(r21278); mpfr_init_set_str(r21279, "-1.4555907476059826e-116", 10, MPFR_RNDN); mpfr_init(r21280); mpfr_init(r21281); mpfr_init(r21282); mpfr_init_set_str(r21283, "1/2", 10, MPFR_RNDN); mpfr_init(r21284); mpfr_init(r21285); mpfr_init(r21286); mpfr_init(r21287); mpfr_init(r21288); mpfr_init(r21289); mpfr_init(r21290); mpfr_init_set_str(r21291, "-1.5275620307579182e-258", 10, MPFR_RNDN); mpfr_init(r21292); mpfr_init(r21293); mpfr_init(r21294); mpfr_init(r21295); mpfr_init(r21296); mpfr_init(r21297); mpfr_init(r21298); mpfr_init(r21299); mpfr_init(r21300); mpfr_init(r21301); mpfr_init(r21302); mpfr_init_set_str(r21303, "3.711468536197641e-292", 10, MPFR_RNDN); mpfr_init(r21304); mpfr_init_set_str(r21305, "1.5724217142242113e-266", 10, MPFR_RNDN); mpfr_init(r21306); mpfr_init_set_str(r21307, "6.762441816806913e-258", 10, MPFR_RNDN); mpfr_init(r21308); mpfr_init_set_str(r21309, "5.85591473108959e-110", 10, MPFR_RNDN); mpfr_init(r21310); mpfr_init_set_str(r21311, "3.594664189067294e-15", 10, MPFR_RNDN); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init(r21315); mpfr_init(r21316); mpfr_init(r21317); mpfr_init(r21318); mpfr_init(r21319); mpfr_init(r21320); mpfr_init(r21321); mpfr_init(r21322); mpfr_init(r21323); mpfr_init(r21324); } double f_fm(double x, double n) { ; mpfr_set_d(r21249, n, MPFR_RNDN); mpfr_div(r21250, r21248, r21249, MPFR_RNDN); ; mpfr_set_si(r21252, mpfr_cmp(r21250, r21251) <= 0, MPFR_RNDN); mpfr_set_d(r21253, x, MPFR_RNDN); mpfr_add(r21254, r21253, r21248, MPFR_RNDN); mpfr_pow(r21255, r21254, r21250, MPFR_RNDN); mpfr_pow(r21256, r21253, r21250, MPFR_RNDN); mpfr_sub(r21257, r21255, r21256, MPFR_RNDN); mpfr_cbrt(r21258, r21257, MPFR_RNDN); mpfr_mul(r21259, r21258, r21258, MPFR_RNDN); mpfr_sqrt(r21260, r21254, MPFR_RNDN); mpfr_pow(r21261, r21260, r21250, MPFR_RNDN); mpfr_sqrt(r21262, r21256, MPFR_RNDN); mpfr_add(r21263, r21261, r21262, MPFR_RNDN); mpfr_sqrt(r21264, r21260, MPFR_RNDN); mpfr_pow(r21265, r21264, r21250, MPFR_RNDN); mpfr_sqrt(r21266, r21262, MPFR_RNDN); mpfr_add(r21267, r21265, r21266, MPFR_RNDN); mpfr_add(r21268, r21248, r21253, MPFR_RNDN); mpfr_sqrt(r21269, r21268, MPFR_RNDN); mpfr_sqrt(r21270, r21269, MPFR_RNDN); mpfr_pow(r21271, r21270, r21250, MPFR_RNDN); mpfr_sub(r21272, r21271, r21266, MPFR_RNDN); mpfr_exp(r21273, r21272, MPFR_RNDN); mpfr_log(r21274, r21273, MPFR_RNDN); mpfr_mul(r21275, r21267, r21274, MPFR_RNDN); mpfr_mul(r21276, r21263, r21275, MPFR_RNDN); mpfr_cbrt(r21277, r21276, MPFR_RNDN); mpfr_mul(r21278, r21259, r21277, MPFR_RNDN); ; mpfr_set_si(r21280, mpfr_cmp(r21250, r21279) <= 0, MPFR_RNDN); mpfr_pow(r21281, r21268, r21250, MPFR_RNDN); mpfr_sub(r21282, r21281, r21248, MPFR_RNDN); ; mpfr_div(r21284, r21283, r21249, MPFR_RNDN); mpfr_div(r21285, r21284, r21249, MPFR_RNDN); mpfr_log(r21286, r21253, MPFR_RNDN); mpfr_mul(r21287, r21286, r21286, MPFR_RNDN); mpfr_div(r21288, r21286, r21249, MPFR_RNDN); mpfr_fma(r21289, r21285, r21287, r21288, MPFR_RNDN); mpfr_sub(r21290, r21282, r21289, MPFR_RNDN); ; mpfr_set_si(r21292, mpfr_cmp(r21250, r21291) <= 0, MPFR_RNDN); mpfr_mul(r21293, r21253, r21249, MPFR_RNDN); mpfr_div(r21294, r21248, r21293, MPFR_RNDN); mpfr_mul(r21295, r21253, r21253, MPFR_RNDN); mpfr_div(r21296, r21284, r21295, MPFR_RNDN); mpfr_sub(r21297, r21294, r21296, MPFR_RNDN); mpfr_neg(r21298, r21286, MPFR_RNDN); mpfr_mul(r21299, r21249, r21249, MPFR_RNDN); mpfr_mul(r21300, r21299, r21253, MPFR_RNDN); mpfr_div(r21301, r21298, r21300, MPFR_RNDN); mpfr_sub(r21302, r21297, r21301, MPFR_RNDN); ; mpfr_set_si(r21304, mpfr_cmp(r21250, r21303) <= 0, MPFR_RNDN); ; mpfr_set_si(r21306, mpfr_cmp(r21250, r21305) <= 0, MPFR_RNDN); ; mpfr_set_si(r21308, mpfr_cmp(r21250, r21307) <= 0, MPFR_RNDN); ; mpfr_set_si(r21310, mpfr_cmp(r21250, r21309) <= 0, MPFR_RNDN); ; mpfr_set_si(r21312, mpfr_cmp(r21250, r21311) <= 0, MPFR_RNDN); mpfr_log1p(r21313, r21253, MPFR_RNDN); mpfr_div(r21314, r21313, r21249, MPFR_RNDN); mpfr_exp(r21315, r21314, MPFR_RNDN); mpfr_sub(r21316, r21315, r21256, MPFR_RNDN); if (mpfr_get_si(r21312, MPFR_RNDN)) { mpfr_set(r21317, r21290, MPFR_RNDN); } else { mpfr_set(r21317, r21316, MPFR_RNDN); }; if (mpfr_get_si(r21310, MPFR_RNDN)) { mpfr_set(r21318, r21302, MPFR_RNDN); } else { mpfr_set(r21318, r21317, MPFR_RNDN); }; if (mpfr_get_si(r21308, MPFR_RNDN)) { mpfr_set(r21319, r21290, MPFR_RNDN); } else { mpfr_set(r21319, r21318, MPFR_RNDN); }; if (mpfr_get_si(r21306, MPFR_RNDN)) { mpfr_set(r21320, r21302, MPFR_RNDN); } else { mpfr_set(r21320, r21319, MPFR_RNDN); }; if (mpfr_get_si(r21304, MPFR_RNDN)) { mpfr_set(r21321, r21290, MPFR_RNDN); } else { mpfr_set(r21321, r21320, MPFR_RNDN); }; if (mpfr_get_si(r21292, MPFR_RNDN)) { mpfr_set(r21322, r21302, MPFR_RNDN); } else { mpfr_set(r21322, r21321, MPFR_RNDN); }; if (mpfr_get_si(r21280, MPFR_RNDN)) { mpfr_set(r21323, r21290, MPFR_RNDN); } else { mpfr_set(r21323, r21322, MPFR_RNDN); }; if (mpfr_get_si(r21252, MPFR_RNDN)) { mpfr_set(r21324, r21278, MPFR_RNDN); } else { mpfr_set(r21324, r21323, MPFR_RNDN); }; return mpfr_get_d(r21324, MPFR_RNDN); } static mpfr_t r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344, r21345, r21346, r21347, r21348, r21349, r21350, r21351, r21352, r21353, r21354, r21355, r21356, r21357, r21358, r21359, r21360, r21361, r21362, r21363, r21364, r21365, r21366, r21367, r21368, r21369, r21370, r21371, r21372, r21373, r21374, r21375, r21376, r21377, r21378, r21379, r21380, r21381, r21382, r21383, r21384, r21385, r21386, r21387, r21388, r21389, r21390, r21391, r21392, r21393, r21394, r21395, r21396, r21397, r21398, r21399, r21400, r21401; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21325, "1", 10, MPFR_RNDN); mpfr_init(r21326); mpfr_init(r21327); mpfr_init_set_str(r21328, "-5.363094132516128e-12", 10, MPFR_RNDN); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init(r21332); mpfr_init(r21333); mpfr_init(r21334); mpfr_init(r21335); mpfr_init(r21336); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); mpfr_init(r21345); mpfr_init(r21346); mpfr_init(r21347); mpfr_init(r21348); mpfr_init(r21349); mpfr_init(r21350); mpfr_init(r21351); mpfr_init(r21352); mpfr_init(r21353); mpfr_init(r21354); mpfr_init(r21355); mpfr_init_set_str(r21356, "-1.4555907476059826e-116", 10, MPFR_RNDN); mpfr_init(r21357); mpfr_init(r21358); mpfr_init(r21359); mpfr_init_set_str(r21360, "1/2", 10, MPFR_RNDN); mpfr_init(r21361); mpfr_init(r21362); mpfr_init(r21363); mpfr_init(r21364); mpfr_init(r21365); mpfr_init(r21366); mpfr_init(r21367); mpfr_init_set_str(r21368, "-1.5275620307579182e-258", 10, MPFR_RNDN); mpfr_init(r21369); mpfr_init(r21370); mpfr_init(r21371); mpfr_init(r21372); mpfr_init(r21373); mpfr_init(r21374); mpfr_init(r21375); mpfr_init(r21376); mpfr_init(r21377); mpfr_init(r21378); mpfr_init(r21379); mpfr_init_set_str(r21380, "3.711468536197641e-292", 10, MPFR_RNDN); mpfr_init(r21381); mpfr_init_set_str(r21382, "1.5724217142242113e-266", 10, MPFR_RNDN); mpfr_init(r21383); mpfr_init_set_str(r21384, "6.762441816806913e-258", 10, MPFR_RNDN); mpfr_init(r21385); mpfr_init_set_str(r21386, "5.85591473108959e-110", 10, MPFR_RNDN); mpfr_init(r21387); mpfr_init_set_str(r21388, "3.594664189067294e-15", 10, MPFR_RNDN); mpfr_init(r21389); mpfr_init(r21390); mpfr_init(r21391); mpfr_init(r21392); mpfr_init(r21393); mpfr_init(r21394); mpfr_init(r21395); mpfr_init(r21396); mpfr_init(r21397); mpfr_init(r21398); mpfr_init(r21399); mpfr_init(r21400); mpfr_init(r21401); } double f_dm(double x, double n) { ; mpfr_set_d(r21326, n, MPFR_RNDN); mpfr_div(r21327, r21325, r21326, MPFR_RNDN); ; mpfr_set_si(r21329, mpfr_cmp(r21327, r21328) <= 0, MPFR_RNDN); mpfr_set_d(r21330, x, MPFR_RNDN); mpfr_add(r21331, r21330, r21325, MPFR_RNDN); mpfr_pow(r21332, r21331, r21327, MPFR_RNDN); mpfr_pow(r21333, r21330, r21327, MPFR_RNDN); mpfr_sub(r21334, r21332, r21333, MPFR_RNDN); mpfr_cbrt(r21335, r21334, MPFR_RNDN); mpfr_mul(r21336, r21335, r21335, MPFR_RNDN); mpfr_sqrt(r21337, r21331, MPFR_RNDN); mpfr_pow(r21338, r21337, r21327, MPFR_RNDN); mpfr_sqrt(r21339, r21333, MPFR_RNDN); mpfr_add(r21340, r21338, r21339, MPFR_RNDN); mpfr_sqrt(r21341, r21337, MPFR_RNDN); mpfr_pow(r21342, r21341, r21327, MPFR_RNDN); mpfr_sqrt(r21343, r21339, MPFR_RNDN); mpfr_add(r21344, r21342, r21343, MPFR_RNDN); mpfr_add(r21345, r21325, r21330, MPFR_RNDN); mpfr_sqrt(r21346, r21345, MPFR_RNDN); mpfr_sqrt(r21347, r21346, MPFR_RNDN); mpfr_pow(r21348, r21347, r21327, MPFR_RNDN); mpfr_sub(r21349, r21348, r21343, MPFR_RNDN); mpfr_exp(r21350, r21349, MPFR_RNDN); mpfr_log(r21351, r21350, MPFR_RNDN); mpfr_mul(r21352, r21344, r21351, MPFR_RNDN); mpfr_mul(r21353, r21340, r21352, MPFR_RNDN); mpfr_cbrt(r21354, r21353, MPFR_RNDN); mpfr_mul(r21355, r21336, r21354, MPFR_RNDN); ; mpfr_set_si(r21357, mpfr_cmp(r21327, r21356) <= 0, MPFR_RNDN); mpfr_pow(r21358, r21345, r21327, MPFR_RNDN); mpfr_sub(r21359, r21358, r21325, MPFR_RNDN); ; mpfr_div(r21361, r21360, r21326, MPFR_RNDN); mpfr_div(r21362, r21361, r21326, MPFR_RNDN); mpfr_log(r21363, r21330, MPFR_RNDN); mpfr_mul(r21364, r21363, r21363, MPFR_RNDN); mpfr_div(r21365, r21363, r21326, MPFR_RNDN); mpfr_fma(r21366, r21362, r21364, r21365, MPFR_RNDN); mpfr_sub(r21367, r21359, r21366, MPFR_RNDN); ; mpfr_set_si(r21369, mpfr_cmp(r21327, r21368) <= 0, MPFR_RNDN); mpfr_mul(r21370, r21330, r21326, MPFR_RNDN); mpfr_div(r21371, r21325, r21370, MPFR_RNDN); mpfr_mul(r21372, r21330, r21330, MPFR_RNDN); mpfr_div(r21373, r21361, r21372, MPFR_RNDN); mpfr_sub(r21374, r21371, r21373, MPFR_RNDN); mpfr_neg(r21375, r21363, MPFR_RNDN); mpfr_mul(r21376, r21326, r21326, MPFR_RNDN); mpfr_mul(r21377, r21376, r21330, MPFR_RNDN); mpfr_div(r21378, r21375, r21377, MPFR_RNDN); mpfr_sub(r21379, r21374, r21378, MPFR_RNDN); ; mpfr_set_si(r21381, mpfr_cmp(r21327, r21380) <= 0, MPFR_RNDN); ; mpfr_set_si(r21383, mpfr_cmp(r21327, r21382) <= 0, MPFR_RNDN); ; mpfr_set_si(r21385, mpfr_cmp(r21327, r21384) <= 0, MPFR_RNDN); ; mpfr_set_si(r21387, mpfr_cmp(r21327, r21386) <= 0, MPFR_RNDN); ; mpfr_set_si(r21389, mpfr_cmp(r21327, r21388) <= 0, MPFR_RNDN); mpfr_log1p(r21390, r21330, MPFR_RNDN); mpfr_div(r21391, r21390, r21326, MPFR_RNDN); mpfr_exp(r21392, r21391, MPFR_RNDN); mpfr_sub(r21393, r21392, r21333, MPFR_RNDN); if (mpfr_get_si(r21389, MPFR_RNDN)) { mpfr_set(r21394, r21367, MPFR_RNDN); } else { mpfr_set(r21394, r21393, MPFR_RNDN); }; if (mpfr_get_si(r21387, MPFR_RNDN)) { mpfr_set(r21395, r21379, MPFR_RNDN); } else { mpfr_set(r21395, r21394, MPFR_RNDN); }; if (mpfr_get_si(r21385, MPFR_RNDN)) { mpfr_set(r21396, r21367, MPFR_RNDN); } else { mpfr_set(r21396, r21395, MPFR_RNDN); }; if (mpfr_get_si(r21383, MPFR_RNDN)) { mpfr_set(r21397, r21379, MPFR_RNDN); } else { mpfr_set(r21397, r21396, MPFR_RNDN); }; if (mpfr_get_si(r21381, MPFR_RNDN)) { mpfr_set(r21398, r21367, MPFR_RNDN); } else { mpfr_set(r21398, r21397, MPFR_RNDN); }; if (mpfr_get_si(r21369, MPFR_RNDN)) { mpfr_set(r21399, r21379, MPFR_RNDN); } else { mpfr_set(r21399, r21398, MPFR_RNDN); }; if (mpfr_get_si(r21357, MPFR_RNDN)) { mpfr_set(r21400, r21367, MPFR_RNDN); } else { mpfr_set(r21400, r21399, MPFR_RNDN); }; if (mpfr_get_si(r21329, MPFR_RNDN)) { mpfr_set(r21401, r21355, MPFR_RNDN); } else { mpfr_set(r21401, r21400, MPFR_RNDN); }; return mpfr_get_d(r21401, MPFR_RNDN); }