#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r27343 = 1; float r27344 = 0.1049934947; float r27345 = x; float r27346 = r27345 * r27345; float r27347 = r27344 * r27346; float r27348 = r27343 + r27347; float r27349 = 0.0424060604; float r27350 = r27346 * r27346; float r27351 = r27349 * r27350; float r27352 = r27348 + r27351; float r27353 = 0.0072644182; float r27354 = r27350 * r27346; float r27355 = r27353 * r27354; float r27356 = r27352 + r27355; float r27357 = 0.0005064034; float r27358 = r27354 * r27346; float r27359 = r27357 * r27358; float r27360 = r27356 + r27359; float r27361 = 0.0001789971; float r27362 = r27358 * r27346; float r27363 = r27361 * r27362; float r27364 = r27360 + r27363; float r27365 = 0.7715471019; float r27366 = r27365 * r27346; float r27367 = r27343 + r27366; float r27368 = 0.2909738639; float r27369 = r27368 * r27350; float r27370 = r27367 + r27369; float r27371 = 0.0694555761; float r27372 = r27371 * r27354; float r27373 = r27370 + r27372; float r27374 = 0.0140005442; float r27375 = r27374 * r27358; float r27376 = r27373 + r27375; float r27377 = 0.0008327945; float r27378 = r27377 * r27362; float r27379 = r27376 + r27378; float r27380 = 2; float r27381 = r27380 * r27361; float r27382 = r27362 * r27346; float r27383 = r27381 * r27382; float r27384 = r27379 + r27383; float r27385 = r27364 / r27384; float r27386 = r27385 * r27345; return r27386; } double f_id(double x) { double r27387 = 1; double r27388 = 0.1049934947; double r27389 = x; double r27390 = r27389 * r27389; double r27391 = r27388 * r27390; double r27392 = r27387 + r27391; double r27393 = 0.0424060604; double r27394 = r27390 * r27390; double r27395 = r27393 * r27394; double r27396 = r27392 + r27395; double r27397 = 0.0072644182; double r27398 = r27394 * r27390; double r27399 = r27397 * r27398; double r27400 = r27396 + r27399; double r27401 = 0.0005064034; double r27402 = r27398 * r27390; double r27403 = r27401 * r27402; double r27404 = r27400 + r27403; double r27405 = 0.0001789971; double r27406 = r27402 * r27390; double r27407 = r27405 * r27406; double r27408 = r27404 + r27407; double r27409 = 0.7715471019; double r27410 = r27409 * r27390; double r27411 = r27387 + r27410; double r27412 = 0.2909738639; double r27413 = r27412 * r27394; double r27414 = r27411 + r27413; double r27415 = 0.0694555761; double r27416 = r27415 * r27398; double r27417 = r27414 + r27416; double r27418 = 0.0140005442; double r27419 = r27418 * r27402; double r27420 = r27417 + r27419; double r27421 = 0.0008327945; double r27422 = r27421 * r27406; double r27423 = r27420 + r27422; double r27424 = 2; double r27425 = r27424 * r27405; double r27426 = r27406 * r27390; double r27427 = r27425 * r27426; double r27428 = r27423 + r27427; double r27429 = r27408 / r27428; double r27430 = r27429 * r27389; return r27430; } double f_of(float x) { float r27431 = 0.2514179000665375; float r27432 = x; float r27433 = 4; float r27434 = pow(r27432, r27433); float r27435 = r27431 / r27434; float r27436 = 0.15298196345929327; float r27437 = 6; float r27438 = pow(r27432, r27437); float r27439 = r27436 / r27438; float r27440 = r27435 + r27439; float r27441 = 0.5; float r27442 = r27441 / r27432; float r27443 = fma(r27440, r27432, r27442); float r27444 = -0.0003066489264684035; bool r27445 = r27443 <= r27444; float r27446 = 0.0001789971; float r27447 = r27432 * r27432; float r27448 = r27446 * r27447; float r27449 = r27447 * r27447; float r27450 = r27449 * r27449; float r27451 = r27448 * r27450; float r27452 = 0.0005064034; float r27453 = 0.0072644182; float r27454 = r27447 * r27453; float r27455 = r27449 * r27454; float r27456 = 0.0424060604; float r27457 = 0.1049934947; float r27458 = 1; float r27459 = fma(r27447, r27457, r27458); float r27460 = fma(r27456, r27449, r27459); float r27461 = r27455 + r27460; float r27462 = fma(r27452, r27450, r27461); float r27463 = r27451 + r27462; float r27464 = 2; float r27465 = r27446 * r27464; float r27466 = r27449 * r27465; float r27467 = r27450 * r27466; float r27468 = 0.2909738639; float r27469 = r27468 * r27447; float r27470 = 0.7715471019; float r27471 = fma(r27470, r27447, r27458); float r27472 = fma(r27469, r27447, r27471); float r27473 = 0.0694555761; float r27474 = r27473 * r27447; float r27475 = r27474 * r27449; float r27476 = r27472 + r27475; float r27477 = r27467 + r27476; float r27478 = 0.0140005442; float r27479 = 0.0008327945; float r27480 = r27479 * r27447; float r27481 = r27480 * r27450; float r27482 = fma(r27478, r27450, r27481); float r27483 = r27477 + r27482; float r27484 = r27463 / r27483; float r27485 = log(r27484); float r27486 = exp(r27485); float r27487 = r27486 * r27432; float r27488 = 2.0832175626341745e-07; bool r27489 = r27443 <= r27488; float r27490 = r27457 * r27447; float r27491 = r27458 + r27490; float r27492 = r27456 * r27449; float r27493 = r27491 + r27492; float r27494 = r27449 * r27447; float r27495 = r27453 * r27494; float r27496 = r27493 + r27495; float r27497 = r27494 * r27447; float r27498 = r27452 * r27497; float r27499 = r27496 + r27498; float r27500 = r27497 * r27447; float r27501 = r27446 * r27500; float r27502 = r27499 + r27501; float r27503 = 3; float r27504 = r27503 + r27458; float r27505 = pow(r27447, r27504); float r27506 = r27465 * r27447; float r27507 = r27479 + r27506; float r27508 = r27505 * r27507; float r27509 = r27432 * r27473; float r27510 = pow(r27432, r27503); float r27511 = r27509 * r27510; float r27512 = r27432 * r27468; float r27513 = r27512 * r27510; float r27514 = fma(r27447, r27511, r27513); float r27515 = r27432 * r27470; float r27516 = fma(r27515, r27432, r27458); float r27517 = fma(r27505, r27478, r27516); float r27518 = r27514 + r27517; float r27519 = fma(r27447, r27508, r27518); float r27520 = r27458 / r27519; float r27521 = r27502 * r27520; float r27522 = r27521 * r27432; float r27523 = r27489 ? r27443 : r27522; float r27524 = r27445 ? r27487 : r27523; return r27524; } double f_od(double x) { double r27525 = 0.2514179000665375; double r27526 = x; double r27527 = 4; double r27528 = pow(r27526, r27527); double r27529 = r27525 / r27528; double r27530 = 0.15298196345929327; double r27531 = 6; double r27532 = pow(r27526, r27531); double r27533 = r27530 / r27532; double r27534 = r27529 + r27533; double r27535 = 0.5; double r27536 = r27535 / r27526; double r27537 = fma(r27534, r27526, r27536); double r27538 = -0.0003066489264684035; bool r27539 = r27537 <= r27538; double r27540 = 0.0001789971; double r27541 = r27526 * r27526; double r27542 = r27540 * r27541; double r27543 = r27541 * r27541; double r27544 = r27543 * r27543; double r27545 = r27542 * r27544; double r27546 = 0.0005064034; double r27547 = 0.0072644182; double r27548 = r27541 * r27547; double r27549 = r27543 * r27548; double r27550 = 0.0424060604; double r27551 = 0.1049934947; double r27552 = 1; double r27553 = fma(r27541, r27551, r27552); double r27554 = fma(r27550, r27543, r27553); double r27555 = r27549 + r27554; double r27556 = fma(r27546, r27544, r27555); double r27557 = r27545 + r27556; double r27558 = 2; double r27559 = r27540 * r27558; double r27560 = r27543 * r27559; double r27561 = r27544 * r27560; double r27562 = 0.2909738639; double r27563 = r27562 * r27541; double r27564 = 0.7715471019; double r27565 = fma(r27564, r27541, r27552); double r27566 = fma(r27563, r27541, r27565); double r27567 = 0.0694555761; double r27568 = r27567 * r27541; double r27569 = r27568 * r27543; double r27570 = r27566 + r27569; double r27571 = r27561 + r27570; double r27572 = 0.0140005442; double r27573 = 0.0008327945; double r27574 = r27573 * r27541; double r27575 = r27574 * r27544; double r27576 = fma(r27572, r27544, r27575); double r27577 = r27571 + r27576; double r27578 = r27557 / r27577; double r27579 = log(r27578); double r27580 = exp(r27579); double r27581 = r27580 * r27526; double r27582 = 2.0832175626341745e-07; bool r27583 = r27537 <= r27582; double r27584 = r27551 * r27541; double r27585 = r27552 + r27584; double r27586 = r27550 * r27543; double r27587 = r27585 + r27586; double r27588 = r27543 * r27541; double r27589 = r27547 * r27588; double r27590 = r27587 + r27589; double r27591 = r27588 * r27541; double r27592 = r27546 * r27591; double r27593 = r27590 + r27592; double r27594 = r27591 * r27541; double r27595 = r27540 * r27594; double r27596 = r27593 + r27595; double r27597 = 3; double r27598 = r27597 + r27552; double r27599 = pow(r27541, r27598); double r27600 = r27559 * r27541; double r27601 = r27573 + r27600; double r27602 = r27599 * r27601; double r27603 = r27526 * r27567; double r27604 = pow(r27526, r27597); double r27605 = r27603 * r27604; double r27606 = r27526 * r27562; double r27607 = r27606 * r27604; double r27608 = fma(r27541, r27605, r27607); double r27609 = r27526 * r27564; double r27610 = fma(r27609, r27526, r27552); double r27611 = fma(r27599, r27572, r27610); double r27612 = r27608 + r27611; double r27613 = fma(r27541, r27602, r27612); double r27614 = r27552 / r27613; double r27615 = r27596 * r27614; double r27616 = r27615 * r27526; double r27617 = r27583 ? r27537 : r27616; double r27618 = r27539 ? r27581 : r27617; return r27618; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27619, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27620, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27621); mpfr_init(r27622); mpfr_init(r27623); mpfr_init(r27624); mpfr_init_set_str(r27625, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27626); mpfr_init(r27627); mpfr_init(r27628); mpfr_init_set_str(r27629, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27630); mpfr_init(r27631); mpfr_init(r27632); mpfr_init_set_str(r27633, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27634); mpfr_init(r27635); mpfr_init(r27636); mpfr_init_set_str(r27637, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27638); mpfr_init(r27639); mpfr_init(r27640); mpfr_init_set_str(r27641, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27642); mpfr_init(r27643); mpfr_init_set_str(r27644, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27645); mpfr_init(r27646); mpfr_init_set_str(r27647, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27648); mpfr_init(r27649); mpfr_init_set_str(r27650, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27651); mpfr_init(r27652); mpfr_init_set_str(r27653, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27654); mpfr_init(r27655); mpfr_init_set_str(r27656, "2", 10, MPFR_RNDN); mpfr_init(r27657); mpfr_init(r27658); mpfr_init(r27659); mpfr_init(r27660); mpfr_init(r27661); mpfr_init(r27662); } double f_im(double x) { ; ; mpfr_set_d(r27621, x, MPFR_RNDN); mpfr_mul(r27622, r27621, r27621, MPFR_RNDN); mpfr_mul(r27623, r27620, r27622, MPFR_RNDN); mpfr_add(r27624, r27619, r27623, MPFR_RNDN); ; mpfr_mul(r27626, r27622, r27622, MPFR_RNDN); mpfr_mul(r27627, r27625, r27626, MPFR_RNDN); mpfr_add(r27628, r27624, r27627, MPFR_RNDN); ; mpfr_mul(r27630, r27626, r27622, MPFR_RNDN); mpfr_mul(r27631, r27629, r27630, MPFR_RNDN); mpfr_add(r27632, r27628, r27631, MPFR_RNDN); ; mpfr_mul(r27634, r27630, r27622, MPFR_RNDN); mpfr_mul(r27635, r27633, r27634, MPFR_RNDN); mpfr_add(r27636, r27632, r27635, MPFR_RNDN); ; mpfr_mul(r27638, r27634, r27622, MPFR_RNDN); mpfr_mul(r27639, r27637, r27638, MPFR_RNDN); mpfr_add(r27640, r27636, r27639, MPFR_RNDN); ; mpfr_mul(r27642, r27641, r27622, MPFR_RNDN); mpfr_add(r27643, r27619, r27642, MPFR_RNDN); ; mpfr_mul(r27645, r27644, r27626, MPFR_RNDN); mpfr_add(r27646, r27643, r27645, MPFR_RNDN); ; mpfr_mul(r27648, r27647, r27630, MPFR_RNDN); mpfr_add(r27649, r27646, r27648, MPFR_RNDN); ; mpfr_mul(r27651, r27650, r27634, MPFR_RNDN); mpfr_add(r27652, r27649, r27651, MPFR_RNDN); ; mpfr_mul(r27654, r27653, r27638, MPFR_RNDN); mpfr_add(r27655, r27652, r27654, MPFR_RNDN); ; mpfr_mul(r27657, r27656, r27637, MPFR_RNDN); mpfr_mul(r27658, r27638, r27622, MPFR_RNDN); mpfr_mul(r27659, r27657, r27658, MPFR_RNDN); mpfr_add(r27660, r27655, r27659, MPFR_RNDN); mpfr_div(r27661, r27640, r27660, MPFR_RNDN); mpfr_mul(r27662, r27661, r27621, MPFR_RNDN); return mpfr_get_d(r27662, MPFR_RNDN); } static mpfr_t r27663, r27664, r27665, r27666, r27667, r27668, r27669, r27670, r27671, r27672, r27673, r27674, r27675, r27676, r27677, r27678, r27679, r27680, r27681, r27682, r27683, r27684, r27685, r27686, r27687, r27688, r27689, r27690, r27691, r27692, r27693, r27694, r27695, r27696, r27697, r27698, r27699, r27700, r27701, r27702, r27703, r27704, r27705, r27706, r27707, r27708, r27709, r27710, r27711, r27712, r27713, r27714, r27715, r27716, r27717, r27718, r27719, r27720, r27721, r27722, r27723, r27724, r27725, r27726, r27727, r27728, r27729, r27730, r27731, r27732, r27733, r27734, r27735, r27736, r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755, r27756; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27663, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r27664); mpfr_init_set_str(r27665, "4", 10, MPFR_RNDN); mpfr_init(r27666); mpfr_init(r27667); mpfr_init_set_str(r27668, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27669, "6", 10, MPFR_RNDN); mpfr_init(r27670); mpfr_init(r27671); mpfr_init(r27672); mpfr_init_set_str(r27673, "0.5", 10, MPFR_RNDN); mpfr_init(r27674); mpfr_init(r27675); mpfr_init_set_str(r27676, "-0.0003066489264684035", 10, MPFR_RNDN); mpfr_init(r27677); mpfr_init_set_str(r27678, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27679); mpfr_init(r27680); mpfr_init(r27681); mpfr_init(r27682); mpfr_init(r27683); mpfr_init_set_str(r27684, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r27685, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27686); mpfr_init(r27687); mpfr_init_set_str(r27688, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r27689, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r27690, "1", 10, MPFR_RNDN); mpfr_init(r27691); mpfr_init(r27692); mpfr_init(r27693); mpfr_init(r27694); mpfr_init(r27695); mpfr_init_set_str(r27696, "2", 10, MPFR_RNDN); mpfr_init(r27697); mpfr_init(r27698); mpfr_init(r27699); mpfr_init_set_str(r27700, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27701); mpfr_init_set_str(r27702, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27703); mpfr_init(r27704); mpfr_init_set_str(r27705, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27706); mpfr_init(r27707); mpfr_init(r27708); mpfr_init(r27709); mpfr_init_set_str(r27710, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r27711, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27712); mpfr_init(r27713); mpfr_init(r27714); mpfr_init(r27715); mpfr_init(r27716); mpfr_init(r27717); mpfr_init(r27718); mpfr_init(r27719); mpfr_init_set_str(r27720, "2.0832175626341745e-07", 10, MPFR_RNDN); mpfr_init(r27721); mpfr_init(r27722); mpfr_init(r27723); mpfr_init(r27724); mpfr_init(r27725); mpfr_init(r27726); mpfr_init(r27727); mpfr_init(r27728); mpfr_init(r27729); mpfr_init(r27730); mpfr_init(r27731); mpfr_init(r27732); mpfr_init(r27733); mpfr_init(r27734); mpfr_init_set_str(r27735, "3", 10, MPFR_RNDN); mpfr_init(r27736); mpfr_init(r27737); mpfr_init(r27738); mpfr_init(r27739); mpfr_init(r27740); mpfr_init(r27741); mpfr_init(r27742); mpfr_init(r27743); mpfr_init(r27744); mpfr_init(r27745); mpfr_init(r27746); mpfr_init(r27747); mpfr_init(r27748); mpfr_init(r27749); mpfr_init(r27750); mpfr_init(r27751); mpfr_init(r27752); mpfr_init(r27753); mpfr_init(r27754); mpfr_init(r27755); mpfr_init(r27756); } double f_fm(double x) { ; mpfr_set_d(r27664, x, MPFR_RNDN); ; mpfr_pow(r27666, r27664, r27665, MPFR_RNDN); mpfr_div(r27667, r27663, r27666, MPFR_RNDN); ; ; mpfr_pow(r27670, r27664, r27669, MPFR_RNDN); mpfr_div(r27671, r27668, r27670, MPFR_RNDN); mpfr_add(r27672, r27667, r27671, MPFR_RNDN); ; mpfr_div(r27674, r27673, r27664, MPFR_RNDN); mpfr_fma(r27675, r27672, r27664, r27674, MPFR_RNDN); ; mpfr_set_si(r27677, mpfr_cmp(r27675, r27676) <= 0, MPFR_RNDN); ; mpfr_mul(r27679, r27664, r27664, MPFR_RNDN); mpfr_mul(r27680, r27678, r27679, MPFR_RNDN); mpfr_mul(r27681, r27679, r27679, MPFR_RNDN); mpfr_mul(r27682, r27681, r27681, MPFR_RNDN); mpfr_mul(r27683, r27680, r27682, MPFR_RNDN); ; ; mpfr_mul(r27686, r27679, r27685, MPFR_RNDN); mpfr_mul(r27687, r27681, r27686, MPFR_RNDN); ; ; ; mpfr_fma(r27691, r27679, r27689, r27690, MPFR_RNDN); mpfr_fma(r27692, r27688, r27681, r27691, MPFR_RNDN); mpfr_add(r27693, r27687, r27692, MPFR_RNDN); mpfr_fma(r27694, r27684, r27682, r27693, MPFR_RNDN); mpfr_add(r27695, r27683, r27694, MPFR_RNDN); ; mpfr_mul(r27697, r27678, r27696, MPFR_RNDN); mpfr_mul(r27698, r27681, r27697, MPFR_RNDN); mpfr_mul(r27699, r27682, r27698, MPFR_RNDN); ; mpfr_mul(r27701, r27700, r27679, MPFR_RNDN); ; mpfr_fma(r27703, r27702, r27679, r27690, MPFR_RNDN); mpfr_fma(r27704, r27701, r27679, r27703, MPFR_RNDN); ; mpfr_mul(r27706, r27705, r27679, MPFR_RNDN); mpfr_mul(r27707, r27706, r27681, MPFR_RNDN); mpfr_add(r27708, r27704, r27707, MPFR_RNDN); mpfr_add(r27709, r27699, r27708, MPFR_RNDN); ; ; mpfr_mul(r27712, r27711, r27679, MPFR_RNDN); mpfr_mul(r27713, r27712, r27682, MPFR_RNDN); mpfr_fma(r27714, r27710, r27682, r27713, MPFR_RNDN); mpfr_add(r27715, r27709, r27714, MPFR_RNDN); mpfr_div(r27716, r27695, r27715, MPFR_RNDN); mpfr_log(r27717, r27716, MPFR_RNDN); mpfr_exp(r27718, r27717, MPFR_RNDN); mpfr_mul(r27719, r27718, r27664, MPFR_RNDN); ; mpfr_set_si(r27721, mpfr_cmp(r27675, r27720) <= 0, MPFR_RNDN); mpfr_mul(r27722, r27689, r27679, MPFR_RNDN); mpfr_add(r27723, r27690, r27722, MPFR_RNDN); mpfr_mul(r27724, r27688, r27681, MPFR_RNDN); mpfr_add(r27725, r27723, r27724, MPFR_RNDN); mpfr_mul(r27726, r27681, r27679, MPFR_RNDN); mpfr_mul(r27727, r27685, r27726, MPFR_RNDN); mpfr_add(r27728, r27725, r27727, MPFR_RNDN); mpfr_mul(r27729, r27726, r27679, MPFR_RNDN); mpfr_mul(r27730, r27684, r27729, MPFR_RNDN); mpfr_add(r27731, r27728, r27730, MPFR_RNDN); mpfr_mul(r27732, r27729, r27679, MPFR_RNDN); mpfr_mul(r27733, r27678, r27732, MPFR_RNDN); mpfr_add(r27734, r27731, r27733, MPFR_RNDN); ; mpfr_add(r27736, r27735, r27690, MPFR_RNDN); mpfr_pow(r27737, r27679, r27736, MPFR_RNDN); mpfr_mul(r27738, r27697, r27679, MPFR_RNDN); mpfr_add(r27739, r27711, r27738, MPFR_RNDN); mpfr_mul(r27740, r27737, r27739, MPFR_RNDN); mpfr_mul(r27741, r27664, r27705, MPFR_RNDN); mpfr_pow(r27742, r27664, r27735, MPFR_RNDN); mpfr_mul(r27743, r27741, r27742, MPFR_RNDN); mpfr_mul(r27744, r27664, r27700, MPFR_RNDN); mpfr_mul(r27745, r27744, r27742, MPFR_RNDN); mpfr_fma(r27746, r27679, r27743, r27745, MPFR_RNDN); mpfr_mul(r27747, r27664, r27702, MPFR_RNDN); mpfr_fma(r27748, r27747, r27664, r27690, MPFR_RNDN); mpfr_fma(r27749, r27737, r27710, r27748, MPFR_RNDN); mpfr_add(r27750, r27746, r27749, MPFR_RNDN); mpfr_fma(r27751, r27679, r27740, r27750, MPFR_RNDN); mpfr_div(r27752, r27690, r27751, MPFR_RNDN); mpfr_mul(r27753, r27734, r27752, MPFR_RNDN); mpfr_mul(r27754, r27753, r27664, MPFR_RNDN); if (mpfr_get_si(r27721, MPFR_RNDN)) { mpfr_set(r27755, r27675, MPFR_RNDN); } else { mpfr_set(r27755, r27754, MPFR_RNDN); }; if (mpfr_get_si(r27677, MPFR_RNDN)) { mpfr_set(r27756, r27719, MPFR_RNDN); } else { mpfr_set(r27756, r27755, MPFR_RNDN); }; return mpfr_get_d(r27756, MPFR_RNDN); } static mpfr_t r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802, r27803, r27804, r27805, r27806, r27807, r27808, r27809, r27810, r27811, r27812, r27813, r27814, r27815, r27816, r27817, r27818, r27819, r27820, r27821, r27822, r27823, r27824, r27825, r27826, r27827, r27828, r27829, r27830, r27831, r27832, r27833, r27834, r27835, r27836, r27837, r27838, r27839, r27840, r27841, r27842, r27843, r27844, r27845, r27846, r27847, r27848, r27849, r27850; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27757, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r27758); mpfr_init_set_str(r27759, "4", 10, MPFR_RNDN); mpfr_init(r27760); mpfr_init(r27761); mpfr_init_set_str(r27762, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27763, "6", 10, MPFR_RNDN); mpfr_init(r27764); mpfr_init(r27765); mpfr_init(r27766); mpfr_init_set_str(r27767, "0.5", 10, MPFR_RNDN); mpfr_init(r27768); mpfr_init(r27769); mpfr_init_set_str(r27770, "-0.0003066489264684035", 10, MPFR_RNDN); mpfr_init(r27771); mpfr_init_set_str(r27772, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27773); mpfr_init(r27774); mpfr_init(r27775); mpfr_init(r27776); mpfr_init(r27777); mpfr_init_set_str(r27778, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r27779, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27780); mpfr_init(r27781); mpfr_init_set_str(r27782, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r27783, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r27784, "1", 10, MPFR_RNDN); mpfr_init(r27785); mpfr_init(r27786); mpfr_init(r27787); mpfr_init(r27788); mpfr_init(r27789); mpfr_init_set_str(r27790, "2", 10, MPFR_RNDN); mpfr_init(r27791); mpfr_init(r27792); mpfr_init(r27793); mpfr_init_set_str(r27794, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27795); mpfr_init_set_str(r27796, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27797); mpfr_init(r27798); mpfr_init_set_str(r27799, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27800); mpfr_init(r27801); mpfr_init(r27802); mpfr_init(r27803); mpfr_init_set_str(r27804, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r27805, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27806); mpfr_init(r27807); mpfr_init(r27808); mpfr_init(r27809); mpfr_init(r27810); mpfr_init(r27811); mpfr_init(r27812); mpfr_init(r27813); mpfr_init_set_str(r27814, "2.0832175626341745e-07", 10, MPFR_RNDN); mpfr_init(r27815); mpfr_init(r27816); mpfr_init(r27817); mpfr_init(r27818); mpfr_init(r27819); mpfr_init(r27820); mpfr_init(r27821); mpfr_init(r27822); mpfr_init(r27823); mpfr_init(r27824); mpfr_init(r27825); mpfr_init(r27826); mpfr_init(r27827); mpfr_init(r27828); mpfr_init_set_str(r27829, "3", 10, MPFR_RNDN); mpfr_init(r27830); mpfr_init(r27831); mpfr_init(r27832); mpfr_init(r27833); mpfr_init(r27834); mpfr_init(r27835); mpfr_init(r27836); mpfr_init(r27837); mpfr_init(r27838); mpfr_init(r27839); mpfr_init(r27840); mpfr_init(r27841); mpfr_init(r27842); mpfr_init(r27843); mpfr_init(r27844); mpfr_init(r27845); mpfr_init(r27846); mpfr_init(r27847); mpfr_init(r27848); mpfr_init(r27849); mpfr_init(r27850); } double f_dm(double x) { ; mpfr_set_d(r27758, x, MPFR_RNDN); ; mpfr_pow(r27760, r27758, r27759, MPFR_RNDN); mpfr_div(r27761, r27757, r27760, MPFR_RNDN); ; ; mpfr_pow(r27764, r27758, r27763, MPFR_RNDN); mpfr_div(r27765, r27762, r27764, MPFR_RNDN); mpfr_add(r27766, r27761, r27765, MPFR_RNDN); ; mpfr_div(r27768, r27767, r27758, MPFR_RNDN); mpfr_fma(r27769, r27766, r27758, r27768, MPFR_RNDN); ; mpfr_set_si(r27771, mpfr_cmp(r27769, r27770) <= 0, MPFR_RNDN); ; mpfr_mul(r27773, r27758, r27758, MPFR_RNDN); mpfr_mul(r27774, r27772, r27773, MPFR_RNDN); mpfr_mul(r27775, r27773, r27773, MPFR_RNDN); mpfr_mul(r27776, r27775, r27775, MPFR_RNDN); mpfr_mul(r27777, r27774, r27776, MPFR_RNDN); ; ; mpfr_mul(r27780, r27773, r27779, MPFR_RNDN); mpfr_mul(r27781, r27775, r27780, MPFR_RNDN); ; ; ; mpfr_fma(r27785, r27773, r27783, r27784, MPFR_RNDN); mpfr_fma(r27786, r27782, r27775, r27785, MPFR_RNDN); mpfr_add(r27787, r27781, r27786, MPFR_RNDN); mpfr_fma(r27788, r27778, r27776, r27787, MPFR_RNDN); mpfr_add(r27789, r27777, r27788, MPFR_RNDN); ; mpfr_mul(r27791, r27772, r27790, MPFR_RNDN); mpfr_mul(r27792, r27775, r27791, MPFR_RNDN); mpfr_mul(r27793, r27776, r27792, MPFR_RNDN); ; mpfr_mul(r27795, r27794, r27773, MPFR_RNDN); ; mpfr_fma(r27797, r27796, r27773, r27784, MPFR_RNDN); mpfr_fma(r27798, r27795, r27773, r27797, MPFR_RNDN); ; mpfr_mul(r27800, r27799, r27773, MPFR_RNDN); mpfr_mul(r27801, r27800, r27775, MPFR_RNDN); mpfr_add(r27802, r27798, r27801, MPFR_RNDN); mpfr_add(r27803, r27793, r27802, MPFR_RNDN); ; ; mpfr_mul(r27806, r27805, r27773, MPFR_RNDN); mpfr_mul(r27807, r27806, r27776, MPFR_RNDN); mpfr_fma(r27808, r27804, r27776, r27807, MPFR_RNDN); mpfr_add(r27809, r27803, r27808, MPFR_RNDN); mpfr_div(r27810, r27789, r27809, MPFR_RNDN); mpfr_log(r27811, r27810, MPFR_RNDN); mpfr_exp(r27812, r27811, MPFR_RNDN); mpfr_mul(r27813, r27812, r27758, MPFR_RNDN); ; mpfr_set_si(r27815, mpfr_cmp(r27769, r27814) <= 0, MPFR_RNDN); mpfr_mul(r27816, r27783, r27773, MPFR_RNDN); mpfr_add(r27817, r27784, r27816, MPFR_RNDN); mpfr_mul(r27818, r27782, r27775, MPFR_RNDN); mpfr_add(r27819, r27817, r27818, MPFR_RNDN); mpfr_mul(r27820, r27775, r27773, MPFR_RNDN); mpfr_mul(r27821, r27779, r27820, MPFR_RNDN); mpfr_add(r27822, r27819, r27821, MPFR_RNDN); mpfr_mul(r27823, r27820, r27773, MPFR_RNDN); mpfr_mul(r27824, r27778, r27823, MPFR_RNDN); mpfr_add(r27825, r27822, r27824, MPFR_RNDN); mpfr_mul(r27826, r27823, r27773, MPFR_RNDN); mpfr_mul(r27827, r27772, r27826, MPFR_RNDN); mpfr_add(r27828, r27825, r27827, MPFR_RNDN); ; mpfr_add(r27830, r27829, r27784, MPFR_RNDN); mpfr_pow(r27831, r27773, r27830, MPFR_RNDN); mpfr_mul(r27832, r27791, r27773, MPFR_RNDN); mpfr_add(r27833, r27805, r27832, MPFR_RNDN); mpfr_mul(r27834, r27831, r27833, MPFR_RNDN); mpfr_mul(r27835, r27758, r27799, MPFR_RNDN); mpfr_pow(r27836, r27758, r27829, MPFR_RNDN); mpfr_mul(r27837, r27835, r27836, MPFR_RNDN); mpfr_mul(r27838, r27758, r27794, MPFR_RNDN); mpfr_mul(r27839, r27838, r27836, MPFR_RNDN); mpfr_fma(r27840, r27773, r27837, r27839, MPFR_RNDN); mpfr_mul(r27841, r27758, r27796, MPFR_RNDN); mpfr_fma(r27842, r27841, r27758, r27784, MPFR_RNDN); mpfr_fma(r27843, r27831, r27804, r27842, MPFR_RNDN); mpfr_add(r27844, r27840, r27843, MPFR_RNDN); mpfr_fma(r27845, r27773, r27834, r27844, MPFR_RNDN); mpfr_div(r27846, r27784, r27845, MPFR_RNDN); mpfr_mul(r27847, r27828, r27846, MPFR_RNDN); mpfr_mul(r27848, r27847, r27758, MPFR_RNDN); if (mpfr_get_si(r27815, MPFR_RNDN)) { mpfr_set(r27849, r27769, MPFR_RNDN); } else { mpfr_set(r27849, r27848, MPFR_RNDN); }; if (mpfr_get_si(r27771, MPFR_RNDN)) { mpfr_set(r27850, r27813, MPFR_RNDN); } else { mpfr_set(r27850, r27849, MPFR_RNDN); }; return mpfr_get_d(r27850, MPFR_RNDN); }