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