#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r18767 = 1.0f; float r18768 = 0.1049934947f; float r18769 = x; float r18770 = r18769 * r18769; float r18771 = r18768 * r18770; float r18772 = r18767 + r18771; float r18773 = 0.0424060604f; float r18774 = r18770 * r18770; float r18775 = r18773 * r18774; float r18776 = r18772 + r18775; float r18777 = 0.0072644182f; float r18778 = r18774 * r18770; float r18779 = r18777 * r18778; float r18780 = r18776 + r18779; float r18781 = 0.0005064034f; float r18782 = r18778 * r18770; float r18783 = r18781 * r18782; float r18784 = r18780 + r18783; float r18785 = 0.0001789971f; float r18786 = r18782 * r18770; float r18787 = r18785 * r18786; float r18788 = r18784 + r18787; float r18789 = 0.7715471019f; float r18790 = r18789 * r18770; float r18791 = r18767 + r18790; float r18792 = 0.2909738639f; float r18793 = r18792 * r18774; float r18794 = r18791 + r18793; float r18795 = 0.0694555761f; float r18796 = r18795 * r18778; float r18797 = r18794 + r18796; float r18798 = 0.0140005442f; float r18799 = r18798 * r18782; float r18800 = r18797 + r18799; float r18801 = 0.0008327945f; float r18802 = r18801 * r18786; float r18803 = r18800 + r18802; float r18804 = 2.0f; float r18805 = r18804 * r18785; float r18806 = r18786 * r18770; float r18807 = r18805 * r18806; float r18808 = r18803 + r18807; float r18809 = r18788 / r18808; float r18810 = r18809 * r18769; return r18810; } double f_id(double x) { double r18811 = 1.0; double r18812 = 0.1049934947; double r18813 = x; double r18814 = r18813 * r18813; double r18815 = r18812 * r18814; double r18816 = r18811 + r18815; double r18817 = 0.0424060604; double r18818 = r18814 * r18814; double r18819 = r18817 * r18818; double r18820 = r18816 + r18819; double r18821 = 0.0072644182; double r18822 = r18818 * r18814; double r18823 = r18821 * r18822; double r18824 = r18820 + r18823; double r18825 = 0.0005064034; double r18826 = r18822 * r18814; double r18827 = r18825 * r18826; double r18828 = r18824 + r18827; double r18829 = 0.0001789971; double r18830 = r18826 * r18814; double r18831 = r18829 * r18830; double r18832 = r18828 + r18831; double r18833 = 0.7715471019; double r18834 = r18833 * r18814; double r18835 = r18811 + r18834; double r18836 = 0.2909738639; double r18837 = r18836 * r18818; double r18838 = r18835 + r18837; double r18839 = 0.0694555761; double r18840 = r18839 * r18822; double r18841 = r18838 + r18840; double r18842 = 0.0140005442; double r18843 = r18842 * r18826; double r18844 = r18841 + r18843; double r18845 = 0.0008327945; double r18846 = r18845 * r18830; double r18847 = r18844 + r18846; double r18848 = 2.0; double r18849 = r18848 * r18829; double r18850 = r18830 * r18814; double r18851 = r18849 * r18850; double r18852 = r18847 + r18851; double r18853 = r18832 / r18852; double r18854 = r18853 * r18813; return r18854; } double f_of(float x) { float r18855 = x; float r18856 = 0.0001789971f; float r18857 = r18855 * (r18855 * r18855); float r18858 = r18855 * r18855; float r18859 = r18857 * r18858; float r18860 = r18859 * r18859; float r18861 = 0.0424060604f; float r18862 = r18858 * r18858; float r18863 = 0.1049934947f; float r18864 = r18855 * r18863; float r18865 = 1.0f; float r18866 = fma(r18864, r18855, r18865); float r18867 = fma(r18861, r18862, r18866); float r18868 = fma(r18856, r18860, r18867); float r18869 = 0.0005064034f; float r18870 = r18862 * r18862; float r18871 = 0.0072644182f; float r18872 = r18857 * r18871; float r18873 = r18872 * r18857; float r18874 = fma(r18869, r18870, r18873); float r18875 = r18868 + r18874; float r18876 = r18855 * r18875; float r18877 = 2.0f; float r18878 = r18856 * r18877; float r18879 = r18857 * r18857; float r18880 = r18855 * r18855; float r18881 = r18880 * r18857; float r18882 = r18879 * r18881; float r18883 = r18855 * r18882; float r18884 = 0.0008327945f; float r18885 = r18857 * (r18857 * r18857); float r18886 = r18885 * r18855; float r18887 = 0.2909738639f; float r18888 = r18880 * r18880; float r18889 = 0.7715471019f; float r18890 = r18855 * r18889; float r18891 = fma(r18890, r18855, r18865); float r18892 = fma(r18887, r18888, r18891); float r18893 = fma(r18884, r18886, r18892); float r18894 = r18888 * r18888; float r18895 = 0.0140005442f; float r18896 = 0.0694555761f; float r18897 = r18879 * r18896; float r18898 = fma(r18894, r18895, r18897); float r18899 = r18893 + r18898; float r18900 = fma(r18878, r18883, r18899); float r18901 = r18876 / r18900; return r18901; } double f_od(double x) { double r18902 = x; double r18903 = 0.0001789971; double r18904 = r18902 * (r18902 * r18902); double r18905 = r18902 * r18902; double r18906 = r18904 * r18905; double r18907 = r18906 * r18906; double r18908 = 0.0424060604; double r18909 = r18905 * r18905; double r18910 = 0.1049934947; double r18911 = r18902 * r18910; double r18912 = 1.0; double r18913 = fma(r18911, r18902, r18912); double r18914 = fma(r18908, r18909, r18913); double r18915 = fma(r18903, r18907, r18914); double r18916 = 0.0005064034; double r18917 = r18909 * r18909; double r18918 = 0.0072644182; double r18919 = r18904 * r18918; double r18920 = r18919 * r18904; double r18921 = fma(r18916, r18917, r18920); double r18922 = r18915 + r18921; double r18923 = r18902 * r18922; double r18924 = 2.0; double r18925 = r18903 * r18924; double r18926 = r18904 * r18904; double r18927 = r18902 * r18902; double r18928 = r18927 * r18904; double r18929 = r18926 * r18928; double r18930 = r18902 * r18929; double r18931 = 0.0008327945; double r18932 = r18904 * (r18904 * r18904); double r18933 = r18932 * r18902; double r18934 = 0.2909738639; double r18935 = r18927 * r18927; double r18936 = 0.7715471019; double r18937 = r18902 * r18936; double r18938 = fma(r18937, r18902, r18912); double r18939 = fma(r18934, r18935, r18938); double r18940 = fma(r18931, r18933, r18939); double r18941 = r18935 * r18935; double r18942 = 0.0140005442; double r18943 = 0.0694555761; double r18944 = r18926 * r18943; double r18945 = fma(r18941, r18942, r18944); double r18946 = r18940 + r18945; double r18947 = fma(r18925, r18930, r18946); double r18948 = r18923 / r18947; return r18948; } 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 r18949, r18950, r18951, r18952, r18953, r18954, r18955, r18956, r18957, r18958, r18959, r18960, r18961, r18962, r18963, r18964, r18965, r18966, r18967, r18968, r18969, r18970, r18971, r18972, r18973, r18974, r18975, r18976, r18977, r18978, r18979, r18980, r18981, r18982, r18983, r18984, r18985, r18986, r18987, r18988, r18989, r18990, r18991, r18992; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18949, "1", 10, MPFR_RNDN); mpfr_init_set_str(r18950, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r18951); mpfr_init(r18952); mpfr_init(r18953); mpfr_init(r18954); mpfr_init_set_str(r18955, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r18956); mpfr_init(r18957); mpfr_init(r18958); mpfr_init_set_str(r18959, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r18960); mpfr_init(r18961); mpfr_init(r18962); mpfr_init_set_str(r18963, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r18964); mpfr_init(r18965); mpfr_init(r18966); mpfr_init_set_str(r18967, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r18968); mpfr_init(r18969); mpfr_init(r18970); mpfr_init_set_str(r18971, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r18972); mpfr_init(r18973); mpfr_init_set_str(r18974, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r18975); mpfr_init(r18976); mpfr_init_set_str(r18977, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r18978); mpfr_init(r18979); mpfr_init_set_str(r18980, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r18981); mpfr_init(r18982); mpfr_init_set_str(r18983, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r18984); mpfr_init(r18985); mpfr_init_set_str(r18986, "2", 10, MPFR_RNDN); mpfr_init(r18987); mpfr_init(r18988); mpfr_init(r18989); mpfr_init(r18990); mpfr_init(r18991); mpfr_init(r18992); } double f_im(double x) { ; ; mpfr_set_d(r18951, x, MPFR_RNDN); mpfr_mul(r18952, r18951, r18951, MPFR_RNDN); mpfr_mul(r18953, r18950, r18952, MPFR_RNDN); mpfr_add(r18954, r18949, r18953, MPFR_RNDN); ; mpfr_mul(r18956, r18952, r18952, MPFR_RNDN); mpfr_mul(r18957, r18955, r18956, MPFR_RNDN); mpfr_add(r18958, r18954, r18957, MPFR_RNDN); ; mpfr_mul(r18960, r18956, r18952, MPFR_RNDN); mpfr_mul(r18961, r18959, r18960, MPFR_RNDN); mpfr_add(r18962, r18958, r18961, MPFR_RNDN); ; mpfr_mul(r18964, r18960, r18952, MPFR_RNDN); mpfr_mul(r18965, r18963, r18964, MPFR_RNDN); mpfr_add(r18966, r18962, r18965, MPFR_RNDN); ; mpfr_mul(r18968, r18964, r18952, MPFR_RNDN); mpfr_mul(r18969, r18967, r18968, MPFR_RNDN); mpfr_add(r18970, r18966, r18969, MPFR_RNDN); ; mpfr_mul(r18972, r18971, r18952, MPFR_RNDN); mpfr_add(r18973, r18949, r18972, MPFR_RNDN); ; mpfr_mul(r18975, r18974, r18956, MPFR_RNDN); mpfr_add(r18976, r18973, r18975, MPFR_RNDN); ; mpfr_mul(r18978, r18977, r18960, MPFR_RNDN); mpfr_add(r18979, r18976, r18978, MPFR_RNDN); ; mpfr_mul(r18981, r18980, r18964, MPFR_RNDN); mpfr_add(r18982, r18979, r18981, MPFR_RNDN); ; mpfr_mul(r18984, r18983, r18968, MPFR_RNDN); mpfr_add(r18985, r18982, r18984, MPFR_RNDN); ; mpfr_mul(r18987, r18986, r18967, MPFR_RNDN); mpfr_mul(r18988, r18968, r18952, MPFR_RNDN); mpfr_mul(r18989, r18987, r18988, MPFR_RNDN); mpfr_add(r18990, r18985, r18989, MPFR_RNDN); mpfr_div(r18991, r18970, r18990, MPFR_RNDN); mpfr_mul(r18992, r18991, r18951, MPFR_RNDN); return mpfr_get_d(r18992, MPFR_RNDN); } static mpfr_t r18993, r18994, r18995, r18996, r18997, r18998, r18999, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18993); mpfr_init_set_str(r18994, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r18995); mpfr_init(r18996); mpfr_init(r18997); mpfr_init(r18998); mpfr_init_set_str(r18999, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19000); mpfr_init_set_str(r19001, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19002); mpfr_init_set_str(r19003, "1", 10, MPFR_RNDN); mpfr_init(r19004); mpfr_init(r19005); mpfr_init(r19006); mpfr_init_set_str(r19007, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19008); mpfr_init_set_str(r19009, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19010); mpfr_init(r19011); mpfr_init(r19012); mpfr_init(r19013); mpfr_init(r19014); mpfr_init_set_str(r19015, "2", 10, MPFR_RNDN); mpfr_init(r19016); mpfr_init(r19017); mpfr_init(r19018); mpfr_init(r19019); mpfr_init(r19020); mpfr_init(r19021); mpfr_init_set_str(r19022, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19023); mpfr_init(r19024); mpfr_init_set_str(r19025, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19026); mpfr_init_set_str(r19027, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19028); mpfr_init(r19029); mpfr_init(r19030); mpfr_init(r19031); mpfr_init(r19032); mpfr_init_set_str(r19033, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r19034, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19035); mpfr_init(r19036); mpfr_init(r19037); mpfr_init(r19038); mpfr_init(r19039); } double f_fm(double x) { mpfr_set_d(r18993, x, MPFR_RNDN); ; mpfr_mul(r18995, r18993, r18993, MPFR_RNDN); mpfr_mul(r18995, r18995, r18993, MPFR_RNDN); mpfr_sqr(r18996, r18993, MPFR_RNDN); mpfr_mul(r18997, r18995, r18996, MPFR_RNDN); mpfr_mul(r18998, r18997, r18997, MPFR_RNDN); ; mpfr_mul(r19000, r18996, r18996, MPFR_RNDN); ; mpfr_mul(r19002, r18993, r19001, MPFR_RNDN); ; mpfr_fma(r19004, r19002, r18993, r19003, MPFR_RNDN); mpfr_fma(r19005, r18999, r19000, r19004, MPFR_RNDN); mpfr_fma(r19006, r18994, r18998, r19005, MPFR_RNDN); ; mpfr_mul(r19008, r19000, r19000, MPFR_RNDN); ; mpfr_mul(r19010, r18995, r19009, MPFR_RNDN); mpfr_mul(r19011, r19010, r18995, MPFR_RNDN); mpfr_fma(r19012, r19007, r19008, r19011, MPFR_RNDN); mpfr_add(r19013, r19006, r19012, MPFR_RNDN); mpfr_mul(r19014, r18993, r19013, MPFR_RNDN); ; mpfr_mul(r19016, r18994, r19015, MPFR_RNDN); mpfr_mul(r19017, r18995, r18995, MPFR_RNDN); mpfr_mul(r19018, r18993, r18993, MPFR_RNDN); mpfr_mul(r19019, r19018, r18995, MPFR_RNDN); mpfr_mul(r19020, r19017, r19019, MPFR_RNDN); mpfr_mul(r19021, r18993, r19020, MPFR_RNDN); ; mpfr_mul(r19023, r18995, r18995, MPFR_RNDN); mpfr_mul(r19023, r19023, r18995, MPFR_RNDN); mpfr_mul(r19024, r19023, r18993, MPFR_RNDN); ; mpfr_mul(r19026, r19018, r19018, MPFR_RNDN); ; mpfr_mul(r19028, r18993, r19027, MPFR_RNDN); mpfr_fma(r19029, r19028, r18993, r19003, MPFR_RNDN); mpfr_fma(r19030, r19025, r19026, r19029, MPFR_RNDN); mpfr_fma(r19031, r19022, r19024, r19030, MPFR_RNDN); mpfr_sqr(r19032, r19026, MPFR_RNDN); ; ; mpfr_mul(r19035, r19017, r19034, MPFR_RNDN); mpfr_fma(r19036, r19032, r19033, r19035, MPFR_RNDN); mpfr_add(r19037, r19031, r19036, MPFR_RNDN); mpfr_fma(r19038, r19016, r19021, r19037, MPFR_RNDN); mpfr_div(r19039, r19014, r19038, MPFR_RNDN); return mpfr_get_d(r19039, MPFR_RNDN); } static mpfr_t 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, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19040); mpfr_init_set_str(r19041, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r19042); mpfr_init(r19043); mpfr_init(r19044); mpfr_init(r19045); mpfr_init_set_str(r19046, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r19047); mpfr_init_set_str(r19048, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r19049); mpfr_init_set_str(r19050, "1", 10, MPFR_RNDN); mpfr_init(r19051); mpfr_init(r19052); mpfr_init(r19053); mpfr_init_set_str(r19054, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r19055); mpfr_init_set_str(r19056, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r19057); mpfr_init(r19058); mpfr_init(r19059); 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); mpfr_init_set_str(r19069, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r19070); mpfr_init(r19071); mpfr_init_set_str(r19072, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r19073); mpfr_init_set_str(r19074, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r19075); mpfr_init(r19076); mpfr_init(r19077); mpfr_init(r19078); mpfr_init(r19079); mpfr_init_set_str(r19080, "0.0140005442", 10, MPFR_RNDN); mpfr_init_set_str(r19081, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r19082); mpfr_init(r19083); mpfr_init(r19084); mpfr_init(r19085); mpfr_init(r19086); } double f_dm(double x) { mpfr_set_d(r19040, x, MPFR_RNDN); ; mpfr_mul(r19042, r19040, r19040, MPFR_RNDN); mpfr_mul(r19042, r19042, r19040, MPFR_RNDN); mpfr_sqr(r19043, r19040, MPFR_RNDN); mpfr_mul(r19044, r19042, r19043, MPFR_RNDN); mpfr_mul(r19045, r19044, r19044, MPFR_RNDN); ; mpfr_mul(r19047, r19043, r19043, MPFR_RNDN); ; mpfr_mul(r19049, r19040, r19048, MPFR_RNDN); ; mpfr_fma(r19051, r19049, r19040, r19050, MPFR_RNDN); mpfr_fma(r19052, r19046, r19047, r19051, MPFR_RNDN); mpfr_fma(r19053, r19041, r19045, r19052, MPFR_RNDN); ; mpfr_mul(r19055, r19047, r19047, MPFR_RNDN); ; mpfr_mul(r19057, r19042, r19056, MPFR_RNDN); mpfr_mul(r19058, r19057, r19042, MPFR_RNDN); mpfr_fma(r19059, r19054, r19055, r19058, MPFR_RNDN); mpfr_add(r19060, r19053, r19059, MPFR_RNDN); mpfr_mul(r19061, r19040, r19060, MPFR_RNDN); ; mpfr_mul(r19063, r19041, r19062, MPFR_RNDN); mpfr_mul(r19064, r19042, r19042, MPFR_RNDN); mpfr_mul(r19065, r19040, r19040, MPFR_RNDN); mpfr_mul(r19066, r19065, r19042, MPFR_RNDN); mpfr_mul(r19067, r19064, r19066, MPFR_RNDN); mpfr_mul(r19068, r19040, r19067, MPFR_RNDN); ; mpfr_mul(r19070, r19042, r19042, MPFR_RNDN); mpfr_mul(r19070, r19070, r19042, MPFR_RNDN); mpfr_mul(r19071, r19070, r19040, MPFR_RNDN); ; mpfr_mul(r19073, r19065, r19065, MPFR_RNDN); ; mpfr_mul(r19075, r19040, r19074, MPFR_RNDN); mpfr_fma(r19076, r19075, r19040, r19050, MPFR_RNDN); mpfr_fma(r19077, r19072, r19073, r19076, MPFR_RNDN); mpfr_fma(r19078, r19069, r19071, r19077, MPFR_RNDN); mpfr_sqr(r19079, r19073, MPFR_RNDN); ; ; mpfr_mul(r19082, r19064, r19081, MPFR_RNDN); mpfr_fma(r19083, r19079, r19080, r19082, MPFR_RNDN); mpfr_add(r19084, r19078, r19083, MPFR_RNDN); mpfr_fma(r19085, r19063, r19068, r19084, MPFR_RNDN); mpfr_div(r19086, r19061, r19085, MPFR_RNDN); return mpfr_get_d(r19086, MPFR_RNDN); }