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