#include #include #include #include #include char *name = "Numeric.SpecFunctions:invIncompleteBetaWorker from math-functions-0.1.5.2, I"; double f_if(float x, float y, float z, float t, float a, float b, float c) { float r40236 = x; float r40237 = y; float r40238 = 2.0f; float r40239 = z; float r40240 = t; float r40241 = a; float r40242 = r40240 + r40241; float r40243 = sqrt(r40242); float r40244 = r40239 * r40243; float r40245 = r40244 / r40240; float r40246 = b; float r40247 = c; float r40248 = r40246 - r40247; float r40249 = 5.0f; float r40250 = 6.0f; float r40251 = r40249 / r40250; float r40252 = r40241 + r40251; float r40253 = 3.0f; float r40254 = r40240 * r40253; float r40255 = r40238 / r40254; float r40256 = r40252 - r40255; float r40257 = r40248 * r40256; float r40258 = r40245 - r40257; float r40259 = r40238 * r40258; float r40260 = exp(r40259); float r40261 = r40237 * r40260; float r40262 = r40236 + r40261; float r40263 = r40236 / r40262; return r40263; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r40264 = x; double r40265 = y; double r40266 = 2.0; double r40267 = z; double r40268 = t; double r40269 = a; double r40270 = r40268 + r40269; double r40271 = sqrt(r40270); double r40272 = r40267 * r40271; double r40273 = r40272 / r40268; double r40274 = b; double r40275 = c; double r40276 = r40274 - r40275; double r40277 = 5.0; double r40278 = 6.0; double r40279 = r40277 / r40278; double r40280 = r40269 + r40279; double r40281 = 3.0; double r40282 = r40268 * r40281; double r40283 = r40266 / r40282; double r40284 = r40280 - r40283; double r40285 = r40276 * r40284; double r40286 = r40273 - r40285; double r40287 = r40266 * r40286; double r40288 = exp(r40287); double r40289 = r40265 * r40288; double r40290 = r40264 + r40289; double r40291 = r40264 / r40290; return r40291; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r40292 = t; float r40293 = -2.4368284028755833e+34f; bool r40294 = r40292 <= r40293; float r40295 = x; float r40296 = y; float r40297 = 2.0f; float r40298 = 0.8333333333333334f; float r40299 = c; float r40300 = r40298 * r40299; float r40301 = a; float r40302 = r40299 * r40301; float r40303 = r40300 + r40302; float r40304 = b; float r40305 = r40304 * r40301; float r40306 = r40303 - r40305; float r40307 = r40297 * r40306; float r40308 = exp(r40307); float r40309 = r40296 * r40308; float r40310 = r40295 + r40309; float r40311 = r40295 / r40310; float r40312 = 2.3351005256458113e-142f; bool r40313 = r40292 <= r40312; float r40314 = r40297 / r40292; float r40315 = 3.0f; float r40316 = r40314 / r40315; float r40317 = 5.0f; float r40318 = 6.0f; float r40319 = r40317 / r40318; float r40320 = r40319 + r40301; float r40321 = r40316 + r40320; float r40322 = z; float r40323 = r40301 + r40292; float r40324 = sqrt(r40323); float r40325 = r40322 * r40324; float r40326 = r40304 - r40299; float r40327 = r40326 * r40292; float r40328 = r40320 - r40316; float r40329 = r40327 * r40328; float r40330 = r40325 - r40329; float r40331 = r40321 * r40330; float r40332 = r40301 + r40319; float r40333 = r40292 * r40315; float r40334 = r40297 / r40333; float r40335 = r40332 + r40334; float r40336 = r40292 * r40335; float r40337 = r40331 / r40336; float r40338 = r40297 * r40337; float r40339 = exp(r40338); float r40340 = r40296 * r40339; float r40341 = r40295 + r40340; float r40342 = r40295 / r40341; float r40343 = r40292 + r40301; float r40344 = sqrt(r40343); float r40345 = r40322 * r40344; float r40346 = r40345 / r40292; float r40347 = r40332 - r40334; float r40348 = r40326 * r40347; float r40349 = r40346 - r40348; float r40350 = r40297 * r40349; float r40351 = exp(r40350); float r40352 = r40296 * r40351; float r40353 = r40295 + r40352; float r40354 = r40295 / r40353; float r40355 = r40313 ? r40342 : r40354; float r40356 = r40294 ? r40311 : r40355; return r40356; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r40357 = t; double r40358 = -2.4368284028755833e+34; bool r40359 = r40357 <= r40358; double r40360 = x; double r40361 = y; double r40362 = 2.0; double r40363 = 0.8333333333333334; double r40364 = c; double r40365 = r40363 * r40364; double r40366 = a; double r40367 = r40364 * r40366; double r40368 = r40365 + r40367; double r40369 = b; double r40370 = r40369 * r40366; double r40371 = r40368 - r40370; double r40372 = r40362 * r40371; double r40373 = exp(r40372); double r40374 = r40361 * r40373; double r40375 = r40360 + r40374; double r40376 = r40360 / r40375; double r40377 = 2.3351005256458113e-142; bool r40378 = r40357 <= r40377; double r40379 = r40362 / r40357; double r40380 = 3.0; double r40381 = r40379 / r40380; double r40382 = 5.0; double r40383 = 6.0; double r40384 = r40382 / r40383; double r40385 = r40384 + r40366; double r40386 = r40381 + r40385; double r40387 = z; double r40388 = r40366 + r40357; double r40389 = sqrt(r40388); double r40390 = r40387 * r40389; double r40391 = r40369 - r40364; double r40392 = r40391 * r40357; double r40393 = r40385 - r40381; double r40394 = r40392 * r40393; double r40395 = r40390 - r40394; double r40396 = r40386 * r40395; double r40397 = r40366 + r40384; double r40398 = r40357 * r40380; double r40399 = r40362 / r40398; double r40400 = r40397 + r40399; double r40401 = r40357 * r40400; double r40402 = r40396 / r40401; double r40403 = r40362 * r40402; double r40404 = exp(r40403); double r40405 = r40361 * r40404; double r40406 = r40360 + r40405; double r40407 = r40360 / r40406; double r40408 = r40357 + r40366; double r40409 = sqrt(r40408); double r40410 = r40387 * r40409; double r40411 = r40410 / r40357; double r40412 = r40397 - r40399; double r40413 = r40391 * r40412; double r40414 = r40411 - r40413; double r40415 = r40362 * r40414; double r40416 = exp(r40415); double r40417 = r40361 * r40416; double r40418 = r40360 + r40417; double r40419 = r40360 / r40418; double r40420 = r40378 ? r40407 : r40419; double r40421 = r40359 ? r40376 : r40420; return r40421; } 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 r40422, r40423, r40424, r40425, r40426, r40427, r40428, r40429, r40430, r40431, r40432, r40433, r40434, r40435, r40436, r40437, r40438, r40439, r40440, r40441, r40442, r40443, r40444, r40445, r40446, r40447, r40448, r40449; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r40422); mpfr_init(r40423); mpfr_init_set_str(r40424, "2.0", 10, MPFR_RNDN); mpfr_init(r40425); mpfr_init(r40426); mpfr_init(r40427); mpfr_init(r40428); mpfr_init(r40429); mpfr_init(r40430); mpfr_init(r40431); mpfr_init(r40432); mpfr_init(r40433); mpfr_init(r40434); mpfr_init_set_str(r40435, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r40436, "6.0", 10, MPFR_RNDN); mpfr_init(r40437); mpfr_init(r40438); mpfr_init_set_str(r40439, "3.0", 10, MPFR_RNDN); mpfr_init(r40440); mpfr_init(r40441); mpfr_init(r40442); mpfr_init(r40443); mpfr_init(r40444); mpfr_init(r40445); mpfr_init(r40446); mpfr_init(r40447); mpfr_init(r40448); mpfr_init(r40449); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r40422, x, MPFR_RNDN); mpfr_set_d(r40423, y, MPFR_RNDN); ; mpfr_set_d(r40425, z, MPFR_RNDN); mpfr_set_d(r40426, t, MPFR_RNDN); mpfr_set_d(r40427, a, MPFR_RNDN); mpfr_add(r40428, r40426, r40427, MPFR_RNDN); mpfr_sqrt(r40429, r40428, MPFR_RNDN); mpfr_mul(r40430, r40425, r40429, MPFR_RNDN); mpfr_div(r40431, r40430, r40426, MPFR_RNDN); mpfr_set_d(r40432, b, MPFR_RNDN); mpfr_set_d(r40433, c, MPFR_RNDN); mpfr_sub(r40434, r40432, r40433, MPFR_RNDN); ; ; mpfr_div(r40437, r40435, r40436, MPFR_RNDN); mpfr_add(r40438, r40427, r40437, MPFR_RNDN); ; mpfr_mul(r40440, r40426, r40439, MPFR_RNDN); mpfr_div(r40441, r40424, r40440, MPFR_RNDN); mpfr_sub(r40442, r40438, r40441, MPFR_RNDN); mpfr_mul(r40443, r40434, r40442, MPFR_RNDN); mpfr_sub(r40444, r40431, r40443, MPFR_RNDN); mpfr_mul(r40445, r40424, r40444, MPFR_RNDN); mpfr_exp(r40446, r40445, MPFR_RNDN); mpfr_mul(r40447, r40423, r40446, MPFR_RNDN); mpfr_add(r40448, r40422, r40447, MPFR_RNDN); mpfr_div(r40449, r40422, r40448, MPFR_RNDN); return mpfr_get_d(r40449, MPFR_RNDN); } static mpfr_t r40450, r40451, r40452, r40453, r40454, r40455, r40456, r40457, r40458, r40459, r40460, r40461, r40462, r40463, r40464, r40465, r40466, r40467, r40468, r40469, r40470, r40471, r40472, r40473, r40474, r40475, r40476, r40477, r40478, r40479, r40480, r40481, r40482, r40483, r40484, r40485, r40486, r40487, r40488, r40489, r40490, r40491, r40492, r40493, r40494, r40495, r40496, r40497, r40498, r40499, r40500, r40501, r40502, r40503, r40504, r40505, r40506, r40507, r40508, r40509, r40510, r40511, r40512, r40513, r40514; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r40450); mpfr_init_set_str(r40451, "-2.4368284028755833e+34", 10, MPFR_RNDN); mpfr_init(r40452); mpfr_init(r40453); mpfr_init(r40454); mpfr_init_set_str(r40455, "2.0", 10, MPFR_RNDN); mpfr_init_set_str(r40456, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r40457); mpfr_init(r40458); mpfr_init(r40459); mpfr_init(r40460); mpfr_init(r40461); mpfr_init(r40462); mpfr_init(r40463); mpfr_init(r40464); mpfr_init(r40465); mpfr_init(r40466); mpfr_init(r40467); mpfr_init(r40468); mpfr_init(r40469); mpfr_init_set_str(r40470, "2.3351005256458113e-142", 10, MPFR_RNDN); mpfr_init(r40471); mpfr_init(r40472); mpfr_init_set_str(r40473, "3.0", 10, MPFR_RNDN); mpfr_init(r40474); mpfr_init_set_str(r40475, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r40476, "6.0", 10, MPFR_RNDN); mpfr_init(r40477); mpfr_init(r40478); mpfr_init(r40479); mpfr_init(r40480); mpfr_init(r40481); mpfr_init(r40482); mpfr_init(r40483); mpfr_init(r40484); mpfr_init(r40485); mpfr_init(r40486); mpfr_init(r40487); mpfr_init(r40488); mpfr_init(r40489); mpfr_init(r40490); mpfr_init(r40491); mpfr_init(r40492); mpfr_init(r40493); mpfr_init(r40494); mpfr_init(r40495); mpfr_init(r40496); mpfr_init(r40497); mpfr_init(r40498); mpfr_init(r40499); mpfr_init(r40500); mpfr_init(r40501); mpfr_init(r40502); mpfr_init(r40503); mpfr_init(r40504); mpfr_init(r40505); mpfr_init(r40506); mpfr_init(r40507); mpfr_init(r40508); mpfr_init(r40509); mpfr_init(r40510); mpfr_init(r40511); mpfr_init(r40512); mpfr_init(r40513); mpfr_init(r40514); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r40450, t, MPFR_RNDN); ; mpfr_set_si(r40452, mpfr_cmp(r40450, r40451) <= 0, MPFR_RNDN); mpfr_set_d(r40453, x, MPFR_RNDN); mpfr_set_d(r40454, y, MPFR_RNDN); ; ; mpfr_set_d(r40457, c, MPFR_RNDN); mpfr_mul(r40458, r40456, r40457, MPFR_RNDN); mpfr_set_d(r40459, a, MPFR_RNDN); mpfr_mul(r40460, r40457, r40459, MPFR_RNDN); mpfr_add(r40461, r40458, r40460, MPFR_RNDN); mpfr_set_d(r40462, b, MPFR_RNDN); mpfr_mul(r40463, r40462, r40459, MPFR_RNDN); mpfr_sub(r40464, r40461, r40463, MPFR_RNDN); mpfr_mul(r40465, r40455, r40464, MPFR_RNDN); mpfr_exp(r40466, r40465, MPFR_RNDN); mpfr_mul(r40467, r40454, r40466, MPFR_RNDN); mpfr_add(r40468, r40453, r40467, MPFR_RNDN); mpfr_div(r40469, r40453, r40468, MPFR_RNDN); ; mpfr_set_si(r40471, mpfr_cmp(r40450, r40470) <= 0, MPFR_RNDN); mpfr_div(r40472, r40455, r40450, MPFR_RNDN); ; mpfr_div(r40474, r40472, r40473, MPFR_RNDN); ; ; mpfr_div(r40477, r40475, r40476, MPFR_RNDN); mpfr_add(r40478, r40477, r40459, MPFR_RNDN); mpfr_add(r40479, r40474, r40478, MPFR_RNDN); mpfr_set_d(r40480, z, MPFR_RNDN); mpfr_add(r40481, r40459, r40450, MPFR_RNDN); mpfr_sqrt(r40482, r40481, MPFR_RNDN); mpfr_mul(r40483, r40480, r40482, MPFR_RNDN); mpfr_sub(r40484, r40462, r40457, MPFR_RNDN); mpfr_mul(r40485, r40484, r40450, MPFR_RNDN); mpfr_sub(r40486, r40478, r40474, MPFR_RNDN); mpfr_mul(r40487, r40485, r40486, MPFR_RNDN); mpfr_sub(r40488, r40483, r40487, MPFR_RNDN); mpfr_mul(r40489, r40479, r40488, MPFR_RNDN); mpfr_add(r40490, r40459, r40477, MPFR_RNDN); mpfr_mul(r40491, r40450, r40473, MPFR_RNDN); mpfr_div(r40492, r40455, r40491, MPFR_RNDN); mpfr_add(r40493, r40490, r40492, MPFR_RNDN); mpfr_mul(r40494, r40450, r40493, MPFR_RNDN); mpfr_div(r40495, r40489, r40494, MPFR_RNDN); mpfr_mul(r40496, r40455, r40495, MPFR_RNDN); mpfr_exp(r40497, r40496, MPFR_RNDN); mpfr_mul(r40498, r40454, r40497, MPFR_RNDN); mpfr_add(r40499, r40453, r40498, MPFR_RNDN); mpfr_div(r40500, r40453, r40499, MPFR_RNDN); mpfr_add(r40501, r40450, r40459, MPFR_RNDN); mpfr_sqrt(r40502, r40501, MPFR_RNDN); mpfr_mul(r40503, r40480, r40502, MPFR_RNDN); mpfr_div(r40504, r40503, r40450, MPFR_RNDN); mpfr_sub(r40505, r40490, r40492, MPFR_RNDN); mpfr_mul(r40506, r40484, r40505, MPFR_RNDN); mpfr_sub(r40507, r40504, r40506, MPFR_RNDN); mpfr_mul(r40508, r40455, r40507, MPFR_RNDN); mpfr_exp(r40509, r40508, MPFR_RNDN); mpfr_mul(r40510, r40454, r40509, MPFR_RNDN); mpfr_add(r40511, r40453, r40510, MPFR_RNDN); mpfr_div(r40512, r40453, r40511, MPFR_RNDN); if (mpfr_get_si(r40471, MPFR_RNDN)) { mpfr_set(r40513, r40500, MPFR_RNDN); } else { mpfr_set(r40513, r40512, MPFR_RNDN); }; if (mpfr_get_si(r40452, MPFR_RNDN)) { mpfr_set(r40514, r40469, MPFR_RNDN); } else { mpfr_set(r40514, r40513, MPFR_RNDN); }; return mpfr_get_d(r40514, MPFR_RNDN); } static mpfr_t r40515, r40516, r40517, r40518, r40519, r40520, r40521, r40522, r40523, r40524, r40525, r40526, r40527, r40528, r40529, r40530, r40531, r40532, r40533, r40534, r40535, r40536, r40537, r40538, r40539, r40540, r40541, r40542, r40543, r40544, r40545, r40546, r40547, r40548, r40549, r40550, r40551, r40552, r40553, r40554, r40555, r40556, r40557, r40558, r40559, r40560, r40561, r40562, r40563, r40564, r40565, r40566, r40567, r40568, r40569, r40570, r40571, r40572, r40573, r40574, r40575, r40576, r40577, r40578, r40579; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r40515); mpfr_init_set_str(r40516, "-2.4368284028755833e+34", 10, MPFR_RNDN); mpfr_init(r40517); mpfr_init(r40518); mpfr_init(r40519); mpfr_init_set_str(r40520, "2.0", 10, MPFR_RNDN); mpfr_init_set_str(r40521, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r40522); mpfr_init(r40523); mpfr_init(r40524); mpfr_init(r40525); mpfr_init(r40526); mpfr_init(r40527); mpfr_init(r40528); mpfr_init(r40529); mpfr_init(r40530); mpfr_init(r40531); mpfr_init(r40532); mpfr_init(r40533); mpfr_init(r40534); mpfr_init_set_str(r40535, "2.3351005256458113e-142", 10, MPFR_RNDN); mpfr_init(r40536); mpfr_init(r40537); mpfr_init_set_str(r40538, "3.0", 10, MPFR_RNDN); mpfr_init(r40539); mpfr_init_set_str(r40540, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r40541, "6.0", 10, MPFR_RNDN); mpfr_init(r40542); mpfr_init(r40543); mpfr_init(r40544); mpfr_init(r40545); mpfr_init(r40546); mpfr_init(r40547); mpfr_init(r40548); mpfr_init(r40549); mpfr_init(r40550); mpfr_init(r40551); mpfr_init(r40552); mpfr_init(r40553); mpfr_init(r40554); mpfr_init(r40555); mpfr_init(r40556); mpfr_init(r40557); mpfr_init(r40558); mpfr_init(r40559); mpfr_init(r40560); mpfr_init(r40561); mpfr_init(r40562); mpfr_init(r40563); mpfr_init(r40564); mpfr_init(r40565); mpfr_init(r40566); mpfr_init(r40567); mpfr_init(r40568); mpfr_init(r40569); mpfr_init(r40570); mpfr_init(r40571); mpfr_init(r40572); mpfr_init(r40573); mpfr_init(r40574); mpfr_init(r40575); mpfr_init(r40576); mpfr_init(r40577); mpfr_init(r40578); mpfr_init(r40579); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r40515, t, MPFR_RNDN); ; mpfr_set_si(r40517, mpfr_cmp(r40515, r40516) <= 0, MPFR_RNDN); mpfr_set_d(r40518, x, MPFR_RNDN); mpfr_set_d(r40519, y, MPFR_RNDN); ; ; mpfr_set_d(r40522, c, MPFR_RNDN); mpfr_mul(r40523, r40521, r40522, MPFR_RNDN); mpfr_set_d(r40524, a, MPFR_RNDN); mpfr_mul(r40525, r40522, r40524, MPFR_RNDN); mpfr_add(r40526, r40523, r40525, MPFR_RNDN); mpfr_set_d(r40527, b, MPFR_RNDN); mpfr_mul(r40528, r40527, r40524, MPFR_RNDN); mpfr_sub(r40529, r40526, r40528, MPFR_RNDN); mpfr_mul(r40530, r40520, r40529, MPFR_RNDN); mpfr_exp(r40531, r40530, MPFR_RNDN); mpfr_mul(r40532, r40519, r40531, MPFR_RNDN); mpfr_add(r40533, r40518, r40532, MPFR_RNDN); mpfr_div(r40534, r40518, r40533, MPFR_RNDN); ; mpfr_set_si(r40536, mpfr_cmp(r40515, r40535) <= 0, MPFR_RNDN); mpfr_div(r40537, r40520, r40515, MPFR_RNDN); ; mpfr_div(r40539, r40537, r40538, MPFR_RNDN); ; ; mpfr_div(r40542, r40540, r40541, MPFR_RNDN); mpfr_add(r40543, r40542, r40524, MPFR_RNDN); mpfr_add(r40544, r40539, r40543, MPFR_RNDN); mpfr_set_d(r40545, z, MPFR_RNDN); mpfr_add(r40546, r40524, r40515, MPFR_RNDN); mpfr_sqrt(r40547, r40546, MPFR_RNDN); mpfr_mul(r40548, r40545, r40547, MPFR_RNDN); mpfr_sub(r40549, r40527, r40522, MPFR_RNDN); mpfr_mul(r40550, r40549, r40515, MPFR_RNDN); mpfr_sub(r40551, r40543, r40539, MPFR_RNDN); mpfr_mul(r40552, r40550, r40551, MPFR_RNDN); mpfr_sub(r40553, r40548, r40552, MPFR_RNDN); mpfr_mul(r40554, r40544, r40553, MPFR_RNDN); mpfr_add(r40555, r40524, r40542, MPFR_RNDN); mpfr_mul(r40556, r40515, r40538, MPFR_RNDN); mpfr_div(r40557, r40520, r40556, MPFR_RNDN); mpfr_add(r40558, r40555, r40557, MPFR_RNDN); mpfr_mul(r40559, r40515, r40558, MPFR_RNDN); mpfr_div(r40560, r40554, r40559, MPFR_RNDN); mpfr_mul(r40561, r40520, r40560, MPFR_RNDN); mpfr_exp(r40562, r40561, MPFR_RNDN); mpfr_mul(r40563, r40519, r40562, MPFR_RNDN); mpfr_add(r40564, r40518, r40563, MPFR_RNDN); mpfr_div(r40565, r40518, r40564, MPFR_RNDN); mpfr_add(r40566, r40515, r40524, MPFR_RNDN); mpfr_sqrt(r40567, r40566, MPFR_RNDN); mpfr_mul(r40568, r40545, r40567, MPFR_RNDN); mpfr_div(r40569, r40568, r40515, MPFR_RNDN); mpfr_sub(r40570, r40555, r40557, MPFR_RNDN); mpfr_mul(r40571, r40549, r40570, MPFR_RNDN); mpfr_sub(r40572, r40569, r40571, MPFR_RNDN); mpfr_mul(r40573, r40520, r40572, MPFR_RNDN); mpfr_exp(r40574, r40573, MPFR_RNDN); mpfr_mul(r40575, r40519, r40574, MPFR_RNDN); mpfr_add(r40576, r40518, r40575, MPFR_RNDN); mpfr_div(r40577, r40518, r40576, MPFR_RNDN); if (mpfr_get_si(r40536, MPFR_RNDN)) { mpfr_set(r40578, r40565, MPFR_RNDN); } else { mpfr_set(r40578, r40577, MPFR_RNDN); }; if (mpfr_get_si(r40517, MPFR_RNDN)) { mpfr_set(r40579, r40534, MPFR_RNDN); } else { mpfr_set(r40579, r40578, MPFR_RNDN); }; return mpfr_get_d(r40579, MPFR_RNDN); }