#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r27296 = 1; float r27297 = 0.1049934947; float r27298 = x; float r27299 = r27298 * r27298; float r27300 = r27297 * r27299; float r27301 = r27296 + r27300; float r27302 = 0.0424060604; float r27303 = r27299 * r27299; float r27304 = r27302 * r27303; float r27305 = r27301 + r27304; float r27306 = 0.0072644182; float r27307 = r27303 * r27299; float r27308 = r27306 * r27307; float r27309 = r27305 + r27308; float r27310 = 0.0005064034; float r27311 = r27307 * r27299; float r27312 = r27310 * r27311; float r27313 = r27309 + r27312; float r27314 = 0.0001789971; float r27315 = r27311 * r27299; float r27316 = r27314 * r27315; float r27317 = r27313 + r27316; float r27318 = 0.7715471019; float r27319 = r27318 * r27299; float r27320 = r27296 + r27319; float r27321 = 0.2909738639; float r27322 = r27321 * r27303; float r27323 = r27320 + r27322; float r27324 = 0.0694555761; float r27325 = r27324 * r27307; float r27326 = r27323 + r27325; float r27327 = 0.0140005442; float r27328 = r27327 * r27311; float r27329 = r27326 + r27328; float r27330 = 0.0008327945; float r27331 = r27330 * r27315; float r27332 = r27329 + r27331; float r27333 = 2; float r27334 = r27333 * r27314; float r27335 = r27315 * r27299; float r27336 = r27334 * r27335; float r27337 = r27332 + r27336; float r27338 = r27317 / r27337; float r27339 = r27338 * r27298; return r27339; } double f_id(double x) { double r27340 = 1; double r27341 = 0.1049934947; double r27342 = x; double r27343 = r27342 * r27342; double r27344 = r27341 * r27343; double r27345 = r27340 + r27344; double r27346 = 0.0424060604; double r27347 = r27343 * r27343; double r27348 = r27346 * r27347; double r27349 = r27345 + r27348; double r27350 = 0.0072644182; double r27351 = r27347 * r27343; double r27352 = r27350 * r27351; double r27353 = r27349 + r27352; double r27354 = 0.0005064034; double r27355 = r27351 * r27343; double r27356 = r27354 * r27355; double r27357 = r27353 + r27356; double r27358 = 0.0001789971; double r27359 = r27355 * r27343; double r27360 = r27358 * r27359; double r27361 = r27357 + r27360; double r27362 = 0.7715471019; double r27363 = r27362 * r27343; double r27364 = r27340 + r27363; double r27365 = 0.2909738639; double r27366 = r27365 * r27347; double r27367 = r27364 + r27366; double r27368 = 0.0694555761; double r27369 = r27368 * r27351; double r27370 = r27367 + r27369; double r27371 = 0.0140005442; double r27372 = r27371 * r27355; double r27373 = r27370 + r27372; double r27374 = 0.0008327945; double r27375 = r27374 * r27359; double r27376 = r27373 + r27375; double r27377 = 2; double r27378 = r27377 * r27358; double r27379 = r27359 * r27343; double r27380 = r27378 * r27379; double r27381 = r27376 + r27380; double r27382 = r27361 / r27381; double r27383 = r27382 * r27342; return r27383; } double f_of(float x) { float r27384 = x; float r27385 = -89136323189.05368; bool r27386 = r27384 <= r27385; float r27387 = 849.5240389471198; bool r27388 = r27384 <= r27387; float r27389 = !r27388; bool r27390 = r27386 || r27389; float r27391 = 0.2514179000665375; float r27392 = 4; float r27393 = pow(r27384, r27392); float r27394 = r27391 / r27393; float r27395 = 0.15298196345929327; float r27396 = 6; float r27397 = pow(r27384, r27396); float r27398 = r27395 / r27397; float r27399 = r27394 + r27398; float r27400 = 0.5; float r27401 = r27400 / r27384; float r27402 = fma(r27399, r27384, r27401); float r27403 = r27384 * r27384; float r27404 = r27403 * r27403; float r27405 = r27404 * r27403; float r27406 = r27403 * r27405; float r27407 = r27403 * r27406; float r27408 = 0.0001789971; float r27409 = r27407 * r27408; float r27410 = 0.0072644182; float r27411 = r27410 * r27405; float r27412 = 0.1049934947; float r27413 = r27403 * r27412; float r27414 = 1; float r27415 = r27413 + r27414; float r27416 = 0.0424060604; float r27417 = r27404 * r27416; float r27418 = r27415 + r27417; float r27419 = r27411 + r27418; float r27420 = 0.0005064034; float r27421 = r27420 * r27406; float r27422 = r27419 + r27421; float r27423 = r27409 + r27422; float r27424 = 0.0008327945; float r27425 = r27407 * r27424; float r27426 = 0.7715471019; float r27427 = r27426 * r27403; float r27428 = r27427 + r27414; float r27429 = 0.2909738639; float r27430 = r27404 * r27429; float r27431 = r27428 + r27430; float r27432 = 0.0694555761; float r27433 = r27405 * r27432; float r27434 = r27431 + r27433; float r27435 = 0.0140005442; float r27436 = r27435 * r27406; float r27437 = r27434 + r27436; float r27438 = r27425 + r27437; float r27439 = 2; float r27440 = r27439 * r27408; float r27441 = r27403 * r27407; float r27442 = r27440 * r27441; float r27443 = r27438 + r27442; float r27444 = r27423 / r27443; float r27445 = r27384 * r27444; float r27446 = r27390 ? r27402 : r27445; return r27446; } double f_od(double x) { double r27447 = x; double r27448 = -89136323189.05368; bool r27449 = r27447 <= r27448; double r27450 = 849.5240389471198; bool r27451 = r27447 <= r27450; double r27452 = !r27451; bool r27453 = r27449 || r27452; double r27454 = 0.2514179000665375; double r27455 = 4; double r27456 = pow(r27447, r27455); double r27457 = r27454 / r27456; double r27458 = 0.15298196345929327; double r27459 = 6; double r27460 = pow(r27447, r27459); double r27461 = r27458 / r27460; double r27462 = r27457 + r27461; double r27463 = 0.5; double r27464 = r27463 / r27447; double r27465 = fma(r27462, r27447, r27464); double r27466 = r27447 * r27447; double r27467 = r27466 * r27466; double r27468 = r27467 * r27466; double r27469 = r27466 * r27468; double r27470 = r27466 * r27469; double r27471 = 0.0001789971; double r27472 = r27470 * r27471; double r27473 = 0.0072644182; double r27474 = r27473 * r27468; double r27475 = 0.1049934947; double r27476 = r27466 * r27475; double r27477 = 1; double r27478 = r27476 + r27477; double r27479 = 0.0424060604; double r27480 = r27467 * r27479; double r27481 = r27478 + r27480; double r27482 = r27474 + r27481; double r27483 = 0.0005064034; double r27484 = r27483 * r27469; double r27485 = r27482 + r27484; double r27486 = r27472 + r27485; double r27487 = 0.0008327945; double r27488 = r27470 * r27487; double r27489 = 0.7715471019; double r27490 = r27489 * r27466; double r27491 = r27490 + r27477; double r27492 = 0.2909738639; double r27493 = r27467 * r27492; double r27494 = r27491 + r27493; double r27495 = 0.0694555761; double r27496 = r27468 * r27495; double r27497 = r27494 + r27496; double r27498 = 0.0140005442; double r27499 = r27498 * r27469; double r27500 = r27497 + r27499; double r27501 = r27488 + r27500; double r27502 = 2; double r27503 = r27502 * r27471; double r27504 = r27466 * r27470; double r27505 = r27503 * r27504; double r27506 = r27501 + r27505; double r27507 = r27486 / r27506; double r27508 = r27447 * r27507; double r27509 = r27453 ? r27465 : r27508; return r27509; } 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 r27510, r27511, r27512, r27513, r27514, r27515, r27516, r27517, r27518, r27519, r27520, r27521, r27522, r27523, r27524, r27525, r27526, r27527, r27528, r27529, r27530, r27531, r27532, r27533, r27534, r27535, r27536, r27537, r27538, r27539, r27540, r27541, r27542, r27543, r27544, r27545, r27546, r27547, r27548, r27549, r27550, r27551, r27552, r27553; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27510, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27511, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27512); mpfr_init(r27513); mpfr_init(r27514); mpfr_init(r27515); mpfr_init_set_str(r27516, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27517); mpfr_init(r27518); mpfr_init(r27519); mpfr_init_set_str(r27520, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27521); mpfr_init(r27522); mpfr_init(r27523); mpfr_init_set_str(r27524, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27525); mpfr_init(r27526); mpfr_init(r27527); mpfr_init_set_str(r27528, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27529); mpfr_init(r27530); mpfr_init(r27531); mpfr_init_set_str(r27532, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27533); mpfr_init(r27534); mpfr_init_set_str(r27535, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27536); mpfr_init(r27537); mpfr_init_set_str(r27538, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27539); mpfr_init(r27540); mpfr_init_set_str(r27541, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27542); mpfr_init(r27543); mpfr_init_set_str(r27544, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27545); mpfr_init(r27546); mpfr_init_set_str(r27547, "2", 10, MPFR_RNDN); mpfr_init(r27548); mpfr_init(r27549); mpfr_init(r27550); mpfr_init(r27551); mpfr_init(r27552); mpfr_init(r27553); } double f_im(double x) { ; ; mpfr_set_d(r27512, x, MPFR_RNDN); mpfr_mul(r27513, r27512, r27512, MPFR_RNDN); mpfr_mul(r27514, r27511, r27513, MPFR_RNDN); mpfr_add(r27515, r27510, r27514, MPFR_RNDN); ; mpfr_mul(r27517, r27513, r27513, MPFR_RNDN); mpfr_mul(r27518, r27516, r27517, MPFR_RNDN); mpfr_add(r27519, r27515, r27518, MPFR_RNDN); ; mpfr_mul(r27521, r27517, r27513, MPFR_RNDN); mpfr_mul(r27522, r27520, r27521, MPFR_RNDN); mpfr_add(r27523, r27519, r27522, MPFR_RNDN); ; mpfr_mul(r27525, r27521, r27513, MPFR_RNDN); mpfr_mul(r27526, r27524, r27525, MPFR_RNDN); mpfr_add(r27527, r27523, r27526, MPFR_RNDN); ; mpfr_mul(r27529, r27525, r27513, MPFR_RNDN); mpfr_mul(r27530, r27528, r27529, MPFR_RNDN); mpfr_add(r27531, r27527, r27530, MPFR_RNDN); ; mpfr_mul(r27533, r27532, r27513, MPFR_RNDN); mpfr_add(r27534, r27510, r27533, MPFR_RNDN); ; mpfr_mul(r27536, r27535, r27517, MPFR_RNDN); mpfr_add(r27537, r27534, r27536, MPFR_RNDN); ; mpfr_mul(r27539, r27538, r27521, MPFR_RNDN); mpfr_add(r27540, r27537, r27539, MPFR_RNDN); ; mpfr_mul(r27542, r27541, r27525, MPFR_RNDN); mpfr_add(r27543, r27540, r27542, MPFR_RNDN); ; mpfr_mul(r27545, r27544, r27529, MPFR_RNDN); mpfr_add(r27546, r27543, r27545, MPFR_RNDN); ; mpfr_mul(r27548, r27547, r27528, MPFR_RNDN); mpfr_mul(r27549, r27529, r27513, MPFR_RNDN); mpfr_mul(r27550, r27548, r27549, MPFR_RNDN); mpfr_add(r27551, r27546, r27550, MPFR_RNDN); mpfr_div(r27552, r27531, r27551, MPFR_RNDN); mpfr_mul(r27553, r27552, r27512, MPFR_RNDN); return mpfr_get_d(r27553, MPFR_RNDN); } static mpfr_t r27554, r27555, r27556, r27557, r27558, r27559, r27560, r27561, r27562, r27563, r27564, r27565, r27566, r27567, r27568, r27569, r27570, r27571, r27572, r27573, r27574, r27575, r27576, r27577, r27578, r27579, r27580, r27581, r27582, r27583, r27584, r27585, r27586, r27587, r27588, r27589, r27590, r27591, r27592, r27593, r27594, r27595, r27596, r27597, r27598, r27599, r27600, r27601, r27602, r27603, r27604, r27605, r27606, r27607, r27608, r27609, r27610, r27611, r27612, r27613, r27614, r27615, r27616; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27554); mpfr_init_set_str(r27555, "-89136323189.05368", 10, MPFR_RNDN); mpfr_init(r27556); mpfr_init_set_str(r27557, "849.5240389471198", 10, MPFR_RNDN); mpfr_init(r27558); mpfr_init(r27559); mpfr_init(r27560); mpfr_init_set_str(r27561, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r27562, "4", 10, MPFR_RNDN); mpfr_init(r27563); mpfr_init(r27564); mpfr_init_set_str(r27565, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27566, "6", 10, MPFR_RNDN); mpfr_init(r27567); mpfr_init(r27568); mpfr_init(r27569); mpfr_init_set_str(r27570, "0.5", 10, MPFR_RNDN); mpfr_init(r27571); mpfr_init(r27572); mpfr_init(r27573); mpfr_init(r27574); mpfr_init(r27575); mpfr_init(r27576); mpfr_init(r27577); mpfr_init_set_str(r27578, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27579); mpfr_init_set_str(r27580, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27581); mpfr_init_set_str(r27582, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27583); mpfr_init_set_str(r27584, "1", 10, MPFR_RNDN); mpfr_init(r27585); mpfr_init_set_str(r27586, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27587); mpfr_init(r27588); mpfr_init(r27589); mpfr_init_set_str(r27590, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27591); mpfr_init(r27592); mpfr_init(r27593); mpfr_init_set_str(r27594, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27595); mpfr_init_set_str(r27596, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27597); mpfr_init(r27598); mpfr_init_set_str(r27599, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27600); mpfr_init(r27601); mpfr_init_set_str(r27602, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27603); mpfr_init(r27604); mpfr_init_set_str(r27605, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27606); mpfr_init(r27607); mpfr_init(r27608); mpfr_init_set_str(r27609, "2", 10, MPFR_RNDN); mpfr_init(r27610); mpfr_init(r27611); mpfr_init(r27612); mpfr_init(r27613); mpfr_init(r27614); mpfr_init(r27615); mpfr_init(r27616); } double f_fm(double x) { mpfr_set_d(r27554, x, MPFR_RNDN); ; mpfr_set_si(r27556, mpfr_cmp(r27554, r27555) <= 0, MPFR_RNDN); ; mpfr_set_si(r27558, mpfr_cmp(r27554, r27557) <= 0, MPFR_RNDN); mpfr_set_si(r27559, !mpfr_get_si(r27558, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27560, mpfr_get_si(r27556, MPFR_RNDN) || mpfr_get_si(r27559, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r27563, r27554, r27562, MPFR_RNDN); mpfr_div(r27564, r27561, r27563, MPFR_RNDN); ; ; mpfr_pow(r27567, r27554, r27566, MPFR_RNDN); mpfr_div(r27568, r27565, r27567, MPFR_RNDN); mpfr_add(r27569, r27564, r27568, MPFR_RNDN); ; mpfr_div(r27571, r27570, r27554, MPFR_RNDN); mpfr_fma(r27572, r27569, r27554, r27571, MPFR_RNDN); mpfr_mul(r27573, r27554, r27554, MPFR_RNDN); mpfr_mul(r27574, r27573, r27573, MPFR_RNDN); mpfr_mul(r27575, r27574, r27573, MPFR_RNDN); mpfr_mul(r27576, r27573, r27575, MPFR_RNDN); mpfr_mul(r27577, r27573, r27576, MPFR_RNDN); ; mpfr_mul(r27579, r27577, r27578, MPFR_RNDN); ; mpfr_mul(r27581, r27580, r27575, MPFR_RNDN); ; mpfr_mul(r27583, r27573, r27582, MPFR_RNDN); ; mpfr_add(r27585, r27583, r27584, MPFR_RNDN); ; mpfr_mul(r27587, r27574, r27586, MPFR_RNDN); mpfr_add(r27588, r27585, r27587, MPFR_RNDN); mpfr_add(r27589, r27581, r27588, MPFR_RNDN); ; mpfr_mul(r27591, r27590, r27576, MPFR_RNDN); mpfr_add(r27592, r27589, r27591, MPFR_RNDN); mpfr_add(r27593, r27579, r27592, MPFR_RNDN); ; mpfr_mul(r27595, r27577, r27594, MPFR_RNDN); ; mpfr_mul(r27597, r27596, r27573, MPFR_RNDN); mpfr_add(r27598, r27597, r27584, MPFR_RNDN); ; mpfr_mul(r27600, r27574, r27599, MPFR_RNDN); mpfr_add(r27601, r27598, r27600, MPFR_RNDN); ; mpfr_mul(r27603, r27575, r27602, MPFR_RNDN); mpfr_add(r27604, r27601, r27603, MPFR_RNDN); ; mpfr_mul(r27606, r27605, r27576, MPFR_RNDN); mpfr_add(r27607, r27604, r27606, MPFR_RNDN); mpfr_add(r27608, r27595, r27607, MPFR_RNDN); ; mpfr_mul(r27610, r27609, r27578, MPFR_RNDN); mpfr_mul(r27611, r27573, r27577, MPFR_RNDN); mpfr_mul(r27612, r27610, r27611, MPFR_RNDN); mpfr_add(r27613, r27608, r27612, MPFR_RNDN); mpfr_div(r27614, r27593, r27613, MPFR_RNDN); mpfr_mul(r27615, r27554, r27614, MPFR_RNDN); if (mpfr_get_si(r27560, MPFR_RNDN)) { mpfr_set(r27616, r27572, MPFR_RNDN); } else { mpfr_set(r27616, r27615, MPFR_RNDN); }; return mpfr_get_d(r27616, MPFR_RNDN); } static mpfr_t r27617, r27618, r27619, r27620, r27621, r27622, r27623, r27624, r27625, r27626, r27627, r27628, r27629, r27630, r27631, r27632, r27633, r27634, r27635, r27636, r27637, r27638, r27639, r27640, r27641, r27642, r27643, r27644, r27645, r27646, r27647, r27648, r27649, r27650, r27651, r27652, r27653, r27654, r27655, r27656, r27657, r27658, r27659, r27660, r27661, r27662, r27663, r27664, r27665, r27666, r27667, r27668, r27669, r27670, r27671, r27672, r27673, r27674, r27675, r27676, r27677, r27678, r27679; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r27617); mpfr_init_set_str(r27618, "-89136323189.05368", 10, MPFR_RNDN); mpfr_init(r27619); mpfr_init_set_str(r27620, "849.5240389471198", 10, MPFR_RNDN); mpfr_init(r27621); mpfr_init(r27622); mpfr_init(r27623); mpfr_init_set_str(r27624, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r27625, "4", 10, MPFR_RNDN); mpfr_init(r27626); mpfr_init(r27627); mpfr_init_set_str(r27628, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27629, "6", 10, MPFR_RNDN); mpfr_init(r27630); mpfr_init(r27631); mpfr_init(r27632); mpfr_init_set_str(r27633, "0.5", 10, MPFR_RNDN); mpfr_init(r27634); mpfr_init(r27635); mpfr_init(r27636); mpfr_init(r27637); mpfr_init(r27638); mpfr_init(r27639); mpfr_init(r27640); mpfr_init_set_str(r27641, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27642); mpfr_init_set_str(r27643, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27644); mpfr_init_set_str(r27645, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27646); mpfr_init_set_str(r27647, "1", 10, MPFR_RNDN); mpfr_init(r27648); mpfr_init_set_str(r27649, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27650); mpfr_init(r27651); mpfr_init(r27652); mpfr_init_set_str(r27653, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27654); mpfr_init(r27655); mpfr_init(r27656); mpfr_init_set_str(r27657, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27658); mpfr_init_set_str(r27659, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27660); mpfr_init(r27661); mpfr_init_set_str(r27662, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27663); mpfr_init(r27664); mpfr_init_set_str(r27665, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27666); mpfr_init(r27667); mpfr_init_set_str(r27668, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27669); mpfr_init(r27670); mpfr_init(r27671); mpfr_init_set_str(r27672, "2", 10, MPFR_RNDN); mpfr_init(r27673); mpfr_init(r27674); mpfr_init(r27675); mpfr_init(r27676); mpfr_init(r27677); mpfr_init(r27678); mpfr_init(r27679); } double f_dm(double x) { mpfr_set_d(r27617, x, MPFR_RNDN); ; mpfr_set_si(r27619, mpfr_cmp(r27617, r27618) <= 0, MPFR_RNDN); ; mpfr_set_si(r27621, mpfr_cmp(r27617, r27620) <= 0, MPFR_RNDN); mpfr_set_si(r27622, !mpfr_get_si(r27621, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r27623, mpfr_get_si(r27619, MPFR_RNDN) || mpfr_get_si(r27622, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r27626, r27617, r27625, MPFR_RNDN); mpfr_div(r27627, r27624, r27626, MPFR_RNDN); ; ; mpfr_pow(r27630, r27617, r27629, MPFR_RNDN); mpfr_div(r27631, r27628, r27630, MPFR_RNDN); mpfr_add(r27632, r27627, r27631, MPFR_RNDN); ; mpfr_div(r27634, r27633, r27617, MPFR_RNDN); mpfr_fma(r27635, r27632, r27617, r27634, MPFR_RNDN); mpfr_mul(r27636, r27617, r27617, MPFR_RNDN); mpfr_mul(r27637, r27636, r27636, MPFR_RNDN); mpfr_mul(r27638, r27637, r27636, MPFR_RNDN); mpfr_mul(r27639, r27636, r27638, MPFR_RNDN); mpfr_mul(r27640, r27636, r27639, MPFR_RNDN); ; mpfr_mul(r27642, r27640, r27641, MPFR_RNDN); ; mpfr_mul(r27644, r27643, r27638, MPFR_RNDN); ; mpfr_mul(r27646, r27636, r27645, MPFR_RNDN); ; mpfr_add(r27648, r27646, r27647, MPFR_RNDN); ; mpfr_mul(r27650, r27637, r27649, MPFR_RNDN); mpfr_add(r27651, r27648, r27650, MPFR_RNDN); mpfr_add(r27652, r27644, r27651, MPFR_RNDN); ; mpfr_mul(r27654, r27653, r27639, MPFR_RNDN); mpfr_add(r27655, r27652, r27654, MPFR_RNDN); mpfr_add(r27656, r27642, r27655, MPFR_RNDN); ; mpfr_mul(r27658, r27640, r27657, MPFR_RNDN); ; mpfr_mul(r27660, r27659, r27636, MPFR_RNDN); mpfr_add(r27661, r27660, r27647, MPFR_RNDN); ; mpfr_mul(r27663, r27637, r27662, MPFR_RNDN); mpfr_add(r27664, r27661, r27663, MPFR_RNDN); ; mpfr_mul(r27666, r27638, r27665, MPFR_RNDN); mpfr_add(r27667, r27664, r27666, MPFR_RNDN); ; mpfr_mul(r27669, r27668, r27639, MPFR_RNDN); mpfr_add(r27670, r27667, r27669, MPFR_RNDN); mpfr_add(r27671, r27658, r27670, MPFR_RNDN); ; mpfr_mul(r27673, r27672, r27641, MPFR_RNDN); mpfr_mul(r27674, r27636, r27640, MPFR_RNDN); mpfr_mul(r27675, r27673, r27674, MPFR_RNDN); mpfr_add(r27676, r27671, r27675, MPFR_RNDN); mpfr_div(r27677, r27656, r27676, MPFR_RNDN); mpfr_mul(r27678, r27617, r27677, MPFR_RNDN); if (mpfr_get_si(r27623, MPFR_RNDN)) { mpfr_set(r27679, r27635, MPFR_RNDN); } else { mpfr_set(r27679, r27678, MPFR_RNDN); }; return mpfr_get_d(r27679, MPFR_RNDN); }