#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25666 = 1; float r25667 = 0.1049934947; float r25668 = x; float r25669 = r25668 * r25668; float r25670 = r25667 * r25669; float r25671 = r25666 + r25670; float r25672 = 0.0424060604; float r25673 = r25669 * r25669; float r25674 = r25672 * r25673; float r25675 = r25671 + r25674; float r25676 = 0.0072644182; float r25677 = r25673 * r25669; float r25678 = r25676 * r25677; float r25679 = r25675 + r25678; float r25680 = 0.0005064034; float r25681 = r25677 * r25669; float r25682 = r25680 * r25681; float r25683 = r25679 + r25682; float r25684 = 0.0001789971; float r25685 = r25681 * r25669; float r25686 = r25684 * r25685; float r25687 = r25683 + r25686; float r25688 = 0.7715471019; float r25689 = r25688 * r25669; float r25690 = r25666 + r25689; float r25691 = 0.2909738639; float r25692 = r25691 * r25673; float r25693 = r25690 + r25692; float r25694 = 0.0694555761; float r25695 = r25694 * r25677; float r25696 = r25693 + r25695; float r25697 = 0.0140005442; float r25698 = r25697 * r25681; float r25699 = r25696 + r25698; float r25700 = 0.0008327945; float r25701 = r25700 * r25685; float r25702 = r25699 + r25701; float r25703 = 2; float r25704 = r25703 * r25684; float r25705 = r25685 * r25669; float r25706 = r25704 * r25705; float r25707 = r25702 + r25706; float r25708 = r25687 / r25707; float r25709 = r25708 * r25668; return r25709; } double f_id(double x) { double r25710 = 1; double r25711 = 0.1049934947; double r25712 = x; double r25713 = r25712 * r25712; double r25714 = r25711 * r25713; double r25715 = r25710 + r25714; double r25716 = 0.0424060604; double r25717 = r25713 * r25713; double r25718 = r25716 * r25717; double r25719 = r25715 + r25718; double r25720 = 0.0072644182; double r25721 = r25717 * r25713; double r25722 = r25720 * r25721; double r25723 = r25719 + r25722; double r25724 = 0.0005064034; double r25725 = r25721 * r25713; double r25726 = r25724 * r25725; double r25727 = r25723 + r25726; double r25728 = 0.0001789971; double r25729 = r25725 * r25713; double r25730 = r25728 * r25729; double r25731 = r25727 + r25730; double r25732 = 0.7715471019; double r25733 = r25732 * r25713; double r25734 = r25710 + r25733; double r25735 = 0.2909738639; double r25736 = r25735 * r25717; double r25737 = r25734 + r25736; double r25738 = 0.0694555761; double r25739 = r25738 * r25721; double r25740 = r25737 + r25739; double r25741 = 0.0140005442; double r25742 = r25741 * r25725; double r25743 = r25740 + r25742; double r25744 = 0.0008327945; double r25745 = r25744 * r25729; double r25746 = r25743 + r25745; double r25747 = 2; double r25748 = r25747 * r25728; double r25749 = r25729 * r25713; double r25750 = r25748 * r25749; double r25751 = r25746 + r25750; double r25752 = r25731 / r25751; double r25753 = r25752 * r25712; return r25753; } double f_of(float x) { float r25754 = 0.5; float r25755 = x; float r25756 = r25754 / r25755; float r25757 = 0.2514179000665375; float r25758 = r25757 * r25755; float r25759 = 4; float r25760 = pow(r25755, r25759); float r25761 = r25758 / r25760; float r25762 = 0.15298196345929327; float r25763 = r25762 * r25755; float r25764 = 6; float r25765 = pow(r25755, r25764); float r25766 = r25763 / r25765; float r25767 = r25761 + r25766; float r25768 = r25756 + r25767; float r25769 = -0.00014612633325219868; bool r25770 = r25768 <= r25769; float r25771 = r25755 * r25755; float r25772 = 0.1049934947; float r25773 = r25771 * r25772; float r25774 = 1; float r25775 = r25773 + r25774; float r25776 = r25771 * r25771; float r25777 = 0.0424060604; float r25778 = 0.0072644182; float r25779 = r25755 * r25778; float r25780 = r25779 * r25755; float r25781 = r25777 + r25780; float r25782 = r25776 * r25781; float r25783 = 3; float r25784 = r25783 + r25774; float r25785 = pow(r25771, r25784); float r25786 = 0.0001789971; float r25787 = r25771 * r25786; float r25788 = 0.0005064034; float r25789 = r25787 + r25788; float r25790 = r25785 * r25789; float r25791 = r25782 + r25790; float r25792 = r25775 + r25791; float r25793 = 0.2909738639; float r25794 = r25793 * r25755; float r25795 = pow(r25755, r25783); float r25796 = r25794 * r25795; float r25797 = r25774 + r25796; float r25798 = 0.7715471019; float r25799 = r25798 * r25755; float r25800 = r25799 * r25755; float r25801 = r25797 + r25800; float r25802 = r25795 * r25771; float r25803 = r25802 * r25802; float r25804 = 0.0008327945; float r25805 = 2; float r25806 = r25786 * r25805; float r25807 = r25806 * r25771; float r25808 = r25804 + r25807; float r25809 = r25803 * r25808; float r25810 = r25795 * r25795; float r25811 = 0.0140005442; float r25812 = r25755 * r25811; float r25813 = r25812 * r25755; float r25814 = 0.0694555761; float r25815 = r25813 + r25814; float r25816 = r25810 * r25815; float r25817 = r25809 + r25816; float r25818 = r25801 + r25817; float r25819 = r25792 / r25818; float r25820 = sqrt(r25819); float r25821 = r25820 * r25820; float r25822 = r25821 * r25755; float r25823 = 6.050122732948135e-08; bool r25824 = r25768 <= r25823; float r25825 = r25788 * r25771; float r25826 = pow(r25771, r25783); float r25827 = r25825 * r25826; float r25828 = r25776 * r25776; float r25829 = r25828 * r25787; float r25830 = r25827 + r25829; float r25831 = r25778 * r25771; float r25832 = r25831 * r25776; float r25833 = r25777 * r25755; float r25834 = r25833 * r25795; float r25835 = r25773 + r25834; float r25836 = r25774 + r25835; float r25837 = r25832 + r25836; float r25838 = r25830 + r25837; float r25839 = pow(r25838, r25783); float r25840 = r25771 * r25811; float r25841 = r25826 * r25840; float r25842 = r25800 + r25774; float r25843 = r25841 + r25842; float r25844 = r25771 * r25793; float r25845 = r25814 * r25755; float r25846 = r25845 * r25795; float r25847 = r25844 + r25846; float r25848 = r25771 * r25847; float r25849 = r25843 + r25848; float r25850 = r25805 * r25786; float r25851 = r25776 * r25850; float r25852 = r25828 * r25851; float r25853 = r25804 * r25771; float r25854 = r25853 * r25828; float r25855 = r25852 + r25854; float r25856 = r25849 + r25855; float r25857 = pow(r25856, r25783); float r25858 = r25839 / r25857; float r25859 = cbrt(r25858); float r25860 = r25859 * r25755; float r25861 = r25824 ? r25768 : r25860; float r25862 = r25770 ? r25822 : r25861; return r25862; } double f_od(double x) { double r25863 = 0.5; double r25864 = x; double r25865 = r25863 / r25864; double r25866 = 0.2514179000665375; double r25867 = r25866 * r25864; double r25868 = 4; double r25869 = pow(r25864, r25868); double r25870 = r25867 / r25869; double r25871 = 0.15298196345929327; double r25872 = r25871 * r25864; double r25873 = 6; double r25874 = pow(r25864, r25873); double r25875 = r25872 / r25874; double r25876 = r25870 + r25875; double r25877 = r25865 + r25876; double r25878 = -0.00014612633325219868; bool r25879 = r25877 <= r25878; double r25880 = r25864 * r25864; double r25881 = 0.1049934947; double r25882 = r25880 * r25881; double r25883 = 1; double r25884 = r25882 + r25883; double r25885 = r25880 * r25880; double r25886 = 0.0424060604; double r25887 = 0.0072644182; double r25888 = r25864 * r25887; double r25889 = r25888 * r25864; double r25890 = r25886 + r25889; double r25891 = r25885 * r25890; double r25892 = 3; double r25893 = r25892 + r25883; double r25894 = pow(r25880, r25893); double r25895 = 0.0001789971; double r25896 = r25880 * r25895; double r25897 = 0.0005064034; double r25898 = r25896 + r25897; double r25899 = r25894 * r25898; double r25900 = r25891 + r25899; double r25901 = r25884 + r25900; double r25902 = 0.2909738639; double r25903 = r25902 * r25864; double r25904 = pow(r25864, r25892); double r25905 = r25903 * r25904; double r25906 = r25883 + r25905; double r25907 = 0.7715471019; double r25908 = r25907 * r25864; double r25909 = r25908 * r25864; double r25910 = r25906 + r25909; double r25911 = r25904 * r25880; double r25912 = r25911 * r25911; double r25913 = 0.0008327945; double r25914 = 2; double r25915 = r25895 * r25914; double r25916 = r25915 * r25880; double r25917 = r25913 + r25916; double r25918 = r25912 * r25917; double r25919 = r25904 * r25904; double r25920 = 0.0140005442; double r25921 = r25864 * r25920; double r25922 = r25921 * r25864; double r25923 = 0.0694555761; double r25924 = r25922 + r25923; double r25925 = r25919 * r25924; double r25926 = r25918 + r25925; double r25927 = r25910 + r25926; double r25928 = r25901 / r25927; double r25929 = sqrt(r25928); double r25930 = r25929 * r25929; double r25931 = r25930 * r25864; double r25932 = 6.050122732948135e-08; bool r25933 = r25877 <= r25932; double r25934 = r25897 * r25880; double r25935 = pow(r25880, r25892); double r25936 = r25934 * r25935; double r25937 = r25885 * r25885; double r25938 = r25937 * r25896; double r25939 = r25936 + r25938; double r25940 = r25887 * r25880; double r25941 = r25940 * r25885; double r25942 = r25886 * r25864; double r25943 = r25942 * r25904; double r25944 = r25882 + r25943; double r25945 = r25883 + r25944; double r25946 = r25941 + r25945; double r25947 = r25939 + r25946; double r25948 = pow(r25947, r25892); double r25949 = r25880 * r25920; double r25950 = r25935 * r25949; double r25951 = r25909 + r25883; double r25952 = r25950 + r25951; double r25953 = r25880 * r25902; double r25954 = r25923 * r25864; double r25955 = r25954 * r25904; double r25956 = r25953 + r25955; double r25957 = r25880 * r25956; double r25958 = r25952 + r25957; double r25959 = r25914 * r25895; double r25960 = r25885 * r25959; double r25961 = r25937 * r25960; double r25962 = r25913 * r25880; double r25963 = r25962 * r25937; double r25964 = r25961 + r25963; double r25965 = r25958 + r25964; double r25966 = pow(r25965, r25892); double r25967 = r25948 / r25966; double r25968 = cbrt(r25967); double r25969 = r25968 * r25864; double r25970 = r25933 ? r25877 : r25969; double r25971 = r25879 ? r25931 : r25970; return r25971; } 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 r25972, r25973, r25974, r25975, r25976, r25977, r25978, r25979, r25980, r25981, r25982, r25983, r25984, r25985, r25986, r25987, r25988, r25989, r25990, r25991, r25992, r25993, r25994, r25995, r25996, r25997, r25998, r25999, r26000, r26001, r26002, r26003, r26004, r26005, r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r25972, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25973, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25974); mpfr_init(r25975); mpfr_init(r25976); mpfr_init(r25977); mpfr_init_set_str(r25978, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25979); mpfr_init(r25980); mpfr_init(r25981); mpfr_init_set_str(r25982, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25983); mpfr_init(r25984); mpfr_init(r25985); mpfr_init_set_str(r25986, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25987); mpfr_init(r25988); mpfr_init(r25989); mpfr_init_set_str(r25990, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25991); mpfr_init(r25992); mpfr_init(r25993); mpfr_init_set_str(r25994, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25995); mpfr_init(r25996); mpfr_init_set_str(r25997, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25998); mpfr_init(r25999); mpfr_init_set_str(r26000, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26001); mpfr_init(r26002); mpfr_init_set_str(r26003, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26004); mpfr_init(r26005); mpfr_init_set_str(r26006, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26007); mpfr_init(r26008); mpfr_init_set_str(r26009, "2", 10, MPFR_RNDN); mpfr_init(r26010); mpfr_init(r26011); mpfr_init(r26012); mpfr_init(r26013); mpfr_init(r26014); mpfr_init(r26015); } double f_im(double x) { ; ; mpfr_set_d(r25974, x, MPFR_RNDN); mpfr_mul(r25975, r25974, r25974, MPFR_RNDN); mpfr_mul(r25976, r25973, r25975, MPFR_RNDN); mpfr_add(r25977, r25972, r25976, MPFR_RNDN); ; mpfr_mul(r25979, r25975, r25975, MPFR_RNDN); mpfr_mul(r25980, r25978, r25979, MPFR_RNDN); mpfr_add(r25981, r25977, r25980, MPFR_RNDN); ; mpfr_mul(r25983, r25979, r25975, MPFR_RNDN); mpfr_mul(r25984, r25982, r25983, MPFR_RNDN); mpfr_add(r25985, r25981, r25984, MPFR_RNDN); ; mpfr_mul(r25987, r25983, r25975, MPFR_RNDN); mpfr_mul(r25988, r25986, r25987, MPFR_RNDN); mpfr_add(r25989, r25985, r25988, MPFR_RNDN); ; mpfr_mul(r25991, r25987, r25975, MPFR_RNDN); mpfr_mul(r25992, r25990, r25991, MPFR_RNDN); mpfr_add(r25993, r25989, r25992, MPFR_RNDN); ; mpfr_mul(r25995, r25994, r25975, MPFR_RNDN); mpfr_add(r25996, r25972, r25995, MPFR_RNDN); ; mpfr_mul(r25998, r25997, r25979, MPFR_RNDN); mpfr_add(r25999, r25996, r25998, MPFR_RNDN); ; mpfr_mul(r26001, r26000, r25983, MPFR_RNDN); mpfr_add(r26002, r25999, r26001, MPFR_RNDN); ; mpfr_mul(r26004, r26003, r25987, MPFR_RNDN); mpfr_add(r26005, r26002, r26004, MPFR_RNDN); ; mpfr_mul(r26007, r26006, r25991, MPFR_RNDN); mpfr_add(r26008, r26005, r26007, MPFR_RNDN); ; mpfr_mul(r26010, r26009, r25990, MPFR_RNDN); mpfr_mul(r26011, r25991, r25975, MPFR_RNDN); mpfr_mul(r26012, r26010, r26011, MPFR_RNDN); mpfr_add(r26013, r26008, r26012, MPFR_RNDN); mpfr_div(r26014, r25993, r26013, MPFR_RNDN); mpfr_mul(r26015, r26014, r25974, MPFR_RNDN); return mpfr_get_d(r26015, MPFR_RNDN); } static mpfr_t r26016, r26017, r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099, r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115, r26116, r26117, r26118, r26119, r26120, r26121, r26122, r26123, r26124; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r26016, "0.5", 10, MPFR_RNDN); mpfr_init(r26017); mpfr_init(r26018); mpfr_init_set_str(r26019, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26020); mpfr_init_set_str(r26021, "4", 10, MPFR_RNDN); mpfr_init(r26022); mpfr_init(r26023); mpfr_init_set_str(r26024, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26025); mpfr_init_set_str(r26026, "6", 10, MPFR_RNDN); mpfr_init(r26027); mpfr_init(r26028); mpfr_init(r26029); mpfr_init(r26030); mpfr_init_set_str(r26031, "-0.00014612633325219868", 10, MPFR_RNDN); mpfr_init(r26032); mpfr_init(r26033); mpfr_init_set_str(r26034, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26035); mpfr_init_set_str(r26036, "1", 10, MPFR_RNDN); mpfr_init(r26037); mpfr_init(r26038); mpfr_init_set_str(r26039, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r26040, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26041); mpfr_init(r26042); mpfr_init(r26043); mpfr_init(r26044); mpfr_init_set_str(r26045, "3", 10, MPFR_RNDN); mpfr_init(r26046); mpfr_init(r26047); mpfr_init_set_str(r26048, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26049); mpfr_init_set_str(r26050, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26051); mpfr_init(r26052); mpfr_init(r26053); mpfr_init(r26054); mpfr_init_set_str(r26055, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init_set_str(r26060, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); mpfr_init(r26065); mpfr_init_set_str(r26066, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r26067, "2", 10, MPFR_RNDN); mpfr_init(r26068); mpfr_init(r26069); mpfr_init(r26070); mpfr_init(r26071); mpfr_init(r26072); mpfr_init_set_str(r26073, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26074); mpfr_init(r26075); mpfr_init_set_str(r26076, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init(r26084); mpfr_init_set_str(r26085, "6.050122732948135e-08", 10, MPFR_RNDN); mpfr_init(r26086); mpfr_init(r26087); mpfr_init(r26088); mpfr_init(r26089); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); mpfr_init(r26093); mpfr_init(r26094); mpfr_init(r26095); mpfr_init(r26096); mpfr_init(r26097); mpfr_init(r26098); mpfr_init(r26099); mpfr_init(r26100); mpfr_init(r26101); mpfr_init(r26102); mpfr_init(r26103); mpfr_init(r26104); mpfr_init(r26105); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init(r26109); mpfr_init(r26110); mpfr_init(r26111); mpfr_init(r26112); mpfr_init(r26113); mpfr_init(r26114); mpfr_init(r26115); mpfr_init(r26116); mpfr_init(r26117); mpfr_init(r26118); mpfr_init(r26119); mpfr_init(r26120); mpfr_init(r26121); mpfr_init(r26122); mpfr_init(r26123); mpfr_init(r26124); } double f_fm(double x) { ; mpfr_set_d(r26017, x, MPFR_RNDN); mpfr_div(r26018, r26016, r26017, MPFR_RNDN); ; mpfr_mul(r26020, r26019, r26017, MPFR_RNDN); ; mpfr_pow(r26022, r26017, r26021, MPFR_RNDN); mpfr_div(r26023, r26020, r26022, MPFR_RNDN); ; mpfr_mul(r26025, r26024, r26017, MPFR_RNDN); ; mpfr_pow(r26027, r26017, r26026, MPFR_RNDN); mpfr_div(r26028, r26025, r26027, MPFR_RNDN); mpfr_add(r26029, r26023, r26028, MPFR_RNDN); mpfr_add(r26030, r26018, r26029, MPFR_RNDN); ; mpfr_set_si(r26032, mpfr_cmp(r26030, r26031) <= 0, MPFR_RNDN); mpfr_mul(r26033, r26017, r26017, MPFR_RNDN); ; mpfr_mul(r26035, r26033, r26034, MPFR_RNDN); ; mpfr_add(r26037, r26035, r26036, MPFR_RNDN); mpfr_mul(r26038, r26033, r26033, MPFR_RNDN); ; ; mpfr_mul(r26041, r26017, r26040, MPFR_RNDN); mpfr_mul(r26042, r26041, r26017, MPFR_RNDN); mpfr_add(r26043, r26039, r26042, MPFR_RNDN); mpfr_mul(r26044, r26038, r26043, MPFR_RNDN); ; mpfr_add(r26046, r26045, r26036, MPFR_RNDN); mpfr_pow(r26047, r26033, r26046, MPFR_RNDN); ; mpfr_mul(r26049, r26033, r26048, MPFR_RNDN); ; mpfr_add(r26051, r26049, r26050, MPFR_RNDN); mpfr_mul(r26052, r26047, r26051, MPFR_RNDN); mpfr_add(r26053, r26044, r26052, MPFR_RNDN); mpfr_add(r26054, r26037, r26053, MPFR_RNDN); ; mpfr_mul(r26056, r26055, r26017, MPFR_RNDN); mpfr_pow(r26057, r26017, r26045, MPFR_RNDN); mpfr_mul(r26058, r26056, r26057, MPFR_RNDN); mpfr_add(r26059, r26036, r26058, MPFR_RNDN); ; mpfr_mul(r26061, r26060, r26017, MPFR_RNDN); mpfr_mul(r26062, r26061, r26017, MPFR_RNDN); mpfr_add(r26063, r26059, r26062, MPFR_RNDN); mpfr_mul(r26064, r26057, r26033, MPFR_RNDN); mpfr_mul(r26065, r26064, r26064, MPFR_RNDN); ; ; mpfr_mul(r26068, r26048, r26067, MPFR_RNDN); mpfr_mul(r26069, r26068, r26033, MPFR_RNDN); mpfr_add(r26070, r26066, r26069, MPFR_RNDN); mpfr_mul(r26071, r26065, r26070, MPFR_RNDN); mpfr_mul(r26072, r26057, r26057, MPFR_RNDN); ; mpfr_mul(r26074, r26017, r26073, MPFR_RNDN); mpfr_mul(r26075, r26074, r26017, MPFR_RNDN); ; mpfr_add(r26077, r26075, r26076, MPFR_RNDN); mpfr_mul(r26078, r26072, r26077, MPFR_RNDN); mpfr_add(r26079, r26071, r26078, MPFR_RNDN); mpfr_add(r26080, r26063, r26079, MPFR_RNDN); mpfr_div(r26081, r26054, r26080, MPFR_RNDN); mpfr_sqrt(r26082, r26081, MPFR_RNDN); mpfr_mul(r26083, r26082, r26082, MPFR_RNDN); mpfr_mul(r26084, r26083, r26017, MPFR_RNDN); ; mpfr_set_si(r26086, mpfr_cmp(r26030, r26085) <= 0, MPFR_RNDN); mpfr_mul(r26087, r26050, r26033, MPFR_RNDN); mpfr_pow(r26088, r26033, r26045, MPFR_RNDN); mpfr_mul(r26089, r26087, r26088, MPFR_RNDN); mpfr_mul(r26090, r26038, r26038, MPFR_RNDN); mpfr_mul(r26091, r26090, r26049, MPFR_RNDN); mpfr_add(r26092, r26089, r26091, MPFR_RNDN); mpfr_mul(r26093, r26040, r26033, MPFR_RNDN); mpfr_mul(r26094, r26093, r26038, MPFR_RNDN); mpfr_mul(r26095, r26039, r26017, MPFR_RNDN); mpfr_mul(r26096, r26095, r26057, MPFR_RNDN); mpfr_add(r26097, r26035, r26096, MPFR_RNDN); mpfr_add(r26098, r26036, r26097, MPFR_RNDN); mpfr_add(r26099, r26094, r26098, MPFR_RNDN); mpfr_add(r26100, r26092, r26099, MPFR_RNDN); mpfr_pow(r26101, r26100, r26045, MPFR_RNDN); mpfr_mul(r26102, r26033, r26073, MPFR_RNDN); mpfr_mul(r26103, r26088, r26102, MPFR_RNDN); mpfr_add(r26104, r26062, r26036, MPFR_RNDN); mpfr_add(r26105, r26103, r26104, MPFR_RNDN); mpfr_mul(r26106, r26033, r26055, MPFR_RNDN); mpfr_mul(r26107, r26076, r26017, MPFR_RNDN); mpfr_mul(r26108, r26107, r26057, MPFR_RNDN); mpfr_add(r26109, r26106, r26108, MPFR_RNDN); mpfr_mul(r26110, r26033, r26109, MPFR_RNDN); mpfr_add(r26111, r26105, r26110, MPFR_RNDN); mpfr_mul(r26112, r26067, r26048, MPFR_RNDN); mpfr_mul(r26113, r26038, r26112, MPFR_RNDN); mpfr_mul(r26114, r26090, r26113, MPFR_RNDN); mpfr_mul(r26115, r26066, r26033, MPFR_RNDN); mpfr_mul(r26116, r26115, r26090, MPFR_RNDN); mpfr_add(r26117, r26114, r26116, MPFR_RNDN); mpfr_add(r26118, r26111, r26117, MPFR_RNDN); mpfr_pow(r26119, r26118, r26045, MPFR_RNDN); mpfr_div(r26120, r26101, r26119, MPFR_RNDN); mpfr_cbrt(r26121, r26120, MPFR_RNDN); mpfr_mul(r26122, r26121, r26017, MPFR_RNDN); if (mpfr_get_si(r26086, MPFR_RNDN)) { mpfr_set(r26123, r26030, MPFR_RNDN); } else { mpfr_set(r26123, r26122, MPFR_RNDN); }; if (mpfr_get_si(r26032, MPFR_RNDN)) { mpfr_set(r26124, r26084, MPFR_RNDN); } else { mpfr_set(r26124, r26123, MPFR_RNDN); }; return mpfr_get_d(r26124, MPFR_RNDN); } static mpfr_t r26125, r26126, r26127, r26128, r26129, r26130, r26131, r26132, r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140, r26141, r26142, r26143, r26144, r26145, r26146, r26147, r26148, r26149, r26150, r26151, r26152, r26153, r26154, r26155, r26156, r26157, r26158, r26159, r26160, r26161, r26162, r26163, r26164, r26165, r26166, r26167, r26168, r26169, r26170, r26171, r26172, r26173, r26174, r26175, r26176, r26177, r26178, r26179, r26180, r26181, r26182, r26183, r26184, r26185, r26186, r26187, r26188, r26189, r26190, r26191, r26192, r26193, r26194, r26195, r26196, r26197, r26198, r26199, r26200, r26201, r26202, r26203, r26204, r26205, r26206, r26207, r26208, r26209, r26210, r26211, r26212, r26213, r26214, r26215, r26216, r26217, r26218, r26219, r26220, r26221, r26222, r26223, r26224, r26225, r26226, r26227, r26228, r26229, r26230, r26231, r26232, r26233; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r26125, "0.5", 10, MPFR_RNDN); mpfr_init(r26126); mpfr_init(r26127); mpfr_init_set_str(r26128, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26129); mpfr_init_set_str(r26130, "4", 10, MPFR_RNDN); mpfr_init(r26131); mpfr_init(r26132); mpfr_init_set_str(r26133, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26134); mpfr_init_set_str(r26135, "6", 10, MPFR_RNDN); mpfr_init(r26136); mpfr_init(r26137); mpfr_init(r26138); mpfr_init(r26139); mpfr_init_set_str(r26140, "-0.00014612633325219868", 10, MPFR_RNDN); mpfr_init(r26141); mpfr_init(r26142); mpfr_init_set_str(r26143, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26144); mpfr_init_set_str(r26145, "1", 10, MPFR_RNDN); mpfr_init(r26146); mpfr_init(r26147); mpfr_init_set_str(r26148, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r26149, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26150); mpfr_init(r26151); mpfr_init(r26152); mpfr_init(r26153); mpfr_init_set_str(r26154, "3", 10, MPFR_RNDN); mpfr_init(r26155); mpfr_init(r26156); mpfr_init_set_str(r26157, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26158); mpfr_init_set_str(r26159, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26160); mpfr_init(r26161); mpfr_init(r26162); mpfr_init(r26163); mpfr_init_set_str(r26164, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26165); mpfr_init(r26166); mpfr_init(r26167); mpfr_init(r26168); mpfr_init_set_str(r26169, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26170); mpfr_init(r26171); mpfr_init(r26172); mpfr_init(r26173); mpfr_init(r26174); mpfr_init_set_str(r26175, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r26176, "2", 10, MPFR_RNDN); mpfr_init(r26177); mpfr_init(r26178); mpfr_init(r26179); mpfr_init(r26180); mpfr_init(r26181); mpfr_init_set_str(r26182, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26183); mpfr_init(r26184); mpfr_init_set_str(r26185, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26186); mpfr_init(r26187); mpfr_init(r26188); mpfr_init(r26189); mpfr_init(r26190); mpfr_init(r26191); mpfr_init(r26192); mpfr_init(r26193); mpfr_init_set_str(r26194, "6.050122732948135e-08", 10, MPFR_RNDN); mpfr_init(r26195); mpfr_init(r26196); mpfr_init(r26197); mpfr_init(r26198); mpfr_init(r26199); mpfr_init(r26200); mpfr_init(r26201); mpfr_init(r26202); mpfr_init(r26203); mpfr_init(r26204); mpfr_init(r26205); mpfr_init(r26206); mpfr_init(r26207); mpfr_init(r26208); mpfr_init(r26209); mpfr_init(r26210); mpfr_init(r26211); mpfr_init(r26212); mpfr_init(r26213); mpfr_init(r26214); mpfr_init(r26215); mpfr_init(r26216); mpfr_init(r26217); mpfr_init(r26218); mpfr_init(r26219); mpfr_init(r26220); mpfr_init(r26221); mpfr_init(r26222); mpfr_init(r26223); mpfr_init(r26224); mpfr_init(r26225); mpfr_init(r26226); mpfr_init(r26227); mpfr_init(r26228); mpfr_init(r26229); mpfr_init(r26230); mpfr_init(r26231); mpfr_init(r26232); mpfr_init(r26233); } double f_dm(double x) { ; mpfr_set_d(r26126, x, MPFR_RNDN); mpfr_div(r26127, r26125, r26126, MPFR_RNDN); ; mpfr_mul(r26129, r26128, r26126, MPFR_RNDN); ; mpfr_pow(r26131, r26126, r26130, MPFR_RNDN); mpfr_div(r26132, r26129, r26131, MPFR_RNDN); ; mpfr_mul(r26134, r26133, r26126, MPFR_RNDN); ; mpfr_pow(r26136, r26126, r26135, MPFR_RNDN); mpfr_div(r26137, r26134, r26136, MPFR_RNDN); mpfr_add(r26138, r26132, r26137, MPFR_RNDN); mpfr_add(r26139, r26127, r26138, MPFR_RNDN); ; mpfr_set_si(r26141, mpfr_cmp(r26139, r26140) <= 0, MPFR_RNDN); mpfr_mul(r26142, r26126, r26126, MPFR_RNDN); ; mpfr_mul(r26144, r26142, r26143, MPFR_RNDN); ; mpfr_add(r26146, r26144, r26145, MPFR_RNDN); mpfr_mul(r26147, r26142, r26142, MPFR_RNDN); ; ; mpfr_mul(r26150, r26126, r26149, MPFR_RNDN); mpfr_mul(r26151, r26150, r26126, MPFR_RNDN); mpfr_add(r26152, r26148, r26151, MPFR_RNDN); mpfr_mul(r26153, r26147, r26152, MPFR_RNDN); ; mpfr_add(r26155, r26154, r26145, MPFR_RNDN); mpfr_pow(r26156, r26142, r26155, MPFR_RNDN); ; mpfr_mul(r26158, r26142, r26157, MPFR_RNDN); ; mpfr_add(r26160, r26158, r26159, MPFR_RNDN); mpfr_mul(r26161, r26156, r26160, MPFR_RNDN); mpfr_add(r26162, r26153, r26161, MPFR_RNDN); mpfr_add(r26163, r26146, r26162, MPFR_RNDN); ; mpfr_mul(r26165, r26164, r26126, MPFR_RNDN); mpfr_pow(r26166, r26126, r26154, MPFR_RNDN); mpfr_mul(r26167, r26165, r26166, MPFR_RNDN); mpfr_add(r26168, r26145, r26167, MPFR_RNDN); ; mpfr_mul(r26170, r26169, r26126, MPFR_RNDN); mpfr_mul(r26171, r26170, r26126, MPFR_RNDN); mpfr_add(r26172, r26168, r26171, MPFR_RNDN); mpfr_mul(r26173, r26166, r26142, MPFR_RNDN); mpfr_mul(r26174, r26173, r26173, MPFR_RNDN); ; ; mpfr_mul(r26177, r26157, r26176, MPFR_RNDN); mpfr_mul(r26178, r26177, r26142, MPFR_RNDN); mpfr_add(r26179, r26175, r26178, MPFR_RNDN); mpfr_mul(r26180, r26174, r26179, MPFR_RNDN); mpfr_mul(r26181, r26166, r26166, MPFR_RNDN); ; mpfr_mul(r26183, r26126, r26182, MPFR_RNDN); mpfr_mul(r26184, r26183, r26126, MPFR_RNDN); ; mpfr_add(r26186, r26184, r26185, MPFR_RNDN); mpfr_mul(r26187, r26181, r26186, MPFR_RNDN); mpfr_add(r26188, r26180, r26187, MPFR_RNDN); mpfr_add(r26189, r26172, r26188, MPFR_RNDN); mpfr_div(r26190, r26163, r26189, MPFR_RNDN); mpfr_sqrt(r26191, r26190, MPFR_RNDN); mpfr_mul(r26192, r26191, r26191, MPFR_RNDN); mpfr_mul(r26193, r26192, r26126, MPFR_RNDN); ; mpfr_set_si(r26195, mpfr_cmp(r26139, r26194) <= 0, MPFR_RNDN); mpfr_mul(r26196, r26159, r26142, MPFR_RNDN); mpfr_pow(r26197, r26142, r26154, MPFR_RNDN); mpfr_mul(r26198, r26196, r26197, MPFR_RNDN); mpfr_mul(r26199, r26147, r26147, MPFR_RNDN); mpfr_mul(r26200, r26199, r26158, MPFR_RNDN); mpfr_add(r26201, r26198, r26200, MPFR_RNDN); mpfr_mul(r26202, r26149, r26142, MPFR_RNDN); mpfr_mul(r26203, r26202, r26147, MPFR_RNDN); mpfr_mul(r26204, r26148, r26126, MPFR_RNDN); mpfr_mul(r26205, r26204, r26166, MPFR_RNDN); mpfr_add(r26206, r26144, r26205, MPFR_RNDN); mpfr_add(r26207, r26145, r26206, MPFR_RNDN); mpfr_add(r26208, r26203, r26207, MPFR_RNDN); mpfr_add(r26209, r26201, r26208, MPFR_RNDN); mpfr_pow(r26210, r26209, r26154, MPFR_RNDN); mpfr_mul(r26211, r26142, r26182, MPFR_RNDN); mpfr_mul(r26212, r26197, r26211, MPFR_RNDN); mpfr_add(r26213, r26171, r26145, MPFR_RNDN); mpfr_add(r26214, r26212, r26213, MPFR_RNDN); mpfr_mul(r26215, r26142, r26164, MPFR_RNDN); mpfr_mul(r26216, r26185, r26126, MPFR_RNDN); mpfr_mul(r26217, r26216, r26166, MPFR_RNDN); mpfr_add(r26218, r26215, r26217, MPFR_RNDN); mpfr_mul(r26219, r26142, r26218, MPFR_RNDN); mpfr_add(r26220, r26214, r26219, MPFR_RNDN); mpfr_mul(r26221, r26176, r26157, MPFR_RNDN); mpfr_mul(r26222, r26147, r26221, MPFR_RNDN); mpfr_mul(r26223, r26199, r26222, MPFR_RNDN); mpfr_mul(r26224, r26175, r26142, MPFR_RNDN); mpfr_mul(r26225, r26224, r26199, MPFR_RNDN); mpfr_add(r26226, r26223, r26225, MPFR_RNDN); mpfr_add(r26227, r26220, r26226, MPFR_RNDN); mpfr_pow(r26228, r26227, r26154, MPFR_RNDN); mpfr_div(r26229, r26210, r26228, MPFR_RNDN); mpfr_cbrt(r26230, r26229, MPFR_RNDN); mpfr_mul(r26231, r26230, r26126, MPFR_RNDN); if (mpfr_get_si(r26195, MPFR_RNDN)) { mpfr_set(r26232, r26139, MPFR_RNDN); } else { mpfr_set(r26232, r26231, MPFR_RNDN); }; if (mpfr_get_si(r26141, MPFR_RNDN)) { mpfr_set(r26233, r26193, MPFR_RNDN); } else { mpfr_set(r26233, r26232, MPFR_RNDN); }; return mpfr_get_d(r26233, MPFR_RNDN); }