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