#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r27655 = 1; float r27656 = 0.1049934947; float r27657 = x; float r27658 = r27657 * r27657; float r27659 = r27656 * r27658; float r27660 = r27655 + r27659; float r27661 = 0.0424060604; float r27662 = r27658 * r27658; float r27663 = r27661 * r27662; float r27664 = r27660 + r27663; float r27665 = 0.0072644182; float r27666 = r27662 * r27658; float r27667 = r27665 * r27666; float r27668 = r27664 + r27667; float r27669 = 0.0005064034; float r27670 = r27666 * r27658; float r27671 = r27669 * r27670; float r27672 = r27668 + r27671; float r27673 = 0.0001789971; float r27674 = r27670 * r27658; float r27675 = r27673 * r27674; float r27676 = r27672 + r27675; float r27677 = 0.7715471019; float r27678 = r27677 * r27658; float r27679 = r27655 + r27678; float r27680 = 0.2909738639; float r27681 = r27680 * r27662; float r27682 = r27679 + r27681; float r27683 = 0.0694555761; float r27684 = r27683 * r27666; float r27685 = r27682 + r27684; float r27686 = 0.0140005442; float r27687 = r27686 * r27670; float r27688 = r27685 + r27687; float r27689 = 0.0008327945; float r27690 = r27689 * r27674; float r27691 = r27688 + r27690; float r27692 = 2; float r27693 = r27692 * r27673; float r27694 = r27674 * r27658; float r27695 = r27693 * r27694; float r27696 = r27691 + r27695; float r27697 = r27676 / r27696; float r27698 = r27697 * r27657; return r27698; } double f_id(double x) { double r27699 = 1; double r27700 = 0.1049934947; double r27701 = x; double r27702 = r27701 * r27701; double r27703 = r27700 * r27702; double r27704 = r27699 + r27703; double r27705 = 0.0424060604; double r27706 = r27702 * r27702; double r27707 = r27705 * r27706; double r27708 = r27704 + r27707; double r27709 = 0.0072644182; double r27710 = r27706 * r27702; double r27711 = r27709 * r27710; double r27712 = r27708 + r27711; double r27713 = 0.0005064034; double r27714 = r27710 * r27702; double r27715 = r27713 * r27714; double r27716 = r27712 + r27715; double r27717 = 0.0001789971; double r27718 = r27714 * r27702; double r27719 = r27717 * r27718; double r27720 = r27716 + r27719; double r27721 = 0.7715471019; double r27722 = r27721 * r27702; double r27723 = r27699 + r27722; double r27724 = 0.2909738639; double r27725 = r27724 * r27706; double r27726 = r27723 + r27725; double r27727 = 0.0694555761; double r27728 = r27727 * r27710; double r27729 = r27726 + r27728; double r27730 = 0.0140005442; double r27731 = r27730 * r27714; double r27732 = r27729 + r27731; double r27733 = 0.0008327945; double r27734 = r27733 * r27718; double r27735 = r27732 + r27734; double r27736 = 2; double r27737 = r27736 * r27717; double r27738 = r27718 * r27702; double r27739 = r27737 * r27738; double r27740 = r27735 + r27739; double r27741 = r27720 / r27740; double r27742 = r27741 * r27701; return r27742; } double f_of(float x) { float r27743 = 0.2514179000665375; float r27744 = x; float r27745 = 4; float r27746 = pow(r27744, r27745); float r27747 = r27743 / r27746; float r27748 = 0.15298196345929327; float r27749 = 6; float r27750 = pow(r27744, r27749); float r27751 = r27748 / r27750; float r27752 = r27747 + r27751; float r27753 = 0.5; float r27754 = r27753 / r27744; float r27755 = fma(r27752, r27744, r27754); float r27756 = -7.879348480694973e-07; bool r27757 = r27755 <= r27756; float r27758 = r27744 * r27744; float r27759 = r27758 * r27758; float r27760 = r27759 * r27759; float r27761 = 0.0001789971; float r27762 = r27758 * r27761; float r27763 = r27760 * r27762; float r27764 = 0.0005064034; float r27765 = 0.0072644182; float r27766 = r27765 * r27758; float r27767 = r27766 * r27759; float r27768 = 0.0424060604; float r27769 = 0.1049934947; float r27770 = 1; float r27771 = fma(r27758, r27769, r27770); float r27772 = fma(r27768, r27759, r27771); float r27773 = r27767 + r27772; float r27774 = fma(r27764, r27760, r27773); float r27775 = r27763 + r27774; float r27776 = 0.0008327945; float r27777 = r27776 * r27758; float r27778 = r27777 * r27760; float r27779 = 2; float r27780 = r27761 * r27779; float r27781 = r27759 * r27780; float r27782 = r27760 * r27781; float r27783 = r27778 + r27782; float r27784 = 0.0140005442; float r27785 = 0.2909738639; float r27786 = 0.7715471019; float r27787 = r27786 * r27744; float r27788 = fma(r27787, r27744, r27770); float r27789 = fma(r27785, r27759, r27788); float r27790 = 0.0694555761; float r27791 = r27758 * r27790; float r27792 = r27759 * r27791; float r27793 = r27789 + r27792; float r27794 = fma(r27760, r27784, r27793); float r27795 = r27783 + r27794; float r27796 = r27775 / r27795; float r27797 = r27796 * r27744; float r27798 = 1.919972694796288e-21; bool r27799 = r27755 <= r27798; float r27800 = r27769 * r27758; float r27801 = r27770 + r27800; float r27802 = r27768 * r27759; float r27803 = r27801 + r27802; float r27804 = r27759 * r27758; float r27805 = r27765 * r27804; float r27806 = r27803 + r27805; float r27807 = r27804 * r27758; float r27808 = r27764 * r27807; float r27809 = r27806 + r27808; float r27810 = r27807 * r27758; float r27811 = r27761 * r27810; float r27812 = r27809 + r27811; float r27813 = r27786 * r27758; float r27814 = r27770 + r27813; float r27815 = r27785 * r27759; float r27816 = r27814 + r27815; float r27817 = r27790 * r27804; float r27818 = r27816 + r27817; float r27819 = r27784 * r27807; float r27820 = r27818 + r27819; float r27821 = r27776 * r27810; float r27822 = r27820 + r27821; float r27823 = r27779 * r27761; float r27824 = r27810 * r27758; float r27825 = r27823 * r27824; float r27826 = r27822 + r27825; float r27827 = r27812 / r27826; float r27828 = r27827 * r27744; float r27829 = r27799 ? r27755 : r27828; float r27830 = r27757 ? r27797 : r27829; return r27830; } double f_od(double x) { double r27831 = 0.2514179000665375; double r27832 = x; double r27833 = 4; double r27834 = pow(r27832, r27833); double r27835 = r27831 / r27834; double r27836 = 0.15298196345929327; double r27837 = 6; double r27838 = pow(r27832, r27837); double r27839 = r27836 / r27838; double r27840 = r27835 + r27839; double r27841 = 0.5; double r27842 = r27841 / r27832; double r27843 = fma(r27840, r27832, r27842); double r27844 = -7.879348480694973e-07; bool r27845 = r27843 <= r27844; double r27846 = r27832 * r27832; double r27847 = r27846 * r27846; double r27848 = r27847 * r27847; double r27849 = 0.0001789971; double r27850 = r27846 * r27849; double r27851 = r27848 * r27850; double r27852 = 0.0005064034; double r27853 = 0.0072644182; double r27854 = r27853 * r27846; double r27855 = r27854 * r27847; double r27856 = 0.0424060604; double r27857 = 0.1049934947; double r27858 = 1; double r27859 = fma(r27846, r27857, r27858); double r27860 = fma(r27856, r27847, r27859); double r27861 = r27855 + r27860; double r27862 = fma(r27852, r27848, r27861); double r27863 = r27851 + r27862; double r27864 = 0.0008327945; double r27865 = r27864 * r27846; double r27866 = r27865 * r27848; double r27867 = 2; double r27868 = r27849 * r27867; double r27869 = r27847 * r27868; double r27870 = r27848 * r27869; double r27871 = r27866 + r27870; double r27872 = 0.0140005442; double r27873 = 0.2909738639; double r27874 = 0.7715471019; double r27875 = r27874 * r27832; double r27876 = fma(r27875, r27832, r27858); double r27877 = fma(r27873, r27847, r27876); double r27878 = 0.0694555761; double r27879 = r27846 * r27878; double r27880 = r27847 * r27879; double r27881 = r27877 + r27880; double r27882 = fma(r27848, r27872, r27881); double r27883 = r27871 + r27882; double r27884 = r27863 / r27883; double r27885 = r27884 * r27832; double r27886 = 1.919972694796288e-21; bool r27887 = r27843 <= r27886; double r27888 = r27857 * r27846; double r27889 = r27858 + r27888; double r27890 = r27856 * r27847; double r27891 = r27889 + r27890; double r27892 = r27847 * r27846; double r27893 = r27853 * r27892; double r27894 = r27891 + r27893; double r27895 = r27892 * r27846; double r27896 = r27852 * r27895; double r27897 = r27894 + r27896; double r27898 = r27895 * r27846; double r27899 = r27849 * r27898; double r27900 = r27897 + r27899; double r27901 = r27874 * r27846; double r27902 = r27858 + r27901; double r27903 = r27873 * r27847; double r27904 = r27902 + r27903; double r27905 = r27878 * r27892; double r27906 = r27904 + r27905; double r27907 = r27872 * r27895; double r27908 = r27906 + r27907; double r27909 = r27864 * r27898; double r27910 = r27908 + r27909; double r27911 = r27867 * r27849; double r27912 = r27898 * r27846; double r27913 = r27911 * r27912; double r27914 = r27910 + r27913; double r27915 = r27900 / r27914; double r27916 = r27915 * r27832; double r27917 = r27887 ? r27843 : r27916; double r27918 = r27845 ? r27885 : r27917; return r27918; } 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 r27919, r27920, r27921, r27922, r27923, r27924, r27925, r27926, r27927, r27928, r27929, r27930, r27931, r27932, r27933, r27934, r27935, r27936, r27937, r27938, r27939, r27940, r27941, r27942, r27943, r27944, r27945, r27946, r27947, r27948, r27949, r27950, r27951, r27952, r27953, r27954, r27955, r27956, r27957, r27958, r27959, r27960, r27961, r27962; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init_set_str(r27919, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27920, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r27921); mpfr_init(r27922); mpfr_init(r27923); mpfr_init(r27924); mpfr_init_set_str(r27925, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r27926); mpfr_init(r27927); mpfr_init(r27928); mpfr_init_set_str(r27929, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27930); mpfr_init(r27931); mpfr_init(r27932); mpfr_init_set_str(r27933, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r27934); mpfr_init(r27935); mpfr_init(r27936); mpfr_init_set_str(r27937, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27938); mpfr_init(r27939); mpfr_init(r27940); mpfr_init_set_str(r27941, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r27942); mpfr_init(r27943); mpfr_init_set_str(r27944, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r27945); mpfr_init(r27946); mpfr_init_set_str(r27947, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r27948); mpfr_init(r27949); mpfr_init_set_str(r27950, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r27951); mpfr_init(r27952); mpfr_init_set_str(r27953, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27954); mpfr_init(r27955); mpfr_init_set_str(r27956, "2", 10, MPFR_RNDN); mpfr_init(r27957); mpfr_init(r27958); mpfr_init(r27959); mpfr_init(r27960); mpfr_init(r27961); mpfr_init(r27962); } double f_im(double x) { ; ; mpfr_set_d(r27921, x, MPFR_RNDN); mpfr_mul(r27922, r27921, r27921, MPFR_RNDN); mpfr_mul(r27923, r27920, r27922, MPFR_RNDN); mpfr_add(r27924, r27919, r27923, MPFR_RNDN); ; mpfr_mul(r27926, r27922, r27922, MPFR_RNDN); mpfr_mul(r27927, r27925, r27926, MPFR_RNDN); mpfr_add(r27928, r27924, r27927, MPFR_RNDN); ; mpfr_mul(r27930, r27926, r27922, MPFR_RNDN); mpfr_mul(r27931, r27929, r27930, MPFR_RNDN); mpfr_add(r27932, r27928, r27931, MPFR_RNDN); ; mpfr_mul(r27934, r27930, r27922, MPFR_RNDN); mpfr_mul(r27935, r27933, r27934, MPFR_RNDN); mpfr_add(r27936, r27932, r27935, MPFR_RNDN); ; mpfr_mul(r27938, r27934, r27922, MPFR_RNDN); mpfr_mul(r27939, r27937, r27938, MPFR_RNDN); mpfr_add(r27940, r27936, r27939, MPFR_RNDN); ; mpfr_mul(r27942, r27941, r27922, MPFR_RNDN); mpfr_add(r27943, r27919, r27942, MPFR_RNDN); ; mpfr_mul(r27945, r27944, r27926, MPFR_RNDN); mpfr_add(r27946, r27943, r27945, MPFR_RNDN); ; mpfr_mul(r27948, r27947, r27930, MPFR_RNDN); mpfr_add(r27949, r27946, r27948, MPFR_RNDN); ; mpfr_mul(r27951, r27950, r27934, MPFR_RNDN); mpfr_add(r27952, r27949, r27951, MPFR_RNDN); ; mpfr_mul(r27954, r27953, r27938, MPFR_RNDN); mpfr_add(r27955, r27952, r27954, MPFR_RNDN); ; mpfr_mul(r27957, r27956, r27937, MPFR_RNDN); mpfr_mul(r27958, r27938, r27922, MPFR_RNDN); mpfr_mul(r27959, r27957, r27958, MPFR_RNDN); mpfr_add(r27960, r27955, r27959, MPFR_RNDN); mpfr_div(r27961, r27940, r27960, MPFR_RNDN); mpfr_mul(r27962, r27961, r27921, MPFR_RNDN); return mpfr_get_d(r27962, MPFR_RNDN); } static mpfr_t r27963, r27964, r27965, r27966, r27967, r27968, r27969, r27970, r27971, r27972, r27973, r27974, r27975, r27976, r27977, r27978, r27979, r27980, r27981, r27982, r27983, r27984, r27985, r27986, r27987, r27988, r27989, r27990, r27991, r27992, r27993, r27994, r27995, r27996, r27997, r27998, r27999, r28000, r28001, r28002, r28003, r28004, r28005, r28006, r28007, r28008, r28009, r28010, r28011, r28012, r28013, r28014, r28015, r28016, r28017, r28018, r28019, r28020, r28021, r28022, r28023, r28024, r28025, r28026, r28027, r28028, r28029, r28030, r28031, r28032, r28033, r28034, r28035, r28036, r28037, r28038, r28039, r28040, r28041, r28042, r28043, r28044, r28045, r28046, r28047, r28048, r28049, r28050; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r27963, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r27964); mpfr_init_set_str(r27965, "4", 10, MPFR_RNDN); mpfr_init(r27966); mpfr_init(r27967); mpfr_init_set_str(r27968, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r27969, "6", 10, MPFR_RNDN); mpfr_init(r27970); mpfr_init(r27971); mpfr_init(r27972); mpfr_init_set_str(r27973, "0.5", 10, MPFR_RNDN); mpfr_init(r27974); mpfr_init(r27975); mpfr_init_set_str(r27976, "-7.879348480694973e-07", 10, MPFR_RNDN); mpfr_init(r27977); mpfr_init(r27978); mpfr_init(r27979); mpfr_init(r27980); mpfr_init_set_str(r27981, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r27982); mpfr_init(r27983); mpfr_init_set_str(r27984, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r27985, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r27986); mpfr_init(r27987); mpfr_init_set_str(r27988, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r27989, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r27990, "1", 10, MPFR_RNDN); mpfr_init(r27991); mpfr_init(r27992); mpfr_init(r27993); mpfr_init(r27994); mpfr_init(r27995); mpfr_init_set_str(r27996, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r27997); mpfr_init(r27998); mpfr_init_set_str(r27999, "2", 10, MPFR_RNDN); mpfr_init(r28000); mpfr_init(r28001); mpfr_init(r28002); mpfr_init(r28003); mpfr_init_set_str(r28004, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r28005, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r28006, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28007); mpfr_init(r28008); mpfr_init(r28009); mpfr_init_set_str(r28010, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28011); mpfr_init(r28012); mpfr_init(r28013); mpfr_init(r28014); mpfr_init(r28015); mpfr_init(r28016); mpfr_init(r28017); mpfr_init_set_str(r28018, "1.919972694796288e-21", 10, MPFR_RNDN); mpfr_init(r28019); mpfr_init(r28020); mpfr_init(r28021); mpfr_init(r28022); mpfr_init(r28023); mpfr_init(r28024); mpfr_init(r28025); mpfr_init(r28026); mpfr_init(r28027); mpfr_init(r28028); mpfr_init(r28029); mpfr_init(r28030); mpfr_init(r28031); mpfr_init(r28032); mpfr_init(r28033); mpfr_init(r28034); mpfr_init(r28035); mpfr_init(r28036); mpfr_init(r28037); mpfr_init(r28038); mpfr_init(r28039); mpfr_init(r28040); mpfr_init(r28041); mpfr_init(r28042); mpfr_init(r28043); mpfr_init(r28044); mpfr_init(r28045); mpfr_init(r28046); mpfr_init(r28047); mpfr_init(r28048); mpfr_init(r28049); mpfr_init(r28050); } double f_fm(double x) { ; mpfr_set_d(r27964, x, MPFR_RNDN); ; mpfr_pow(r27966, r27964, r27965, MPFR_RNDN); mpfr_div(r27967, r27963, r27966, MPFR_RNDN); ; ; mpfr_pow(r27970, r27964, r27969, MPFR_RNDN); mpfr_div(r27971, r27968, r27970, MPFR_RNDN); mpfr_add(r27972, r27967, r27971, MPFR_RNDN); ; mpfr_div(r27974, r27973, r27964, MPFR_RNDN); mpfr_fma(r27975, r27972, r27964, r27974, MPFR_RNDN); ; mpfr_set_si(r27977, mpfr_cmp(r27975, r27976) <= 0, MPFR_RNDN); mpfr_mul(r27978, r27964, r27964, MPFR_RNDN); mpfr_mul(r27979, r27978, r27978, MPFR_RNDN); mpfr_mul(r27980, r27979, r27979, MPFR_RNDN); ; mpfr_mul(r27982, r27978, r27981, MPFR_RNDN); mpfr_mul(r27983, r27980, r27982, MPFR_RNDN); ; ; mpfr_mul(r27986, r27985, r27978, MPFR_RNDN); mpfr_mul(r27987, r27986, r27979, MPFR_RNDN); ; ; ; mpfr_fma(r27991, r27978, r27989, r27990, MPFR_RNDN); mpfr_fma(r27992, r27988, r27979, r27991, MPFR_RNDN); mpfr_add(r27993, r27987, r27992, MPFR_RNDN); mpfr_fma(r27994, r27984, r27980, r27993, MPFR_RNDN); mpfr_add(r27995, r27983, r27994, MPFR_RNDN); ; mpfr_mul(r27997, r27996, r27978, MPFR_RNDN); mpfr_mul(r27998, r27997, r27980, MPFR_RNDN); ; mpfr_mul(r28000, r27981, r27999, MPFR_RNDN); mpfr_mul(r28001, r27979, r28000, MPFR_RNDN); mpfr_mul(r28002, r27980, r28001, MPFR_RNDN); mpfr_add(r28003, r27998, r28002, MPFR_RNDN); ; ; ; mpfr_mul(r28007, r28006, r27964, MPFR_RNDN); mpfr_fma(r28008, r28007, r27964, r27990, MPFR_RNDN); mpfr_fma(r28009, r28005, r27979, r28008, MPFR_RNDN); ; mpfr_mul(r28011, r27978, r28010, MPFR_RNDN); mpfr_mul(r28012, r27979, r28011, MPFR_RNDN); mpfr_add(r28013, r28009, r28012, MPFR_RNDN); mpfr_fma(r28014, r27980, r28004, r28013, MPFR_RNDN); mpfr_add(r28015, r28003, r28014, MPFR_RNDN); mpfr_div(r28016, r27995, r28015, MPFR_RNDN); mpfr_mul(r28017, r28016, r27964, MPFR_RNDN); ; mpfr_set_si(r28019, mpfr_cmp(r27975, r28018) <= 0, MPFR_RNDN); mpfr_mul(r28020, r27989, r27978, MPFR_RNDN); mpfr_add(r28021, r27990, r28020, MPFR_RNDN); mpfr_mul(r28022, r27988, r27979, MPFR_RNDN); mpfr_add(r28023, r28021, r28022, MPFR_RNDN); mpfr_mul(r28024, r27979, r27978, MPFR_RNDN); mpfr_mul(r28025, r27985, r28024, MPFR_RNDN); mpfr_add(r28026, r28023, r28025, MPFR_RNDN); mpfr_mul(r28027, r28024, r27978, MPFR_RNDN); mpfr_mul(r28028, r27984, r28027, MPFR_RNDN); mpfr_add(r28029, r28026, r28028, MPFR_RNDN); mpfr_mul(r28030, r28027, r27978, MPFR_RNDN); mpfr_mul(r28031, r27981, r28030, MPFR_RNDN); mpfr_add(r28032, r28029, r28031, MPFR_RNDN); mpfr_mul(r28033, r28006, r27978, MPFR_RNDN); mpfr_add(r28034, r27990, r28033, MPFR_RNDN); mpfr_mul(r28035, r28005, r27979, MPFR_RNDN); mpfr_add(r28036, r28034, r28035, MPFR_RNDN); mpfr_mul(r28037, r28010, r28024, MPFR_RNDN); mpfr_add(r28038, r28036, r28037, MPFR_RNDN); mpfr_mul(r28039, r28004, r28027, MPFR_RNDN); mpfr_add(r28040, r28038, r28039, MPFR_RNDN); mpfr_mul(r28041, r27996, r28030, MPFR_RNDN); mpfr_add(r28042, r28040, r28041, MPFR_RNDN); mpfr_mul(r28043, r27999, r27981, MPFR_RNDN); mpfr_mul(r28044, r28030, r27978, MPFR_RNDN); mpfr_mul(r28045, r28043, r28044, MPFR_RNDN); mpfr_add(r28046, r28042, r28045, MPFR_RNDN); mpfr_div(r28047, r28032, r28046, MPFR_RNDN); mpfr_mul(r28048, r28047, r27964, MPFR_RNDN); if (mpfr_get_si(r28019, MPFR_RNDN)) { mpfr_set(r28049, r27975, MPFR_RNDN); } else { mpfr_set(r28049, r28048, MPFR_RNDN); }; if (mpfr_get_si(r27977, MPFR_RNDN)) { mpfr_set(r28050, r28017, MPFR_RNDN); } else { mpfr_set(r28050, r28049, MPFR_RNDN); }; return mpfr_get_d(r28050, MPFR_RNDN); } static mpfr_t r28051, r28052, r28053, r28054, r28055, r28056, r28057, r28058, r28059, r28060, r28061, r28062, r28063, r28064, r28065, r28066, r28067, r28068, r28069, r28070, r28071, r28072, r28073, r28074, r28075, r28076, r28077, r28078, r28079, r28080, r28081, r28082, r28083, r28084, r28085, r28086, r28087, r28088, r28089, r28090, r28091, r28092, r28093, r28094, r28095, r28096, r28097, r28098, r28099, r28100, r28101, r28102, r28103, r28104, r28105, r28106, r28107, r28108, r28109, r28110, r28111, r28112, r28113, r28114, r28115, r28116, r28117, r28118, r28119, r28120, r28121, r28122, r28123, r28124, r28125, r28126, r28127, r28128, r28129, r28130, r28131, r28132, r28133, r28134, r28135, r28136, r28137, r28138; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r28051, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r28052); mpfr_init_set_str(r28053, "4", 10, MPFR_RNDN); mpfr_init(r28054); mpfr_init(r28055); mpfr_init_set_str(r28056, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r28057, "6", 10, MPFR_RNDN); mpfr_init(r28058); mpfr_init(r28059); mpfr_init(r28060); mpfr_init_set_str(r28061, "0.5", 10, MPFR_RNDN); mpfr_init(r28062); mpfr_init(r28063); mpfr_init_set_str(r28064, "-7.879348480694973e-07", 10, MPFR_RNDN); mpfr_init(r28065); mpfr_init(r28066); mpfr_init(r28067); mpfr_init(r28068); mpfr_init_set_str(r28069, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28070); mpfr_init(r28071); mpfr_init_set_str(r28072, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r28073, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r28074); mpfr_init(r28075); mpfr_init_set_str(r28076, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r28077, "0.1049934947", 10, MPFR_RNDN); mpfr_init_set_str(r28078, "1", 10, MPFR_RNDN); mpfr_init(r28079); mpfr_init(r28080); mpfr_init(r28081); mpfr_init(r28082); mpfr_init(r28083); mpfr_init_set_str(r28084, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28085); mpfr_init(r28086); mpfr_init_set_str(r28087, "2", 10, MPFR_RNDN); mpfr_init(r28088); mpfr_init(r28089); mpfr_init(r28090); mpfr_init(r28091); mpfr_init_set_str(r28092, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r28093, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r28094, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28095); mpfr_init(r28096); mpfr_init(r28097); mpfr_init_set_str(r28098, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28099); mpfr_init(r28100); mpfr_init(r28101); mpfr_init(r28102); mpfr_init(r28103); mpfr_init(r28104); mpfr_init(r28105); mpfr_init_set_str(r28106, "1.919972694796288e-21", 10, MPFR_RNDN); mpfr_init(r28107); mpfr_init(r28108); mpfr_init(r28109); mpfr_init(r28110); mpfr_init(r28111); mpfr_init(r28112); mpfr_init(r28113); mpfr_init(r28114); mpfr_init(r28115); mpfr_init(r28116); mpfr_init(r28117); mpfr_init(r28118); mpfr_init(r28119); mpfr_init(r28120); mpfr_init(r28121); mpfr_init(r28122); mpfr_init(r28123); mpfr_init(r28124); mpfr_init(r28125); mpfr_init(r28126); mpfr_init(r28127); mpfr_init(r28128); mpfr_init(r28129); mpfr_init(r28130); mpfr_init(r28131); mpfr_init(r28132); mpfr_init(r28133); mpfr_init(r28134); mpfr_init(r28135); mpfr_init(r28136); mpfr_init(r28137); mpfr_init(r28138); } double f_dm(double x) { ; mpfr_set_d(r28052, x, MPFR_RNDN); ; mpfr_pow(r28054, r28052, r28053, MPFR_RNDN); mpfr_div(r28055, r28051, r28054, MPFR_RNDN); ; ; mpfr_pow(r28058, r28052, r28057, MPFR_RNDN); mpfr_div(r28059, r28056, r28058, MPFR_RNDN); mpfr_add(r28060, r28055, r28059, MPFR_RNDN); ; mpfr_div(r28062, r28061, r28052, MPFR_RNDN); mpfr_fma(r28063, r28060, r28052, r28062, MPFR_RNDN); ; mpfr_set_si(r28065, mpfr_cmp(r28063, r28064) <= 0, MPFR_RNDN); mpfr_mul(r28066, r28052, r28052, MPFR_RNDN); mpfr_mul(r28067, r28066, r28066, MPFR_RNDN); mpfr_mul(r28068, r28067, r28067, MPFR_RNDN); ; mpfr_mul(r28070, r28066, r28069, MPFR_RNDN); mpfr_mul(r28071, r28068, r28070, MPFR_RNDN); ; ; mpfr_mul(r28074, r28073, r28066, MPFR_RNDN); mpfr_mul(r28075, r28074, r28067, MPFR_RNDN); ; ; ; mpfr_fma(r28079, r28066, r28077, r28078, MPFR_RNDN); mpfr_fma(r28080, r28076, r28067, r28079, MPFR_RNDN); mpfr_add(r28081, r28075, r28080, MPFR_RNDN); mpfr_fma(r28082, r28072, r28068, r28081, MPFR_RNDN); mpfr_add(r28083, r28071, r28082, MPFR_RNDN); ; mpfr_mul(r28085, r28084, r28066, MPFR_RNDN); mpfr_mul(r28086, r28085, r28068, MPFR_RNDN); ; mpfr_mul(r28088, r28069, r28087, MPFR_RNDN); mpfr_mul(r28089, r28067, r28088, MPFR_RNDN); mpfr_mul(r28090, r28068, r28089, MPFR_RNDN); mpfr_add(r28091, r28086, r28090, MPFR_RNDN); ; ; ; mpfr_mul(r28095, r28094, r28052, MPFR_RNDN); mpfr_fma(r28096, r28095, r28052, r28078, MPFR_RNDN); mpfr_fma(r28097, r28093, r28067, r28096, MPFR_RNDN); ; mpfr_mul(r28099, r28066, r28098, MPFR_RNDN); mpfr_mul(r28100, r28067, r28099, MPFR_RNDN); mpfr_add(r28101, r28097, r28100, MPFR_RNDN); mpfr_fma(r28102, r28068, r28092, r28101, MPFR_RNDN); mpfr_add(r28103, r28091, r28102, MPFR_RNDN); mpfr_div(r28104, r28083, r28103, MPFR_RNDN); mpfr_mul(r28105, r28104, r28052, MPFR_RNDN); ; mpfr_set_si(r28107, mpfr_cmp(r28063, r28106) <= 0, MPFR_RNDN); mpfr_mul(r28108, r28077, r28066, MPFR_RNDN); mpfr_add(r28109, r28078, r28108, MPFR_RNDN); mpfr_mul(r28110, r28076, r28067, MPFR_RNDN); mpfr_add(r28111, r28109, r28110, MPFR_RNDN); mpfr_mul(r28112, r28067, r28066, MPFR_RNDN); mpfr_mul(r28113, r28073, r28112, MPFR_RNDN); mpfr_add(r28114, r28111, r28113, MPFR_RNDN); mpfr_mul(r28115, r28112, r28066, MPFR_RNDN); mpfr_mul(r28116, r28072, r28115, MPFR_RNDN); mpfr_add(r28117, r28114, r28116, MPFR_RNDN); mpfr_mul(r28118, r28115, r28066, MPFR_RNDN); mpfr_mul(r28119, r28069, r28118, MPFR_RNDN); mpfr_add(r28120, r28117, r28119, MPFR_RNDN); mpfr_mul(r28121, r28094, r28066, MPFR_RNDN); mpfr_add(r28122, r28078, r28121, MPFR_RNDN); mpfr_mul(r28123, r28093, r28067, MPFR_RNDN); mpfr_add(r28124, r28122, r28123, MPFR_RNDN); mpfr_mul(r28125, r28098, r28112, MPFR_RNDN); mpfr_add(r28126, r28124, r28125, MPFR_RNDN); mpfr_mul(r28127, r28092, r28115, MPFR_RNDN); mpfr_add(r28128, r28126, r28127, MPFR_RNDN); mpfr_mul(r28129, r28084, r28118, MPFR_RNDN); mpfr_add(r28130, r28128, r28129, MPFR_RNDN); mpfr_mul(r28131, r28087, r28069, MPFR_RNDN); mpfr_mul(r28132, r28118, r28066, MPFR_RNDN); mpfr_mul(r28133, r28131, r28132, MPFR_RNDN); mpfr_add(r28134, r28130, r28133, MPFR_RNDN); mpfr_div(r28135, r28120, r28134, MPFR_RNDN); mpfr_mul(r28136, r28135, r28052, MPFR_RNDN); if (mpfr_get_si(r28107, MPFR_RNDN)) { mpfr_set(r28137, r28063, MPFR_RNDN); } else { mpfr_set(r28137, r28136, MPFR_RNDN); }; if (mpfr_get_si(r28065, MPFR_RNDN)) { mpfr_set(r28138, r28105, MPFR_RNDN); } else { mpfr_set(r28138, r28137, MPFR_RNDN); }; return mpfr_get_d(r28138, MPFR_RNDN); }