#include #include #include #include #include char *name = "jeff quadratic root 2"; double f_if(float a, float b, float c) { float r8099 = b; float r8100 = 0; bool r8101 = r8099 >= r8100; float r8102 = 2; float r8103 = c; float r8104 = r8102 * r8103; float r8105 = -r8099; float r8106 = r8099 * r8099; float r8107 = 4; float r8108 = a; float r8109 = r8107 * r8108; float r8110 = r8109 * r8103; float r8111 = r8106 - r8110; float r8112 = sqrt(r8111); float r8113 = r8105 - r8112; float r8114 = r8104 / r8113; float r8115 = r8105 + r8112; float r8116 = r8102 * r8108; float r8117 = r8115 / r8116; float r8118 = r8101 ? r8114 : r8117; return r8118; } double f_id(double a, double b, double c) { double r8119 = b; double r8120 = 0; bool r8121 = r8119 >= r8120; double r8122 = 2; double r8123 = c; double r8124 = r8122 * r8123; double r8125 = -r8119; double r8126 = r8119 * r8119; double r8127 = 4; double r8128 = a; double r8129 = r8127 * r8128; double r8130 = r8129 * r8123; double r8131 = r8126 - r8130; double r8132 = sqrt(r8131); double r8133 = r8125 - r8132; double r8134 = r8124 / r8133; double r8135 = r8125 + r8132; double r8136 = r8122 * r8128; double r8137 = r8135 / r8136; double r8138 = r8121 ? r8134 : r8137; return r8138; } double f_of(float a, float b, float c) { float r8139 = b; float r8140 = 0; bool r8141 = r8139 >= r8140; float r8142 = 2; float r8143 = c; float r8144 = r8142 * r8143; float r8145 = -r8139; float r8146 = r8139 * r8139; float r8147 = a; float r8148 = 4; float r8149 = r8147 * r8148; float r8150 = r8143 * r8149; float r8151 = r8146 - r8150; float r8152 = sqrt(r8151); float r8153 = r8145 - r8152; float r8154 = r8144 / r8153; float r8155 = r8145 + r8152; float r8156 = r8142 * r8147; float r8157 = r8155 / r8156; float r8158 = r8141 ? r8154 : r8157; float r8159 = -1.0439296684993796e+241; bool r8160 = r8158 <= r8159; float r8161 = r8143 * r8147; float r8162 = r8161 * r8148; float r8163 = r8146 - r8162; float r8164 = sqrt(r8163); float r8165 = r8145 - r8164; float r8166 = r8144 / r8165; float r8167 = r8143 / r8139; float r8168 = r8139 / r8147; float r8169 = r8167 - r8168; float r8170 = r8141 ? r8166 : r8169; float r8171 = -1.9180879887997911e-295; bool r8172 = r8158 <= r8171; float r8173 = r8143 * r8148; float r8174 = r8173 * r8147; float r8175 = r8146 - r8174; float r8176 = sqrt(r8175); float r8177 = r8145 - r8176; float r8178 = cbrt(r8177); float r8179 = r8178 * r8178; float r8180 = r8179 * r8178; float r8181 = r8144 / r8180; float r8182 = r8176 - r8139; float r8183 = r8182 / r8156; float r8184 = r8141 ? r8181 : r8183; float r8185 = 2.3689218286716864e-293; bool r8186 = r8158 <= r8185; float r8187 = r8147 / r8139; float r8188 = cbrt(r8187); float r8189 = cbrt(r8188); float r8190 = r8189 * r8189; float r8191 = r8190 * r8189; float r8192 = r8191 * r8144; float r8193 = r8188 * r8188; float r8194 = r8192 * r8193; float r8195 = r8139 + r8139; float r8196 = r8194 - r8195; float r8197 = r8144 / r8196; float r8198 = r8164 - r8139; float r8199 = r8198 / r8156; float r8200 = r8141 ? r8197 : r8199; float r8201 = 1.1619586500950238e+284; bool r8202 = r8158 <= r8201; float r8203 = r8202 ? r8184 : r8170; float r8204 = r8186 ? r8200 : r8203; float r8205 = r8172 ? r8184 : r8204; float r8206 = r8160 ? r8170 : r8205; return r8206; } double f_od(double a, double b, double c) { double r8207 = b; double r8208 = 0; bool r8209 = r8207 >= r8208; double r8210 = 2; double r8211 = c; double r8212 = r8210 * r8211; double r8213 = -r8207; double r8214 = r8207 * r8207; double r8215 = a; double r8216 = 4; double r8217 = r8215 * r8216; double r8218 = r8211 * r8217; double r8219 = r8214 - r8218; double r8220 = sqrt(r8219); double r8221 = r8213 - r8220; double r8222 = r8212 / r8221; double r8223 = r8213 + r8220; double r8224 = r8210 * r8215; double r8225 = r8223 / r8224; double r8226 = r8209 ? r8222 : r8225; double r8227 = -1.0439296684993796e+241; bool r8228 = r8226 <= r8227; double r8229 = r8211 * r8215; double r8230 = r8229 * r8216; double r8231 = r8214 - r8230; double r8232 = sqrt(r8231); double r8233 = r8213 - r8232; double r8234 = r8212 / r8233; double r8235 = r8211 / r8207; double r8236 = r8207 / r8215; double r8237 = r8235 - r8236; double r8238 = r8209 ? r8234 : r8237; double r8239 = -1.9180879887997911e-295; bool r8240 = r8226 <= r8239; double r8241 = r8211 * r8216; double r8242 = r8241 * r8215; double r8243 = r8214 - r8242; double r8244 = sqrt(r8243); double r8245 = r8213 - r8244; double r8246 = cbrt(r8245); double r8247 = r8246 * r8246; double r8248 = r8247 * r8246; double r8249 = r8212 / r8248; double r8250 = r8244 - r8207; double r8251 = r8250 / r8224; double r8252 = r8209 ? r8249 : r8251; double r8253 = 2.3689218286716864e-293; bool r8254 = r8226 <= r8253; double r8255 = r8215 / r8207; double r8256 = cbrt(r8255); double r8257 = cbrt(r8256); double r8258 = r8257 * r8257; double r8259 = r8258 * r8257; double r8260 = r8259 * r8212; double r8261 = r8256 * r8256; double r8262 = r8260 * r8261; double r8263 = r8207 + r8207; double r8264 = r8262 - r8263; double r8265 = r8212 / r8264; double r8266 = r8232 - r8207; double r8267 = r8266 / r8224; double r8268 = r8209 ? r8265 : r8267; double r8269 = 1.1619586500950238e+284; bool r8270 = r8226 <= r8269; double r8271 = r8270 ? r8252 : r8238; double r8272 = r8254 ? r8268 : r8271; double r8273 = r8240 ? r8252 : r8272; double r8274 = r8228 ? r8238 : r8273; return r8274; } 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 r8275, r8276, r8277, r8278, r8279, r8280, r8281, r8282, r8283, r8284, r8285, r8286, r8287, r8288, r8289, r8290, r8291, r8292, r8293, r8294; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r8275); mpfr_init_set_str(r8276, "0", 10, MPFR_RNDN); mpfr_init(r8277); mpfr_init_set_str(r8278, "2", 10, MPFR_RNDN); mpfr_init(r8279); mpfr_init(r8280); mpfr_init(r8281); mpfr_init(r8282); mpfr_init_set_str(r8283, "4", 10, MPFR_RNDN); mpfr_init(r8284); mpfr_init(r8285); mpfr_init(r8286); mpfr_init(r8287); mpfr_init(r8288); mpfr_init(r8289); mpfr_init(r8290); mpfr_init(r8291); mpfr_init(r8292); mpfr_init(r8293); mpfr_init(r8294); } double f_im(double a, double b, double c) { mpfr_set_d(r8275, b, MPFR_RNDN); ; mpfr_set_si(r8277, mpfr_cmp(r8275, r8276) >= 0, MPFR_RNDN); ; mpfr_set_d(r8279, c, MPFR_RNDN); mpfr_mul(r8280, r8278, r8279, MPFR_RNDN); mpfr_neg(r8281, r8275, MPFR_RNDN); mpfr_mul(r8282, r8275, r8275, MPFR_RNDN); ; mpfr_set_d(r8284, a, MPFR_RNDN); mpfr_mul(r8285, r8283, r8284, MPFR_RNDN); mpfr_mul(r8286, r8285, r8279, MPFR_RNDN); mpfr_sub(r8287, r8282, r8286, MPFR_RNDN); mpfr_sqrt(r8288, r8287, MPFR_RNDN); mpfr_sub(r8289, r8281, r8288, MPFR_RNDN); mpfr_div(r8290, r8280, r8289, MPFR_RNDN); mpfr_add(r8291, r8281, r8288, MPFR_RNDN); mpfr_mul(r8292, r8278, r8284, MPFR_RNDN); mpfr_div(r8293, r8291, r8292, MPFR_RNDN); if (mpfr_get_si(r8277, MPFR_RNDN)) { mpfr_set(r8294, r8290, MPFR_RNDN); } else { mpfr_set(r8294, r8293, MPFR_RNDN); }; return mpfr_get_d(r8294, MPFR_RNDN); } static mpfr_t r8295, r8296, r8297, r8298, r8299, r8300, r8301, r8302, r8303, r8304, r8305, r8306, r8307, r8308, r8309, r8310, r8311, r8312, r8313, r8314, r8315, r8316, r8317, r8318, r8319, r8320, r8321, r8322, r8323, r8324, r8325, r8326, r8327, r8328, r8329, r8330, r8331, r8332, r8333, r8334, r8335, r8336, r8337, r8338, r8339, r8340, r8341, r8342, r8343, r8344, r8345, r8346, r8347, r8348, r8349, r8350, r8351, r8352, r8353, r8354, r8355, r8356, r8357, r8358, r8359, r8360, r8361, r8362; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r8295); mpfr_init_set_str(r8296, "0", 10, MPFR_RNDN); mpfr_init(r8297); mpfr_init_set_str(r8298, "2", 10, MPFR_RNDN); mpfr_init(r8299); mpfr_init(r8300); mpfr_init(r8301); mpfr_init(r8302); mpfr_init(r8303); mpfr_init_set_str(r8304, "4", 10, MPFR_RNDN); mpfr_init(r8305); mpfr_init(r8306); mpfr_init(r8307); mpfr_init(r8308); mpfr_init(r8309); mpfr_init(r8310); mpfr_init(r8311); mpfr_init(r8312); mpfr_init(r8313); mpfr_init(r8314); mpfr_init_set_str(r8315, "-1.0439296684993796e+241", 10, MPFR_RNDN); mpfr_init(r8316); mpfr_init(r8317); mpfr_init(r8318); mpfr_init(r8319); mpfr_init(r8320); mpfr_init(r8321); mpfr_init(r8322); mpfr_init(r8323); mpfr_init(r8324); mpfr_init(r8325); mpfr_init(r8326); mpfr_init_set_str(r8327, "-1.9180879887997911e-295", 10, MPFR_RNDN); mpfr_init(r8328); mpfr_init(r8329); mpfr_init(r8330); mpfr_init(r8331); mpfr_init(r8332); mpfr_init(r8333); mpfr_init(r8334); mpfr_init(r8335); mpfr_init(r8336); mpfr_init(r8337); mpfr_init(r8338); mpfr_init(r8339); mpfr_init(r8340); mpfr_init_set_str(r8341, "2.3689218286716864e-293", 10, MPFR_RNDN); mpfr_init(r8342); mpfr_init(r8343); mpfr_init(r8344); mpfr_init(r8345); mpfr_init(r8346); mpfr_init(r8347); mpfr_init(r8348); mpfr_init(r8349); mpfr_init(r8350); mpfr_init(r8351); mpfr_init(r8352); mpfr_init(r8353); mpfr_init(r8354); mpfr_init(r8355); mpfr_init(r8356); mpfr_init_set_str(r8357, "1.1619586500950238e+284", 10, MPFR_RNDN); mpfr_init(r8358); mpfr_init(r8359); mpfr_init(r8360); mpfr_init(r8361); mpfr_init(r8362); } double f_fm(double a, double b, double c) { mpfr_set_d(r8295, b, MPFR_RNDN); ; mpfr_set_si(r8297, mpfr_cmp(r8295, r8296) >= 0, MPFR_RNDN); ; mpfr_set_d(r8299, c, MPFR_RNDN); mpfr_mul(r8300, r8298, r8299, MPFR_RNDN); mpfr_neg(r8301, r8295, MPFR_RNDN); mpfr_mul(r8302, r8295, r8295, MPFR_RNDN); mpfr_set_d(r8303, a, MPFR_RNDN); ; mpfr_mul(r8305, r8303, r8304, MPFR_RNDN); mpfr_mul(r8306, r8299, r8305, MPFR_RNDN); mpfr_sub(r8307, r8302, r8306, MPFR_RNDN); mpfr_sqrt(r8308, r8307, MPFR_RNDN); mpfr_sub(r8309, r8301, r8308, MPFR_RNDN); mpfr_div(r8310, r8300, r8309, MPFR_RNDN); mpfr_add(r8311, r8301, r8308, MPFR_RNDN); mpfr_mul(r8312, r8298, r8303, MPFR_RNDN); mpfr_div(r8313, r8311, r8312, MPFR_RNDN); if (mpfr_get_si(r8297, MPFR_RNDN)) { mpfr_set(r8314, r8310, MPFR_RNDN); } else { mpfr_set(r8314, r8313, MPFR_RNDN); }; ; mpfr_set_si(r8316, mpfr_cmp(r8314, r8315) <= 0, MPFR_RNDN); mpfr_mul(r8317, r8299, r8303, MPFR_RNDN); mpfr_mul(r8318, r8317, r8304, MPFR_RNDN); mpfr_sub(r8319, r8302, r8318, MPFR_RNDN); mpfr_sqrt(r8320, r8319, MPFR_RNDN); mpfr_sub(r8321, r8301, r8320, MPFR_RNDN); mpfr_div(r8322, r8300, r8321, MPFR_RNDN); mpfr_div(r8323, r8299, r8295, MPFR_RNDN); mpfr_div(r8324, r8295, r8303, MPFR_RNDN); mpfr_sub(r8325, r8323, r8324, MPFR_RNDN); if (mpfr_get_si(r8297, MPFR_RNDN)) { mpfr_set(r8326, r8322, MPFR_RNDN); } else { mpfr_set(r8326, r8325, MPFR_RNDN); }; ; mpfr_set_si(r8328, mpfr_cmp(r8314, r8327) <= 0, MPFR_RNDN); mpfr_mul(r8329, r8299, r8304, MPFR_RNDN); mpfr_mul(r8330, r8329, r8303, MPFR_RNDN); mpfr_sub(r8331, r8302, r8330, MPFR_RNDN); mpfr_sqrt(r8332, r8331, MPFR_RNDN); mpfr_sub(r8333, r8301, r8332, MPFR_RNDN); mpfr_cbrt(r8334, r8333, MPFR_RNDN); mpfr_mul(r8335, r8334, r8334, MPFR_RNDN); mpfr_mul(r8336, r8335, r8334, MPFR_RNDN); mpfr_div(r8337, r8300, r8336, MPFR_RNDN); mpfr_sub(r8338, r8332, r8295, MPFR_RNDN); mpfr_div(r8339, r8338, r8312, MPFR_RNDN); if (mpfr_get_si(r8297, MPFR_RNDN)) { mpfr_set(r8340, r8337, MPFR_RNDN); } else { mpfr_set(r8340, r8339, MPFR_RNDN); }; ; mpfr_set_si(r8342, mpfr_cmp(r8314, r8341) <= 0, MPFR_RNDN); mpfr_div(r8343, r8303, r8295, MPFR_RNDN); mpfr_cbrt(r8344, r8343, MPFR_RNDN); mpfr_cbrt(r8345, r8344, MPFR_RNDN); mpfr_mul(r8346, r8345, r8345, MPFR_RNDN); mpfr_mul(r8347, r8346, r8345, MPFR_RNDN); mpfr_mul(r8348, r8347, r8300, MPFR_RNDN); mpfr_mul(r8349, r8344, r8344, MPFR_RNDN); mpfr_mul(r8350, r8348, r8349, MPFR_RNDN); mpfr_add(r8351, r8295, r8295, MPFR_RNDN); mpfr_sub(r8352, r8350, r8351, MPFR_RNDN); mpfr_div(r8353, r8300, r8352, MPFR_RNDN); mpfr_sub(r8354, r8320, r8295, MPFR_RNDN); mpfr_div(r8355, r8354, r8312, MPFR_RNDN); if (mpfr_get_si(r8297, MPFR_RNDN)) { mpfr_set(r8356, r8353, MPFR_RNDN); } else { mpfr_set(r8356, r8355, MPFR_RNDN); }; ; mpfr_set_si(r8358, mpfr_cmp(r8314, r8357) <= 0, MPFR_RNDN); if (mpfr_get_si(r8358, MPFR_RNDN)) { mpfr_set(r8359, r8340, MPFR_RNDN); } else { mpfr_set(r8359, r8326, MPFR_RNDN); }; if (mpfr_get_si(r8342, MPFR_RNDN)) { mpfr_set(r8360, r8356, MPFR_RNDN); } else { mpfr_set(r8360, r8359, MPFR_RNDN); }; if (mpfr_get_si(r8328, MPFR_RNDN)) { mpfr_set(r8361, r8340, MPFR_RNDN); } else { mpfr_set(r8361, r8360, MPFR_RNDN); }; if (mpfr_get_si(r8316, MPFR_RNDN)) { mpfr_set(r8362, r8326, MPFR_RNDN); } else { mpfr_set(r8362, r8361, MPFR_RNDN); }; return mpfr_get_d(r8362, MPFR_RNDN); } static mpfr_t r8363, r8364, r8365, r8366, r8367, r8368, r8369, r8370, r8371, r8372, r8373, r8374, r8375, r8376, r8377, r8378, r8379, r8380, r8381, r8382, r8383, r8384, r8385, r8386, r8387, r8388, r8389, r8390, r8391, r8392, r8393, r8394, r8395, r8396, r8397, r8398, r8399, r8400, r8401, r8402, r8403, r8404, r8405, r8406, r8407, r8408, r8409, r8410, r8411, r8412, r8413, r8414, r8415, r8416, r8417, r8418, r8419, r8420, r8421, r8422, r8423, r8424, r8425, r8426, r8427, r8428, r8429, r8430; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r8363); mpfr_init_set_str(r8364, "0", 10, MPFR_RNDN); mpfr_init(r8365); mpfr_init_set_str(r8366, "2", 10, MPFR_RNDN); mpfr_init(r8367); mpfr_init(r8368); mpfr_init(r8369); mpfr_init(r8370); mpfr_init(r8371); mpfr_init_set_str(r8372, "4", 10, MPFR_RNDN); mpfr_init(r8373); mpfr_init(r8374); mpfr_init(r8375); mpfr_init(r8376); mpfr_init(r8377); mpfr_init(r8378); mpfr_init(r8379); mpfr_init(r8380); mpfr_init(r8381); mpfr_init(r8382); mpfr_init_set_str(r8383, "-1.0439296684993796e+241", 10, MPFR_RNDN); mpfr_init(r8384); mpfr_init(r8385); mpfr_init(r8386); mpfr_init(r8387); mpfr_init(r8388); mpfr_init(r8389); mpfr_init(r8390); mpfr_init(r8391); mpfr_init(r8392); mpfr_init(r8393); mpfr_init(r8394); mpfr_init_set_str(r8395, "-1.9180879887997911e-295", 10, MPFR_RNDN); mpfr_init(r8396); mpfr_init(r8397); mpfr_init(r8398); mpfr_init(r8399); mpfr_init(r8400); mpfr_init(r8401); mpfr_init(r8402); mpfr_init(r8403); mpfr_init(r8404); mpfr_init(r8405); mpfr_init(r8406); mpfr_init(r8407); mpfr_init(r8408); mpfr_init_set_str(r8409, "2.3689218286716864e-293", 10, MPFR_RNDN); mpfr_init(r8410); mpfr_init(r8411); mpfr_init(r8412); mpfr_init(r8413); mpfr_init(r8414); mpfr_init(r8415); mpfr_init(r8416); mpfr_init(r8417); mpfr_init(r8418); mpfr_init(r8419); mpfr_init(r8420); mpfr_init(r8421); mpfr_init(r8422); mpfr_init(r8423); mpfr_init(r8424); mpfr_init_set_str(r8425, "1.1619586500950238e+284", 10, MPFR_RNDN); mpfr_init(r8426); mpfr_init(r8427); mpfr_init(r8428); mpfr_init(r8429); mpfr_init(r8430); } double f_dm(double a, double b, double c) { mpfr_set_d(r8363, b, MPFR_RNDN); ; mpfr_set_si(r8365, mpfr_cmp(r8363, r8364) >= 0, MPFR_RNDN); ; mpfr_set_d(r8367, c, MPFR_RNDN); mpfr_mul(r8368, r8366, r8367, MPFR_RNDN); mpfr_neg(r8369, r8363, MPFR_RNDN); mpfr_mul(r8370, r8363, r8363, MPFR_RNDN); mpfr_set_d(r8371, a, MPFR_RNDN); ; mpfr_mul(r8373, r8371, r8372, MPFR_RNDN); mpfr_mul(r8374, r8367, r8373, MPFR_RNDN); mpfr_sub(r8375, r8370, r8374, MPFR_RNDN); mpfr_sqrt(r8376, r8375, MPFR_RNDN); mpfr_sub(r8377, r8369, r8376, MPFR_RNDN); mpfr_div(r8378, r8368, r8377, MPFR_RNDN); mpfr_add(r8379, r8369, r8376, MPFR_RNDN); mpfr_mul(r8380, r8366, r8371, MPFR_RNDN); mpfr_div(r8381, r8379, r8380, MPFR_RNDN); if (mpfr_get_si(r8365, MPFR_RNDN)) { mpfr_set(r8382, r8378, MPFR_RNDN); } else { mpfr_set(r8382, r8381, MPFR_RNDN); }; ; mpfr_set_si(r8384, mpfr_cmp(r8382, r8383) <= 0, MPFR_RNDN); mpfr_mul(r8385, r8367, r8371, MPFR_RNDN); mpfr_mul(r8386, r8385, r8372, MPFR_RNDN); mpfr_sub(r8387, r8370, r8386, MPFR_RNDN); mpfr_sqrt(r8388, r8387, MPFR_RNDN); mpfr_sub(r8389, r8369, r8388, MPFR_RNDN); mpfr_div(r8390, r8368, r8389, MPFR_RNDN); mpfr_div(r8391, r8367, r8363, MPFR_RNDN); mpfr_div(r8392, r8363, r8371, MPFR_RNDN); mpfr_sub(r8393, r8391, r8392, MPFR_RNDN); if (mpfr_get_si(r8365, MPFR_RNDN)) { mpfr_set(r8394, r8390, MPFR_RNDN); } else { mpfr_set(r8394, r8393, MPFR_RNDN); }; ; mpfr_set_si(r8396, mpfr_cmp(r8382, r8395) <= 0, MPFR_RNDN); mpfr_mul(r8397, r8367, r8372, MPFR_RNDN); mpfr_mul(r8398, r8397, r8371, MPFR_RNDN); mpfr_sub(r8399, r8370, r8398, MPFR_RNDN); mpfr_sqrt(r8400, r8399, MPFR_RNDN); mpfr_sub(r8401, r8369, r8400, MPFR_RNDN); mpfr_cbrt(r8402, r8401, MPFR_RNDN); mpfr_mul(r8403, r8402, r8402, MPFR_RNDN); mpfr_mul(r8404, r8403, r8402, MPFR_RNDN); mpfr_div(r8405, r8368, r8404, MPFR_RNDN); mpfr_sub(r8406, r8400, r8363, MPFR_RNDN); mpfr_div(r8407, r8406, r8380, MPFR_RNDN); if (mpfr_get_si(r8365, MPFR_RNDN)) { mpfr_set(r8408, r8405, MPFR_RNDN); } else { mpfr_set(r8408, r8407, MPFR_RNDN); }; ; mpfr_set_si(r8410, mpfr_cmp(r8382, r8409) <= 0, MPFR_RNDN); mpfr_div(r8411, r8371, r8363, MPFR_RNDN); mpfr_cbrt(r8412, r8411, MPFR_RNDN); mpfr_cbrt(r8413, r8412, MPFR_RNDN); mpfr_mul(r8414, r8413, r8413, MPFR_RNDN); mpfr_mul(r8415, r8414, r8413, MPFR_RNDN); mpfr_mul(r8416, r8415, r8368, MPFR_RNDN); mpfr_mul(r8417, r8412, r8412, MPFR_RNDN); mpfr_mul(r8418, r8416, r8417, MPFR_RNDN); mpfr_add(r8419, r8363, r8363, MPFR_RNDN); mpfr_sub(r8420, r8418, r8419, MPFR_RNDN); mpfr_div(r8421, r8368, r8420, MPFR_RNDN); mpfr_sub(r8422, r8388, r8363, MPFR_RNDN); mpfr_div(r8423, r8422, r8380, MPFR_RNDN); if (mpfr_get_si(r8365, MPFR_RNDN)) { mpfr_set(r8424, r8421, MPFR_RNDN); } else { mpfr_set(r8424, r8423, MPFR_RNDN); }; ; mpfr_set_si(r8426, mpfr_cmp(r8382, r8425) <= 0, MPFR_RNDN); if (mpfr_get_si(r8426, MPFR_RNDN)) { mpfr_set(r8427, r8408, MPFR_RNDN); } else { mpfr_set(r8427, r8394, MPFR_RNDN); }; if (mpfr_get_si(r8410, MPFR_RNDN)) { mpfr_set(r8428, r8424, MPFR_RNDN); } else { mpfr_set(r8428, r8427, MPFR_RNDN); }; if (mpfr_get_si(r8396, MPFR_RNDN)) { mpfr_set(r8429, r8408, MPFR_RNDN); } else { mpfr_set(r8429, r8428, MPFR_RNDN); }; if (mpfr_get_si(r8384, MPFR_RNDN)) { mpfr_set(r8430, r8394, MPFR_RNDN); } else { mpfr_set(r8430, r8429, MPFR_RNDN); }; return mpfr_get_d(r8430, MPFR_RNDN); }