#include #include #include #include #include char *name = "From Warwick Tucker's Validated Numerics"; double f_if() { float r8122 = 333.75f; float r8123 = 33096.0f; float r8124 = 6.0f; float r8125 = pow(r8123, r8124); float r8126 = r8122 * r8125; float r8127 = 77617.0f; float r8128 = r8127 * r8127; float r8129 = 11.0f; float r8130 = r8129 * r8128; float r8131 = r8123 * r8123; float r8132 = r8130 * r8131; float r8133 = -r8125; float r8134 = r8132 + r8133; float r8135 = -121.0f; float r8136 = 4.0f; float r8137 = pow(r8123, r8136); float r8138 = r8135 * r8137; float r8139 = r8134 + r8138; float r8140 = -2.0f; float r8141 = r8139 + r8140; float r8142 = r8128 * r8141; float r8143 = r8126 + r8142; float r8144 = 5.5f; float r8145 = 8.0f; float r8146 = pow(r8123, r8145); float r8147 = r8144 * r8146; float r8148 = r8143 + r8147; float r8149 = 2.0f; float r8150 = r8149 * r8123; float r8151 = r8127 / r8150; float r8152 = r8148 + r8151; return r8152; } double f_id() { double r8153 = 333.75; double r8154 = 33096.0; double r8155 = 6.0; double r8156 = pow(r8154, r8155); double r8157 = r8153 * r8156; double r8158 = 77617.0; double r8159 = r8158 * r8158; double r8160 = 11.0; double r8161 = r8160 * r8159; double r8162 = r8154 * r8154; double r8163 = r8161 * r8162; double r8164 = -r8156; double r8165 = r8163 + r8164; double r8166 = -121.0; double r8167 = 4.0; double r8168 = pow(r8154, r8167); double r8169 = r8166 * r8168; double r8170 = r8165 + r8169; double r8171 = -2.0; double r8172 = r8170 + r8171; double r8173 = r8159 * r8172; double r8174 = r8157 + r8173; double r8175 = 5.5; double r8176 = 8.0; double r8177 = pow(r8154, r8176); double r8178 = r8175 * r8177; double r8179 = r8174 + r8178; double r8180 = 2.0; double r8181 = r8180 * r8154; double r8182 = r8158 / r8181; double r8183 = r8179 + r8182; return r8183; } double f_of() { float r8184 = 77617.0f; float r8185 = r8184 * r8184; float r8186 = 33096.0f; float r8187 = r8184 * r8186; float r8188 = r8187 * r8184; float r8189 = 11.0f; float r8190 = r8186 * r8189; float r8191 = r8188 * r8190; float r8192 = r8185 * r8191; float r8193 = r8192 * r8192; float r8194 = 5.5f; float r8195 = 8.0f; float r8196 = pow(r8186, r8195); float r8197 = r8194 * r8196; float r8198 = 2.0f; float r8199 = r8184 / r8198; float r8200 = r8199 / r8186; float r8201 = r8197 + r8200; float r8202 = 333.75f; float r8203 = 6.0f; float r8204 = pow(r8186, r8203); float r8205 = r8202 * r8204; float r8206 = r8201 + r8205; float r8207 = 3.0f; float r8208 = pow(r8206, r8207); float r8209 = -r8184; float r8210 = r8209 * r8184; float r8211 = -2.0f; float r8212 = r8204 - r8211; float r8213 = -121.0f; float r8214 = 4.0f; float r8215 = pow(r8186, r8214); float r8216 = r8213 * r8215; float r8217 = r8212 - r8216; float r8218 = r8210 * r8217; float r8219 = pow(r8218, r8207); float r8220 = r8208 + r8219; float r8221 = r8206 * r8206; float r8222 = r8218 * r8206; float r8223 = r8221 - r8222; float r8224 = r8185 * r8217; float r8225 = r8224 * r8224; float r8226 = r8223 + r8225; float r8227 = r8220 / r8226; float r8228 = r8227 * r8227; float r8229 = r8193 - r8228; float r8230 = r8216 + r8212; float r8231 = pow(r8230, r8207); float r8232 = r8205 * r8205; float r8233 = r8201 * r8201; float r8234 = r8232 - r8233; float r8235 = pow(r8234, r8207); float r8236 = r8231 * r8235; float r8237 = r8212 * r8212; float r8238 = r8216 * r8216; float r8239 = r8237 - r8238; float r8240 = r8210 * r8239; float r8241 = pow(r8240, r8207); float r8242 = r8205 - r8201; float r8243 = pow(r8242, r8207); float r8244 = r8241 * r8243; float r8245 = r8236 + r8244; float r8246 = r8243 * r8231; float r8247 = r8245 / r8246; float r8248 = r8247 / r8226; float r8249 = r8192 - r8248; float r8250 = r8229 / r8249; return r8250; } double f_od() { double r8251 = 77617.0; double r8252 = r8251 * r8251; double r8253 = 33096.0; double r8254 = r8251 * r8253; double r8255 = r8254 * r8251; double r8256 = 11.0; double r8257 = r8253 * r8256; double r8258 = r8255 * r8257; double r8259 = r8252 * r8258; double r8260 = r8259 * r8259; double r8261 = 5.5; double r8262 = 8.0; double r8263 = pow(r8253, r8262); double r8264 = r8261 * r8263; double r8265 = 2.0; double r8266 = r8251 / r8265; double r8267 = r8266 / r8253; double r8268 = r8264 + r8267; double r8269 = 333.75; double r8270 = 6.0; double r8271 = pow(r8253, r8270); double r8272 = r8269 * r8271; double r8273 = r8268 + r8272; double r8274 = 3.0; double r8275 = pow(r8273, r8274); double r8276 = -r8251; double r8277 = r8276 * r8251; double r8278 = -2.0; double r8279 = r8271 - r8278; double r8280 = -121.0; double r8281 = 4.0; double r8282 = pow(r8253, r8281); double r8283 = r8280 * r8282; double r8284 = r8279 - r8283; double r8285 = r8277 * r8284; double r8286 = pow(r8285, r8274); double r8287 = r8275 + r8286; double r8288 = r8273 * r8273; double r8289 = r8285 * r8273; double r8290 = r8288 - r8289; double r8291 = r8252 * r8284; double r8292 = r8291 * r8291; double r8293 = r8290 + r8292; double r8294 = r8287 / r8293; double r8295 = r8294 * r8294; double r8296 = r8260 - r8295; double r8297 = r8283 + r8279; double r8298 = pow(r8297, r8274); double r8299 = r8272 * r8272; double r8300 = r8268 * r8268; double r8301 = r8299 - r8300; double r8302 = pow(r8301, r8274); double r8303 = r8298 * r8302; double r8304 = r8279 * r8279; double r8305 = r8283 * r8283; double r8306 = r8304 - r8305; double r8307 = r8277 * r8306; double r8308 = pow(r8307, r8274); double r8309 = r8272 - r8268; double r8310 = pow(r8309, r8274); double r8311 = r8308 * r8310; double r8312 = r8303 + r8311; double r8313 = r8310 * r8298; double r8314 = r8312 / r8313; double r8315 = r8314 / r8293; double r8316 = r8259 - r8315; double r8317 = r8296 / r8316; return r8317; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r8318, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r8319, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r8320, "6", 10, MPFR_RNDN); mpfr_init(r8321); mpfr_init(r8322); mpfr_init_set_str(r8323, "77617", 10, MPFR_RNDN); mpfr_init(r8324); mpfr_init_set_str(r8325, "11", 10, MPFR_RNDN); mpfr_init(r8326); mpfr_init(r8327); mpfr_init(r8328); mpfr_init(r8329); mpfr_init(r8330); mpfr_init_set_str(r8331, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r8332, "4", 10, MPFR_RNDN); mpfr_init(r8333); mpfr_init(r8334); mpfr_init(r8335); mpfr_init_set_str(r8336, "-2", 10, MPFR_RNDN); mpfr_init(r8337); mpfr_init(r8338); mpfr_init(r8339); mpfr_init_set_str(r8340, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r8341, "8", 10, MPFR_RNDN); mpfr_init(r8342); mpfr_init(r8343); mpfr_init(r8344); mpfr_init_set_str(r8345, "2", 10, MPFR_RNDN); mpfr_init(r8346); mpfr_init(r8347); mpfr_init(r8348); } double f_im() { ; ; ; mpfr_pow(r8321, r8319, r8320, MPFR_RNDN); mpfr_mul(r8322, r8318, r8321, MPFR_RNDN); ; mpfr_mul(r8324, r8323, r8323, MPFR_RNDN); ; mpfr_mul(r8326, r8325, r8324, MPFR_RNDN); mpfr_mul(r8327, r8319, r8319, MPFR_RNDN); mpfr_mul(r8328, r8326, r8327, MPFR_RNDN); mpfr_neg(r8329, r8321, MPFR_RNDN); mpfr_add(r8330, r8328, r8329, MPFR_RNDN); ; ; mpfr_pow(r8333, r8319, r8332, MPFR_RNDN); mpfr_mul(r8334, r8331, r8333, MPFR_RNDN); mpfr_add(r8335, r8330, r8334, MPFR_RNDN); ; mpfr_add(r8337, r8335, r8336, MPFR_RNDN); mpfr_mul(r8338, r8324, r8337, MPFR_RNDN); mpfr_add(r8339, r8322, r8338, MPFR_RNDN); ; ; mpfr_pow(r8342, r8319, r8341, MPFR_RNDN); mpfr_mul(r8343, r8340, r8342, MPFR_RNDN); mpfr_add(r8344, r8339, r8343, MPFR_RNDN); ; mpfr_mul(r8346, r8345, r8319, MPFR_RNDN); mpfr_div(r8347, r8323, r8346, MPFR_RNDN); mpfr_add(r8348, r8344, r8347, MPFR_RNDN); return mpfr_get_d(r8348, MPFR_RNDN); } static mpfr_t r8349, r8350, r8351, r8352, r8353, r8354, r8355, r8356, r8357, r8358, r8359, r8360, r8361, r8362, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8349, "77617", 10, MPFR_RNDN); mpfr_init(r8350); mpfr_init_set_str(r8351, "33096", 10, MPFR_RNDN); mpfr_init(r8352); mpfr_init(r8353); mpfr_init_set_str(r8354, "11", 10, MPFR_RNDN); mpfr_init(r8355); mpfr_init(r8356); mpfr_init(r8357); mpfr_init(r8358); mpfr_init_set_str(r8359, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r8360, "8", 10, MPFR_RNDN); mpfr_init(r8361); mpfr_init(r8362); mpfr_init_set_str(r8363, "2", 10, MPFR_RNDN); mpfr_init(r8364); mpfr_init(r8365); mpfr_init(r8366); mpfr_init_set_str(r8367, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r8368, "6", 10, MPFR_RNDN); mpfr_init(r8369); mpfr_init(r8370); mpfr_init(r8371); mpfr_init_set_str(r8372, "3", 10, MPFR_RNDN); mpfr_init(r8373); mpfr_init(r8374); mpfr_init(r8375); mpfr_init_set_str(r8376, "-2", 10, MPFR_RNDN); mpfr_init(r8377); mpfr_init_set_str(r8378, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r8379, "4", 10, MPFR_RNDN); mpfr_init(r8380); mpfr_init(r8381); mpfr_init(r8382); mpfr_init(r8383); 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(r8395); 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(r8409); mpfr_init(r8410); mpfr_init(r8411); mpfr_init(r8412); mpfr_init(r8413); mpfr_init(r8414); mpfr_init(r8415); } double f_fm() { ; mpfr_mul(r8350, r8349, r8349, MPFR_RNDN); ; mpfr_mul(r8352, r8349, r8351, MPFR_RNDN); mpfr_mul(r8353, r8352, r8349, MPFR_RNDN); ; mpfr_mul(r8355, r8351, r8354, MPFR_RNDN); mpfr_mul(r8356, r8353, r8355, MPFR_RNDN); mpfr_mul(r8357, r8350, r8356, MPFR_RNDN); mpfr_mul(r8358, r8357, r8357, MPFR_RNDN); ; ; mpfr_pow(r8361, r8351, r8360, MPFR_RNDN); mpfr_mul(r8362, r8359, r8361, MPFR_RNDN); ; mpfr_div(r8364, r8349, r8363, MPFR_RNDN); mpfr_div(r8365, r8364, r8351, MPFR_RNDN); mpfr_add(r8366, r8362, r8365, MPFR_RNDN); ; ; mpfr_pow(r8369, r8351, r8368, MPFR_RNDN); mpfr_mul(r8370, r8367, r8369, MPFR_RNDN); mpfr_add(r8371, r8366, r8370, MPFR_RNDN); ; mpfr_pow(r8373, r8371, r8372, MPFR_RNDN); mpfr_neg(r8374, r8349, MPFR_RNDN); mpfr_mul(r8375, r8374, r8349, MPFR_RNDN); ; mpfr_sub(r8377, r8369, r8376, MPFR_RNDN); ; ; mpfr_pow(r8380, r8351, r8379, MPFR_RNDN); mpfr_mul(r8381, r8378, r8380, MPFR_RNDN); mpfr_sub(r8382, r8377, r8381, MPFR_RNDN); mpfr_mul(r8383, r8375, r8382, MPFR_RNDN); mpfr_pow(r8384, r8383, r8372, MPFR_RNDN); mpfr_add(r8385, r8373, r8384, MPFR_RNDN); mpfr_mul(r8386, r8371, r8371, MPFR_RNDN); mpfr_mul(r8387, r8383, r8371, MPFR_RNDN); mpfr_sub(r8388, r8386, r8387, MPFR_RNDN); mpfr_mul(r8389, r8350, r8382, MPFR_RNDN); mpfr_mul(r8390, r8389, r8389, MPFR_RNDN); mpfr_add(r8391, r8388, r8390, MPFR_RNDN); mpfr_div(r8392, r8385, r8391, MPFR_RNDN); mpfr_mul(r8393, r8392, r8392, MPFR_RNDN); mpfr_sub(r8394, r8358, r8393, MPFR_RNDN); mpfr_add(r8395, r8381, r8377, MPFR_RNDN); mpfr_pow(r8396, r8395, r8372, MPFR_RNDN); mpfr_mul(r8397, r8370, r8370, MPFR_RNDN); mpfr_mul(r8398, r8366, r8366, MPFR_RNDN); mpfr_sub(r8399, r8397, r8398, MPFR_RNDN); mpfr_pow(r8400, r8399, r8372, MPFR_RNDN); mpfr_mul(r8401, r8396, r8400, MPFR_RNDN); mpfr_mul(r8402, r8377, r8377, MPFR_RNDN); mpfr_mul(r8403, r8381, r8381, MPFR_RNDN); mpfr_sub(r8404, r8402, r8403, MPFR_RNDN); mpfr_mul(r8405, r8375, r8404, MPFR_RNDN); mpfr_pow(r8406, r8405, r8372, MPFR_RNDN); mpfr_sub(r8407, r8370, r8366, MPFR_RNDN); mpfr_pow(r8408, r8407, r8372, MPFR_RNDN); mpfr_mul(r8409, r8406, r8408, MPFR_RNDN); mpfr_add(r8410, r8401, r8409, MPFR_RNDN); mpfr_mul(r8411, r8408, r8396, MPFR_RNDN); mpfr_div(r8412, r8410, r8411, MPFR_RNDN); mpfr_div(r8413, r8412, r8391, MPFR_RNDN); mpfr_sub(r8414, r8357, r8413, MPFR_RNDN); mpfr_div(r8415, r8394, r8414, MPFR_RNDN); return mpfr_get_d(r8415, MPFR_RNDN); } static mpfr_t r8416, r8417, r8418, r8419, r8420, r8421, r8422, r8423, r8424, r8425, r8426, r8427, r8428, r8429, r8430, r8431, r8432, r8433, r8434, r8435, r8436, r8437, r8438, r8439, r8440, r8441, r8442, r8443, r8444, r8445, r8446, r8447, r8448, r8449, r8450, r8451, r8452, r8453, r8454, r8455, r8456, r8457, r8458, r8459, r8460, r8461, r8462, r8463, r8464, r8465, r8466, r8467, r8468, r8469, r8470, r8471, r8472, r8473, r8474, r8475, r8476, r8477, r8478, r8479, r8480, r8481, r8482; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8416, "77617", 10, MPFR_RNDN); mpfr_init(r8417); mpfr_init_set_str(r8418, "33096", 10, MPFR_RNDN); mpfr_init(r8419); mpfr_init(r8420); mpfr_init_set_str(r8421, "11", 10, MPFR_RNDN); mpfr_init(r8422); mpfr_init(r8423); mpfr_init(r8424); mpfr_init(r8425); mpfr_init_set_str(r8426, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r8427, "8", 10, MPFR_RNDN); mpfr_init(r8428); mpfr_init(r8429); mpfr_init_set_str(r8430, "2", 10, MPFR_RNDN); mpfr_init(r8431); mpfr_init(r8432); mpfr_init(r8433); mpfr_init_set_str(r8434, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r8435, "6", 10, MPFR_RNDN); mpfr_init(r8436); mpfr_init(r8437); mpfr_init(r8438); mpfr_init_set_str(r8439, "3", 10, MPFR_RNDN); mpfr_init(r8440); mpfr_init(r8441); mpfr_init(r8442); mpfr_init_set_str(r8443, "-2", 10, MPFR_RNDN); mpfr_init(r8444); mpfr_init_set_str(r8445, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r8446, "4", 10, MPFR_RNDN); mpfr_init(r8447); mpfr_init(r8448); mpfr_init(r8449); mpfr_init(r8450); mpfr_init(r8451); mpfr_init(r8452); mpfr_init(r8453); mpfr_init(r8454); mpfr_init(r8455); mpfr_init(r8456); mpfr_init(r8457); mpfr_init(r8458); mpfr_init(r8459); mpfr_init(r8460); mpfr_init(r8461); mpfr_init(r8462); mpfr_init(r8463); mpfr_init(r8464); mpfr_init(r8465); mpfr_init(r8466); mpfr_init(r8467); mpfr_init(r8468); mpfr_init(r8469); mpfr_init(r8470); mpfr_init(r8471); mpfr_init(r8472); mpfr_init(r8473); mpfr_init(r8474); mpfr_init(r8475); mpfr_init(r8476); mpfr_init(r8477); mpfr_init(r8478); mpfr_init(r8479); mpfr_init(r8480); mpfr_init(r8481); mpfr_init(r8482); } double f_dm() { ; mpfr_mul(r8417, r8416, r8416, MPFR_RNDN); ; mpfr_mul(r8419, r8416, r8418, MPFR_RNDN); mpfr_mul(r8420, r8419, r8416, MPFR_RNDN); ; mpfr_mul(r8422, r8418, r8421, MPFR_RNDN); mpfr_mul(r8423, r8420, r8422, MPFR_RNDN); mpfr_mul(r8424, r8417, r8423, MPFR_RNDN); mpfr_mul(r8425, r8424, r8424, MPFR_RNDN); ; ; mpfr_pow(r8428, r8418, r8427, MPFR_RNDN); mpfr_mul(r8429, r8426, r8428, MPFR_RNDN); ; mpfr_div(r8431, r8416, r8430, MPFR_RNDN); mpfr_div(r8432, r8431, r8418, MPFR_RNDN); mpfr_add(r8433, r8429, r8432, MPFR_RNDN); ; ; mpfr_pow(r8436, r8418, r8435, MPFR_RNDN); mpfr_mul(r8437, r8434, r8436, MPFR_RNDN); mpfr_add(r8438, r8433, r8437, MPFR_RNDN); ; mpfr_pow(r8440, r8438, r8439, MPFR_RNDN); mpfr_neg(r8441, r8416, MPFR_RNDN); mpfr_mul(r8442, r8441, r8416, MPFR_RNDN); ; mpfr_sub(r8444, r8436, r8443, MPFR_RNDN); ; ; mpfr_pow(r8447, r8418, r8446, MPFR_RNDN); mpfr_mul(r8448, r8445, r8447, MPFR_RNDN); mpfr_sub(r8449, r8444, r8448, MPFR_RNDN); mpfr_mul(r8450, r8442, r8449, MPFR_RNDN); mpfr_pow(r8451, r8450, r8439, MPFR_RNDN); mpfr_add(r8452, r8440, r8451, MPFR_RNDN); mpfr_mul(r8453, r8438, r8438, MPFR_RNDN); mpfr_mul(r8454, r8450, r8438, MPFR_RNDN); mpfr_sub(r8455, r8453, r8454, MPFR_RNDN); mpfr_mul(r8456, r8417, r8449, MPFR_RNDN); mpfr_mul(r8457, r8456, r8456, MPFR_RNDN); mpfr_add(r8458, r8455, r8457, MPFR_RNDN); mpfr_div(r8459, r8452, r8458, MPFR_RNDN); mpfr_mul(r8460, r8459, r8459, MPFR_RNDN); mpfr_sub(r8461, r8425, r8460, MPFR_RNDN); mpfr_add(r8462, r8448, r8444, MPFR_RNDN); mpfr_pow(r8463, r8462, r8439, MPFR_RNDN); mpfr_mul(r8464, r8437, r8437, MPFR_RNDN); mpfr_mul(r8465, r8433, r8433, MPFR_RNDN); mpfr_sub(r8466, r8464, r8465, MPFR_RNDN); mpfr_pow(r8467, r8466, r8439, MPFR_RNDN); mpfr_mul(r8468, r8463, r8467, MPFR_RNDN); mpfr_mul(r8469, r8444, r8444, MPFR_RNDN); mpfr_mul(r8470, r8448, r8448, MPFR_RNDN); mpfr_sub(r8471, r8469, r8470, MPFR_RNDN); mpfr_mul(r8472, r8442, r8471, MPFR_RNDN); mpfr_pow(r8473, r8472, r8439, MPFR_RNDN); mpfr_sub(r8474, r8437, r8433, MPFR_RNDN); mpfr_pow(r8475, r8474, r8439, MPFR_RNDN); mpfr_mul(r8476, r8473, r8475, MPFR_RNDN); mpfr_add(r8477, r8468, r8476, MPFR_RNDN); mpfr_mul(r8478, r8475, r8463, MPFR_RNDN); mpfr_div(r8479, r8477, r8478, MPFR_RNDN); mpfr_div(r8480, r8479, r8458, MPFR_RNDN); mpfr_sub(r8481, r8424, r8480, MPFR_RNDN); mpfr_div(r8482, r8461, r8481, MPFR_RNDN); return mpfr_get_d(r8482, MPFR_RNDN); }