#include #include #include #include #include char *name = "Jmat.Real.erfi, branch x less than or equal to 0.5"; double f_if(float x) { float r17836 = 1.0f; float r17837 = atan2(1.0, 0.0); float r17838 = sqrt(r17837); float r17839 = r17836 / r17838; float r17840 = 2.0f; float r17841 = x; float r17842 = fabs(r17841); float r17843 = r17840 * r17842; float r17844 = 3.0f; float r17845 = r17840 / r17844; float r17846 = r17842 * r17842; float r17847 = r17846 * r17842; float r17848 = r17845 * r17847; float r17849 = r17843 + r17848; float r17850 = 5.0f; float r17851 = r17836 / r17850; float r17852 = r17847 * r17842; float r17853 = r17852 * r17842; float r17854 = r17851 * r17853; float r17855 = r17849 + r17854; float r17856 = 21.0f; float r17857 = r17836 / r17856; float r17858 = r17853 * r17842; float r17859 = r17858 * r17842; float r17860 = r17857 * r17859; float r17861 = r17855 + r17860; float r17862 = r17839 * r17861; float r17863 = fabs(r17862); return r17863; } double f_id(double x) { double r17864 = 1.0; double r17865 = atan2(1.0, 0.0); double r17866 = sqrt(r17865); double r17867 = r17864 / r17866; double r17868 = 2.0; double r17869 = x; double r17870 = fabs(r17869); double r17871 = r17868 * r17870; double r17872 = 3.0; double r17873 = r17868 / r17872; double r17874 = r17870 * r17870; double r17875 = r17874 * r17870; double r17876 = r17873 * r17875; double r17877 = r17871 + r17876; double r17878 = 5.0; double r17879 = r17864 / r17878; double r17880 = r17875 * r17870; double r17881 = r17880 * r17870; double r17882 = r17879 * r17881; double r17883 = r17877 + r17882; double r17884 = 21.0; double r17885 = r17864 / r17884; double r17886 = r17881 * r17870; double r17887 = r17886 * r17870; double r17888 = r17885 * r17887; double r17889 = r17883 + r17888; double r17890 = r17867 * r17889; double r17891 = fabs(r17890); return r17891; } double f_of(float x) { float r17892 = 1.0f; float r17893 = atan2(1.0, 0.0); float r17894 = r17892 / r17893; float r17895 = sqrt(r17894); float r17896 = 0.047619047619047616f; float r17897 = x; float r17898 = fabs(r17897); float r17899 = r17898 * (r17898 * r17898); float r17900 = r17899 * r17899; float r17901 = r17898 * r17900; float r17902 = r17896 * r17901; float r17903 = 0.6666666666666666f; float r17904 = 3.0f; float r17905 = pow(r17898, r17904); float r17906 = r17903 * r17905; float r17907 = 2.0f; float r17908 = r17907 * r17898; float r17909 = 0.2f; float r17910 = r17898 * r17898; float r17911 = r17910 * r17899; float r17912 = r17909 * r17911; float r17913 = r17908 + r17912; float r17914 = r17906 + r17913; float r17915 = r17902 + r17914; float r17916 = r17895 * r17915; float r17917 = fabs(r17916); return r17917; } double f_od(double x) { double r17918 = 1.0; double r17919 = atan2(1.0, 0.0); double r17920 = r17918 / r17919; double r17921 = sqrt(r17920); double r17922 = 0.047619047619047616; double r17923 = x; double r17924 = fabs(r17923); double r17925 = r17924 * (r17924 * r17924); double r17926 = r17925 * r17925; double r17927 = r17924 * r17926; double r17928 = r17922 * r17927; double r17929 = 0.6666666666666666; double r17930 = 3.0; double r17931 = pow(r17924, r17930); double r17932 = r17929 * r17931; double r17933 = 2.0; double r17934 = r17933 * r17924; double r17935 = 0.2; double r17936 = r17924 * r17924; double r17937 = r17936 * r17925; double r17938 = r17935 * r17937; double r17939 = r17934 + r17938; double r17940 = r17932 + r17939; double r17941 = r17928 + r17940; double r17942 = r17921 * r17941; double r17943 = fabs(r17942); return r17943; } 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 r17944, r17945, r17946, r17947, r17948, r17949, r17950, r17951, r17952, r17953, r17954, r17955, r17956, r17957, r17958, r17959, r17960, r17961, r17962, r17963, r17964, r17965, r17966, r17967, r17968, r17969, r17970, r17971; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17944, "1", 10, MPFR_RNDN); mpfr_init(r17945); mpfr_init(r17946); mpfr_init(r17947); mpfr_init_set_str(r17948, "2", 10, MPFR_RNDN); mpfr_init(r17949); mpfr_init(r17950); mpfr_init(r17951); mpfr_init_set_str(r17952, "3", 10, MPFR_RNDN); mpfr_init(r17953); mpfr_init(r17954); mpfr_init(r17955); mpfr_init(r17956); mpfr_init(r17957); mpfr_init_set_str(r17958, "5", 10, MPFR_RNDN); mpfr_init(r17959); mpfr_init(r17960); mpfr_init(r17961); mpfr_init(r17962); mpfr_init(r17963); mpfr_init_set_str(r17964, "21", 10, MPFR_RNDN); mpfr_init(r17965); mpfr_init(r17966); mpfr_init(r17967); mpfr_init(r17968); mpfr_init(r17969); mpfr_init(r17970); mpfr_init(r17971); } double f_im(double x) { ; mpfr_const_pi(r17945, MPFR_RNDN); mpfr_sqrt(r17946, r17945, MPFR_RNDN); mpfr_div(r17947, r17944, r17946, MPFR_RNDN); ; mpfr_set_d(r17949, x, MPFR_RNDN); mpfr_abs(r17950, r17949, MPFR_RNDN); mpfr_mul(r17951, r17948, r17950, MPFR_RNDN); ; mpfr_div(r17953, r17948, r17952, MPFR_RNDN); mpfr_mul(r17954, r17950, r17950, MPFR_RNDN); mpfr_mul(r17955, r17954, r17950, MPFR_RNDN); mpfr_mul(r17956, r17953, r17955, MPFR_RNDN); mpfr_add(r17957, r17951, r17956, MPFR_RNDN); ; mpfr_div(r17959, r17944, r17958, MPFR_RNDN); mpfr_mul(r17960, r17955, r17950, MPFR_RNDN); mpfr_mul(r17961, r17960, r17950, MPFR_RNDN); mpfr_mul(r17962, r17959, r17961, MPFR_RNDN); mpfr_add(r17963, r17957, r17962, MPFR_RNDN); ; mpfr_div(r17965, r17944, r17964, MPFR_RNDN); mpfr_mul(r17966, r17961, r17950, MPFR_RNDN); mpfr_mul(r17967, r17966, r17950, MPFR_RNDN); mpfr_mul(r17968, r17965, r17967, MPFR_RNDN); mpfr_add(r17969, r17963, r17968, MPFR_RNDN); mpfr_mul(r17970, r17947, r17969, MPFR_RNDN); mpfr_abs(r17971, r17970, MPFR_RNDN); return mpfr_get_d(r17971, MPFR_RNDN); } static mpfr_t r17972, r17973, r17974, r17975, r17976, r17977, r17978, r17979, r17980, r17981, r17982, r17983, r17984, r17985, r17986, r17987, r17988, r17989, r17990, r17991, r17992, r17993, r17994, r17995, r17996, r17997; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17972, "1", 10, MPFR_RNDN); mpfr_init(r17973); mpfr_init(r17974); mpfr_init(r17975); mpfr_init_set_str(r17976, "1/21", 10, MPFR_RNDN); mpfr_init(r17977); mpfr_init(r17978); mpfr_init(r17979); mpfr_init(r17980); mpfr_init(r17981); mpfr_init(r17982); mpfr_init_set_str(r17983, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r17984, "3", 10, MPFR_RNDN); mpfr_init(r17985); mpfr_init(r17986); mpfr_init_set_str(r17987, "2", 10, MPFR_RNDN); mpfr_init(r17988); mpfr_init_set_str(r17989, "1/5", 10, MPFR_RNDN); mpfr_init(r17990); mpfr_init(r17991); mpfr_init(r17992); mpfr_init(r17993); mpfr_init(r17994); mpfr_init(r17995); mpfr_init(r17996); mpfr_init(r17997); } double f_fm(double x) { ; mpfr_const_pi(r17973, MPFR_RNDN); mpfr_div(r17974, r17972, r17973, MPFR_RNDN); mpfr_sqrt(r17975, r17974, MPFR_RNDN); ; mpfr_set_d(r17977, x, MPFR_RNDN); mpfr_abs(r17978, r17977, MPFR_RNDN); mpfr_mul(r17979, r17978, r17978, MPFR_RNDN); mpfr_mul(r17979, r17979, r17978, MPFR_RNDN); mpfr_sqr(r17980, r17979, MPFR_RNDN); mpfr_mul(r17981, r17978, r17980, MPFR_RNDN); mpfr_mul(r17982, r17976, r17981, MPFR_RNDN); ; ; mpfr_pow(r17985, r17978, r17984, MPFR_RNDN); mpfr_mul(r17986, r17983, r17985, MPFR_RNDN); ; mpfr_mul(r17988, r17987, r17978, MPFR_RNDN); ; mpfr_sqr(r17990, r17978, MPFR_RNDN); mpfr_mul(r17991, r17990, r17979, MPFR_RNDN); mpfr_mul(r17992, r17989, r17991, MPFR_RNDN); mpfr_add(r17993, r17988, r17992, MPFR_RNDN); mpfr_add(r17994, r17986, r17993, MPFR_RNDN); mpfr_add(r17995, r17982, r17994, MPFR_RNDN); mpfr_mul(r17996, r17975, r17995, MPFR_RNDN); mpfr_abs(r17997, r17996, MPFR_RNDN); return mpfr_get_d(r17997, MPFR_RNDN); } static mpfr_t r17998, r17999, r18000, r18001, r18002, r18003, r18004, r18005, r18006, r18007, r18008, r18009, r18010, r18011, r18012, r18013, r18014, r18015, r18016, r18017, r18018, r18019, r18020, r18021, r18022, r18023; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17998, "1", 10, MPFR_RNDN); mpfr_init(r17999); mpfr_init(r18000); mpfr_init(r18001); mpfr_init_set_str(r18002, "1/21", 10, MPFR_RNDN); mpfr_init(r18003); mpfr_init(r18004); mpfr_init(r18005); mpfr_init(r18006); mpfr_init(r18007); mpfr_init(r18008); mpfr_init_set_str(r18009, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18010, "3", 10, MPFR_RNDN); mpfr_init(r18011); mpfr_init(r18012); mpfr_init_set_str(r18013, "2", 10, MPFR_RNDN); mpfr_init(r18014); mpfr_init_set_str(r18015, "1/5", 10, MPFR_RNDN); mpfr_init(r18016); mpfr_init(r18017); mpfr_init(r18018); mpfr_init(r18019); mpfr_init(r18020); mpfr_init(r18021); mpfr_init(r18022); mpfr_init(r18023); } double f_dm(double x) { ; mpfr_const_pi(r17999, MPFR_RNDN); mpfr_div(r18000, r17998, r17999, MPFR_RNDN); mpfr_sqrt(r18001, r18000, MPFR_RNDN); ; mpfr_set_d(r18003, x, MPFR_RNDN); mpfr_abs(r18004, r18003, MPFR_RNDN); mpfr_mul(r18005, r18004, r18004, MPFR_RNDN); mpfr_mul(r18005, r18005, r18004, MPFR_RNDN); mpfr_sqr(r18006, r18005, MPFR_RNDN); mpfr_mul(r18007, r18004, r18006, MPFR_RNDN); mpfr_mul(r18008, r18002, r18007, MPFR_RNDN); ; ; mpfr_pow(r18011, r18004, r18010, MPFR_RNDN); mpfr_mul(r18012, r18009, r18011, MPFR_RNDN); ; mpfr_mul(r18014, r18013, r18004, MPFR_RNDN); ; mpfr_sqr(r18016, r18004, MPFR_RNDN); mpfr_mul(r18017, r18016, r18005, MPFR_RNDN); mpfr_mul(r18018, r18015, r18017, MPFR_RNDN); mpfr_add(r18019, r18014, r18018, MPFR_RNDN); mpfr_add(r18020, r18012, r18019, MPFR_RNDN); mpfr_add(r18021, r18008, r18020, MPFR_RNDN); mpfr_mul(r18022, r18001, r18021, MPFR_RNDN); mpfr_abs(r18023, r18022, MPFR_RNDN); return mpfr_get_d(r18023, MPFR_RNDN); }