#include #include #include #include #include char *name = "Jmat.Real.dawson"; double f_if(float x) { float r28727 = 1.0f; float r28728 = 0.1049934947f; float r28729 = x; float r28730 = r28729 * r28729; float r28731 = r28728 * r28730; float r28732 = r28727 + r28731; float r28733 = 0.0424060604f; float r28734 = r28730 * r28730; float r28735 = r28733 * r28734; float r28736 = r28732 + r28735; float r28737 = 0.0072644182f; float r28738 = r28734 * r28730; float r28739 = r28737 * r28738; float r28740 = r28736 + r28739; float r28741 = 0.0005064034f; float r28742 = r28738 * r28730; float r28743 = r28741 * r28742; float r28744 = r28740 + r28743; float r28745 = 0.0001789971f; float r28746 = r28742 * r28730; float r28747 = r28745 * r28746; float r28748 = r28744 + r28747; float r28749 = 0.7715471019f; float r28750 = r28749 * r28730; float r28751 = r28727 + r28750; float r28752 = 0.2909738639f; float r28753 = r28752 * r28734; float r28754 = r28751 + r28753; float r28755 = 0.0694555761f; float r28756 = r28755 * r28738; float r28757 = r28754 + r28756; float r28758 = 0.0140005442f; float r28759 = r28758 * r28742; float r28760 = r28757 + r28759; float r28761 = 0.0008327945f; float r28762 = r28761 * r28746; float r28763 = r28760 + r28762; float r28764 = 2.0f; float r28765 = r28764 * r28745; float r28766 = r28746 * r28730; float r28767 = r28765 * r28766; float r28768 = r28763 + r28767; float r28769 = r28748 / r28768; float r28770 = r28769 * r28729; return r28770; } double f_id(double x) { double r28771 = 1.0; double r28772 = 0.1049934947; double r28773 = x; double r28774 = r28773 * r28773; double r28775 = r28772 * r28774; double r28776 = r28771 + r28775; double r28777 = 0.0424060604; double r28778 = r28774 * r28774; double r28779 = r28777 * r28778; double r28780 = r28776 + r28779; double r28781 = 0.0072644182; double r28782 = r28778 * r28774; double r28783 = r28781 * r28782; double r28784 = r28780 + r28783; double r28785 = 0.0005064034; double r28786 = r28782 * r28774; double r28787 = r28785 * r28786; double r28788 = r28784 + r28787; double r28789 = 0.0001789971; double r28790 = r28786 * r28774; double r28791 = r28789 * r28790; double r28792 = r28788 + r28791; double r28793 = 0.7715471019; double r28794 = r28793 * r28774; double r28795 = r28771 + r28794; double r28796 = 0.2909738639; double r28797 = r28796 * r28778; double r28798 = r28795 + r28797; double r28799 = 0.0694555761; double r28800 = r28799 * r28782; double r28801 = r28798 + r28800; double r28802 = 0.0140005442; double r28803 = r28802 * r28786; double r28804 = r28801 + r28803; double r28805 = 0.0008327945; double r28806 = r28805 * r28790; double r28807 = r28804 + r28806; double r28808 = 2.0; double r28809 = r28808 * r28789; double r28810 = r28790 * r28774; double r28811 = r28809 * r28810; double r28812 = r28807 + r28811; double r28813 = r28792 / r28812; double r28814 = r28813 * r28773; return r28814; } double f_of(float x) { float r28815 = 0.0001789971f; float r28816 = x; float r28817 = r28816 * (r28816 * r28816); float r28818 = r28817 * r28817; float r28819 = r28816 * r28816; float r28820 = r28819 * r28817; float r28821 = r28818 * r28820; float r28822 = r28815 * r28821; float r28823 = r28819 * r28819; float r28824 = r28818 * r28823; float r28825 = r28816 * r28815; float r28826 = r28816 + r28816; float r28827 = r28825 * r28826; float r28828 = 0.0008327945f; float r28829 = r28827 + r28828; float r28830 = r28824 * r28829; float r28831 = 0.0140005442f; float r28832 = r28831 * r28819; float r28833 = 0.0694555761f; float r28834 = r28832 + r28833; float r28835 = r28818 * r28834; float r28836 = r28830 + r28835; float r28837 = 0.7715471019f; float r28838 = r28837 * r28819; float r28839 = 1.0f; float r28840 = 0.2909738639f; float r28841 = r28816 * r28840; float r28842 = r28841 * r28817; float r28843 = r28839 + r28842; float r28844 = r28838 + r28843; float r28845 = r28836 + r28844; float r28846 = r28822 / r28845; float r28847 = 0.1049934947f; float r28848 = r28819 * r28847; float r28849 = 0.0424060604f; float r28850 = r28849 * r28816; float r28851 = r28850 * r28817; float r28852 = r28851 + r28839; float r28853 = r28848 + r28852; float r28854 = 0.0005064034f; float r28855 = r28819 * r28854; float r28856 = 0.0072644182f; float r28857 = r28855 + r28856; float r28858 = r28818 * r28857; float r28859 = r28853 + r28858; float r28860 = r28816 * r28859; float r28861 = r28815 + r28815; float r28862 = r28819 * r28861; float r28863 = r28862 + r28828; float r28864 = r28824 * r28863; float r28865 = r28834 * r28818; float r28866 = r28844 + r28865; float r28867 = r28864 + r28866; float r28868 = r28860 / r28867; float r28869 = r28846 + r28868; return r28869; } double f_od(double x) { double r28870 = 0.0001789971; double r28871 = x; double r28872 = r28871 * (r28871 * r28871); double r28873 = r28872 * r28872; double r28874 = r28871 * r28871; double r28875 = r28874 * r28872; double r28876 = r28873 * r28875; double r28877 = r28870 * r28876; double r28878 = r28874 * r28874; double r28879 = r28873 * r28878; double r28880 = r28871 * r28870; double r28881 = r28871 + r28871; double r28882 = r28880 * r28881; double r28883 = 0.0008327945; double r28884 = r28882 + r28883; double r28885 = r28879 * r28884; double r28886 = 0.0140005442; double r28887 = r28886 * r28874; double r28888 = 0.0694555761; double r28889 = r28887 + r28888; double r28890 = r28873 * r28889; double r28891 = r28885 + r28890; double r28892 = 0.7715471019; double r28893 = r28892 * r28874; double r28894 = 1.0; double r28895 = 0.2909738639; double r28896 = r28871 * r28895; double r28897 = r28896 * r28872; double r28898 = r28894 + r28897; double r28899 = r28893 + r28898; double r28900 = r28891 + r28899; double r28901 = r28877 / r28900; double r28902 = 0.1049934947; double r28903 = r28874 * r28902; double r28904 = 0.0424060604; double r28905 = r28904 * r28871; double r28906 = r28905 * r28872; double r28907 = r28906 + r28894; double r28908 = r28903 + r28907; double r28909 = 0.0005064034; double r28910 = r28874 * r28909; double r28911 = 0.0072644182; double r28912 = r28910 + r28911; double r28913 = r28873 * r28912; double r28914 = r28908 + r28913; double r28915 = r28871 * r28914; double r28916 = r28870 + r28870; double r28917 = r28874 * r28916; double r28918 = r28917 + r28883; double r28919 = r28879 * r28918; double r28920 = r28889 * r28873; double r28921 = r28899 + r28920; double r28922 = r28919 + r28921; double r28923 = r28915 / r28922; double r28924 = r28901 + r28923; return r28924; } 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 r28925, r28926, r28927, r28928, r28929, r28930, r28931, r28932, r28933, r28934, r28935, r28936, r28937, r28938, r28939, r28940, r28941, r28942, r28943, r28944, r28945, r28946, r28947, r28948, r28949, r28950, r28951, r28952, r28953, r28954, r28955, r28956, r28957, r28958, r28959, r28960, r28961, r28962, r28963, r28964, r28965, r28966, r28967, r28968; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r28925, "1", 10, MPFR_RNDN); mpfr_init_set_str(r28926, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r28927); mpfr_init(r28928); mpfr_init(r28929); mpfr_init(r28930); mpfr_init_set_str(r28931, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r28932); mpfr_init(r28933); mpfr_init(r28934); mpfr_init_set_str(r28935, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r28936); mpfr_init(r28937); mpfr_init(r28938); mpfr_init_set_str(r28939, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r28940); mpfr_init(r28941); mpfr_init(r28942); mpfr_init_set_str(r28943, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28944); mpfr_init(r28945); mpfr_init(r28946); mpfr_init_set_str(r28947, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28948); mpfr_init(r28949); mpfr_init_set_str(r28950, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r28951); mpfr_init(r28952); mpfr_init_set_str(r28953, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28954); mpfr_init(r28955); mpfr_init_set_str(r28956, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r28957); mpfr_init(r28958); mpfr_init_set_str(r28959, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28960); mpfr_init(r28961); mpfr_init_set_str(r28962, "2", 10, MPFR_RNDN); mpfr_init(r28963); mpfr_init(r28964); mpfr_init(r28965); mpfr_init(r28966); mpfr_init(r28967); mpfr_init(r28968); } double f_im(double x) { ; ; mpfr_set_d(r28927, x, MPFR_RNDN); mpfr_mul(r28928, r28927, r28927, MPFR_RNDN); mpfr_mul(r28929, r28926, r28928, MPFR_RNDN); mpfr_add(r28930, r28925, r28929, MPFR_RNDN); ; mpfr_mul(r28932, r28928, r28928, MPFR_RNDN); mpfr_mul(r28933, r28931, r28932, MPFR_RNDN); mpfr_add(r28934, r28930, r28933, MPFR_RNDN); ; mpfr_mul(r28936, r28932, r28928, MPFR_RNDN); mpfr_mul(r28937, r28935, r28936, MPFR_RNDN); mpfr_add(r28938, r28934, r28937, MPFR_RNDN); ; mpfr_mul(r28940, r28936, r28928, MPFR_RNDN); mpfr_mul(r28941, r28939, r28940, MPFR_RNDN); mpfr_add(r28942, r28938, r28941, MPFR_RNDN); ; mpfr_mul(r28944, r28940, r28928, MPFR_RNDN); mpfr_mul(r28945, r28943, r28944, MPFR_RNDN); mpfr_add(r28946, r28942, r28945, MPFR_RNDN); ; mpfr_mul(r28948, r28947, r28928, MPFR_RNDN); mpfr_add(r28949, r28925, r28948, MPFR_RNDN); ; mpfr_mul(r28951, r28950, r28932, MPFR_RNDN); mpfr_add(r28952, r28949, r28951, MPFR_RNDN); ; mpfr_mul(r28954, r28953, r28936, MPFR_RNDN); mpfr_add(r28955, r28952, r28954, MPFR_RNDN); ; mpfr_mul(r28957, r28956, r28940, MPFR_RNDN); mpfr_add(r28958, r28955, r28957, MPFR_RNDN); ; mpfr_mul(r28960, r28959, r28944, MPFR_RNDN); mpfr_add(r28961, r28958, r28960, MPFR_RNDN); ; mpfr_mul(r28963, r28962, r28943, MPFR_RNDN); mpfr_mul(r28964, r28944, r28928, MPFR_RNDN); mpfr_mul(r28965, r28963, r28964, MPFR_RNDN); mpfr_add(r28966, r28961, r28965, MPFR_RNDN); mpfr_div(r28967, r28946, r28966, MPFR_RNDN); mpfr_mul(r28968, r28967, r28927, MPFR_RNDN); return mpfr_get_d(r28968, MPFR_RNDN); } static mpfr_t r28969, r28970, r28971, r28972, r28973, r28974, r28975, r28976, r28977, r28978, r28979, r28980, r28981, r28982, r28983, r28984, r28985, r28986, r28987, r28988, r28989, r28990, r28991, r28992, r28993, r28994, r28995, r28996, r28997, r28998, r28999, r29000, r29001, r29002, r29003, r29004, r29005, r29006, r29007, r29008, r29009, r29010, r29011, r29012, r29013, r29014, r29015, r29016, r29017, r29018, r29019, r29020, r29021, r29022, r29023; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r28969, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r28970); mpfr_init(r28971); mpfr_init(r28972); mpfr_init(r28973); mpfr_init(r28974); mpfr_init(r28975); mpfr_init(r28976); mpfr_init(r28977); mpfr_init(r28978); mpfr_init(r28979); mpfr_init(r28980); mpfr_init(r28981); mpfr_init_set_str(r28982, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r28983); mpfr_init(r28984); mpfr_init_set_str(r28985, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r28986); mpfr_init_set_str(r28987, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r28988); mpfr_init(r28989); mpfr_init(r28990); mpfr_init_set_str(r28991, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r28992); mpfr_init_set_str(r28993, "1", 10, MPFR_RNDN); mpfr_init_set_str(r28994, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r28995); mpfr_init(r28996); mpfr_init(r28997); mpfr_init(r28998); mpfr_init(r28999); mpfr_init(r29000); mpfr_init_set_str(r29001, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r29002); mpfr_init_set_str(r29003, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r29004); mpfr_init(r29005); mpfr_init(r29006); mpfr_init(r29007); mpfr_init_set_str(r29008, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r29009); mpfr_init_set_str(r29010, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r29011); mpfr_init(r29012); mpfr_init(r29013); mpfr_init(r29014); mpfr_init(r29015); mpfr_init(r29016); mpfr_init(r29017); mpfr_init(r29018); mpfr_init(r29019); mpfr_init(r29020); mpfr_init(r29021); mpfr_init(r29022); mpfr_init(r29023); } double f_fm(double x) { ; mpfr_set_d(r28970, x, MPFR_RNDN); mpfr_mul(r28971, r28970, r28970, MPFR_RNDN); mpfr_mul(r28971, r28971, r28970, MPFR_RNDN); mpfr_sqr(r28972, r28971, MPFR_RNDN); mpfr_sqr(r28973, r28970, MPFR_RNDN); mpfr_mul(r28974, r28973, r28971, MPFR_RNDN); mpfr_mul(r28975, r28972, r28974, MPFR_RNDN); mpfr_mul(r28976, r28969, r28975, MPFR_RNDN); mpfr_sqr(r28977, r28973, MPFR_RNDN); mpfr_mul(r28978, r28972, r28977, MPFR_RNDN); mpfr_mul(r28979, r28970, r28969, MPFR_RNDN); mpfr_add(r28980, r28970, r28970, MPFR_RNDN); mpfr_mul(r28981, r28979, r28980, MPFR_RNDN); ; mpfr_add(r28983, r28981, r28982, MPFR_RNDN); mpfr_mul(r28984, r28978, r28983, MPFR_RNDN); ; mpfr_mul(r28986, r28985, r28973, MPFR_RNDN); ; mpfr_add(r28988, r28986, r28987, MPFR_RNDN); mpfr_mul(r28989, r28972, r28988, MPFR_RNDN); mpfr_add(r28990, r28984, r28989, MPFR_RNDN); ; mpfr_mul(r28992, r28991, r28973, MPFR_RNDN); ; ; mpfr_mul(r28995, r28970, r28994, MPFR_RNDN); mpfr_mul(r28996, r28995, r28971, MPFR_RNDN); mpfr_add(r28997, r28993, r28996, MPFR_RNDN); mpfr_add(r28998, r28992, r28997, MPFR_RNDN); mpfr_add(r28999, r28990, r28998, MPFR_RNDN); mpfr_div(r29000, r28976, r28999, MPFR_RNDN); ; mpfr_mul(r29002, r28973, r29001, MPFR_RNDN); ; mpfr_mul(r29004, r29003, r28970, MPFR_RNDN); mpfr_mul(r29005, r29004, r28971, MPFR_RNDN); mpfr_add(r29006, r29005, r28993, MPFR_RNDN); mpfr_add(r29007, r29002, r29006, MPFR_RNDN); ; mpfr_mul(r29009, r28973, r29008, MPFR_RNDN); ; mpfr_add(r29011, r29009, r29010, MPFR_RNDN); mpfr_mul(r29012, r28972, r29011, MPFR_RNDN); mpfr_add(r29013, r29007, r29012, MPFR_RNDN); mpfr_mul(r29014, r28970, r29013, MPFR_RNDN); mpfr_add(r29015, r28969, r28969, MPFR_RNDN); mpfr_mul(r29016, r28973, r29015, MPFR_RNDN); mpfr_add(r29017, r29016, r28982, MPFR_RNDN); mpfr_mul(r29018, r28978, r29017, MPFR_RNDN); mpfr_mul(r29019, r28988, r28972, MPFR_RNDN); mpfr_add(r29020, r28998, r29019, MPFR_RNDN); mpfr_add(r29021, r29018, r29020, MPFR_RNDN); mpfr_div(r29022, r29014, r29021, MPFR_RNDN); mpfr_add(r29023, r29000, r29022, MPFR_RNDN); return mpfr_get_d(r29023, MPFR_RNDN); } static mpfr_t r29024, r29025, r29026, r29027, r29028, r29029, r29030, r29031, r29032, r29033, r29034, r29035, r29036, r29037, r29038, r29039, r29040, r29041, r29042, r29043, r29044, r29045, r29046, r29047, r29048, r29049, r29050, r29051, r29052, r29053, r29054, r29055, r29056, r29057, r29058, r29059, r29060, r29061, r29062, r29063, r29064, r29065, r29066, r29067, r29068, r29069, r29070, r29071, r29072, r29073, r29074, r29075, r29076, r29077, r29078; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r29024, "0.0001789971", 10, MPFR_RNDN); mpfr_init(r29025); mpfr_init(r29026); mpfr_init(r29027); mpfr_init(r29028); mpfr_init(r29029); mpfr_init(r29030); mpfr_init(r29031); mpfr_init(r29032); mpfr_init(r29033); mpfr_init(r29034); mpfr_init(r29035); mpfr_init(r29036); mpfr_init_set_str(r29037, "0.0008327945", 10, MPFR_RNDN); mpfr_init(r29038); mpfr_init(r29039); mpfr_init_set_str(r29040, "0.0140005442", 10, MPFR_RNDN); mpfr_init(r29041); mpfr_init_set_str(r29042, "0.0694555761", 10, MPFR_RNDN); mpfr_init(r29043); mpfr_init(r29044); mpfr_init(r29045); mpfr_init_set_str(r29046, "0.7715471019", 10, MPFR_RNDN); mpfr_init(r29047); mpfr_init_set_str(r29048, "1", 10, MPFR_RNDN); mpfr_init_set_str(r29049, "0.2909738639", 10, MPFR_RNDN); mpfr_init(r29050); mpfr_init(r29051); mpfr_init(r29052); mpfr_init(r29053); mpfr_init(r29054); mpfr_init(r29055); mpfr_init_set_str(r29056, "0.1049934947", 10, MPFR_RNDN); mpfr_init(r29057); mpfr_init_set_str(r29058, "0.0424060604", 10, MPFR_RNDN); mpfr_init(r29059); mpfr_init(r29060); mpfr_init(r29061); mpfr_init(r29062); mpfr_init_set_str(r29063, "0.0005064034", 10, MPFR_RNDN); mpfr_init(r29064); mpfr_init_set_str(r29065, "0.0072644182", 10, MPFR_RNDN); mpfr_init(r29066); mpfr_init(r29067); mpfr_init(r29068); mpfr_init(r29069); mpfr_init(r29070); mpfr_init(r29071); mpfr_init(r29072); mpfr_init(r29073); mpfr_init(r29074); mpfr_init(r29075); mpfr_init(r29076); mpfr_init(r29077); mpfr_init(r29078); } double f_dm(double x) { ; mpfr_set_d(r29025, x, MPFR_RNDN); mpfr_mul(r29026, r29025, r29025, MPFR_RNDN); mpfr_mul(r29026, r29026, r29025, MPFR_RNDN); mpfr_sqr(r29027, r29026, MPFR_RNDN); mpfr_sqr(r29028, r29025, MPFR_RNDN); mpfr_mul(r29029, r29028, r29026, MPFR_RNDN); mpfr_mul(r29030, r29027, r29029, MPFR_RNDN); mpfr_mul(r29031, r29024, r29030, MPFR_RNDN); mpfr_sqr(r29032, r29028, MPFR_RNDN); mpfr_mul(r29033, r29027, r29032, MPFR_RNDN); mpfr_mul(r29034, r29025, r29024, MPFR_RNDN); mpfr_add(r29035, r29025, r29025, MPFR_RNDN); mpfr_mul(r29036, r29034, r29035, MPFR_RNDN); ; mpfr_add(r29038, r29036, r29037, MPFR_RNDN); mpfr_mul(r29039, r29033, r29038, MPFR_RNDN); ; mpfr_mul(r29041, r29040, r29028, MPFR_RNDN); ; mpfr_add(r29043, r29041, r29042, MPFR_RNDN); mpfr_mul(r29044, r29027, r29043, MPFR_RNDN); mpfr_add(r29045, r29039, r29044, MPFR_RNDN); ; mpfr_mul(r29047, r29046, r29028, MPFR_RNDN); ; ; mpfr_mul(r29050, r29025, r29049, MPFR_RNDN); mpfr_mul(r29051, r29050, r29026, MPFR_RNDN); mpfr_add(r29052, r29048, r29051, MPFR_RNDN); mpfr_add(r29053, r29047, r29052, MPFR_RNDN); mpfr_add(r29054, r29045, r29053, MPFR_RNDN); mpfr_div(r29055, r29031, r29054, MPFR_RNDN); ; mpfr_mul(r29057, r29028, r29056, MPFR_RNDN); ; mpfr_mul(r29059, r29058, r29025, MPFR_RNDN); mpfr_mul(r29060, r29059, r29026, MPFR_RNDN); mpfr_add(r29061, r29060, r29048, MPFR_RNDN); mpfr_add(r29062, r29057, r29061, MPFR_RNDN); ; mpfr_mul(r29064, r29028, r29063, MPFR_RNDN); ; mpfr_add(r29066, r29064, r29065, MPFR_RNDN); mpfr_mul(r29067, r29027, r29066, MPFR_RNDN); mpfr_add(r29068, r29062, r29067, MPFR_RNDN); mpfr_mul(r29069, r29025, r29068, MPFR_RNDN); mpfr_add(r29070, r29024, r29024, MPFR_RNDN); mpfr_mul(r29071, r29028, r29070, MPFR_RNDN); mpfr_add(r29072, r29071, r29037, MPFR_RNDN); mpfr_mul(r29073, r29033, r29072, MPFR_RNDN); mpfr_mul(r29074, r29043, r29027, MPFR_RNDN); mpfr_add(r29075, r29053, r29074, MPFR_RNDN); mpfr_add(r29076, r29073, r29075, MPFR_RNDN); mpfr_div(r29077, r29069, r29076, MPFR_RNDN); mpfr_add(r29078, r29055, r29077, MPFR_RNDN); return mpfr_get_d(r29078, MPFR_RNDN); }