#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18812 = 1.0f; float r18813 = 0.10499349236488342f; float r18814 = x; float r18815 = r18814 * r18814; float r18816 = r18813 * r18815; float r18817 = r18812 + r18816; float r18818 = 0.04240605980157852f; float r18819 = r18815 * r18815; float r18820 = r18818 * r18819; float r18821 = r18817 + r18820; float r18822 = 0.007264418061822653f; float r18823 = r18819 * r18815; float r18824 = r18822 * r18823; float r18825 = r18821 + r18824; float r18826 = 0.0005064033903181553f; float r18827 = r18823 * r18815; float r18828 = r18826 * r18827; float r18829 = r18825 + r18828; float r18830 = 0.00017899709928315133f; float r18831 = r18827 * r18815; float r18832 = r18830 * r18831; float r18833 = r18829 + r18832; float r18834 = 0.7715470790863037f; float r18835 = r18834 * r18815; float r18836 = r18812 + r18835; float r18837 = 0.29097387194633484f; float r18838 = r18837 * r18819; float r18839 = r18836 + r18838; float r18840 = 0.0694555789232254f; float r18841 = r18840 * r18823; float r18842 = r18839 + r18841; float r18843 = 0.01400054432451725f; float r18844 = r18843 * r18827; float r18845 = r18842 + r18844; float r18846 = 0.0008327945251949131f; float r18847 = r18846 * r18831; float r18848 = r18845 + r18847; float r18849 = 2.0f; float r18850 = r18849 * r18830; float r18851 = r18831 * r18815; float r18852 = r18850 * r18851; float r18853 = r18848 + r18852; float r18854 = r18833 / r18853; float r18855 = r18854 * r18814; return r18855; } double f_id(double x) { double r18856 = 1.0; double r18857 = 0.10499349236488342; double r18858 = x; double r18859 = r18858 * r18858; double r18860 = r18857 * r18859; double r18861 = r18856 + r18860; double r18862 = 0.04240605980157852; double r18863 = r18859 * r18859; double r18864 = r18862 * r18863; double r18865 = r18861 + r18864; double r18866 = 0.007264418061822653; double r18867 = r18863 * r18859; double r18868 = r18866 * r18867; double r18869 = r18865 + r18868; double r18870 = 0.0005064033903181553; double r18871 = r18867 * r18859; double r18872 = r18870 * r18871; double r18873 = r18869 + r18872; double r18874 = 0.00017899709928315133; double r18875 = r18871 * r18859; double r18876 = r18874 * r18875; double r18877 = r18873 + r18876; double r18878 = 0.7715470790863037; double r18879 = r18878 * r18859; double r18880 = r18856 + r18879; double r18881 = 0.29097387194633484; double r18882 = r18881 * r18863; double r18883 = r18880 + r18882; double r18884 = 0.0694555789232254; double r18885 = r18884 * r18867; double r18886 = r18883 + r18885; double r18887 = 0.01400054432451725; double r18888 = r18887 * r18871; double r18889 = r18886 + r18888; double r18890 = 0.0008327945251949131; double r18891 = r18890 * r18875; double r18892 = r18889 + r18891; double r18893 = 2.0; double r18894 = r18893 * r18874; double r18895 = r18875 * r18859; double r18896 = r18894 * r18895; double r18897 = r18892 + r18896; double r18898 = r18877 / r18897; double r18899 = r18898 * r18858; return r18899; } double f_of(float x) { float r18900 = x; float r18901 = r18900 * r18900; float r18902 = r18901 * (r18901 * r18901); float r18903 = 0.0694555789232254f; float r18904 = 0.01400054432451725f; float r18905 = r18901 * r18904; float r18906 = r18903 + r18905; float r18907 = r18902 * r18906; float r18908 = 0.7715470790863037f; float r18909 = r18908 * r18900; float r18910 = r18909 * r18900; float r18911 = 1.0f; float r18912 = r18910 + r18911; float r18913 = 0.29097387194633484f; float r18914 = r18900 * r18913; float r18915 = r18900 * (r18900 * r18900); float r18916 = r18914 * r18915; float r18917 = r18912 + r18916; float r18918 = r18907 + r18917; float r18919 = 2.0f; float r18920 = r18900 * r18919; float r18921 = 0.00017899709928315133f; float r18922 = r18900 * r18921; float r18923 = r18920 * r18922; float r18924 = 0.0008327945251949131f; float r18925 = r18923 + r18924; float r18926 = r18901 * r18901; float r18927 = r18926 * r18902; float r18928 = r18925 * r18927; float r18929 = r18918 + r18928; float r18930 = 0.0005064033903181553f; float r18931 = r18901 * r18930; float r18932 = 0.007264418061822653f; float r18933 = r18931 + r18932; float r18934 = r18902 * r18933; float r18935 = r18926 * r18926; float r18936 = r18922 * r18900; float r18937 = r18935 * r18936; float r18938 = 0.10499349236488342f; float r18939 = r18900 * r18938; float r18940 = r18900 * r18939; float r18941 = 0.04240605980157852f; float r18942 = r18900 * r18941; float r18943 = r18915 * r18942; float r18944 = r18940 + r18943; float r18945 = r18944 + r18911; float r18946 = r18937 + r18945; float r18947 = r18934 + r18946; float r18948 = r18929 / r18947; float r18949 = r18900 / r18948; return r18949; } double f_od(double x) { double r18950 = x; double r18951 = r18950 * r18950; double r18952 = r18951 * (r18951 * r18951); double r18953 = 0.0694555789232254; double r18954 = 0.01400054432451725; double r18955 = r18951 * r18954; double r18956 = r18953 + r18955; double r18957 = r18952 * r18956; double r18958 = 0.7715470790863037; double r18959 = r18958 * r18950; double r18960 = r18959 * r18950; double r18961 = 1.0; double r18962 = r18960 + r18961; double r18963 = 0.29097387194633484; double r18964 = r18950 * r18963; double r18965 = r18950 * (r18950 * r18950); double r18966 = r18964 * r18965; double r18967 = r18962 + r18966; double r18968 = r18957 + r18967; double r18969 = 2.0; double r18970 = r18950 * r18969; double r18971 = 0.00017899709928315133; double r18972 = r18950 * r18971; double r18973 = r18970 * r18972; double r18974 = 0.0008327945251949131; double r18975 = r18973 + r18974; double r18976 = r18951 * r18951; double r18977 = r18976 * r18952; double r18978 = r18975 * r18977; double r18979 = r18968 + r18978; double r18980 = 0.0005064033903181553; double r18981 = r18951 * r18980; double r18982 = 0.007264418061822653; double r18983 = r18981 + r18982; double r18984 = r18952 * r18983; double r18985 = r18976 * r18976; double r18986 = r18972 * r18950; double r18987 = r18985 * r18986; double r18988 = 0.10499349236488342; double r18989 = r18950 * r18988; double r18990 = r18950 * r18989; double r18991 = 0.04240605980157852; double r18992 = r18950 * r18991; double r18993 = r18965 * r18992; double r18994 = r18990 + r18993; double r18995 = r18994 + r18961; double r18996 = r18987 + r18995; double r18997 = r18984 + r18996; double r18998 = r18979 / r18997; double r18999 = r18950 / r18998; return r18999; } 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 r19000, r19001, r19002, r19003, r19004, r19005, r19006, r19007, r19008, r19009, r19010, r19011, r19012, r19013, r19014, r19015, r19016, r19017, r19018, r19019, r19020, r19021, r19022, r19023, r19024, r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033, r19034, r19035, r19036, r19037, r19038, r19039, r19040, r19041, r19042, r19043; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19000, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19001, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19002); mpfr_init(r19003); mpfr_init(r19004); mpfr_init(r19005); mpfr_init_set_str(r19006, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19007); mpfr_init(r19008); mpfr_init(r19009); mpfr_init_set_str(r19010, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19011); mpfr_init(r19012); mpfr_init(r19013); mpfr_init_set_str(r19014, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19015); mpfr_init(r19016); mpfr_init(r19017); mpfr_init_set_str(r19018, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19019); mpfr_init(r19020); mpfr_init(r19021); mpfr_init_set_str(r19022, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19023); mpfr_init(r19024); mpfr_init_set_str(r19025, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19026); mpfr_init(r19027); mpfr_init_set_str(r19028, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19029); mpfr_init(r19030); mpfr_init_set_str(r19031, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19032); mpfr_init(r19033); mpfr_init_set_str(r19034, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19035); mpfr_init(r19036); mpfr_init_set_str(r19037, "2", 10, MPFR_RNDN); mpfr_init(r19038); mpfr_init(r19039); mpfr_init(r19040); mpfr_init(r19041); mpfr_init(r19042); mpfr_init(r19043); } double f_im(double x) { ; ; mpfr_set_d(r19002, x, MPFR_RNDN); mpfr_mul(r19003, r19002, r19002, MPFR_RNDN); mpfr_mul(r19004, r19001, r19003, MPFR_RNDN); mpfr_add(r19005, r19000, r19004, MPFR_RNDN); ; mpfr_mul(r19007, r19003, r19003, MPFR_RNDN); mpfr_mul(r19008, r19006, r19007, MPFR_RNDN); mpfr_add(r19009, r19005, r19008, MPFR_RNDN); ; mpfr_mul(r19011, r19007, r19003, MPFR_RNDN); mpfr_mul(r19012, r19010, r19011, MPFR_RNDN); mpfr_add(r19013, r19009, r19012, MPFR_RNDN); ; mpfr_mul(r19015, r19011, r19003, MPFR_RNDN); mpfr_mul(r19016, r19014, r19015, MPFR_RNDN); mpfr_add(r19017, r19013, r19016, MPFR_RNDN); ; mpfr_mul(r19019, r19015, r19003, MPFR_RNDN); mpfr_mul(r19020, r19018, r19019, MPFR_RNDN); mpfr_add(r19021, r19017, r19020, MPFR_RNDN); ; mpfr_mul(r19023, r19022, r19003, MPFR_RNDN); mpfr_add(r19024, r19000, r19023, MPFR_RNDN); ; mpfr_mul(r19026, r19025, r19007, MPFR_RNDN); mpfr_add(r19027, r19024, r19026, MPFR_RNDN); ; mpfr_mul(r19029, r19028, r19011, MPFR_RNDN); mpfr_add(r19030, r19027, r19029, MPFR_RNDN); ; mpfr_mul(r19032, r19031, r19015, MPFR_RNDN); mpfr_add(r19033, r19030, r19032, MPFR_RNDN); ; mpfr_mul(r19035, r19034, r19019, MPFR_RNDN); mpfr_add(r19036, r19033, r19035, MPFR_RNDN); ; mpfr_mul(r19038, r19037, r19018, MPFR_RNDN); mpfr_mul(r19039, r19019, r19003, MPFR_RNDN); mpfr_mul(r19040, r19038, r19039, MPFR_RNDN); mpfr_add(r19041, r19036, r19040, MPFR_RNDN); mpfr_div(r19042, r19021, r19041, MPFR_RNDN); mpfr_mul(r19043, r19042, r19002, MPFR_RNDN); return mpfr_get_d(r19043, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19044); mpfr_init(r19045); mpfr_init(r19046); mpfr_init_set_str(r19047, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19048, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19049); mpfr_init(r19050); mpfr_init(r19051); mpfr_init_set_str(r19052, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19053); mpfr_init(r19054); mpfr_init_set_str(r19055, "1", 10, MPFR_RNDN); mpfr_init(r19056); mpfr_init_set_str(r19057, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19058); mpfr_init(r19059); mpfr_init(r19060); mpfr_init(r19061); mpfr_init(r19062); mpfr_init_set_str(r19063, "2", 10, MPFR_RNDN); mpfr_init(r19064); mpfr_init_set_str(r19065, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19066); mpfr_init(r19067); mpfr_init_set_str(r19068, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19069); mpfr_init(r19070); mpfr_init(r19071); mpfr_init(r19072); mpfr_init(r19073); mpfr_init_set_str(r19074, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19075); mpfr_init_set_str(r19076, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19077); mpfr_init(r19078); mpfr_init(r19079); mpfr_init(r19080); mpfr_init(r19081); mpfr_init_set_str(r19082, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19083); mpfr_init(r19084); mpfr_init_set_str(r19085, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19086); mpfr_init(r19087); mpfr_init(r19088); mpfr_init(r19089); mpfr_init(r19090); mpfr_init(r19091); mpfr_init(r19092); mpfr_init(r19093); } double f_fm(double x) { mpfr_set_d(r19044, x, MPFR_RNDN); mpfr_sqr(r19045, r19044, MPFR_RNDN); mpfr_mul(r19046, r19045, r19045, MPFR_RNDN); mpfr_mul(r19046, r19046, r19045, MPFR_RNDN); ; ; mpfr_mul(r19049, r19045, r19048, MPFR_RNDN); mpfr_add(r19050, r19047, r19049, MPFR_RNDN); mpfr_mul(r19051, r19046, r19050, MPFR_RNDN); ; mpfr_mul(r19053, r19052, r19044, MPFR_RNDN); mpfr_mul(r19054, r19053, r19044, MPFR_RNDN); ; mpfr_add(r19056, r19054, r19055, MPFR_RNDN); ; mpfr_mul(r19058, r19044, r19057, MPFR_RNDN); mpfr_mul(r19059, r19044, r19044, MPFR_RNDN); mpfr_mul(r19059, r19059, r19044, MPFR_RNDN); mpfr_mul(r19060, r19058, r19059, MPFR_RNDN); mpfr_add(r19061, r19056, r19060, MPFR_RNDN); mpfr_add(r19062, r19051, r19061, MPFR_RNDN); ; mpfr_mul(r19064, r19044, r19063, MPFR_RNDN); ; mpfr_mul(r19066, r19044, r19065, MPFR_RNDN); mpfr_mul(r19067, r19064, r19066, MPFR_RNDN); ; mpfr_add(r19069, r19067, r19068, MPFR_RNDN); mpfr_mul(r19070, r19045, r19045, MPFR_RNDN); mpfr_mul(r19071, r19070, r19046, MPFR_RNDN); mpfr_mul(r19072, r19069, r19071, MPFR_RNDN); mpfr_add(r19073, r19062, r19072, MPFR_RNDN); ; mpfr_mul(r19075, r19045, r19074, MPFR_RNDN); ; mpfr_add(r19077, r19075, r19076, MPFR_RNDN); mpfr_mul(r19078, r19046, r19077, MPFR_RNDN); mpfr_mul(r19079, r19070, r19070, MPFR_RNDN); mpfr_mul(r19080, r19066, r19044, MPFR_RNDN); mpfr_mul(r19081, r19079, r19080, MPFR_RNDN); ; mpfr_mul(r19083, r19044, r19082, MPFR_RNDN); mpfr_mul(r19084, r19044, r19083, MPFR_RNDN); ; mpfr_mul(r19086, r19044, r19085, MPFR_RNDN); mpfr_mul(r19087, r19059, r19086, MPFR_RNDN); mpfr_add(r19088, r19084, r19087, MPFR_RNDN); mpfr_add(r19089, r19088, r19055, MPFR_RNDN); mpfr_add(r19090, r19081, r19089, MPFR_RNDN); mpfr_add(r19091, r19078, r19090, MPFR_RNDN); mpfr_div(r19092, r19073, r19091, MPFR_RNDN); mpfr_div(r19093, r19044, r19092, MPFR_RNDN); return mpfr_get_d(r19093, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); mpfr_init_set_str(r19097, "0.0694555761", 10, MPFR_RNDN); mpfr_init_set_str(r19098, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r19099); mpfr_init(r19100); mpfr_init(r19101); mpfr_init_set_str(r19102, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19103); mpfr_init(r19104); mpfr_init_set_str(r19105, "1", 10, MPFR_RNDN); mpfr_init(r19106); mpfr_init_set_str(r19107, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19108); mpfr_init(r19109); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init_set_str(r19113, "2", 10, MPFR_RNDN); mpfr_init(r19114); mpfr_init_set_str(r19115, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19116); mpfr_init(r19117); mpfr_init_set_str(r19118, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19119); mpfr_init(r19120); mpfr_init(r19121); mpfr_init(r19122); mpfr_init(r19123); mpfr_init_set_str(r19124, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19125); mpfr_init_set_str(r19126, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19127); mpfr_init(r19128); mpfr_init(r19129); mpfr_init(r19130); mpfr_init(r19131); mpfr_init_set_str(r19132, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19133); mpfr_init(r19134); mpfr_init_set_str(r19135, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init(r19142); mpfr_init(r19143); } double f_dm(double x) { mpfr_set_d(r19094, x, MPFR_RNDN); mpfr_sqr(r19095, r19094, MPFR_RNDN); mpfr_mul(r19096, r19095, r19095, MPFR_RNDN); mpfr_mul(r19096, r19096, r19095, MPFR_RNDN); ; ; mpfr_mul(r19099, r19095, r19098, MPFR_RNDN); mpfr_add(r19100, r19097, r19099, MPFR_RNDN); mpfr_mul(r19101, r19096, r19100, MPFR_RNDN); ; mpfr_mul(r19103, r19102, r19094, MPFR_RNDN); mpfr_mul(r19104, r19103, r19094, MPFR_RNDN); ; mpfr_add(r19106, r19104, r19105, MPFR_RNDN); ; mpfr_mul(r19108, r19094, r19107, MPFR_RNDN); mpfr_mul(r19109, r19094, r19094, MPFR_RNDN); mpfr_mul(r19109, r19109, r19094, MPFR_RNDN); mpfr_mul(r19110, r19108, r19109, MPFR_RNDN); mpfr_add(r19111, r19106, r19110, MPFR_RNDN); mpfr_add(r19112, r19101, r19111, MPFR_RNDN); ; mpfr_mul(r19114, r19094, r19113, MPFR_RNDN); ; mpfr_mul(r19116, r19094, r19115, MPFR_RNDN); mpfr_mul(r19117, r19114, r19116, MPFR_RNDN); ; mpfr_add(r19119, r19117, r19118, MPFR_RNDN); mpfr_mul(r19120, r19095, r19095, MPFR_RNDN); mpfr_mul(r19121, r19120, r19096, MPFR_RNDN); mpfr_mul(r19122, r19119, r19121, MPFR_RNDN); mpfr_add(r19123, r19112, r19122, MPFR_RNDN); ; mpfr_mul(r19125, r19095, r19124, MPFR_RNDN); ; mpfr_add(r19127, r19125, r19126, MPFR_RNDN); mpfr_mul(r19128, r19096, r19127, MPFR_RNDN); mpfr_mul(r19129, r19120, r19120, MPFR_RNDN); mpfr_mul(r19130, r19116, r19094, MPFR_RNDN); mpfr_mul(r19131, r19129, r19130, MPFR_RNDN); ; mpfr_mul(r19133, r19094, r19132, MPFR_RNDN); mpfr_mul(r19134, r19094, r19133, MPFR_RNDN); ; mpfr_mul(r19136, r19094, r19135, MPFR_RNDN); mpfr_mul(r19137, r19109, r19136, MPFR_RNDN); mpfr_add(r19138, r19134, r19137, MPFR_RNDN); mpfr_add(r19139, r19138, r19105, MPFR_RNDN); mpfr_add(r19140, r19131, r19139, MPFR_RNDN); mpfr_add(r19141, r19128, r19140, MPFR_RNDN); mpfr_div(r19142, r19123, r19141, MPFR_RNDN); mpfr_div(r19143, r19094, r19142, MPFR_RNDN); return mpfr_get_d(r19143, MPFR_RNDN); }