#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18743 = 1.0f; float r18744 = 0.1049934947f; float r18745 = x; float r18746 = r18745 * r18745; float r18747 = r18744 * r18746; float r18748 = r18743 + r18747; float r18749 = 0.0424060604f; float r18750 = r18746 * r18746; float r18751 = r18749 * r18750; float r18752 = r18748 + r18751; float r18753 = 0.0072644182f; float r18754 = r18750 * r18746; float r18755 = r18753 * r18754; float r18756 = r18752 + r18755; float r18757 = 0.0005064034f; float r18758 = r18754 * r18746; float r18759 = r18757 * r18758; float r18760 = r18756 + r18759; float r18761 = 0.0001789971f; float r18762 = r18758 * r18746; float r18763 = r18761 * r18762; float r18764 = r18760 + r18763; float r18765 = 0.7715471019f; float r18766 = r18765 * r18746; float r18767 = r18743 + r18766; float r18768 = 0.2909738639f; float r18769 = r18768 * r18750; float r18770 = r18767 + r18769; float r18771 = 0.0694555761f; float r18772 = r18771 * r18754; float r18773 = r18770 + r18772; float r18774 = 0.0140005442f; float r18775 = r18774 * r18758; float r18776 = r18773 + r18775; float r18777 = 0.0008327945f; float r18778 = r18777 * r18762; float r18779 = r18776 + r18778; float r18780 = 2.0f; float r18781 = r18780 * r18761; float r18782 = r18762 * r18746; float r18783 = r18781 * r18782; float r18784 = r18779 + r18783; float r18785 = r18764 / r18784; float r18786 = r18785 * r18745; return r18786; } double f_id(double x) { double r18787 = 1.0; double r18788 = 0.1049934947; double r18789 = x; double r18790 = r18789 * r18789; double r18791 = r18788 * r18790; double r18792 = r18787 + r18791; double r18793 = 0.0424060604; double r18794 = r18790 * r18790; double r18795 = r18793 * r18794; double r18796 = r18792 + r18795; double r18797 = 0.0072644182; double r18798 = r18794 * r18790; double r18799 = r18797 * r18798; double r18800 = r18796 + r18799; double r18801 = 0.0005064034; double r18802 = r18798 * r18790; double r18803 = r18801 * r18802; double r18804 = r18800 + r18803; double r18805 = 0.0001789971; double r18806 = r18802 * r18790; double r18807 = r18805 * r18806; double r18808 = r18804 + r18807; double r18809 = 0.7715471019; double r18810 = r18809 * r18790; double r18811 = r18787 + r18810; double r18812 = 0.2909738639; double r18813 = r18812 * r18794; double r18814 = r18811 + r18813; double r18815 = 0.0694555761; double r18816 = r18815 * r18798; double r18817 = r18814 + r18816; double r18818 = 0.0140005442; double r18819 = r18818 * r18802; double r18820 = r18817 + r18819; double r18821 = 0.0008327945; double r18822 = r18821 * r18806; double r18823 = r18820 + r18822; double r18824 = 2.0; double r18825 = r18824 * r18805; double r18826 = r18806 * r18790; double r18827 = r18825 * r18826; double r18828 = r18823 + r18827; double r18829 = r18808 / r18828; double r18830 = r18829 * r18789; return r18830; } double f_of(float x) { float r18831 = x; float r18832 = -1.1392497874101738e+24f; bool r18833 = r18831 <= r18832; float r18834 = 1.0f; float r18835 = r18834 / r18831; float r18836 = r18835 * (r18835 * r18835); float r18837 = r18836 * (r18836 * r18836); float r18838 = r18831 * r18831; float r18839 = 0.0003579942f; float r18840 = r18834 / r18838; float r18841 = r18840 * (r18840 * r18840); float r18842 = 6.0f; float r18843 = pow(r18831, r18842); float r18844 = r18841 / r18843; float r18845 = 8.0f; float r18846 = pow(r18835, r18845); float r18847 = 0.0140005442f; float r18848 = 0.0694555761f; float r18849 = r18848 * r18841; float r18850 = fma(r18846, r18847, r18849); float r18851 = 0.0008327945f; float r18852 = 10.0f; float r18853 = pow(r18831, r18852); float r18854 = r18834 / r18853; float r18855 = 0.2909738639f; float r18856 = 4.0f; float r18857 = pow(r18835, r18856); float r18858 = 0.7715471019f; float r18859 = r18858 / r18831; float r18860 = fma(r18859, r18835, r18834); float r18861 = fma(r18855, r18857, r18860); float r18862 = fma(r18851, r18854, r18861); float r18863 = r18850 + r18862; float r18864 = fma(r18839, r18844, r18863); float r18865 = r18838 * r18864; float r18866 = r18837 / r18865; float r18867 = 0.0001789971f; float r18868 = 0.0005064034f; float r18869 = 0.0072644182f; float r18870 = r18869 * r18841; float r18871 = fma(r18868, r18846, r18870); float r18872 = 0.0424060604f; float r18873 = 0.1049934947f; float r18874 = r18873 / r18831; float r18875 = fma(r18874, r18835, r18834); float r18876 = fma(r18872, r18857, r18875); float r18877 = r18871 + r18876; float r18878 = r18831 * r18864; float r18879 = r18877 / r18878; float r18880 = fma(r18866, r18867, r18879); float r18881 = 1.653196896254347e+29f; bool r18882 = r18831 <= r18881; float r18883 = r18831 * (r18831 * r18831); float r18884 = r18883 * (r18883 * r18883); float r18885 = r18884 * r18838; float r18886 = r18883 * r18843; float r18887 = r18883 * r18886; float r18888 = r18838 * (r18838 * r18838); float r18889 = pow(r18831, r18856); float r18890 = r18888 * r18889; float r18891 = r18831 * r18858; float r18892 = fma(r18891, r18831, r18834); float r18893 = fma(r18855, r18889, r18892); float r18894 = fma(r18851, r18890, r18893); float r18895 = pow(r18831, r18845); float r18896 = r18888 * r18848; float r18897 = fma(r18895, r18847, r18896); float r18898 = r18894 + r18897; float r18899 = fma(r18839, r18887, r18898); float r18900 = r18885 / r18899; float r18901 = r18831 * r18831; float r18902 = r18901 * (r18901 * r18901); float r18903 = r18843 * r18902; float r18904 = r18901 * r18901; float r18905 = r18901 * r18889; float r18906 = r18904 * r18905; float r18907 = r18858 * r18831; float r18908 = fma(r18907, r18831, r18834); float r18909 = fma(r18855, r18889, r18908); float r18910 = fma(r18851, r18906, r18909); float r18911 = r18848 * r18902; float r18912 = fma(r18895, r18847, r18911); float r18913 = r18910 + r18912; float r18914 = fma(r18839, r18903, r18913); float r18915 = r18831 / r18914; float r18916 = r18869 * r18902; float r18917 = fma(r18868, r18895, r18916); float r18918 = r18873 * r18831; float r18919 = fma(r18918, r18831, r18834); float r18920 = fma(r18872, r18889, r18919); float r18921 = r18917 + r18920; float r18922 = r18915 * r18921; float r18923 = fma(r18900, r18867, r18922); float r18924 = r18882 ? r18923 : r18880; float r18925 = r18833 ? r18880 : r18924; return r18925; } double f_od(double x) { double r18926 = x; double r18927 = -1.1392497874101738e+24; bool r18928 = r18926 <= r18927; double r18929 = 1.0; double r18930 = r18929 / r18926; double r18931 = r18930 * (r18930 * r18930); double r18932 = r18931 * (r18931 * r18931); double r18933 = r18926 * r18926; double r18934 = 0.0003579942; double r18935 = r18929 / r18933; double r18936 = r18935 * (r18935 * r18935); double r18937 = 6.0; double r18938 = pow(r18926, r18937); double r18939 = r18936 / r18938; double r18940 = 8.0; double r18941 = pow(r18930, r18940); double r18942 = 0.0140005442; double r18943 = 0.0694555761; double r18944 = r18943 * r18936; double r18945 = fma(r18941, r18942, r18944); double r18946 = 0.0008327945; double r18947 = 10.0; double r18948 = pow(r18926, r18947); double r18949 = r18929 / r18948; double r18950 = 0.2909738639; double r18951 = 4.0; double r18952 = pow(r18930, r18951); double r18953 = 0.7715471019; double r18954 = r18953 / r18926; double r18955 = fma(r18954, r18930, r18929); double r18956 = fma(r18950, r18952, r18955); double r18957 = fma(r18946, r18949, r18956); double r18958 = r18945 + r18957; double r18959 = fma(r18934, r18939, r18958); double r18960 = r18933 * r18959; double r18961 = r18932 / r18960; double r18962 = 0.0001789971; double r18963 = 0.0005064034; double r18964 = 0.0072644182; double r18965 = r18964 * r18936; double r18966 = fma(r18963, r18941, r18965); double r18967 = 0.0424060604; double r18968 = 0.1049934947; double r18969 = r18968 / r18926; double r18970 = fma(r18969, r18930, r18929); double r18971 = fma(r18967, r18952, r18970); double r18972 = r18966 + r18971; double r18973 = r18926 * r18959; double r18974 = r18972 / r18973; double r18975 = fma(r18961, r18962, r18974); double r18976 = 1.653196896254347e+29; bool r18977 = r18926 <= r18976; double r18978 = r18926 * (r18926 * r18926); double r18979 = r18978 * (r18978 * r18978); double r18980 = r18979 * r18933; double r18981 = r18978 * r18938; double r18982 = r18978 * r18981; double r18983 = r18933 * (r18933 * r18933); double r18984 = pow(r18926, r18951); double r18985 = r18983 * r18984; double r18986 = r18926 * r18953; double r18987 = fma(r18986, r18926, r18929); double r18988 = fma(r18950, r18984, r18987); double r18989 = fma(r18946, r18985, r18988); double r18990 = pow(r18926, r18940); double r18991 = r18983 * r18943; double r18992 = fma(r18990, r18942, r18991); double r18993 = r18989 + r18992; double r18994 = fma(r18934, r18982, r18993); double r18995 = r18980 / r18994; double r18996 = r18926 * r18926; double r18997 = r18996 * (r18996 * r18996); double r18998 = r18938 * r18997; double r18999 = r18996 * r18996; double r19000 = r18996 * r18984; double r19001 = r18999 * r19000; double r19002 = r18953 * r18926; double r19003 = fma(r19002, r18926, r18929); double r19004 = fma(r18950, r18984, r19003); double r19005 = fma(r18946, r19001, r19004); double r19006 = r18943 * r18997; double r19007 = fma(r18990, r18942, r19006); double r19008 = r19005 + r19007; double r19009 = fma(r18934, r18998, r19008); double r19010 = r18926 / r19009; double r19011 = r18964 * r18997; double r19012 = fma(r18963, r18990, r19011); double r19013 = r18968 * r18926; double r19014 = fma(r19013, r18926, r18929); double r19015 = fma(r18967, r18984, r19014); double r19016 = r19012 + r19015; double r19017 = r19010 * r19016; double r19018 = fma(r18995, r18962, r19017); double r19019 = r18977 ? r19018 : r18975; double r19020 = r18928 ? r18975 : r19019; return r19020; } 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 r19021, r19022, r19023, r19024, r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033, r19034, r19035, r19036, r19037, r19038, r19039, r19040, r19041, r19042, r19043, r19044, r19045, r19046, r19047, r19048, r19049, r19050, r19051, r19052, r19053, r19054, r19055, r19056, r19057, r19058, r19059, r19060, r19061, r19062, r19063, r19064; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19021, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19022, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19023); mpfr_init(r19024); mpfr_init(r19025); mpfr_init(r19026); mpfr_init_set_str(r19027, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19028); mpfr_init(r19029); mpfr_init(r19030); mpfr_init_set_str(r19031, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19032); mpfr_init(r19033); mpfr_init(r19034); mpfr_init_set_str(r19035, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19036); mpfr_init(r19037); mpfr_init(r19038); mpfr_init_set_str(r19039, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19040); mpfr_init(r19041); mpfr_init(r19042); mpfr_init_set_str(r19043, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19044); mpfr_init(r19045); mpfr_init_set_str(r19046, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19047); mpfr_init(r19048); mpfr_init_set_str(r19049, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19050); mpfr_init(r19051); mpfr_init_set_str(r19052, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19053); mpfr_init(r19054); mpfr_init_set_str(r19055, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19056); mpfr_init(r19057); mpfr_init_set_str(r19058, "2", 10, MPFR_RNDN); mpfr_init(r19059); mpfr_init(r19060); mpfr_init(r19061); mpfr_init(r19062); mpfr_init(r19063); mpfr_init(r19064); } double f_im(double x) { ; ; mpfr_set_d(r19023, x, MPFR_RNDN); mpfr_mul(r19024, r19023, r19023, MPFR_RNDN); mpfr_mul(r19025, r19022, r19024, MPFR_RNDN); mpfr_add(r19026, r19021, r19025, MPFR_RNDN); ; mpfr_mul(r19028, r19024, r19024, MPFR_RNDN); mpfr_mul(r19029, r19027, r19028, MPFR_RNDN); mpfr_add(r19030, r19026, r19029, MPFR_RNDN); ; mpfr_mul(r19032, r19028, r19024, MPFR_RNDN); mpfr_mul(r19033, r19031, r19032, MPFR_RNDN); mpfr_add(r19034, r19030, r19033, MPFR_RNDN); ; mpfr_mul(r19036, r19032, r19024, MPFR_RNDN); mpfr_mul(r19037, r19035, r19036, MPFR_RNDN); mpfr_add(r19038, r19034, r19037, MPFR_RNDN); ; mpfr_mul(r19040, r19036, r19024, MPFR_RNDN); mpfr_mul(r19041, r19039, r19040, MPFR_RNDN); mpfr_add(r19042, r19038, r19041, MPFR_RNDN); ; mpfr_mul(r19044, r19043, r19024, MPFR_RNDN); mpfr_add(r19045, r19021, r19044, MPFR_RNDN); ; mpfr_mul(r19047, r19046, r19028, MPFR_RNDN); mpfr_add(r19048, r19045, r19047, MPFR_RNDN); ; mpfr_mul(r19050, r19049, r19032, MPFR_RNDN); mpfr_add(r19051, r19048, r19050, MPFR_RNDN); ; mpfr_mul(r19053, r19052, r19036, MPFR_RNDN); mpfr_add(r19054, r19051, r19053, MPFR_RNDN); ; mpfr_mul(r19056, r19055, r19040, MPFR_RNDN); mpfr_add(r19057, r19054, r19056, MPFR_RNDN); ; mpfr_mul(r19059, r19058, r19039, MPFR_RNDN); mpfr_mul(r19060, r19040, r19024, MPFR_RNDN); mpfr_mul(r19061, r19059, r19060, MPFR_RNDN); mpfr_add(r19062, r19057, r19061, MPFR_RNDN); mpfr_div(r19063, r19042, r19062, MPFR_RNDN); mpfr_mul(r19064, r19063, r19023, MPFR_RNDN); return mpfr_get_d(r19064, MPFR_RNDN); } static mpfr_t r19065, r19066, r19067, r19068, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086, r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117, r19118, r19119, r19120, r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138, r19139, r19140, r19141, r19142, r19143, r19144, r19145, r19146, r19147, r19148, r19149, r19150, r19151, r19152, r19153, r19154, r19155, r19156, r19157, r19158, r19159; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19065); mpfr_init_set_str(r19066, "-1.1392497874101738e+24", 10, MPFR_RNDN); mpfr_init(r19067); mpfr_init_set_str(r19068, "1", 10, MPFR_RNDN); mpfr_init(r19069); mpfr_init(r19070); mpfr_init(r19071); mpfr_init(r19072); mpfr_init_set_str(r19073, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19074); mpfr_init(r19075); mpfr_init_set_str(r19076, "6", 10, MPFR_RNDN); mpfr_init(r19077); mpfr_init(r19078); mpfr_init_set_str(r19079, "8", 10, MPFR_RNDN); mpfr_init(r19080); mpfr_init_set_str(r19081, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r19082, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19083); mpfr_init(r19084); mpfr_init_set_str(r19085, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r19086, "10", 10, MPFR_RNDN); mpfr_init(r19087); mpfr_init(r19088); mpfr_init_set_str(r19089, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r19090, "4", 10, MPFR_RNDN); mpfr_init(r19091); mpfr_init_set_str(r19092, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19093); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); mpfr_init(r19097); mpfr_init(r19098); mpfr_init(r19099); mpfr_init(r19100); mpfr_init_set_str(r19101, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19102, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r19103, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19104); mpfr_init(r19105); mpfr_init_set_str(r19106, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r19107, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19108); mpfr_init(r19109); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init(r19113); mpfr_init(r19114); mpfr_init_set_str(r19115, "1.653196896254347e+29", 10, MPFR_RNDN); mpfr_init(r19116); mpfr_init(r19117); mpfr_init(r19118); mpfr_init(r19119); mpfr_init(r19120); mpfr_init(r19121); mpfr_init(r19122); mpfr_init(r19123); mpfr_init(r19124); mpfr_init(r19125); mpfr_init(r19126); mpfr_init(r19127); mpfr_init(r19128); mpfr_init(r19129); mpfr_init(r19130); mpfr_init(r19131); mpfr_init(r19132); mpfr_init(r19133); mpfr_init(r19134); mpfr_init(r19135); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init(r19142); mpfr_init(r19143); mpfr_init(r19144); mpfr_init(r19145); mpfr_init(r19146); mpfr_init(r19147); mpfr_init(r19148); mpfr_init(r19149); mpfr_init(r19150); mpfr_init(r19151); mpfr_init(r19152); mpfr_init(r19153); mpfr_init(r19154); mpfr_init(r19155); mpfr_init(r19156); mpfr_init(r19157); mpfr_init(r19158); mpfr_init(r19159); } double f_fm(double x) { mpfr_set_d(r19065, x, MPFR_RNDN); ; mpfr_set_si(r19067, mpfr_cmp(r19065, r19066) <= 0, MPFR_RNDN); ; mpfr_div(r19069, r19068, r19065, MPFR_RNDN); mpfr_mul(r19070, r19069, r19069, MPFR_RNDN); mpfr_mul(r19070, r19070, r19069, MPFR_RNDN); mpfr_mul(r19071, r19070, r19070, MPFR_RNDN); mpfr_mul(r19071, r19071, r19070, MPFR_RNDN); mpfr_sqr(r19072, r19065, MPFR_RNDN); ; mpfr_div(r19074, r19068, r19072, MPFR_RNDN); mpfr_mul(r19075, r19074, r19074, MPFR_RNDN); mpfr_mul(r19075, r19075, r19074, MPFR_RNDN); ; mpfr_pow(r19077, r19065, r19076, MPFR_RNDN); mpfr_div(r19078, r19075, r19077, MPFR_RNDN); ; mpfr_pow(r19080, r19069, r19079, MPFR_RNDN); ; ; mpfr_mul(r19083, r19082, r19075, MPFR_RNDN); mpfr_fma(r19084, r19080, r19081, r19083, MPFR_RNDN); ; ; mpfr_pow(r19087, r19065, r19086, MPFR_RNDN); mpfr_div(r19088, r19068, r19087, MPFR_RNDN); ; ; mpfr_pow(r19091, r19069, r19090, MPFR_RNDN); ; mpfr_div(r19093, r19092, r19065, MPFR_RNDN); mpfr_fma(r19094, r19093, r19069, r19068, MPFR_RNDN); mpfr_fma(r19095, r19089, r19091, r19094, MPFR_RNDN); mpfr_fma(r19096, r19085, r19088, r19095, MPFR_RNDN); mpfr_add(r19097, r19084, r19096, MPFR_RNDN); mpfr_fma(r19098, r19073, r19078, r19097, MPFR_RNDN); mpfr_mul(r19099, r19072, r19098, MPFR_RNDN); mpfr_div(r19100, r19071, r19099, MPFR_RNDN); ; ; ; mpfr_mul(r19104, r19103, r19075, MPFR_RNDN); mpfr_fma(r19105, r19102, r19080, r19104, MPFR_RNDN); ; ; mpfr_div(r19108, r19107, r19065, MPFR_RNDN); mpfr_fma(r19109, r19108, r19069, r19068, MPFR_RNDN); mpfr_fma(r19110, r19106, r19091, r19109, MPFR_RNDN); mpfr_add(r19111, r19105, r19110, MPFR_RNDN); mpfr_mul(r19112, r19065, r19098, MPFR_RNDN); mpfr_div(r19113, r19111, r19112, MPFR_RNDN); mpfr_fma(r19114, r19100, r19101, r19113, MPFR_RNDN); ; mpfr_set_si(r19116, mpfr_cmp(r19065, r19115) <= 0, MPFR_RNDN); mpfr_mul(r19117, r19065, r19065, MPFR_RNDN); mpfr_mul(r19117, r19117, r19065, MPFR_RNDN); mpfr_mul(r19118, r19117, r19117, MPFR_RNDN); mpfr_mul(r19118, r19118, r19117, MPFR_RNDN); mpfr_mul(r19119, r19118, r19072, MPFR_RNDN); mpfr_mul(r19120, r19117, r19077, MPFR_RNDN); mpfr_mul(r19121, r19117, r19120, MPFR_RNDN); mpfr_mul(r19122, r19072, r19072, MPFR_RNDN); mpfr_mul(r19122, r19122, r19072, MPFR_RNDN); mpfr_pow(r19123, r19065, r19090, MPFR_RNDN); mpfr_mul(r19124, r19122, r19123, MPFR_RNDN); mpfr_mul(r19125, r19065, r19092, MPFR_RNDN); mpfr_fma(r19126, r19125, r19065, r19068, MPFR_RNDN); mpfr_fma(r19127, r19089, r19123, r19126, MPFR_RNDN); mpfr_fma(r19128, r19085, r19124, r19127, MPFR_RNDN); mpfr_pow(r19129, r19065, r19079, MPFR_RNDN); mpfr_mul(r19130, r19122, r19082, MPFR_RNDN); mpfr_fma(r19131, r19129, r19081, r19130, MPFR_RNDN); mpfr_add(r19132, r19128, r19131, MPFR_RNDN); mpfr_fma(r19133, r19073, r19121, r19132, MPFR_RNDN); mpfr_div(r19134, r19119, r19133, MPFR_RNDN); mpfr_mul(r19135, r19065, r19065, MPFR_RNDN); mpfr_mul(r19136, r19135, r19135, MPFR_RNDN); mpfr_mul(r19136, r19136, r19135, MPFR_RNDN); mpfr_mul(r19137, r19077, r19136, MPFR_RNDN); mpfr_mul(r19138, r19135, r19135, MPFR_RNDN); mpfr_mul(r19139, r19135, r19123, MPFR_RNDN); mpfr_mul(r19140, r19138, r19139, MPFR_RNDN); mpfr_mul(r19141, r19092, r19065, MPFR_RNDN); mpfr_fma(r19142, r19141, r19065, r19068, MPFR_RNDN); mpfr_fma(r19143, r19089, r19123, r19142, MPFR_RNDN); mpfr_fma(r19144, r19085, r19140, r19143, MPFR_RNDN); mpfr_mul(r19145, r19082, r19136, MPFR_RNDN); mpfr_fma(r19146, r19129, r19081, r19145, MPFR_RNDN); mpfr_add(r19147, r19144, r19146, MPFR_RNDN); mpfr_fma(r19148, r19073, r19137, r19147, MPFR_RNDN); mpfr_div(r19149, r19065, r19148, MPFR_RNDN); mpfr_mul(r19150, r19103, r19136, MPFR_RNDN); mpfr_fma(r19151, r19102, r19129, r19150, MPFR_RNDN); mpfr_mul(r19152, r19107, r19065, MPFR_RNDN); mpfr_fma(r19153, r19152, r19065, r19068, MPFR_RNDN); mpfr_fma(r19154, r19106, r19123, r19153, MPFR_RNDN); mpfr_add(r19155, r19151, r19154, MPFR_RNDN); mpfr_mul(r19156, r19149, r19155, MPFR_RNDN); mpfr_fma(r19157, r19134, r19101, r19156, MPFR_RNDN); if (mpfr_get_si(r19116, MPFR_RNDN)) { mpfr_set(r19158, r19157, MPFR_RNDN); } else { mpfr_set(r19158, r19114, MPFR_RNDN); }; if (mpfr_get_si(r19067, MPFR_RNDN)) { mpfr_set(r19159, r19114, MPFR_RNDN); } else { mpfr_set(r19159, r19158, MPFR_RNDN); }; return mpfr_get_d(r19159, MPFR_RNDN); } static mpfr_t r19160, r19161, r19162, r19163, r19164, r19165, r19166, r19167, r19168, r19169, r19170, r19171, r19172, r19173, r19174, r19175, r19176, r19177, r19178, r19179, r19180, r19181, r19182, r19183, r19184, r19185, r19186, r19187, r19188, r19189, r19190, r19191, r19192, r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201, r19202, r19203, r19204, r19205, r19206, r19207, r19208, r19209, r19210, r19211, r19212, r19213, r19214, r19215, r19216, r19217, r19218, r19219, r19220, r19221, r19222, r19223, r19224, r19225, r19226, r19227, r19228, r19229, r19230, r19231, r19232, r19233, r19234, r19235, r19236, r19237, r19238, r19239, r19240, r19241, r19242, r19243, r19244, r19245, r19246, r19247, r19248, r19249, r19250, r19251, r19252, r19253, r19254; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19160); mpfr_init_set_str(r19161, "-1.1392497874101738e+24", 10, MPFR_RNDN); mpfr_init(r19162); mpfr_init_set_str(r19163, "1", 10, MPFR_RNDN); mpfr_init(r19164); mpfr_init(r19165); mpfr_init(r19166); mpfr_init(r19167); mpfr_init_set_str(r19168, "0.0003579942", 10, MPFR_RNDN); mpfr_init(r19169); mpfr_init(r19170); mpfr_init_set_str(r19171, "6", 10, MPFR_RNDN); mpfr_init(r19172); mpfr_init(r19173); mpfr_init_set_str(r19174, "8", 10, MPFR_RNDN); mpfr_init(r19175); mpfr_init_set_str(r19176, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r19177, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19178); mpfr_init(r19179); mpfr_init_set_str(r19180, "0.0008327945", 10, MPFR_RNDN); mpfr_init_set_str(r19181, "10", 10, MPFR_RNDN); mpfr_init(r19182); mpfr_init(r19183); mpfr_init_set_str(r19184, "0.2909738639", 10, MPFR_RNDN); mpfr_init_set_str(r19185, "4", 10, MPFR_RNDN); mpfr_init(r19186); mpfr_init_set_str(r19187, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19188); mpfr_init(r19189); mpfr_init(r19190); mpfr_init(r19191); mpfr_init(r19192); mpfr_init(r19193); mpfr_init(r19194); mpfr_init(r19195); mpfr_init_set_str(r19196, "0.0001789971", 10, MPFR_RNDN); mpfr_init_set_str(r19197, "0.0005064034", 10, MPFR_RNDN); mpfr_init_set_str(r19198, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19199); mpfr_init(r19200); mpfr_init_set_str(r19201, "0.0424060604", 10, MPFR_RNDN); mpfr_init_set_str(r19202, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19203); mpfr_init(r19204); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); mpfr_init(r19208); mpfr_init(r19209); mpfr_init_set_str(r19210, "1.653196896254347e+29", 10, MPFR_RNDN); mpfr_init(r19211); mpfr_init(r19212); mpfr_init(r19213); mpfr_init(r19214); mpfr_init(r19215); mpfr_init(r19216); mpfr_init(r19217); mpfr_init(r19218); mpfr_init(r19219); mpfr_init(r19220); mpfr_init(r19221); mpfr_init(r19222); mpfr_init(r19223); mpfr_init(r19224); mpfr_init(r19225); mpfr_init(r19226); mpfr_init(r19227); mpfr_init(r19228); mpfr_init(r19229); mpfr_init(r19230); mpfr_init(r19231); mpfr_init(r19232); mpfr_init(r19233); mpfr_init(r19234); mpfr_init(r19235); mpfr_init(r19236); mpfr_init(r19237); mpfr_init(r19238); mpfr_init(r19239); mpfr_init(r19240); mpfr_init(r19241); mpfr_init(r19242); mpfr_init(r19243); mpfr_init(r19244); mpfr_init(r19245); mpfr_init(r19246); mpfr_init(r19247); mpfr_init(r19248); mpfr_init(r19249); mpfr_init(r19250); mpfr_init(r19251); mpfr_init(r19252); mpfr_init(r19253); mpfr_init(r19254); } double f_dm(double x) { mpfr_set_d(r19160, x, MPFR_RNDN); ; mpfr_set_si(r19162, mpfr_cmp(r19160, r19161) <= 0, MPFR_RNDN); ; mpfr_div(r19164, r19163, r19160, MPFR_RNDN); mpfr_mul(r19165, r19164, r19164, MPFR_RNDN); mpfr_mul(r19165, r19165, r19164, MPFR_RNDN); mpfr_mul(r19166, r19165, r19165, MPFR_RNDN); mpfr_mul(r19166, r19166, r19165, MPFR_RNDN); mpfr_sqr(r19167, r19160, MPFR_RNDN); ; mpfr_div(r19169, r19163, r19167, MPFR_RNDN); mpfr_mul(r19170, r19169, r19169, MPFR_RNDN); mpfr_mul(r19170, r19170, r19169, MPFR_RNDN); ; mpfr_pow(r19172, r19160, r19171, MPFR_RNDN); mpfr_div(r19173, r19170, r19172, MPFR_RNDN); ; mpfr_pow(r19175, r19164, r19174, MPFR_RNDN); ; ; mpfr_mul(r19178, r19177, r19170, MPFR_RNDN); mpfr_fma(r19179, r19175, r19176, r19178, MPFR_RNDN); ; ; mpfr_pow(r19182, r19160, r19181, MPFR_RNDN); mpfr_div(r19183, r19163, r19182, MPFR_RNDN); ; ; mpfr_pow(r19186, r19164, r19185, MPFR_RNDN); ; mpfr_div(r19188, r19187, r19160, MPFR_RNDN); mpfr_fma(r19189, r19188, r19164, r19163, MPFR_RNDN); mpfr_fma(r19190, r19184, r19186, r19189, MPFR_RNDN); mpfr_fma(r19191, r19180, r19183, r19190, MPFR_RNDN); mpfr_add(r19192, r19179, r19191, MPFR_RNDN); mpfr_fma(r19193, r19168, r19173, r19192, MPFR_RNDN); mpfr_mul(r19194, r19167, r19193, MPFR_RNDN); mpfr_div(r19195, r19166, r19194, MPFR_RNDN); ; ; ; mpfr_mul(r19199, r19198, r19170, MPFR_RNDN); mpfr_fma(r19200, r19197, r19175, r19199, MPFR_RNDN); ; ; mpfr_div(r19203, r19202, r19160, MPFR_RNDN); mpfr_fma(r19204, r19203, r19164, r19163, MPFR_RNDN); mpfr_fma(r19205, r19201, r19186, r19204, MPFR_RNDN); mpfr_add(r19206, r19200, r19205, MPFR_RNDN); mpfr_mul(r19207, r19160, r19193, MPFR_RNDN); mpfr_div(r19208, r19206, r19207, MPFR_RNDN); mpfr_fma(r19209, r19195, r19196, r19208, MPFR_RNDN); ; mpfr_set_si(r19211, mpfr_cmp(r19160, r19210) <= 0, MPFR_RNDN); mpfr_mul(r19212, r19160, r19160, MPFR_RNDN); mpfr_mul(r19212, r19212, r19160, MPFR_RNDN); mpfr_mul(r19213, r19212, r19212, MPFR_RNDN); mpfr_mul(r19213, r19213, r19212, MPFR_RNDN); mpfr_mul(r19214, r19213, r19167, MPFR_RNDN); mpfr_mul(r19215, r19212, r19172, MPFR_RNDN); mpfr_mul(r19216, r19212, r19215, MPFR_RNDN); mpfr_mul(r19217, r19167, r19167, MPFR_RNDN); mpfr_mul(r19217, r19217, r19167, MPFR_RNDN); mpfr_pow(r19218, r19160, r19185, MPFR_RNDN); mpfr_mul(r19219, r19217, r19218, MPFR_RNDN); mpfr_mul(r19220, r19160, r19187, MPFR_RNDN); mpfr_fma(r19221, r19220, r19160, r19163, MPFR_RNDN); mpfr_fma(r19222, r19184, r19218, r19221, MPFR_RNDN); mpfr_fma(r19223, r19180, r19219, r19222, MPFR_RNDN); mpfr_pow(r19224, r19160, r19174, MPFR_RNDN); mpfr_mul(r19225, r19217, r19177, MPFR_RNDN); mpfr_fma(r19226, r19224, r19176, r19225, MPFR_RNDN); mpfr_add(r19227, r19223, r19226, MPFR_RNDN); mpfr_fma(r19228, r19168, r19216, r19227, MPFR_RNDN); mpfr_div(r19229, r19214, r19228, MPFR_RNDN); mpfr_mul(r19230, r19160, r19160, MPFR_RNDN); mpfr_mul(r19231, r19230, r19230, MPFR_RNDN); mpfr_mul(r19231, r19231, r19230, MPFR_RNDN); mpfr_mul(r19232, r19172, r19231, MPFR_RNDN); mpfr_mul(r19233, r19230, r19230, MPFR_RNDN); mpfr_mul(r19234, r19230, r19218, MPFR_RNDN); mpfr_mul(r19235, r19233, r19234, MPFR_RNDN); mpfr_mul(r19236, r19187, r19160, MPFR_RNDN); mpfr_fma(r19237, r19236, r19160, r19163, MPFR_RNDN); mpfr_fma(r19238, r19184, r19218, r19237, MPFR_RNDN); mpfr_fma(r19239, r19180, r19235, r19238, MPFR_RNDN); mpfr_mul(r19240, r19177, r19231, MPFR_RNDN); mpfr_fma(r19241, r19224, r19176, r19240, MPFR_RNDN); mpfr_add(r19242, r19239, r19241, MPFR_RNDN); mpfr_fma(r19243, r19168, r19232, r19242, MPFR_RNDN); mpfr_div(r19244, r19160, r19243, MPFR_RNDN); mpfr_mul(r19245, r19198, r19231, MPFR_RNDN); mpfr_fma(r19246, r19197, r19224, r19245, MPFR_RNDN); mpfr_mul(r19247, r19202, r19160, MPFR_RNDN); mpfr_fma(r19248, r19247, r19160, r19163, MPFR_RNDN); mpfr_fma(r19249, r19201, r19218, r19248, MPFR_RNDN); mpfr_add(r19250, r19246, r19249, MPFR_RNDN); mpfr_mul(r19251, r19244, r19250, MPFR_RNDN); mpfr_fma(r19252, r19229, r19196, r19251, MPFR_RNDN); if (mpfr_get_si(r19211, MPFR_RNDN)) { mpfr_set(r19253, r19252, MPFR_RNDN); } else { mpfr_set(r19253, r19209, MPFR_RNDN); }; if (mpfr_get_si(r19162, MPFR_RNDN)) { mpfr_set(r19254, r19209, MPFR_RNDN); } else { mpfr_set(r19254, r19253, MPFR_RNDN); }; return mpfr_get_d(r19254, MPFR_RNDN); }