#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25793 = 1; float r25794 = 0.1049934947; float r25795 = x; float r25796 = r25795 * r25795; float r25797 = r25794 * r25796; float r25798 = r25793 + r25797; float r25799 = 0.0424060604; float r25800 = r25796 * r25796; float r25801 = r25799 * r25800; float r25802 = r25798 + r25801; float r25803 = 0.0072644182; float r25804 = r25800 * r25796; float r25805 = r25803 * r25804; float r25806 = r25802 + r25805; float r25807 = 0.0005064034; float r25808 = r25804 * r25796; float r25809 = r25807 * r25808; float r25810 = r25806 + r25809; float r25811 = 0.0001789971; float r25812 = r25808 * r25796; float r25813 = r25811 * r25812; float r25814 = r25810 + r25813; float r25815 = 0.7715471019; float r25816 = r25815 * r25796; float r25817 = r25793 + r25816; float r25818 = 0.2909738639; float r25819 = r25818 * r25800; float r25820 = r25817 + r25819; float r25821 = 0.0694555761; float r25822 = r25821 * r25804; float r25823 = r25820 + r25822; float r25824 = 0.0140005442; float r25825 = r25824 * r25808; float r25826 = r25823 + r25825; float r25827 = 0.0008327945; float r25828 = r25827 * r25812; float r25829 = r25826 + r25828; float r25830 = 2; float r25831 = r25830 * r25811; float r25832 = r25812 * r25796; float r25833 = r25831 * r25832; float r25834 = r25829 + r25833; float r25835 = r25814 / r25834; float r25836 = r25835 * r25795; return r25836; } double f_id(double x) { double r25837 = 1; double r25838 = 0.1049934947; double r25839 = x; double r25840 = r25839 * r25839; double r25841 = r25838 * r25840; double r25842 = r25837 + r25841; double r25843 = 0.0424060604; double r25844 = r25840 * r25840; double r25845 = r25843 * r25844; double r25846 = r25842 + r25845; double r25847 = 0.0072644182; double r25848 = r25844 * r25840; double r25849 = r25847 * r25848; double r25850 = r25846 + r25849; double r25851 = 0.0005064034; double r25852 = r25848 * r25840; double r25853 = r25851 * r25852; double r25854 = r25850 + r25853; double r25855 = 0.0001789971; double r25856 = r25852 * r25840; double r25857 = r25855 * r25856; double r25858 = r25854 + r25857; double r25859 = 0.7715471019; double r25860 = r25859 * r25840; double r25861 = r25837 + r25860; double r25862 = 0.2909738639; double r25863 = r25862 * r25844; double r25864 = r25861 + r25863; double r25865 = 0.0694555761; double r25866 = r25865 * r25848; double r25867 = r25864 + r25866; double r25868 = 0.0140005442; double r25869 = r25868 * r25852; double r25870 = r25867 + r25869; double r25871 = 0.0008327945; double r25872 = r25871 * r25856; double r25873 = r25870 + r25872; double r25874 = 2; double r25875 = r25874 * r25855; double r25876 = r25856 * r25840; double r25877 = r25875 * r25876; double r25878 = r25873 + r25877; double r25879 = r25858 / r25878; double r25880 = r25879 * r25839; return r25880; } double f_of(float x) { float r25881 = x; float r25882 = -5.161343281950335e+18; bool r25883 = r25881 <= r25882; float r25884 = 0.5; float r25885 = r25884 / r25881; float r25886 = 0.2514179000665375; float r25887 = r25886 * r25881; float r25888 = 4; float r25889 = pow(r25881, r25888); float r25890 = r25887 / r25889; float r25891 = 0.15298196345929327; float r25892 = r25891 * r25881; float r25893 = 6; float r25894 = pow(r25881, r25893); float r25895 = r25892 / r25894; float r25896 = r25890 + r25895; float r25897 = r25885 + r25896; float r25898 = 510048.78143239743; bool r25899 = r25881 <= r25898; float r25900 = 0.1049934947; float r25901 = r25881 * r25900; float r25902 = r25901 * r25881; float r25903 = 1; float r25904 = 0.0424060604; float r25905 = r25904 * r25881; float r25906 = 3; float r25907 = pow(r25881, r25906); float r25908 = r25905 * r25907; float r25909 = r25903 + r25908; float r25910 = r25902 + r25909; float r25911 = r25881 * r25881; float r25912 = r25911 * r25911; float r25913 = 0.0072644182; float r25914 = r25911 * r25913; float r25915 = r25912 * r25914; float r25916 = 0.0005064034; float r25917 = r25916 * r25911; float r25918 = r25907 * r25907; float r25919 = r25917 * r25918; float r25920 = r25915 + r25919; float r25921 = r25912 * r25912; float r25922 = 0.0001789971; float r25923 = r25911 * r25922; float r25924 = r25921 * r25923; float r25925 = r25920 + r25924; float r25926 = r25910 + r25925; float r25927 = 0.0008327945; float r25928 = r25911 * r25927; float r25929 = r25921 * r25928; float r25930 = 2; float r25931 = r25922 * r25930; float r25932 = r25912 * r25931; float r25933 = r25921 * r25932; float r25934 = r25929 + r25933; float r25935 = 0.0140005442; float r25936 = r25935 * r25911; float r25937 = r25936 * r25918; float r25938 = 0.7715471019; float r25939 = r25938 * r25911; float r25940 = r25939 + r25903; float r25941 = r25937 + r25940; float r25942 = 0.2909738639; float r25943 = r25942 * r25911; float r25944 = 0.0694555761; float r25945 = r25944 * r25881; float r25946 = r25945 * r25907; float r25947 = r25943 + r25946; float r25948 = r25911 * r25947; float r25949 = r25941 + r25948; float r25950 = r25934 + r25949; float r25951 = r25926 / r25950; float r25952 = pow(r25951, r25906); float r25953 = cbrt(r25952); float r25954 = r25953 * r25881; float r25955 = r25899 ? r25954 : r25897; float r25956 = r25883 ? r25897 : r25955; return r25956; } double f_od(double x) { double r25957 = x; double r25958 = -5.161343281950335e+18; bool r25959 = r25957 <= r25958; double r25960 = 0.5; double r25961 = r25960 / r25957; double r25962 = 0.2514179000665375; double r25963 = r25962 * r25957; double r25964 = 4; double r25965 = pow(r25957, r25964); double r25966 = r25963 / r25965; double r25967 = 0.15298196345929327; double r25968 = r25967 * r25957; double r25969 = 6; double r25970 = pow(r25957, r25969); double r25971 = r25968 / r25970; double r25972 = r25966 + r25971; double r25973 = r25961 + r25972; double r25974 = 510048.78143239743; bool r25975 = r25957 <= r25974; double r25976 = 0.1049934947; double r25977 = r25957 * r25976; double r25978 = r25977 * r25957; double r25979 = 1; double r25980 = 0.0424060604; double r25981 = r25980 * r25957; double r25982 = 3; double r25983 = pow(r25957, r25982); double r25984 = r25981 * r25983; double r25985 = r25979 + r25984; double r25986 = r25978 + r25985; double r25987 = r25957 * r25957; double r25988 = r25987 * r25987; double r25989 = 0.0072644182; double r25990 = r25987 * r25989; double r25991 = r25988 * r25990; double r25992 = 0.0005064034; double r25993 = r25992 * r25987; double r25994 = r25983 * r25983; double r25995 = r25993 * r25994; double r25996 = r25991 + r25995; double r25997 = r25988 * r25988; double r25998 = 0.0001789971; double r25999 = r25987 * r25998; double r26000 = r25997 * r25999; double r26001 = r25996 + r26000; double r26002 = r25986 + r26001; double r26003 = 0.0008327945; double r26004 = r25987 * r26003; double r26005 = r25997 * r26004; double r26006 = 2; double r26007 = r25998 * r26006; double r26008 = r25988 * r26007; double r26009 = r25997 * r26008; double r26010 = r26005 + r26009; double r26011 = 0.0140005442; double r26012 = r26011 * r25987; double r26013 = r26012 * r25994; double r26014 = 0.7715471019; double r26015 = r26014 * r25987; double r26016 = r26015 + r25979; double r26017 = r26013 + r26016; double r26018 = 0.2909738639; double r26019 = r26018 * r25987; double r26020 = 0.0694555761; double r26021 = r26020 * r25957; double r26022 = r26021 * r25983; double r26023 = r26019 + r26022; double r26024 = r25987 * r26023; double r26025 = r26017 + r26024; double r26026 = r26010 + r26025; double r26027 = r26002 / r26026; double r26028 = pow(r26027, r25982); double r26029 = cbrt(r26028); double r26030 = r26029 * r25957; double r26031 = r25975 ? r26030 : r25973; double r26032 = r25959 ? r25973 : r26031; return r26032; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r26033, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26034, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init(r26038); mpfr_init_set_str(r26039, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26040); mpfr_init(r26041); mpfr_init(r26042); mpfr_init_set_str(r26043, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); mpfr_init_set_str(r26047, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26048); mpfr_init(r26049); mpfr_init(r26050); mpfr_init_set_str(r26051, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26052); mpfr_init(r26053); mpfr_init(r26054); mpfr_init_set_str(r26055, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26056); mpfr_init(r26057); mpfr_init_set_str(r26058, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26059); mpfr_init(r26060); mpfr_init_set_str(r26061, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26062); mpfr_init(r26063); mpfr_init_set_str(r26064, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26065); mpfr_init(r26066); mpfr_init_set_str(r26067, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26068); mpfr_init(r26069); mpfr_init_set_str(r26070, "2", 10, MPFR_RNDN); mpfr_init(r26071); mpfr_init(r26072); mpfr_init(r26073); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); } double f_im(double x) { ; ; mpfr_set_d(r26035, x, MPFR_RNDN); mpfr_mul(r26036, r26035, r26035, MPFR_RNDN); mpfr_mul(r26037, r26034, r26036, MPFR_RNDN); mpfr_add(r26038, r26033, r26037, MPFR_RNDN); ; mpfr_mul(r26040, r26036, r26036, MPFR_RNDN); mpfr_mul(r26041, r26039, r26040, MPFR_RNDN); mpfr_add(r26042, r26038, r26041, MPFR_RNDN); ; mpfr_mul(r26044, r26040, r26036, MPFR_RNDN); mpfr_mul(r26045, r26043, r26044, MPFR_RNDN); mpfr_add(r26046, r26042, r26045, MPFR_RNDN); ; mpfr_mul(r26048, r26044, r26036, MPFR_RNDN); mpfr_mul(r26049, r26047, r26048, MPFR_RNDN); mpfr_add(r26050, r26046, r26049, MPFR_RNDN); ; mpfr_mul(r26052, r26048, r26036, MPFR_RNDN); mpfr_mul(r26053, r26051, r26052, MPFR_RNDN); mpfr_add(r26054, r26050, r26053, MPFR_RNDN); ; mpfr_mul(r26056, r26055, r26036, MPFR_RNDN); mpfr_add(r26057, r26033, r26056, MPFR_RNDN); ; mpfr_mul(r26059, r26058, r26040, MPFR_RNDN); mpfr_add(r26060, r26057, r26059, MPFR_RNDN); ; mpfr_mul(r26062, r26061, r26044, MPFR_RNDN); mpfr_add(r26063, r26060, r26062, MPFR_RNDN); ; mpfr_mul(r26065, r26064, r26048, MPFR_RNDN); mpfr_add(r26066, r26063, r26065, MPFR_RNDN); ; mpfr_mul(r26068, r26067, r26052, MPFR_RNDN); mpfr_add(r26069, r26066, r26068, MPFR_RNDN); ; mpfr_mul(r26071, r26070, r26051, MPFR_RNDN); mpfr_mul(r26072, r26052, r26036, MPFR_RNDN); mpfr_mul(r26073, r26071, r26072, MPFR_RNDN); mpfr_add(r26074, r26069, r26073, MPFR_RNDN); mpfr_div(r26075, r26054, r26074, MPFR_RNDN); mpfr_mul(r26076, r26075, r26035, MPFR_RNDN); return mpfr_get_d(r26076, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r26077); mpfr_init_set_str(r26078, "-5.161343281950335e+18", 10, MPFR_RNDN); mpfr_init(r26079); mpfr_init_set_str(r26080, "0.5", 10, MPFR_RNDN); mpfr_init(r26081); mpfr_init_set_str(r26082, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26083); mpfr_init_set_str(r26084, "4", 10, MPFR_RNDN); mpfr_init(r26085); mpfr_init(r26086); mpfr_init_set_str(r26087, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26088); mpfr_init_set_str(r26089, "6", 10, MPFR_RNDN); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); mpfr_init(r26093); mpfr_init_set_str(r26094, "510048.78143239743", 10, MPFR_RNDN); mpfr_init(r26095); mpfr_init_set_str(r26096, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26097); mpfr_init(r26098); mpfr_init_set_str(r26099, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26100, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26101); mpfr_init_set_str(r26102, "3", 10, MPFR_RNDN); mpfr_init(r26103); mpfr_init(r26104); mpfr_init(r26105); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init_set_str(r26109, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26110); mpfr_init(r26111); mpfr_init_set_str(r26112, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26113); mpfr_init(r26114); mpfr_init(r26115); mpfr_init(r26116); mpfr_init(r26117); mpfr_init_set_str(r26118, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26119); mpfr_init(r26120); mpfr_init(r26121); mpfr_init(r26122); mpfr_init_set_str(r26123, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26124); mpfr_init(r26125); mpfr_init_set_str(r26126, "2", 10, MPFR_RNDN); mpfr_init(r26127); mpfr_init(r26128); mpfr_init(r26129); mpfr_init(r26130); mpfr_init_set_str(r26131, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26132); mpfr_init(r26133); mpfr_init_set_str(r26134, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26135); mpfr_init(r26136); mpfr_init(r26137); mpfr_init_set_str(r26138, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26139); mpfr_init_set_str(r26140, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26141); mpfr_init(r26142); mpfr_init(r26143); mpfr_init(r26144); mpfr_init(r26145); mpfr_init(r26146); mpfr_init(r26147); mpfr_init(r26148); mpfr_init(r26149); mpfr_init(r26150); mpfr_init(r26151); mpfr_init(r26152); } double f_fm(double x) { mpfr_set_d(r26077, x, MPFR_RNDN); ; mpfr_set_si(r26079, mpfr_cmp(r26077, r26078) <= 0, MPFR_RNDN); ; mpfr_div(r26081, r26080, r26077, MPFR_RNDN); ; mpfr_mul(r26083, r26082, r26077, MPFR_RNDN); ; mpfr_pow(r26085, r26077, r26084, MPFR_RNDN); mpfr_div(r26086, r26083, r26085, MPFR_RNDN); ; mpfr_mul(r26088, r26087, r26077, MPFR_RNDN); ; mpfr_pow(r26090, r26077, r26089, MPFR_RNDN); mpfr_div(r26091, r26088, r26090, MPFR_RNDN); mpfr_add(r26092, r26086, r26091, MPFR_RNDN); mpfr_add(r26093, r26081, r26092, MPFR_RNDN); ; mpfr_set_si(r26095, mpfr_cmp(r26077, r26094) <= 0, MPFR_RNDN); ; mpfr_mul(r26097, r26077, r26096, MPFR_RNDN); mpfr_mul(r26098, r26097, r26077, MPFR_RNDN); ; ; mpfr_mul(r26101, r26100, r26077, MPFR_RNDN); ; mpfr_pow(r26103, r26077, r26102, MPFR_RNDN); mpfr_mul(r26104, r26101, r26103, MPFR_RNDN); mpfr_add(r26105, r26099, r26104, MPFR_RNDN); mpfr_add(r26106, r26098, r26105, MPFR_RNDN); mpfr_mul(r26107, r26077, r26077, MPFR_RNDN); mpfr_mul(r26108, r26107, r26107, MPFR_RNDN); ; mpfr_mul(r26110, r26107, r26109, MPFR_RNDN); mpfr_mul(r26111, r26108, r26110, MPFR_RNDN); ; mpfr_mul(r26113, r26112, r26107, MPFR_RNDN); mpfr_mul(r26114, r26103, r26103, MPFR_RNDN); mpfr_mul(r26115, r26113, r26114, MPFR_RNDN); mpfr_add(r26116, r26111, r26115, MPFR_RNDN); mpfr_mul(r26117, r26108, r26108, MPFR_RNDN); ; mpfr_mul(r26119, r26107, r26118, MPFR_RNDN); mpfr_mul(r26120, r26117, r26119, MPFR_RNDN); mpfr_add(r26121, r26116, r26120, MPFR_RNDN); mpfr_add(r26122, r26106, r26121, MPFR_RNDN); ; mpfr_mul(r26124, r26107, r26123, MPFR_RNDN); mpfr_mul(r26125, r26117, r26124, MPFR_RNDN); ; mpfr_mul(r26127, r26118, r26126, MPFR_RNDN); mpfr_mul(r26128, r26108, r26127, MPFR_RNDN); mpfr_mul(r26129, r26117, r26128, MPFR_RNDN); mpfr_add(r26130, r26125, r26129, MPFR_RNDN); ; mpfr_mul(r26132, r26131, r26107, MPFR_RNDN); mpfr_mul(r26133, r26132, r26114, MPFR_RNDN); ; mpfr_mul(r26135, r26134, r26107, MPFR_RNDN); mpfr_add(r26136, r26135, r26099, MPFR_RNDN); mpfr_add(r26137, r26133, r26136, MPFR_RNDN); ; mpfr_mul(r26139, r26138, r26107, MPFR_RNDN); ; mpfr_mul(r26141, r26140, r26077, MPFR_RNDN); mpfr_mul(r26142, r26141, r26103, MPFR_RNDN); mpfr_add(r26143, r26139, r26142, MPFR_RNDN); mpfr_mul(r26144, r26107, r26143, MPFR_RNDN); mpfr_add(r26145, r26137, r26144, MPFR_RNDN); mpfr_add(r26146, r26130, r26145, MPFR_RNDN); mpfr_div(r26147, r26122, r26146, MPFR_RNDN); mpfr_pow(r26148, r26147, r26102, MPFR_RNDN); mpfr_cbrt(r26149, r26148, MPFR_RNDN); mpfr_mul(r26150, r26149, r26077, MPFR_RNDN); if (mpfr_get_si(r26095, MPFR_RNDN)) { mpfr_set(r26151, r26150, MPFR_RNDN); } else { mpfr_set(r26151, r26093, MPFR_RNDN); }; if (mpfr_get_si(r26079, MPFR_RNDN)) { mpfr_set(r26152, r26093, MPFR_RNDN); } else { mpfr_set(r26152, r26151, MPFR_RNDN); }; return mpfr_get_d(r26152, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26153); mpfr_init_set_str(r26154, "-5.161343281950335e+18", 10, MPFR_RNDN); mpfr_init(r26155); mpfr_init_set_str(r26156, "0.5", 10, MPFR_RNDN); mpfr_init(r26157); mpfr_init_set_str(r26158, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init(r26159); mpfr_init_set_str(r26160, "4", 10, MPFR_RNDN); mpfr_init(r26161); mpfr_init(r26162); mpfr_init_set_str(r26163, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init(r26164); mpfr_init_set_str(r26165, "6", 10, MPFR_RNDN); mpfr_init(r26166); mpfr_init(r26167); mpfr_init(r26168); mpfr_init(r26169); mpfr_init_set_str(r26170, "510048.78143239743", 10, MPFR_RNDN); mpfr_init(r26171); mpfr_init_set_str(r26172, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26173); mpfr_init(r26174); mpfr_init_set_str(r26175, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26176, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26177); mpfr_init_set_str(r26178, "3", 10, MPFR_RNDN); mpfr_init(r26179); mpfr_init(r26180); mpfr_init(r26181); mpfr_init(r26182); mpfr_init(r26183); mpfr_init(r26184); mpfr_init_set_str(r26185, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26186); mpfr_init(r26187); mpfr_init_set_str(r26188, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26189); mpfr_init(r26190); mpfr_init(r26191); mpfr_init(r26192); mpfr_init(r26193); mpfr_init_set_str(r26194, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26195); mpfr_init(r26196); mpfr_init(r26197); mpfr_init(r26198); mpfr_init_set_str(r26199, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26200); mpfr_init(r26201); mpfr_init_set_str(r26202, "2", 10, MPFR_RNDN); mpfr_init(r26203); mpfr_init(r26204); mpfr_init(r26205); mpfr_init(r26206); mpfr_init_set_str(r26207, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26208); mpfr_init(r26209); mpfr_init_set_str(r26210, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26211); mpfr_init(r26212); mpfr_init(r26213); mpfr_init_set_str(r26214, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26215); mpfr_init_set_str(r26216, "0.0694555761", 10, MPFR_RNDN); 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); } double f_dm(double x) { mpfr_set_d(r26153, x, MPFR_RNDN); ; mpfr_set_si(r26155, mpfr_cmp(r26153, r26154) <= 0, MPFR_RNDN); ; mpfr_div(r26157, r26156, r26153, MPFR_RNDN); ; mpfr_mul(r26159, r26158, r26153, MPFR_RNDN); ; mpfr_pow(r26161, r26153, r26160, MPFR_RNDN); mpfr_div(r26162, r26159, r26161, MPFR_RNDN); ; mpfr_mul(r26164, r26163, r26153, MPFR_RNDN); ; mpfr_pow(r26166, r26153, r26165, MPFR_RNDN); mpfr_div(r26167, r26164, r26166, MPFR_RNDN); mpfr_add(r26168, r26162, r26167, MPFR_RNDN); mpfr_add(r26169, r26157, r26168, MPFR_RNDN); ; mpfr_set_si(r26171, mpfr_cmp(r26153, r26170) <= 0, MPFR_RNDN); ; mpfr_mul(r26173, r26153, r26172, MPFR_RNDN); mpfr_mul(r26174, r26173, r26153, MPFR_RNDN); ; ; mpfr_mul(r26177, r26176, r26153, MPFR_RNDN); ; mpfr_pow(r26179, r26153, r26178, MPFR_RNDN); mpfr_mul(r26180, r26177, r26179, MPFR_RNDN); mpfr_add(r26181, r26175, r26180, MPFR_RNDN); mpfr_add(r26182, r26174, r26181, MPFR_RNDN); mpfr_mul(r26183, r26153, r26153, MPFR_RNDN); mpfr_mul(r26184, r26183, r26183, MPFR_RNDN); ; mpfr_mul(r26186, r26183, r26185, MPFR_RNDN); mpfr_mul(r26187, r26184, r26186, MPFR_RNDN); ; mpfr_mul(r26189, r26188, r26183, MPFR_RNDN); mpfr_mul(r26190, r26179, r26179, MPFR_RNDN); mpfr_mul(r26191, r26189, r26190, MPFR_RNDN); mpfr_add(r26192, r26187, r26191, MPFR_RNDN); mpfr_mul(r26193, r26184, r26184, MPFR_RNDN); ; mpfr_mul(r26195, r26183, r26194, MPFR_RNDN); mpfr_mul(r26196, r26193, r26195, MPFR_RNDN); mpfr_add(r26197, r26192, r26196, MPFR_RNDN); mpfr_add(r26198, r26182, r26197, MPFR_RNDN); ; mpfr_mul(r26200, r26183, r26199, MPFR_RNDN); mpfr_mul(r26201, r26193, r26200, MPFR_RNDN); ; mpfr_mul(r26203, r26194, r26202, MPFR_RNDN); mpfr_mul(r26204, r26184, r26203, MPFR_RNDN); mpfr_mul(r26205, r26193, r26204, MPFR_RNDN); mpfr_add(r26206, r26201, r26205, MPFR_RNDN); ; mpfr_mul(r26208, r26207, r26183, MPFR_RNDN); mpfr_mul(r26209, r26208, r26190, MPFR_RNDN); ; mpfr_mul(r26211, r26210, r26183, MPFR_RNDN); mpfr_add(r26212, r26211, r26175, MPFR_RNDN); mpfr_add(r26213, r26209, r26212, MPFR_RNDN); ; mpfr_mul(r26215, r26214, r26183, MPFR_RNDN); ; mpfr_mul(r26217, r26216, r26153, MPFR_RNDN); mpfr_mul(r26218, r26217, r26179, MPFR_RNDN); mpfr_add(r26219, r26215, r26218, MPFR_RNDN); mpfr_mul(r26220, r26183, r26219, MPFR_RNDN); mpfr_add(r26221, r26213, r26220, MPFR_RNDN); mpfr_add(r26222, r26206, r26221, MPFR_RNDN); mpfr_div(r26223, r26198, r26222, MPFR_RNDN); mpfr_pow(r26224, r26223, r26178, MPFR_RNDN); mpfr_cbrt(r26225, r26224, MPFR_RNDN); mpfr_mul(r26226, r26225, r26153, MPFR_RNDN); if (mpfr_get_si(r26171, MPFR_RNDN)) { mpfr_set(r26227, r26226, MPFR_RNDN); } else { mpfr_set(r26227, r26169, MPFR_RNDN); }; if (mpfr_get_si(r26155, MPFR_RNDN)) { mpfr_set(r26228, r26169, MPFR_RNDN); } else { mpfr_set(r26228, r26227, MPFR_RNDN); }; return mpfr_get_d(r26228, MPFR_RNDN); }