#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r25720 = 1; float r25721 = 0.1049934947; float r25722 = x; float r25723 = r25722 * r25722; float r25724 = r25721 * r25723; float r25725 = r25720 + r25724; float r25726 = 0.0424060604; float r25727 = r25723 * r25723; float r25728 = r25726 * r25727; float r25729 = r25725 + r25728; float r25730 = 0.0072644182; float r25731 = r25727 * r25723; float r25732 = r25730 * r25731; float r25733 = r25729 + r25732; float r25734 = 0.0005064034; float r25735 = r25731 * r25723; float r25736 = r25734 * r25735; float r25737 = r25733 + r25736; float r25738 = 0.0001789971; float r25739 = r25735 * r25723; float r25740 = r25738 * r25739; float r25741 = r25737 + r25740; float r25742 = 0.7715471019; float r25743 = r25742 * r25723; float r25744 = r25720 + r25743; float r25745 = 0.2909738639; float r25746 = r25745 * r25727; float r25747 = r25744 + r25746; float r25748 = 0.0694555761; float r25749 = r25748 * r25731; float r25750 = r25747 + r25749; float r25751 = 0.0140005442; float r25752 = r25751 * r25735; float r25753 = r25750 + r25752; float r25754 = 0.0008327945; float r25755 = r25754 * r25739; float r25756 = r25753 + r25755; float r25757 = 2; float r25758 = r25757 * r25738; float r25759 = r25739 * r25723; float r25760 = r25758 * r25759; float r25761 = r25756 + r25760; float r25762 = r25741 / r25761; float r25763 = r25762 * r25722; return r25763; } double f_id(double x) { double r25764 = 1; double r25765 = 0.1049934947; double r25766 = x; double r25767 = r25766 * r25766; double r25768 = r25765 * r25767; double r25769 = r25764 + r25768; double r25770 = 0.0424060604; double r25771 = r25767 * r25767; double r25772 = r25770 * r25771; double r25773 = r25769 + r25772; double r25774 = 0.0072644182; double r25775 = r25771 * r25767; double r25776 = r25774 * r25775; double r25777 = r25773 + r25776; double r25778 = 0.0005064034; double r25779 = r25775 * r25767; double r25780 = r25778 * r25779; double r25781 = r25777 + r25780; double r25782 = 0.0001789971; double r25783 = r25779 * r25767; double r25784 = r25782 * r25783; double r25785 = r25781 + r25784; double r25786 = 0.7715471019; double r25787 = r25786 * r25767; double r25788 = r25764 + r25787; double r25789 = 0.2909738639; double r25790 = r25789 * r25771; double r25791 = r25788 + r25790; double r25792 = 0.0694555761; double r25793 = r25792 * r25775; double r25794 = r25791 + r25793; double r25795 = 0.0140005442; double r25796 = r25795 * r25779; double r25797 = r25794 + r25796; double r25798 = 0.0008327945; double r25799 = r25798 * r25783; double r25800 = r25797 + r25799; double r25801 = 2; double r25802 = r25801 * r25782; double r25803 = r25783 * r25767; double r25804 = r25802 * r25803; double r25805 = r25800 + r25804; double r25806 = r25785 / r25805; double r25807 = r25806 * r25766; return r25807; } double f_of(float x) { float r25808 = x; float r25809 = -70779069.15005009; bool r25810 = r25808 <= r25809; float r25811 = 106836641.40103465; bool r25812 = r25808 <= r25811; float r25813 = !r25812; bool r25814 = r25810 || r25813; float r25815 = 0.2514179000665375; float r25816 = 4; float r25817 = pow(r25808, r25816); float r25818 = r25815 / r25817; float r25819 = 0.15298196345929327; float r25820 = 6; float r25821 = pow(r25808, r25820); float r25822 = r25819 / r25821; float r25823 = r25818 + r25822; float r25824 = 0.5; float r25825 = r25824 / r25808; float r25826 = fma(r25823, r25808, r25825); float r25827 = r25808 * r25808; float r25828 = r25827 * r25827; float r25829 = r25828 * r25827; float r25830 = r25827 * r25829; float r25831 = r25827 * r25830; float r25832 = 0.0001789971; float r25833 = r25831 * r25832; float r25834 = 0.0072644182; float r25835 = r25834 * r25829; float r25836 = 0.1049934947; float r25837 = r25827 * r25836; float r25838 = 1; float r25839 = r25837 + r25838; float r25840 = 0.0424060604; float r25841 = r25828 * r25840; float r25842 = r25839 + r25841; float r25843 = r25835 + r25842; float r25844 = 0.0005064034; float r25845 = r25844 * r25830; float r25846 = r25843 + r25845; float r25847 = r25833 + r25846; float r25848 = 0.0008327945; float r25849 = r25831 * r25848; float r25850 = 0.7715471019; float r25851 = r25850 * r25827; float r25852 = r25851 + r25838; float r25853 = 0.2909738639; float r25854 = r25828 * r25853; float r25855 = r25852 + r25854; float r25856 = 0.0694555761; float r25857 = r25829 * r25856; float r25858 = r25855 + r25857; float r25859 = 0.0140005442; float r25860 = r25859 * r25830; float r25861 = r25858 + r25860; float r25862 = r25849 + r25861; float r25863 = 2; float r25864 = r25863 * r25832; float r25865 = r25827 * r25831; float r25866 = r25864 * r25865; float r25867 = r25862 + r25866; float r25868 = r25847 / r25867; float r25869 = r25808 * r25868; float r25870 = r25814 ? r25826 : r25869; return r25870; } double f_od(double x) { double r25871 = x; double r25872 = -70779069.15005009; bool r25873 = r25871 <= r25872; double r25874 = 106836641.40103465; bool r25875 = r25871 <= r25874; double r25876 = !r25875; bool r25877 = r25873 || r25876; double r25878 = 0.2514179000665375; double r25879 = 4; double r25880 = pow(r25871, r25879); double r25881 = r25878 / r25880; double r25882 = 0.15298196345929327; double r25883 = 6; double r25884 = pow(r25871, r25883); double r25885 = r25882 / r25884; double r25886 = r25881 + r25885; double r25887 = 0.5; double r25888 = r25887 / r25871; double r25889 = fma(r25886, r25871, r25888); double r25890 = r25871 * r25871; double r25891 = r25890 * r25890; double r25892 = r25891 * r25890; double r25893 = r25890 * r25892; double r25894 = r25890 * r25893; double r25895 = 0.0001789971; double r25896 = r25894 * r25895; double r25897 = 0.0072644182; double r25898 = r25897 * r25892; double r25899 = 0.1049934947; double r25900 = r25890 * r25899; double r25901 = 1; double r25902 = r25900 + r25901; double r25903 = 0.0424060604; double r25904 = r25891 * r25903; double r25905 = r25902 + r25904; double r25906 = r25898 + r25905; double r25907 = 0.0005064034; double r25908 = r25907 * r25893; double r25909 = r25906 + r25908; double r25910 = r25896 + r25909; double r25911 = 0.0008327945; double r25912 = r25894 * r25911; double r25913 = 0.7715471019; double r25914 = r25913 * r25890; double r25915 = r25914 + r25901; double r25916 = 0.2909738639; double r25917 = r25891 * r25916; double r25918 = r25915 + r25917; double r25919 = 0.0694555761; double r25920 = r25892 * r25919; double r25921 = r25918 + r25920; double r25922 = 0.0140005442; double r25923 = r25922 * r25893; double r25924 = r25921 + r25923; double r25925 = r25912 + r25924; double r25926 = 2; double r25927 = r25926 * r25895; double r25928 = r25890 * r25894; double r25929 = r25927 * r25928; double r25930 = r25925 + r25929; double r25931 = r25910 / r25930; double r25932 = r25871 * r25931; double r25933 = r25877 ? r25889 : r25932; return r25933; } 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 r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945, r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966, r25967, r25968, r25969, r25970, r25971, r25972, r25973, r25974, r25975, r25976, r25977; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r25934, "1", 10, MPFR_RNDN); mpfr_init_set_str(r25935, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r25936); mpfr_init(r25937); mpfr_init(r25938); mpfr_init(r25939); mpfr_init_set_str(r25940, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r25941); mpfr_init(r25942); mpfr_init(r25943); mpfr_init_set_str(r25944, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r25945); mpfr_init(r25946); mpfr_init(r25947); mpfr_init_set_str(r25948, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r25949); mpfr_init(r25950); mpfr_init(r25951); mpfr_init_set_str(r25952, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r25953); mpfr_init(r25954); mpfr_init(r25955); mpfr_init_set_str(r25956, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r25957); mpfr_init(r25958); mpfr_init_set_str(r25959, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r25960); mpfr_init(r25961); mpfr_init_set_str(r25962, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r25963); mpfr_init(r25964); mpfr_init_set_str(r25965, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r25966); mpfr_init(r25967); mpfr_init_set_str(r25968, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r25969); mpfr_init(r25970); mpfr_init_set_str(r25971, "2", 10, MPFR_RNDN); mpfr_init(r25972); mpfr_init(r25973); mpfr_init(r25974); mpfr_init(r25975); mpfr_init(r25976); mpfr_init(r25977); } double f_im(double x) { ; ; mpfr_set_d(r25936, x, MPFR_RNDN); mpfr_mul(r25937, r25936, r25936, MPFR_RNDN); mpfr_mul(r25938, r25935, r25937, MPFR_RNDN); mpfr_add(r25939, r25934, r25938, MPFR_RNDN); ; mpfr_mul(r25941, r25937, r25937, MPFR_RNDN); mpfr_mul(r25942, r25940, r25941, MPFR_RNDN); mpfr_add(r25943, r25939, r25942, MPFR_RNDN); ; mpfr_mul(r25945, r25941, r25937, MPFR_RNDN); mpfr_mul(r25946, r25944, r25945, MPFR_RNDN); mpfr_add(r25947, r25943, r25946, MPFR_RNDN); ; mpfr_mul(r25949, r25945, r25937, MPFR_RNDN); mpfr_mul(r25950, r25948, r25949, MPFR_RNDN); mpfr_add(r25951, r25947, r25950, MPFR_RNDN); ; mpfr_mul(r25953, r25949, r25937, MPFR_RNDN); mpfr_mul(r25954, r25952, r25953, MPFR_RNDN); mpfr_add(r25955, r25951, r25954, MPFR_RNDN); ; mpfr_mul(r25957, r25956, r25937, MPFR_RNDN); mpfr_add(r25958, r25934, r25957, MPFR_RNDN); ; mpfr_mul(r25960, r25959, r25941, MPFR_RNDN); mpfr_add(r25961, r25958, r25960, MPFR_RNDN); ; mpfr_mul(r25963, r25962, r25945, MPFR_RNDN); mpfr_add(r25964, r25961, r25963, MPFR_RNDN); ; mpfr_mul(r25966, r25965, r25949, MPFR_RNDN); mpfr_add(r25967, r25964, r25966, MPFR_RNDN); ; mpfr_mul(r25969, r25968, r25953, MPFR_RNDN); mpfr_add(r25970, r25967, r25969, MPFR_RNDN); ; mpfr_mul(r25972, r25971, r25952, MPFR_RNDN); mpfr_mul(r25973, r25953, r25937, MPFR_RNDN); mpfr_mul(r25974, r25972, r25973, MPFR_RNDN); mpfr_add(r25975, r25970, r25974, MPFR_RNDN); mpfr_div(r25976, r25955, r25975, MPFR_RNDN); mpfr_mul(r25977, r25976, r25936, MPFR_RNDN); return mpfr_get_d(r25977, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25978); mpfr_init_set_str(r25979, "-70779069.15005009", 10, MPFR_RNDN); mpfr_init(r25980); mpfr_init_set_str(r25981, "106836641.40103465", 10, MPFR_RNDN); mpfr_init(r25982); mpfr_init(r25983); mpfr_init(r25984); mpfr_init_set_str(r25985, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r25986, "4", 10, MPFR_RNDN); mpfr_init(r25987); mpfr_init(r25988); mpfr_init_set_str(r25989, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r25990, "6", 10, MPFR_RNDN); mpfr_init(r25991); mpfr_init(r25992); mpfr_init(r25993); mpfr_init_set_str(r25994, "0.5", 10, MPFR_RNDN); mpfr_init(r25995); mpfr_init(r25996); mpfr_init(r25997); mpfr_init(r25998); mpfr_init(r25999); mpfr_init(r26000); mpfr_init(r26001); mpfr_init_set_str(r26002, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26003); mpfr_init_set_str(r26004, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26005); mpfr_init_set_str(r26006, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26007); mpfr_init_set_str(r26008, "1", 10, MPFR_RNDN); mpfr_init(r26009); mpfr_init_set_str(r26010, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26011); mpfr_init(r26012); mpfr_init(r26013); mpfr_init_set_str(r26014, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26015); mpfr_init(r26016); mpfr_init(r26017); mpfr_init_set_str(r26018, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26019); mpfr_init_set_str(r26020, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26021); mpfr_init(r26022); mpfr_init_set_str(r26023, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26024); mpfr_init(r26025); mpfr_init_set_str(r26026, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26027); mpfr_init(r26028); mpfr_init_set_str(r26029, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26030); mpfr_init(r26031); mpfr_init(r26032); mpfr_init_set_str(r26033, "2", 10, MPFR_RNDN); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init(r26038); mpfr_init(r26039); mpfr_init(r26040); } double f_fm(double x) { mpfr_set_d(r25978, x, MPFR_RNDN); ; mpfr_set_si(r25980, mpfr_cmp(r25978, r25979) <= 0, MPFR_RNDN); ; mpfr_set_si(r25982, mpfr_cmp(r25978, r25981) <= 0, MPFR_RNDN); mpfr_set_si(r25983, !mpfr_get_si(r25982, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r25984, mpfr_get_si(r25980, MPFR_RNDN) || mpfr_get_si(r25983, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r25987, r25978, r25986, MPFR_RNDN); mpfr_div(r25988, r25985, r25987, MPFR_RNDN); ; ; mpfr_pow(r25991, r25978, r25990, MPFR_RNDN); mpfr_div(r25992, r25989, r25991, MPFR_RNDN); mpfr_add(r25993, r25988, r25992, MPFR_RNDN); ; mpfr_div(r25995, r25994, r25978, MPFR_RNDN); mpfr_fma(r25996, r25993, r25978, r25995, MPFR_RNDN); mpfr_mul(r25997, r25978, r25978, MPFR_RNDN); mpfr_mul(r25998, r25997, r25997, MPFR_RNDN); mpfr_mul(r25999, r25998, r25997, MPFR_RNDN); mpfr_mul(r26000, r25997, r25999, MPFR_RNDN); mpfr_mul(r26001, r25997, r26000, MPFR_RNDN); ; mpfr_mul(r26003, r26001, r26002, MPFR_RNDN); ; mpfr_mul(r26005, r26004, r25999, MPFR_RNDN); ; mpfr_mul(r26007, r25997, r26006, MPFR_RNDN); ; mpfr_add(r26009, r26007, r26008, MPFR_RNDN); ; mpfr_mul(r26011, r25998, r26010, MPFR_RNDN); mpfr_add(r26012, r26009, r26011, MPFR_RNDN); mpfr_add(r26013, r26005, r26012, MPFR_RNDN); ; mpfr_mul(r26015, r26014, r26000, MPFR_RNDN); mpfr_add(r26016, r26013, r26015, MPFR_RNDN); mpfr_add(r26017, r26003, r26016, MPFR_RNDN); ; mpfr_mul(r26019, r26001, r26018, MPFR_RNDN); ; mpfr_mul(r26021, r26020, r25997, MPFR_RNDN); mpfr_add(r26022, r26021, r26008, MPFR_RNDN); ; mpfr_mul(r26024, r25998, r26023, MPFR_RNDN); mpfr_add(r26025, r26022, r26024, MPFR_RNDN); ; mpfr_mul(r26027, r25999, r26026, MPFR_RNDN); mpfr_add(r26028, r26025, r26027, MPFR_RNDN); ; mpfr_mul(r26030, r26029, r26000, MPFR_RNDN); mpfr_add(r26031, r26028, r26030, MPFR_RNDN); mpfr_add(r26032, r26019, r26031, MPFR_RNDN); ; mpfr_mul(r26034, r26033, r26002, MPFR_RNDN); mpfr_mul(r26035, r25997, r26001, MPFR_RNDN); mpfr_mul(r26036, r26034, r26035, MPFR_RNDN); mpfr_add(r26037, r26032, r26036, MPFR_RNDN); mpfr_div(r26038, r26017, r26037, MPFR_RNDN); mpfr_mul(r26039, r25978, r26038, MPFR_RNDN); if (mpfr_get_si(r25984, MPFR_RNDN)) { mpfr_set(r26040, r25996, MPFR_RNDN); } else { mpfr_set(r26040, r26039, MPFR_RNDN); }; return mpfr_get_d(r26040, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26041); mpfr_init_set_str(r26042, "-70779069.15005009", 10, MPFR_RNDN); mpfr_init(r26043); mpfr_init_set_str(r26044, "106836641.40103465", 10, MPFR_RNDN); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); mpfr_init_set_str(r26048, "0.2514179000665375", 10, MPFR_RNDN); mpfr_init_set_str(r26049, "4", 10, MPFR_RNDN); mpfr_init(r26050); mpfr_init(r26051); mpfr_init_set_str(r26052, "0.15298196345929327", 10, MPFR_RNDN); mpfr_init_set_str(r26053, "6", 10, MPFR_RNDN); mpfr_init(r26054); mpfr_init(r26055); mpfr_init(r26056); mpfr_init_set_str(r26057, "0.5", 10, MPFR_RNDN); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); mpfr_init_set_str(r26065, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r26066); mpfr_init_set_str(r26067, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r26068); mpfr_init_set_str(r26069, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r26070); mpfr_init_set_str(r26071, "1", 10, MPFR_RNDN); mpfr_init(r26072); mpfr_init_set_str(r26073, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init_set_str(r26077, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init_set_str(r26081, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r26082); mpfr_init_set_str(r26083, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r26084); mpfr_init(r26085); mpfr_init_set_str(r26086, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r26087); mpfr_init(r26088); mpfr_init_set_str(r26089, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r26090); mpfr_init(r26091); mpfr_init_set_str(r26092, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r26093); mpfr_init(r26094); mpfr_init(r26095); mpfr_init_set_str(r26096, "2", 10, MPFR_RNDN); mpfr_init(r26097); mpfr_init(r26098); mpfr_init(r26099); mpfr_init(r26100); mpfr_init(r26101); mpfr_init(r26102); mpfr_init(r26103); } double f_dm(double x) { mpfr_set_d(r26041, x, MPFR_RNDN); ; mpfr_set_si(r26043, mpfr_cmp(r26041, r26042) <= 0, MPFR_RNDN); ; mpfr_set_si(r26045, mpfr_cmp(r26041, r26044) <= 0, MPFR_RNDN); mpfr_set_si(r26046, !mpfr_get_si(r26045, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r26047, mpfr_get_si(r26043, MPFR_RNDN) || mpfr_get_si(r26046, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r26050, r26041, r26049, MPFR_RNDN); mpfr_div(r26051, r26048, r26050, MPFR_RNDN); ; ; mpfr_pow(r26054, r26041, r26053, MPFR_RNDN); mpfr_div(r26055, r26052, r26054, MPFR_RNDN); mpfr_add(r26056, r26051, r26055, MPFR_RNDN); ; mpfr_div(r26058, r26057, r26041, MPFR_RNDN); mpfr_fma(r26059, r26056, r26041, r26058, MPFR_RNDN); mpfr_mul(r26060, r26041, r26041, MPFR_RNDN); mpfr_mul(r26061, r26060, r26060, MPFR_RNDN); mpfr_mul(r26062, r26061, r26060, MPFR_RNDN); mpfr_mul(r26063, r26060, r26062, MPFR_RNDN); mpfr_mul(r26064, r26060, r26063, MPFR_RNDN); ; mpfr_mul(r26066, r26064, r26065, MPFR_RNDN); ; mpfr_mul(r26068, r26067, r26062, MPFR_RNDN); ; mpfr_mul(r26070, r26060, r26069, MPFR_RNDN); ; mpfr_add(r26072, r26070, r26071, MPFR_RNDN); ; mpfr_mul(r26074, r26061, r26073, MPFR_RNDN); mpfr_add(r26075, r26072, r26074, MPFR_RNDN); mpfr_add(r26076, r26068, r26075, MPFR_RNDN); ; mpfr_mul(r26078, r26077, r26063, MPFR_RNDN); mpfr_add(r26079, r26076, r26078, MPFR_RNDN); mpfr_add(r26080, r26066, r26079, MPFR_RNDN); ; mpfr_mul(r26082, r26064, r26081, MPFR_RNDN); ; mpfr_mul(r26084, r26083, r26060, MPFR_RNDN); mpfr_add(r26085, r26084, r26071, MPFR_RNDN); ; mpfr_mul(r26087, r26061, r26086, MPFR_RNDN); mpfr_add(r26088, r26085, r26087, MPFR_RNDN); ; mpfr_mul(r26090, r26062, r26089, MPFR_RNDN); mpfr_add(r26091, r26088, r26090, MPFR_RNDN); ; mpfr_mul(r26093, r26092, r26063, MPFR_RNDN); mpfr_add(r26094, r26091, r26093, MPFR_RNDN); mpfr_add(r26095, r26082, r26094, MPFR_RNDN); ; mpfr_mul(r26097, r26096, r26065, MPFR_RNDN); mpfr_mul(r26098, r26060, r26064, MPFR_RNDN); mpfr_mul(r26099, r26097, r26098, MPFR_RNDN); mpfr_add(r26100, r26095, r26099, MPFR_RNDN); mpfr_div(r26101, r26080, r26100, MPFR_RNDN); mpfr_mul(r26102, r26041, r26101, MPFR_RNDN); if (mpfr_get_si(r26047, MPFR_RNDN)) { mpfr_set(r26103, r26059, MPFR_RNDN); } else { mpfr_set(r26103, r26102, MPFR_RNDN); }; return mpfr_get_d(r26103, MPFR_RNDN); }