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