#include #include #include #include #include char *name = "Toniolo and Linder, Equation (7)"; double f_if(float x, float l, float t) { float r21003 = 2; float r21004 = sqrt(r21003); float r21005 = t; float r21006 = r21004 * r21005; float r21007 = x; float r21008 = 1; float r21009 = r21007 + r21008; float r21010 = r21007 - r21008; float r21011 = r21009 / r21010; float r21012 = l; float r21013 = r21012 * r21012; float r21014 = r21005 * r21005; float r21015 = r21003 * r21014; float r21016 = r21013 + r21015; float r21017 = r21011 * r21016; float r21018 = r21017 - r21013; float r21019 = sqrt(r21018); float r21020 = r21006 / r21019; return r21020; } double f_id(double x, double l, double t) { double r21021 = 2; double r21022 = sqrt(r21021); double r21023 = t; double r21024 = r21022 * r21023; double r21025 = x; double r21026 = 1; double r21027 = r21025 + r21026; double r21028 = r21025 - r21026; double r21029 = r21027 / r21028; double r21030 = l; double r21031 = r21030 * r21030; double r21032 = r21023 * r21023; double r21033 = r21021 * r21032; double r21034 = r21031 + r21033; double r21035 = r21029 * r21034; double r21036 = r21035 - r21031; double r21037 = sqrt(r21036); double r21038 = r21024 / r21037; return r21038; } double f_of(float x, float l, float t) { float r21039 = t; float r21040 = -1.42998038934624e+69; bool r21041 = r21039 <= r21040; float r21042 = 2; float r21043 = sqrt(r21042); float r21044 = r21039 * r21043; float r21045 = x; float r21046 = r21039 / r21045; float r21047 = r21046 / r21045; float r21048 = r21047 / r21043; float r21049 = r21048 - r21044; float r21050 = r21039 + r21039; float r21051 = r21050 / r21043; float r21052 = r21051 / r21045; float r21053 = r21045 * r21045; float r21054 = r21051 / r21053; float r21055 = r21052 + r21054; float r21056 = r21049 - r21055; float r21057 = r21044 / r21056; float r21058 = -2.5401539047568627e-163; bool r21059 = r21039 <= r21058; float r21060 = cbrt(r21043); float r21061 = r21060 * r21060; float r21062 = sqrt(r21061); float r21063 = sqrt(r21060); float r21064 = sqrt(r21043); float r21065 = r21064 * r21039; float r21066 = r21063 * r21065; float r21067 = r21062 * r21066; float r21068 = pow(r21039, r21042); float r21069 = r21042 * r21068; float r21070 = l; float r21071 = r21045 / r21070; float r21072 = r21070 / r21071; float r21073 = r21042 * r21072; float r21074 = 4; float r21075 = r21068 / r21045; float r21076 = r21074 * r21075; float r21077 = r21073 + r21076; float r21078 = r21069 + r21077; float r21079 = sqrt(r21078); float r21080 = r21067 / r21079; float r21081 = -5.553561555528792e-295; bool r21082 = r21039 <= r21081; float r21083 = 1.2691519409648163e-243; bool r21084 = r21039 <= r21083; float r21085 = r21064 * r21065; float r21086 = r21085 / r21079; float r21087 = 7.679813441312647e-161; bool r21088 = r21039 <= r21087; float r21089 = r21042 / r21045; float r21090 = r21089 / r21043; float r21091 = r21043 + r21090; float r21092 = r21039 * r21091; float r21093 = r21042 / r21043; float r21094 = 1; float r21095 = r21094 / r21043; float r21096 = r21093 - r21095; float r21097 = r21047 * r21096; float r21098 = r21092 + r21097; float r21099 = r21044 / r21098; float r21100 = 4.634093238500572e+140; bool r21101 = r21039 <= r21100; float r21102 = +inf.0; bool r21103 = r21039 <= r21102; float r21104 = r21103 ? r21099 : r21099; float r21105 = r21101 ? r21080 : r21104; float r21106 = r21088 ? r21099 : r21105; float r21107 = r21084 ? r21086 : r21106; float r21108 = r21082 ? r21057 : r21107; float r21109 = r21059 ? r21080 : r21108; float r21110 = r21041 ? r21057 : r21109; return r21110; } double f_od(double x, double l, double t) { double r21111 = t; double r21112 = -1.42998038934624e+69; bool r21113 = r21111 <= r21112; double r21114 = 2; double r21115 = sqrt(r21114); double r21116 = r21111 * r21115; double r21117 = x; double r21118 = r21111 / r21117; double r21119 = r21118 / r21117; double r21120 = r21119 / r21115; double r21121 = r21120 - r21116; double r21122 = r21111 + r21111; double r21123 = r21122 / r21115; double r21124 = r21123 / r21117; double r21125 = r21117 * r21117; double r21126 = r21123 / r21125; double r21127 = r21124 + r21126; double r21128 = r21121 - r21127; double r21129 = r21116 / r21128; double r21130 = -2.5401539047568627e-163; bool r21131 = r21111 <= r21130; double r21132 = cbrt(r21115); double r21133 = r21132 * r21132; double r21134 = sqrt(r21133); double r21135 = sqrt(r21132); double r21136 = sqrt(r21115); double r21137 = r21136 * r21111; double r21138 = r21135 * r21137; double r21139 = r21134 * r21138; double r21140 = pow(r21111, r21114); double r21141 = r21114 * r21140; double r21142 = l; double r21143 = r21117 / r21142; double r21144 = r21142 / r21143; double r21145 = r21114 * r21144; double r21146 = 4; double r21147 = r21140 / r21117; double r21148 = r21146 * r21147; double r21149 = r21145 + r21148; double r21150 = r21141 + r21149; double r21151 = sqrt(r21150); double r21152 = r21139 / r21151; double r21153 = -5.553561555528792e-295; bool r21154 = r21111 <= r21153; double r21155 = 1.2691519409648163e-243; bool r21156 = r21111 <= r21155; double r21157 = r21136 * r21137; double r21158 = r21157 / r21151; double r21159 = 7.679813441312647e-161; bool r21160 = r21111 <= r21159; double r21161 = r21114 / r21117; double r21162 = r21161 / r21115; double r21163 = r21115 + r21162; double r21164 = r21111 * r21163; double r21165 = r21114 / r21115; double r21166 = 1; double r21167 = r21166 / r21115; double r21168 = r21165 - r21167; double r21169 = r21119 * r21168; double r21170 = r21164 + r21169; double r21171 = r21116 / r21170; double r21172 = 4.634093238500572e+140; bool r21173 = r21111 <= r21172; double r21174 = +inf.0; bool r21175 = r21111 <= r21174; double r21176 = r21175 ? r21171 : r21171; double r21177 = r21173 ? r21152 : r21176; double r21178 = r21160 ? r21171 : r21177; double r21179 = r21156 ? r21158 : r21178; double r21180 = r21154 ? r21129 : r21179; double r21181 = r21131 ? r21152 : r21180; double r21182 = r21113 ? r21129 : r21181; return r21182; } 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 r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190, r21191, r21192, r21193, r21194, r21195, r21196, r21197, r21198, r21199, r21200; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r21183, "2", 10, MPFR_RNDN); mpfr_init(r21184); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init_set_str(r21188, "1", 10, MPFR_RNDN); mpfr_init(r21189); mpfr_init(r21190); mpfr_init(r21191); mpfr_init(r21192); mpfr_init(r21193); mpfr_init(r21194); mpfr_init(r21195); mpfr_init(r21196); mpfr_init(r21197); mpfr_init(r21198); mpfr_init(r21199); mpfr_init(r21200); } double f_im(double x, double l, double t) { ; mpfr_sqrt(r21184, r21183, MPFR_RNDN); mpfr_set_d(r21185, t, MPFR_RNDN); mpfr_mul(r21186, r21184, r21185, MPFR_RNDN); mpfr_set_d(r21187, x, MPFR_RNDN); ; mpfr_add(r21189, r21187, r21188, MPFR_RNDN); mpfr_sub(r21190, r21187, r21188, MPFR_RNDN); mpfr_div(r21191, r21189, r21190, MPFR_RNDN); mpfr_set_d(r21192, l, MPFR_RNDN); mpfr_mul(r21193, r21192, r21192, MPFR_RNDN); mpfr_mul(r21194, r21185, r21185, MPFR_RNDN); mpfr_mul(r21195, r21183, r21194, MPFR_RNDN); mpfr_add(r21196, r21193, r21195, MPFR_RNDN); mpfr_mul(r21197, r21191, r21196, MPFR_RNDN); mpfr_sub(r21198, r21197, r21193, MPFR_RNDN); mpfr_sqrt(r21199, r21198, MPFR_RNDN); mpfr_div(r21200, r21186, r21199, MPFR_RNDN); return mpfr_get_d(r21200, MPFR_RNDN); } static mpfr_t r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216, r21217, r21218, r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229, r21230, r21231, r21232, r21233, r21234, r21235, r21236, r21237, r21238, r21239, r21240, r21241, r21242, r21243, r21244, r21245, r21246, r21247, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r21201); mpfr_init_set_str(r21202, "-1.42998038934624e+69", 10, MPFR_RNDN); mpfr_init(r21203); mpfr_init_set_str(r21204, "2", 10, MPFR_RNDN); mpfr_init(r21205); mpfr_init(r21206); mpfr_init(r21207); mpfr_init(r21208); mpfr_init(r21209); mpfr_init(r21210); mpfr_init(r21211); mpfr_init(r21212); mpfr_init(r21213); mpfr_init(r21214); mpfr_init(r21215); mpfr_init(r21216); mpfr_init(r21217); mpfr_init(r21218); mpfr_init(r21219); mpfr_init_set_str(r21220, "-2.5401539047568627e-163", 10, MPFR_RNDN); mpfr_init(r21221); mpfr_init(r21222); mpfr_init(r21223); mpfr_init(r21224); mpfr_init(r21225); mpfr_init(r21226); mpfr_init(r21227); mpfr_init(r21228); mpfr_init(r21229); mpfr_init(r21230); mpfr_init(r21231); mpfr_init(r21232); mpfr_init(r21233); mpfr_init(r21234); mpfr_init(r21235); mpfr_init_set_str(r21236, "4", 10, MPFR_RNDN); mpfr_init(r21237); mpfr_init(r21238); mpfr_init(r21239); mpfr_init(r21240); mpfr_init(r21241); mpfr_init(r21242); mpfr_init_set_str(r21243, "-5.553561555528792e-295", 10, MPFR_RNDN); mpfr_init(r21244); mpfr_init_set_str(r21245, "1.2691519409648163e-243", 10, MPFR_RNDN); mpfr_init(r21246); mpfr_init(r21247); mpfr_init(r21248); mpfr_init_set_str(r21249, "7.679813441312647e-161", 10, MPFR_RNDN); mpfr_init(r21250); mpfr_init(r21251); mpfr_init(r21252); mpfr_init(r21253); mpfr_init(r21254); mpfr_init(r21255); mpfr_init_set_str(r21256, "1", 10, MPFR_RNDN); mpfr_init(r21257); mpfr_init(r21258); mpfr_init(r21259); mpfr_init(r21260); mpfr_init(r21261); mpfr_init_set_str(r21262, "4.634093238500572e+140", 10, MPFR_RNDN); mpfr_init(r21263); mpfr_init_set_str(r21264, "+inf.0", 10, MPFR_RNDN); mpfr_init(r21265); mpfr_init(r21266); mpfr_init(r21267); mpfr_init(r21268); mpfr_init(r21269); mpfr_init(r21270); mpfr_init(r21271); mpfr_init(r21272); } double f_fm(double x, double l, double t) { mpfr_set_d(r21201, t, MPFR_RNDN); ; mpfr_set_si(r21203, mpfr_cmp(r21201, r21202) <= 0, MPFR_RNDN); ; mpfr_sqrt(r21205, r21204, MPFR_RNDN); mpfr_mul(r21206, r21201, r21205, MPFR_RNDN); mpfr_set_d(r21207, x, MPFR_RNDN); mpfr_div(r21208, r21201, r21207, MPFR_RNDN); mpfr_div(r21209, r21208, r21207, MPFR_RNDN); mpfr_div(r21210, r21209, r21205, MPFR_RNDN); mpfr_sub(r21211, r21210, r21206, MPFR_RNDN); mpfr_add(r21212, r21201, r21201, MPFR_RNDN); mpfr_div(r21213, r21212, r21205, MPFR_RNDN); mpfr_div(r21214, r21213, r21207, MPFR_RNDN); mpfr_mul(r21215, r21207, r21207, MPFR_RNDN); mpfr_div(r21216, r21213, r21215, MPFR_RNDN); mpfr_add(r21217, r21214, r21216, MPFR_RNDN); mpfr_sub(r21218, r21211, r21217, MPFR_RNDN); mpfr_div(r21219, r21206, r21218, MPFR_RNDN); ; mpfr_set_si(r21221, mpfr_cmp(r21201, r21220) <= 0, MPFR_RNDN); mpfr_cbrt(r21222, r21205, MPFR_RNDN); mpfr_mul(r21223, r21222, r21222, MPFR_RNDN); mpfr_sqrt(r21224, r21223, MPFR_RNDN); mpfr_sqrt(r21225, r21222, MPFR_RNDN); mpfr_sqrt(r21226, r21205, MPFR_RNDN); mpfr_mul(r21227, r21226, r21201, MPFR_RNDN); mpfr_mul(r21228, r21225, r21227, MPFR_RNDN); mpfr_mul(r21229, r21224, r21228, MPFR_RNDN); mpfr_pow(r21230, r21201, r21204, MPFR_RNDN); mpfr_mul(r21231, r21204, r21230, MPFR_RNDN); mpfr_set_d(r21232, l, MPFR_RNDN); mpfr_div(r21233, r21207, r21232, MPFR_RNDN); mpfr_div(r21234, r21232, r21233, MPFR_RNDN); mpfr_mul(r21235, r21204, r21234, MPFR_RNDN); ; mpfr_div(r21237, r21230, r21207, MPFR_RNDN); mpfr_mul(r21238, r21236, r21237, MPFR_RNDN); mpfr_add(r21239, r21235, r21238, MPFR_RNDN); mpfr_add(r21240, r21231, r21239, MPFR_RNDN); mpfr_sqrt(r21241, r21240, MPFR_RNDN); mpfr_div(r21242, r21229, r21241, MPFR_RNDN); ; mpfr_set_si(r21244, mpfr_cmp(r21201, r21243) <= 0, MPFR_RNDN); ; mpfr_set_si(r21246, mpfr_cmp(r21201, r21245) <= 0, MPFR_RNDN); mpfr_mul(r21247, r21226, r21227, MPFR_RNDN); mpfr_div(r21248, r21247, r21241, MPFR_RNDN); ; mpfr_set_si(r21250, mpfr_cmp(r21201, r21249) <= 0, MPFR_RNDN); mpfr_div(r21251, r21204, r21207, MPFR_RNDN); mpfr_div(r21252, r21251, r21205, MPFR_RNDN); mpfr_add(r21253, r21205, r21252, MPFR_RNDN); mpfr_mul(r21254, r21201, r21253, MPFR_RNDN); mpfr_div(r21255, r21204, r21205, MPFR_RNDN); ; mpfr_div(r21257, r21256, r21205, MPFR_RNDN); mpfr_sub(r21258, r21255, r21257, MPFR_RNDN); mpfr_mul(r21259, r21209, r21258, MPFR_RNDN); mpfr_add(r21260, r21254, r21259, MPFR_RNDN); mpfr_div(r21261, r21206, r21260, MPFR_RNDN); ; mpfr_set_si(r21263, mpfr_cmp(r21201, r21262) <= 0, MPFR_RNDN); ; mpfr_set_si(r21265, mpfr_cmp(r21201, r21264) <= 0, MPFR_RNDN); if (mpfr_get_si(r21265, MPFR_RNDN)) { mpfr_set(r21266, r21261, MPFR_RNDN); } else { mpfr_set(r21266, r21261, MPFR_RNDN); }; if (mpfr_get_si(r21263, MPFR_RNDN)) { mpfr_set(r21267, r21242, MPFR_RNDN); } else { mpfr_set(r21267, r21266, MPFR_RNDN); }; if (mpfr_get_si(r21250, MPFR_RNDN)) { mpfr_set(r21268, r21261, MPFR_RNDN); } else { mpfr_set(r21268, r21267, MPFR_RNDN); }; if (mpfr_get_si(r21246, MPFR_RNDN)) { mpfr_set(r21269, r21248, MPFR_RNDN); } else { mpfr_set(r21269, r21268, MPFR_RNDN); }; if (mpfr_get_si(r21244, MPFR_RNDN)) { mpfr_set(r21270, r21219, MPFR_RNDN); } else { mpfr_set(r21270, r21269, MPFR_RNDN); }; if (mpfr_get_si(r21221, MPFR_RNDN)) { mpfr_set(r21271, r21242, MPFR_RNDN); } else { mpfr_set(r21271, r21270, MPFR_RNDN); }; if (mpfr_get_si(r21203, MPFR_RNDN)) { mpfr_set(r21272, r21219, MPFR_RNDN); } else { mpfr_set(r21272, r21271, MPFR_RNDN); }; return mpfr_get_d(r21272, MPFR_RNDN); } static mpfr_t 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, r21325, r21326, r21327, r21328, r21329, r21330, r21331, r21332, r21333, r21334, r21335, r21336, r21337, r21338, r21339, r21340, r21341, r21342, r21343, r21344; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r21273); mpfr_init_set_str(r21274, "-1.42998038934624e+69", 10, MPFR_RNDN); mpfr_init(r21275); mpfr_init_set_str(r21276, "2", 10, MPFR_RNDN); mpfr_init(r21277); mpfr_init(r21278); mpfr_init(r21279); mpfr_init(r21280); mpfr_init(r21281); mpfr_init(r21282); mpfr_init(r21283); mpfr_init(r21284); mpfr_init(r21285); mpfr_init(r21286); mpfr_init(r21287); mpfr_init(r21288); mpfr_init(r21289); mpfr_init(r21290); mpfr_init(r21291); mpfr_init_set_str(r21292, "-2.5401539047568627e-163", 10, MPFR_RNDN); 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(r21303); mpfr_init(r21304); mpfr_init(r21305); mpfr_init(r21306); mpfr_init(r21307); mpfr_init_set_str(r21308, "4", 10, MPFR_RNDN); mpfr_init(r21309); mpfr_init(r21310); mpfr_init(r21311); mpfr_init(r21312); mpfr_init(r21313); mpfr_init(r21314); mpfr_init_set_str(r21315, "-5.553561555528792e-295", 10, MPFR_RNDN); mpfr_init(r21316); mpfr_init_set_str(r21317, "1.2691519409648163e-243", 10, MPFR_RNDN); mpfr_init(r21318); mpfr_init(r21319); mpfr_init(r21320); mpfr_init_set_str(r21321, "7.679813441312647e-161", 10, MPFR_RNDN); mpfr_init(r21322); mpfr_init(r21323); mpfr_init(r21324); mpfr_init(r21325); mpfr_init(r21326); mpfr_init(r21327); mpfr_init_set_str(r21328, "1", 10, MPFR_RNDN); mpfr_init(r21329); mpfr_init(r21330); mpfr_init(r21331); mpfr_init(r21332); mpfr_init(r21333); mpfr_init_set_str(r21334, "4.634093238500572e+140", 10, MPFR_RNDN); mpfr_init(r21335); mpfr_init_set_str(r21336, "+inf.0", 10, MPFR_RNDN); mpfr_init(r21337); mpfr_init(r21338); mpfr_init(r21339); mpfr_init(r21340); mpfr_init(r21341); mpfr_init(r21342); mpfr_init(r21343); mpfr_init(r21344); } double f_dm(double x, double l, double t) { mpfr_set_d(r21273, t, MPFR_RNDN); ; mpfr_set_si(r21275, mpfr_cmp(r21273, r21274) <= 0, MPFR_RNDN); ; mpfr_sqrt(r21277, r21276, MPFR_RNDN); mpfr_mul(r21278, r21273, r21277, MPFR_RNDN); mpfr_set_d(r21279, x, MPFR_RNDN); mpfr_div(r21280, r21273, r21279, MPFR_RNDN); mpfr_div(r21281, r21280, r21279, MPFR_RNDN); mpfr_div(r21282, r21281, r21277, MPFR_RNDN); mpfr_sub(r21283, r21282, r21278, MPFR_RNDN); mpfr_add(r21284, r21273, r21273, MPFR_RNDN); mpfr_div(r21285, r21284, r21277, MPFR_RNDN); mpfr_div(r21286, r21285, r21279, MPFR_RNDN); mpfr_mul(r21287, r21279, r21279, MPFR_RNDN); mpfr_div(r21288, r21285, r21287, MPFR_RNDN); mpfr_add(r21289, r21286, r21288, MPFR_RNDN); mpfr_sub(r21290, r21283, r21289, MPFR_RNDN); mpfr_div(r21291, r21278, r21290, MPFR_RNDN); ; mpfr_set_si(r21293, mpfr_cmp(r21273, r21292) <= 0, MPFR_RNDN); mpfr_cbrt(r21294, r21277, MPFR_RNDN); mpfr_mul(r21295, r21294, r21294, MPFR_RNDN); mpfr_sqrt(r21296, r21295, MPFR_RNDN); mpfr_sqrt(r21297, r21294, MPFR_RNDN); mpfr_sqrt(r21298, r21277, MPFR_RNDN); mpfr_mul(r21299, r21298, r21273, MPFR_RNDN); mpfr_mul(r21300, r21297, r21299, MPFR_RNDN); mpfr_mul(r21301, r21296, r21300, MPFR_RNDN); mpfr_pow(r21302, r21273, r21276, MPFR_RNDN); mpfr_mul(r21303, r21276, r21302, MPFR_RNDN); mpfr_set_d(r21304, l, MPFR_RNDN); mpfr_div(r21305, r21279, r21304, MPFR_RNDN); mpfr_div(r21306, r21304, r21305, MPFR_RNDN); mpfr_mul(r21307, r21276, r21306, MPFR_RNDN); ; mpfr_div(r21309, r21302, r21279, MPFR_RNDN); mpfr_mul(r21310, r21308, r21309, MPFR_RNDN); mpfr_add(r21311, r21307, r21310, MPFR_RNDN); mpfr_add(r21312, r21303, r21311, MPFR_RNDN); mpfr_sqrt(r21313, r21312, MPFR_RNDN); mpfr_div(r21314, r21301, r21313, MPFR_RNDN); ; mpfr_set_si(r21316, mpfr_cmp(r21273, r21315) <= 0, MPFR_RNDN); ; mpfr_set_si(r21318, mpfr_cmp(r21273, r21317) <= 0, MPFR_RNDN); mpfr_mul(r21319, r21298, r21299, MPFR_RNDN); mpfr_div(r21320, r21319, r21313, MPFR_RNDN); ; mpfr_set_si(r21322, mpfr_cmp(r21273, r21321) <= 0, MPFR_RNDN); mpfr_div(r21323, r21276, r21279, MPFR_RNDN); mpfr_div(r21324, r21323, r21277, MPFR_RNDN); mpfr_add(r21325, r21277, r21324, MPFR_RNDN); mpfr_mul(r21326, r21273, r21325, MPFR_RNDN); mpfr_div(r21327, r21276, r21277, MPFR_RNDN); ; mpfr_div(r21329, r21328, r21277, MPFR_RNDN); mpfr_sub(r21330, r21327, r21329, MPFR_RNDN); mpfr_mul(r21331, r21281, r21330, MPFR_RNDN); mpfr_add(r21332, r21326, r21331, MPFR_RNDN); mpfr_div(r21333, r21278, r21332, MPFR_RNDN); ; mpfr_set_si(r21335, mpfr_cmp(r21273, r21334) <= 0, MPFR_RNDN); ; mpfr_set_si(r21337, mpfr_cmp(r21273, r21336) <= 0, MPFR_RNDN); if (mpfr_get_si(r21337, MPFR_RNDN)) { mpfr_set(r21338, r21333, MPFR_RNDN); } else { mpfr_set(r21338, r21333, MPFR_RNDN); }; if (mpfr_get_si(r21335, MPFR_RNDN)) { mpfr_set(r21339, r21314, MPFR_RNDN); } else { mpfr_set(r21339, r21338, MPFR_RNDN); }; if (mpfr_get_si(r21322, MPFR_RNDN)) { mpfr_set(r21340, r21333, MPFR_RNDN); } else { mpfr_set(r21340, r21339, MPFR_RNDN); }; if (mpfr_get_si(r21318, MPFR_RNDN)) { mpfr_set(r21341, r21320, MPFR_RNDN); } else { mpfr_set(r21341, r21340, MPFR_RNDN); }; if (mpfr_get_si(r21316, MPFR_RNDN)) { mpfr_set(r21342, r21291, MPFR_RNDN); } else { mpfr_set(r21342, r21341, MPFR_RNDN); }; if (mpfr_get_si(r21293, MPFR_RNDN)) { mpfr_set(r21343, r21314, MPFR_RNDN); } else { mpfr_set(r21343, r21342, MPFR_RNDN); }; if (mpfr_get_si(r21275, MPFR_RNDN)) { mpfr_set(r21344, r21291, MPFR_RNDN); } else { mpfr_set(r21344, r21343, MPFR_RNDN); }; return mpfr_get_d(r21344, MPFR_RNDN); }