#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r19331 = 1.0f; float r19332 = 0.1049934947f; float r19333 = x; float r19334 = r19333 * r19333; float r19335 = r19332 * r19334; float r19336 = r19331 + r19335; float r19337 = 0.0424060604f; float r19338 = r19334 * r19334; float r19339 = r19337 * r19338; float r19340 = r19336 + r19339; float r19341 = 0.0072644182f; float r19342 = r19338 * r19334; float r19343 = r19341 * r19342; float r19344 = r19340 + r19343; float r19345 = 0.0005064034f; float r19346 = r19342 * r19334; float r19347 = r19345 * r19346; float r19348 = r19344 + r19347; float r19349 = 0.0001789971f; float r19350 = r19346 * r19334; float r19351 = r19349 * r19350; float r19352 = r19348 + r19351; float r19353 = 0.7715471019f; float r19354 = r19353 * r19334; float r19355 = r19331 + r19354; float r19356 = 0.2909738639f; float r19357 = r19356 * r19338; float r19358 = r19355 + r19357; float r19359 = 0.0694555761f; float r19360 = r19359 * r19342; float r19361 = r19358 + r19360; float r19362 = 0.0140005442f; float r19363 = r19362 * r19346; float r19364 = r19361 + r19363; float r19365 = 0.0008327945f; float r19366 = r19365 * r19350; float r19367 = r19364 + r19366; float r19368 = 2.0f; float r19369 = r19368 * r19349; float r19370 = r19350 * r19334; float r19371 = r19369 * r19370; float r19372 = r19367 + r19371; float r19373 = r19352 / r19372; float r19374 = r19373 * r19333; return r19374; } double f_id(double x) { double r19375 = 1.0; double r19376 = 0.1049934947; double r19377 = x; double r19378 = r19377 * r19377; double r19379 = r19376 * r19378; double r19380 = r19375 + r19379; double r19381 = 0.0424060604; double r19382 = r19378 * r19378; double r19383 = r19381 * r19382; double r19384 = r19380 + r19383; double r19385 = 0.0072644182; double r19386 = r19382 * r19378; double r19387 = r19385 * r19386; double r19388 = r19384 + r19387; double r19389 = 0.0005064034; double r19390 = r19386 * r19378; double r19391 = r19389 * r19390; double r19392 = r19388 + r19391; double r19393 = 0.0001789971; double r19394 = r19390 * r19378; double r19395 = r19393 * r19394; double r19396 = r19392 + r19395; double r19397 = 0.7715471019; double r19398 = r19397 * r19378; double r19399 = r19375 + r19398; double r19400 = 0.2909738639; double r19401 = r19400 * r19382; double r19402 = r19399 + r19401; double r19403 = 0.0694555761; double r19404 = r19403 * r19386; double r19405 = r19402 + r19404; double r19406 = 0.0140005442; double r19407 = r19406 * r19390; double r19408 = r19405 + r19407; double r19409 = 0.0008327945; double r19410 = r19409 * r19394; double r19411 = r19408 + r19410; double r19412 = 2.0; double r19413 = r19412 * r19393; double r19414 = r19394 * r19378; double r19415 = r19413 * r19414; double r19416 = r19411 + r19415; double r19417 = r19396 / r19416; double r19418 = r19417 * r19377; return r19418; } double f_of(float x) { float r19419 = x; float r19420 = -2.0717181681580965e+22f; bool r19421 = r19419 <= r19420; float r19422 = 1.0f; float r19423 = 0.1049934947f; float r19424 = r19419 * r19419; float r19425 = r19423 / r19424; float r19426 = r19422 + r19425; float r19427 = 0.0072644182f; float r19428 = r19427 * r19422; float r19429 = r19419 * (r19419 * r19419); float r19430 = r19429 * r19429; float r19431 = r19428 / r19430; float r19432 = r19426 + r19431; float r19433 = 0.0424060604f; float r19434 = r19433 * r19422; float r19435 = r19434 / r19429; float r19436 = r19435 / r19419; float r19437 = 0.0001789971f; float r19438 = 10.0f; float r19439 = pow(r19419, r19438); float r19440 = r19437 / r19439; float r19441 = r19436 + r19440; float r19442 = 0.0005064034f; float r19443 = r19422 * r19442; float r19444 = r19443 / r19430; float r19445 = r19444 / r19424; float r19446 = r19441 + r19445; float r19447 = r19432 + r19446; float r19448 = 0.0694555761f; float r19449 = r19422 * r19448; float r19450 = r19449 / r19430; float r19451 = r19422 + r19450; float r19452 = 0.2909738639f; float r19453 = r19422 / r19419; float r19454 = r19453 / r19429; float r19455 = r19452 * r19454; float r19456 = r19451 + r19455; float r19457 = 0.7715471019f; float r19458 = r19457 / r19424; float r19459 = r19422 / r19429; float r19460 = r19459 * r19459; float r19461 = 6.0f; float r19462 = pow(r19419, r19461); float r19463 = 0.0003579942f; float r19464 = r19462 / r19463; float r19465 = r19460 / r19464; float r19466 = r19458 + r19465; float r19467 = r19456 + r19466; float r19468 = 0.0140005442f; float r19469 = r19468 * r19422; float r19470 = r19469 / r19430; float r19471 = r19470 / r19424; float r19472 = 0.0008327945f; float r19473 = r19472 * r19422; float r19474 = r19473 / r19430; float r19475 = 4.0f; float r19476 = pow(r19419, r19475); float r19477 = r19474 / r19476; float r19478 = r19471 + r19477; float r19479 = r19467 + r19478; float r19480 = r19479 * r19419; float r19481 = r19447 / r19480; float r19482 = 1.8438098325834933e+28f; bool r19483 = r19419 <= r19482; float r19484 = r19419 * r19419; float r19485 = r19429 * r19484; float r19486 = r19485 * r19485; float r19487 = 2.0f; float r19488 = r19487 * r19437; float r19489 = r19484 * r19488; float r19490 = r19472 + r19489; float r19491 = r19486 * r19490; float r19492 = r19419 * r19452; float r19493 = r19492 * r19429; float r19494 = r19457 * r19419; float r19495 = r19419 * r19494; float r19496 = r19422 + r19495; float r19497 = r19493 + r19496; float r19498 = r19491 + r19497; float r19499 = r19484 * r19468; float r19500 = r19499 + r19448; float r19501 = r19484 * (r19484 * r19484); float r19502 = r19500 * r19501; float r19503 = r19498 + r19502; float r19504 = r19419 / r19503; float r19505 = r19419 * r19433; float r19506 = r19505 * r19429; float r19507 = r19424 * r19423; float r19508 = r19507 + r19422; float r19509 = r19506 + r19508; float r19510 = r19424 * r19437; float r19511 = r19424 * r19424; float r19512 = r19511 * r19511; float r19513 = r19510 * r19512; float r19514 = r19509 + r19513; float r19515 = r19442 * r19419; float r19516 = r19515 * r19419; float r19517 = r19427 + r19516; float r19518 = r19511 * r19517; float r19519 = r19424 * r19518; float r19520 = r19514 + r19519; float r19521 = r19504 * r19520; float r19522 = r19483 ? r19521 : r19481; float r19523 = r19421 ? r19481 : r19522; return r19523; } double f_od(double x) { double r19524 = x; double r19525 = -2.0717181681580965e+22; bool r19526 = r19524 <= r19525; double r19527 = 1.0; double r19528 = 0.1049934947; double r19529 = r19524 * r19524; double r19530 = r19528 / r19529; double r19531 = r19527 + r19530; double r19532 = 0.0072644182; double r19533 = r19532 * r19527; double r19534 = r19524 * (r19524 * r19524); double r19535 = r19534 * r19534; double r19536 = r19533 / r19535; double r19537 = r19531 + r19536; double r19538 = 0.0424060604; double r19539 = r19538 * r19527; double r19540 = r19539 / r19534; double r19541 = r19540 / r19524; double r19542 = 0.0001789971; double r19543 = 10.0; double r19544 = pow(r19524, r19543); double r19545 = r19542 / r19544; double r19546 = r19541 + r19545; double r19547 = 0.0005064034; double r19548 = r19527 * r19547; double r19549 = r19548 / r19535; double r19550 = r19549 / r19529; double r19551 = r19546 + r19550; double r19552 = r19537 + r19551; double r19553 = 0.0694555761; double r19554 = r19527 * r19553; double r19555 = r19554 / r19535; double r19556 = r19527 + r19555; double r19557 = 0.2909738639; double r19558 = r19527 / r19524; double r19559 = r19558 / r19534; double r19560 = r19557 * r19559; double r19561 = r19556 + r19560; double r19562 = 0.7715471019; double r19563 = r19562 / r19529; double r19564 = r19527 / r19534; double r19565 = r19564 * r19564; double r19566 = 6.0; double r19567 = pow(r19524, r19566); double r19568 = 0.0003579942; double r19569 = r19567 / r19568; double r19570 = r19565 / r19569; double r19571 = r19563 + r19570; double r19572 = r19561 + r19571; double r19573 = 0.0140005442; double r19574 = r19573 * r19527; double r19575 = r19574 / r19535; double r19576 = r19575 / r19529; double r19577 = 0.0008327945; double r19578 = r19577 * r19527; double r19579 = r19578 / r19535; double r19580 = 4.0; double r19581 = pow(r19524, r19580); double r19582 = r19579 / r19581; double r19583 = r19576 + r19582; double r19584 = r19572 + r19583; double r19585 = r19584 * r19524; double r19586 = r19552 / r19585; double r19587 = 1.8438098325834933e+28; bool r19588 = r19524 <= r19587; double r19589 = r19524 * r19524; double r19590 = r19534 * r19589; double r19591 = r19590 * r19590; double r19592 = 2.0; double r19593 = r19592 * r19542; double r19594 = r19589 * r19593; double r19595 = r19577 + r19594; double r19596 = r19591 * r19595; double r19597 = r19524 * r19557; double r19598 = r19597 * r19534; double r19599 = r19562 * r19524; double r19600 = r19524 * r19599; double r19601 = r19527 + r19600; double r19602 = r19598 + r19601; double r19603 = r19596 + r19602; double r19604 = r19589 * r19573; double r19605 = r19604 + r19553; double r19606 = r19589 * (r19589 * r19589); double r19607 = r19605 * r19606; double r19608 = r19603 + r19607; double r19609 = r19524 / r19608; double r19610 = r19524 * r19538; double r19611 = r19610 * r19534; double r19612 = r19529 * r19528; double r19613 = r19612 + r19527; double r19614 = r19611 + r19613; double r19615 = r19529 * r19542; double r19616 = r19529 * r19529; double r19617 = r19616 * r19616; double r19618 = r19615 * r19617; double r19619 = r19614 + r19618; double r19620 = r19547 * r19524; double r19621 = r19620 * r19524; double r19622 = r19532 + r19621; double r19623 = r19616 * r19622; double r19624 = r19529 * r19623; double r19625 = r19619 + r19624; double r19626 = r19609 * r19625; double r19627 = r19588 ? r19626 : r19586; double r19628 = r19526 ? r19586 : r19627; return r19628; } 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 r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650, r19651, r19652, r19653, r19654, r19655, r19656, r19657, r19658, r19659, r19660, r19661, r19662, r19663, r19664, r19665, r19666, r19667, r19668, r19669, r19670, r19671, r19672; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19629, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19630, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19631); mpfr_init(r19632); mpfr_init(r19633); mpfr_init(r19634); mpfr_init_set_str(r19635, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19636); mpfr_init(r19637); mpfr_init(r19638); mpfr_init_set_str(r19639, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19640); mpfr_init(r19641); mpfr_init(r19642); mpfr_init_set_str(r19643, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19644); mpfr_init(r19645); mpfr_init(r19646); mpfr_init_set_str(r19647, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19648); mpfr_init(r19649); mpfr_init(r19650); mpfr_init_set_str(r19651, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19652); mpfr_init(r19653); mpfr_init_set_str(r19654, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19655); mpfr_init(r19656); mpfr_init_set_str(r19657, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19658); mpfr_init(r19659); mpfr_init_set_str(r19660, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19661); mpfr_init(r19662); mpfr_init_set_str(r19663, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19664); mpfr_init(r19665); mpfr_init_set_str(r19666, "2", 10, MPFR_RNDN); mpfr_init(r19667); mpfr_init(r19668); mpfr_init(r19669); mpfr_init(r19670); mpfr_init(r19671); mpfr_init(r19672); } double f_im(double x) { ; ; mpfr_set_d(r19631, x, MPFR_RNDN); mpfr_mul(r19632, r19631, r19631, MPFR_RNDN); mpfr_mul(r19633, r19630, r19632, MPFR_RNDN); mpfr_add(r19634, r19629, r19633, MPFR_RNDN); ; mpfr_mul(r19636, r19632, r19632, MPFR_RNDN); mpfr_mul(r19637, r19635, r19636, MPFR_RNDN); mpfr_add(r19638, r19634, r19637, MPFR_RNDN); ; mpfr_mul(r19640, r19636, r19632, MPFR_RNDN); mpfr_mul(r19641, r19639, r19640, MPFR_RNDN); mpfr_add(r19642, r19638, r19641, MPFR_RNDN); ; mpfr_mul(r19644, r19640, r19632, MPFR_RNDN); mpfr_mul(r19645, r19643, r19644, MPFR_RNDN); mpfr_add(r19646, r19642, r19645, MPFR_RNDN); ; mpfr_mul(r19648, r19644, r19632, MPFR_RNDN); mpfr_mul(r19649, r19647, r19648, MPFR_RNDN); mpfr_add(r19650, r19646, r19649, MPFR_RNDN); ; mpfr_mul(r19652, r19651, r19632, MPFR_RNDN); mpfr_add(r19653, r19629, r19652, MPFR_RNDN); ; mpfr_mul(r19655, r19654, r19636, MPFR_RNDN); mpfr_add(r19656, r19653, r19655, MPFR_RNDN); ; mpfr_mul(r19658, r19657, r19640, MPFR_RNDN); mpfr_add(r19659, r19656, r19658, MPFR_RNDN); ; mpfr_mul(r19661, r19660, r19644, MPFR_RNDN); mpfr_add(r19662, r19659, r19661, MPFR_RNDN); ; mpfr_mul(r19664, r19663, r19648, MPFR_RNDN); mpfr_add(r19665, r19662, r19664, MPFR_RNDN); ; mpfr_mul(r19667, r19666, r19647, MPFR_RNDN); mpfr_mul(r19668, r19648, r19632, MPFR_RNDN); mpfr_mul(r19669, r19667, r19668, MPFR_RNDN); mpfr_add(r19670, r19665, r19669, MPFR_RNDN); mpfr_div(r19671, r19650, r19670, MPFR_RNDN); mpfr_mul(r19672, r19671, r19631, MPFR_RNDN); return mpfr_get_d(r19672, MPFR_RNDN); } static mpfr_t r19673, r19674, r19675, r19676, r19677, r19678, r19679, r19680, r19681, r19682, r19683, r19684, r19685, r19686, r19687, r19688, r19689, r19690, r19691, r19692, r19693, r19694, r19695, r19696, r19697, r19698, r19699, r19700, r19701, r19702, r19703, r19704, r19705, r19706, r19707, r19708, r19709, r19710, r19711, r19712, r19713, r19714, r19715, r19716, r19717, r19718, r19719, r19720, r19721, r19722, r19723, r19724, r19725, r19726, r19727, r19728, r19729, r19730, r19731, r19732, r19733, r19734, r19735, r19736, r19737, r19738, r19739, r19740, r19741, r19742, r19743, r19744, r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752, r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19673); mpfr_init_set_str(r19674, "-2.0717181681580965e+22", 10, MPFR_RNDN); mpfr_init(r19675); mpfr_init_set_str(r19676, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19677, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19678); mpfr_init(r19679); mpfr_init(r19680); mpfr_init_set_str(r19681, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19682); mpfr_init(r19683); mpfr_init(r19684); mpfr_init(r19685); mpfr_init(r19686); mpfr_init_set_str(r19687, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19688); mpfr_init(r19689); mpfr_init(r19690); mpfr_init_set_str(r19691, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19692, "10", 10, MPFR_RNDN); mpfr_init(r19693); mpfr_init(r19694); mpfr_init(r19695); mpfr_init_set_str(r19696, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19697); mpfr_init(r19698); mpfr_init(r19699); mpfr_init(r19700); mpfr_init(r19701); mpfr_init_set_str(r19702, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19703); mpfr_init(r19704); mpfr_init(r19705); mpfr_init_set_str(r19706, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19707); mpfr_init(r19708); mpfr_init(r19709); mpfr_init(r19710); mpfr_init_set_str(r19711, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19712); mpfr_init(r19713); mpfr_init(r19714); mpfr_init_set_str(r19715, "6", 10, MPFR_RNDN); mpfr_init(r19716); mpfr_init_set_str(r19717, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19718); mpfr_init(r19719); mpfr_init(r19720); mpfr_init(r19721); mpfr_init_set_str(r19722, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19723); mpfr_init(r19724); mpfr_init(r19725); mpfr_init_set_str(r19726, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19727); mpfr_init(r19728); mpfr_init_set_str(r19729, "4", 10, MPFR_RNDN); mpfr_init(r19730); mpfr_init(r19731); mpfr_init(r19732); mpfr_init(r19733); mpfr_init(r19734); mpfr_init(r19735); mpfr_init_set_str(r19736, "1.8438098325834933e+28", 10, MPFR_RNDN); mpfr_init(r19737); mpfr_init(r19738); mpfr_init(r19739); mpfr_init(r19740); mpfr_init_set_str(r19741, "2", 10, MPFR_RNDN); mpfr_init(r19742); mpfr_init(r19743); mpfr_init(r19744); mpfr_init(r19745); mpfr_init(r19746); mpfr_init(r19747); mpfr_init(r19748); mpfr_init(r19749); mpfr_init(r19750); mpfr_init(r19751); mpfr_init(r19752); mpfr_init(r19753); mpfr_init(r19754); mpfr_init(r19755); mpfr_init(r19756); mpfr_init(r19757); mpfr_init(r19758); mpfr_init(r19759); mpfr_init(r19760); mpfr_init(r19761); mpfr_init(r19762); mpfr_init(r19763); mpfr_init(r19764); mpfr_init(r19765); mpfr_init(r19766); mpfr_init(r19767); mpfr_init(r19768); mpfr_init(r19769); mpfr_init(r19770); mpfr_init(r19771); mpfr_init(r19772); mpfr_init(r19773); mpfr_init(r19774); mpfr_init(r19775); mpfr_init(r19776); mpfr_init(r19777); } double f_fm(double x) { mpfr_set_d(r19673, x, MPFR_RNDN); ; mpfr_set_si(r19675, mpfr_cmp(r19673, r19674) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19678, r19673, r19673, MPFR_RNDN); mpfr_div(r19679, r19677, r19678, MPFR_RNDN); mpfr_add(r19680, r19676, r19679, MPFR_RNDN); ; mpfr_mul(r19682, r19681, r19676, MPFR_RNDN); mpfr_mul(r19683, r19673, r19673, MPFR_RNDN); mpfr_mul(r19683, r19683, r19673, MPFR_RNDN); mpfr_sqr(r19684, r19683, MPFR_RNDN); mpfr_div(r19685, r19682, r19684, MPFR_RNDN); mpfr_add(r19686, r19680, r19685, MPFR_RNDN); ; mpfr_mul(r19688, r19687, r19676, MPFR_RNDN); mpfr_div(r19689, r19688, r19683, MPFR_RNDN); mpfr_div(r19690, r19689, r19673, MPFR_RNDN); ; ; mpfr_pow(r19693, r19673, r19692, MPFR_RNDN); mpfr_div(r19694, r19691, r19693, MPFR_RNDN); mpfr_add(r19695, r19690, r19694, MPFR_RNDN); ; mpfr_mul(r19697, r19676, r19696, MPFR_RNDN); mpfr_div(r19698, r19697, r19684, MPFR_RNDN); mpfr_div(r19699, r19698, r19678, MPFR_RNDN); mpfr_add(r19700, r19695, r19699, MPFR_RNDN); mpfr_add(r19701, r19686, r19700, MPFR_RNDN); ; mpfr_mul(r19703, r19676, r19702, MPFR_RNDN); mpfr_div(r19704, r19703, r19684, MPFR_RNDN); mpfr_add(r19705, r19676, r19704, MPFR_RNDN); ; mpfr_div(r19707, r19676, r19673, MPFR_RNDN); mpfr_div(r19708, r19707, r19683, MPFR_RNDN); mpfr_mul(r19709, r19706, r19708, MPFR_RNDN); mpfr_add(r19710, r19705, r19709, MPFR_RNDN); ; mpfr_div(r19712, r19711, r19678, MPFR_RNDN); mpfr_div(r19713, r19676, r19683, MPFR_RNDN); mpfr_mul(r19714, r19713, r19713, MPFR_RNDN); ; mpfr_pow(r19716, r19673, r19715, MPFR_RNDN); ; mpfr_div(r19718, r19716, r19717, MPFR_RNDN); mpfr_div(r19719, r19714, r19718, MPFR_RNDN); mpfr_add(r19720, r19712, r19719, MPFR_RNDN); mpfr_add(r19721, r19710, r19720, MPFR_RNDN); ; mpfr_mul(r19723, r19722, r19676, MPFR_RNDN); mpfr_div(r19724, r19723, r19684, MPFR_RNDN); mpfr_div(r19725, r19724, r19678, MPFR_RNDN); ; mpfr_mul(r19727, r19726, r19676, MPFR_RNDN); mpfr_div(r19728, r19727, r19684, MPFR_RNDN); ; mpfr_pow(r19730, r19673, r19729, MPFR_RNDN); mpfr_div(r19731, r19728, r19730, MPFR_RNDN); mpfr_add(r19732, r19725, r19731, MPFR_RNDN); mpfr_add(r19733, r19721, r19732, MPFR_RNDN); mpfr_mul(r19734, r19733, r19673, MPFR_RNDN); mpfr_div(r19735, r19701, r19734, MPFR_RNDN); ; mpfr_set_si(r19737, mpfr_cmp(r19673, r19736) <= 0, MPFR_RNDN); mpfr_sqr(r19738, r19673, MPFR_RNDN); mpfr_mul(r19739, r19683, r19738, MPFR_RNDN); mpfr_mul(r19740, r19739, r19739, MPFR_RNDN); ; mpfr_mul(r19742, r19741, r19691, MPFR_RNDN); mpfr_mul(r19743, r19738, r19742, MPFR_RNDN); mpfr_add(r19744, r19726, r19743, MPFR_RNDN); mpfr_mul(r19745, r19740, r19744, MPFR_RNDN); mpfr_mul(r19746, r19673, r19706, MPFR_RNDN); mpfr_mul(r19747, r19746, r19683, MPFR_RNDN); mpfr_mul(r19748, r19711, r19673, MPFR_RNDN); mpfr_mul(r19749, r19673, r19748, MPFR_RNDN); mpfr_add(r19750, r19676, r19749, MPFR_RNDN); mpfr_add(r19751, r19747, r19750, MPFR_RNDN); mpfr_add(r19752, r19745, r19751, MPFR_RNDN); mpfr_mul(r19753, r19738, r19722, MPFR_RNDN); mpfr_add(r19754, r19753, r19702, MPFR_RNDN); mpfr_mul(r19755, r19738, r19738, MPFR_RNDN); mpfr_mul(r19755, r19755, r19738, MPFR_RNDN); mpfr_mul(r19756, r19754, r19755, MPFR_RNDN); mpfr_add(r19757, r19752, r19756, MPFR_RNDN); mpfr_div(r19758, r19673, r19757, MPFR_RNDN); mpfr_mul(r19759, r19673, r19687, MPFR_RNDN); mpfr_mul(r19760, r19759, r19683, MPFR_RNDN); mpfr_mul(r19761, r19678, r19677, MPFR_RNDN); mpfr_add(r19762, r19761, r19676, MPFR_RNDN); mpfr_add(r19763, r19760, r19762, MPFR_RNDN); mpfr_mul(r19764, r19678, r19691, MPFR_RNDN); mpfr_mul(r19765, r19678, r19678, MPFR_RNDN); mpfr_mul(r19766, r19765, r19765, MPFR_RNDN); mpfr_mul(r19767, r19764, r19766, MPFR_RNDN); mpfr_add(r19768, r19763, r19767, MPFR_RNDN); mpfr_mul(r19769, r19696, r19673, MPFR_RNDN); mpfr_mul(r19770, r19769, r19673, MPFR_RNDN); mpfr_add(r19771, r19681, r19770, MPFR_RNDN); mpfr_mul(r19772, r19765, r19771, MPFR_RNDN); mpfr_mul(r19773, r19678, r19772, MPFR_RNDN); mpfr_add(r19774, r19768, r19773, MPFR_RNDN); mpfr_mul(r19775, r19758, r19774, MPFR_RNDN); if (mpfr_get_si(r19737, MPFR_RNDN)) { mpfr_set(r19776, r19775, MPFR_RNDN); } else { mpfr_set(r19776, r19735, MPFR_RNDN); }; if (mpfr_get_si(r19675, MPFR_RNDN)) { mpfr_set(r19777, r19735, MPFR_RNDN); } else { mpfr_set(r19777, r19776, MPFR_RNDN); }; return mpfr_get_d(r19777, MPFR_RNDN); } static mpfr_t r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803, r19804, r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824, r19825, r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843, r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861, r19862, r19863, r19864, r19865, r19866, r19867, r19868, r19869, r19870, r19871, r19872, r19873, r19874, r19875, r19876, r19877, r19878, r19879, r19880, r19881, r19882; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19778); mpfr_init_set_str(r19779, "-2.0717181681580965e+22", 10, MPFR_RNDN); mpfr_init(r19780); mpfr_init_set_str(r19781, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19782, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19783); mpfr_init(r19784); mpfr_init(r19785); mpfr_init_set_str(r19786, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19787); mpfr_init(r19788); mpfr_init(r19789); mpfr_init(r19790); mpfr_init(r19791); mpfr_init_set_str(r19792, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19793); mpfr_init(r19794); mpfr_init(r19795); mpfr_init_set_str(r19796, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19797, "10", 10, MPFR_RNDN); mpfr_init(r19798); mpfr_init(r19799); mpfr_init(r19800); mpfr_init_set_str(r19801, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19802); mpfr_init(r19803); mpfr_init(r19804); mpfr_init(r19805); mpfr_init(r19806); mpfr_init_set_str(r19807, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19808); mpfr_init(r19809); mpfr_init(r19810); mpfr_init_set_str(r19811, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19812); mpfr_init(r19813); mpfr_init(r19814); mpfr_init(r19815); mpfr_init_set_str(r19816, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19817); mpfr_init(r19818); mpfr_init(r19819); mpfr_init_set_str(r19820, "6", 10, MPFR_RNDN); mpfr_init(r19821); mpfr_init_set_str(r19822, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19823); mpfr_init(r19824); mpfr_init(r19825); mpfr_init(r19826); mpfr_init_set_str(r19827, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19828); mpfr_init(r19829); mpfr_init(r19830); mpfr_init_set_str(r19831, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19832); mpfr_init(r19833); mpfr_init_set_str(r19834, "4", 10, MPFR_RNDN); mpfr_init(r19835); mpfr_init(r19836); mpfr_init(r19837); mpfr_init(r19838); mpfr_init(r19839); mpfr_init(r19840); mpfr_init_set_str(r19841, "1.8438098325834933e+28", 10, MPFR_RNDN); mpfr_init(r19842); mpfr_init(r19843); mpfr_init(r19844); mpfr_init(r19845); mpfr_init_set_str(r19846, "2", 10, MPFR_RNDN); mpfr_init(r19847); mpfr_init(r19848); mpfr_init(r19849); mpfr_init(r19850); mpfr_init(r19851); mpfr_init(r19852); mpfr_init(r19853); mpfr_init(r19854); mpfr_init(r19855); mpfr_init(r19856); mpfr_init(r19857); mpfr_init(r19858); mpfr_init(r19859); mpfr_init(r19860); mpfr_init(r19861); mpfr_init(r19862); mpfr_init(r19863); mpfr_init(r19864); mpfr_init(r19865); mpfr_init(r19866); mpfr_init(r19867); mpfr_init(r19868); mpfr_init(r19869); mpfr_init(r19870); mpfr_init(r19871); mpfr_init(r19872); mpfr_init(r19873); mpfr_init(r19874); mpfr_init(r19875); mpfr_init(r19876); mpfr_init(r19877); mpfr_init(r19878); mpfr_init(r19879); mpfr_init(r19880); mpfr_init(r19881); mpfr_init(r19882); } double f_dm(double x) { mpfr_set_d(r19778, x, MPFR_RNDN); ; mpfr_set_si(r19780, mpfr_cmp(r19778, r19779) <= 0, MPFR_RNDN); ; ; mpfr_mul(r19783, r19778, r19778, MPFR_RNDN); mpfr_div(r19784, r19782, r19783, MPFR_RNDN); mpfr_add(r19785, r19781, r19784, MPFR_RNDN); ; mpfr_mul(r19787, r19786, r19781, MPFR_RNDN); mpfr_mul(r19788, r19778, r19778, MPFR_RNDN); mpfr_mul(r19788, r19788, r19778, MPFR_RNDN); mpfr_sqr(r19789, r19788, MPFR_RNDN); mpfr_div(r19790, r19787, r19789, MPFR_RNDN); mpfr_add(r19791, r19785, r19790, MPFR_RNDN); ; mpfr_mul(r19793, r19792, r19781, MPFR_RNDN); mpfr_div(r19794, r19793, r19788, MPFR_RNDN); mpfr_div(r19795, r19794, r19778, MPFR_RNDN); ; ; mpfr_pow(r19798, r19778, r19797, MPFR_RNDN); mpfr_div(r19799, r19796, r19798, MPFR_RNDN); mpfr_add(r19800, r19795, r19799, MPFR_RNDN); ; mpfr_mul(r19802, r19781, r19801, MPFR_RNDN); mpfr_div(r19803, r19802, r19789, MPFR_RNDN); mpfr_div(r19804, r19803, r19783, MPFR_RNDN); mpfr_add(r19805, r19800, r19804, MPFR_RNDN); mpfr_add(r19806, r19791, r19805, MPFR_RNDN); ; mpfr_mul(r19808, r19781, r19807, MPFR_RNDN); mpfr_div(r19809, r19808, r19789, MPFR_RNDN); mpfr_add(r19810, r19781, r19809, MPFR_RNDN); ; mpfr_div(r19812, r19781, r19778, MPFR_RNDN); mpfr_div(r19813, r19812, r19788, MPFR_RNDN); mpfr_mul(r19814, r19811, r19813, MPFR_RNDN); mpfr_add(r19815, r19810, r19814, MPFR_RNDN); ; mpfr_div(r19817, r19816, r19783, MPFR_RNDN); mpfr_div(r19818, r19781, r19788, MPFR_RNDN); mpfr_mul(r19819, r19818, r19818, MPFR_RNDN); ; mpfr_pow(r19821, r19778, r19820, MPFR_RNDN); ; mpfr_div(r19823, r19821, r19822, MPFR_RNDN); mpfr_div(r19824, r19819, r19823, MPFR_RNDN); mpfr_add(r19825, r19817, r19824, MPFR_RNDN); mpfr_add(r19826, r19815, r19825, MPFR_RNDN); ; mpfr_mul(r19828, r19827, r19781, MPFR_RNDN); mpfr_div(r19829, r19828, r19789, MPFR_RNDN); mpfr_div(r19830, r19829, r19783, MPFR_RNDN); ; mpfr_mul(r19832, r19831, r19781, MPFR_RNDN); mpfr_div(r19833, r19832, r19789, MPFR_RNDN); ; mpfr_pow(r19835, r19778, r19834, MPFR_RNDN); mpfr_div(r19836, r19833, r19835, MPFR_RNDN); mpfr_add(r19837, r19830, r19836, MPFR_RNDN); mpfr_add(r19838, r19826, r19837, MPFR_RNDN); mpfr_mul(r19839, r19838, r19778, MPFR_RNDN); mpfr_div(r19840, r19806, r19839, MPFR_RNDN); ; mpfr_set_si(r19842, mpfr_cmp(r19778, r19841) <= 0, MPFR_RNDN); mpfr_sqr(r19843, r19778, MPFR_RNDN); mpfr_mul(r19844, r19788, r19843, MPFR_RNDN); mpfr_mul(r19845, r19844, r19844, MPFR_RNDN); ; mpfr_mul(r19847, r19846, r19796, MPFR_RNDN); mpfr_mul(r19848, r19843, r19847, MPFR_RNDN); mpfr_add(r19849, r19831, r19848, MPFR_RNDN); mpfr_mul(r19850, r19845, r19849, MPFR_RNDN); mpfr_mul(r19851, r19778, r19811, MPFR_RNDN); mpfr_mul(r19852, r19851, r19788, MPFR_RNDN); mpfr_mul(r19853, r19816, r19778, MPFR_RNDN); mpfr_mul(r19854, r19778, r19853, MPFR_RNDN); mpfr_add(r19855, r19781, r19854, MPFR_RNDN); mpfr_add(r19856, r19852, r19855, MPFR_RNDN); mpfr_add(r19857, r19850, r19856, MPFR_RNDN); mpfr_mul(r19858, r19843, r19827, MPFR_RNDN); mpfr_add(r19859, r19858, r19807, MPFR_RNDN); mpfr_mul(r19860, r19843, r19843, MPFR_RNDN); mpfr_mul(r19860, r19860, r19843, MPFR_RNDN); mpfr_mul(r19861, r19859, r19860, MPFR_RNDN); mpfr_add(r19862, r19857, r19861, MPFR_RNDN); mpfr_div(r19863, r19778, r19862, MPFR_RNDN); mpfr_mul(r19864, r19778, r19792, MPFR_RNDN); mpfr_mul(r19865, r19864, r19788, MPFR_RNDN); mpfr_mul(r19866, r19783, r19782, MPFR_RNDN); mpfr_add(r19867, r19866, r19781, MPFR_RNDN); mpfr_add(r19868, r19865, r19867, MPFR_RNDN); mpfr_mul(r19869, r19783, r19796, MPFR_RNDN); mpfr_mul(r19870, r19783, r19783, MPFR_RNDN); mpfr_mul(r19871, r19870, r19870, MPFR_RNDN); mpfr_mul(r19872, r19869, r19871, MPFR_RNDN); mpfr_add(r19873, r19868, r19872, MPFR_RNDN); mpfr_mul(r19874, r19801, r19778, MPFR_RNDN); mpfr_mul(r19875, r19874, r19778, MPFR_RNDN); mpfr_add(r19876, r19786, r19875, MPFR_RNDN); mpfr_mul(r19877, r19870, r19876, MPFR_RNDN); mpfr_mul(r19878, r19783, r19877, MPFR_RNDN); mpfr_add(r19879, r19873, r19878, MPFR_RNDN); mpfr_mul(r19880, r19863, r19879, MPFR_RNDN); if (mpfr_get_si(r19842, MPFR_RNDN)) { mpfr_set(r19881, r19880, MPFR_RNDN); } else { mpfr_set(r19881, r19840, MPFR_RNDN); }; if (mpfr_get_si(r19780, MPFR_RNDN)) { mpfr_set(r19882, r19840, MPFR_RNDN); } else { mpfr_set(r19882, r19881, MPFR_RNDN); }; return mpfr_get_d(r19882, MPFR_RNDN); }