#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 r18864 = 1.0f; float r18865 = atan2(1.0, 0.0); float r18866 = sqrt(r18865); float r18867 = r18864 / r18866; float r18868 = 2.0f; float r18869 = x; float r18870 = fabs(r18869); float r18871 = r18868 * r18870; float r18872 = 3.0f; float r18873 = r18868 / r18872; float r18874 = r18870 * r18870; float r18875 = r18874 * r18870; float r18876 = r18873 * r18875; float r18877 = r18871 + r18876; float r18878 = 5.0f; float r18879 = r18864 / r18878; float r18880 = r18875 * r18870; float r18881 = r18880 * r18870; float r18882 = r18879 * r18881; float r18883 = r18877 + r18882; float r18884 = 21.0f; float r18885 = r18864 / r18884; float r18886 = r18881 * r18870; float r18887 = r18886 * r18870; float r18888 = r18885 * r18887; float r18889 = r18883 + r18888; float r18890 = r18867 * r18889; float r18891 = fabs(r18890); return r18891; } double f_id(double x) { double r18892 = 1.0; double r18893 = atan2(1.0, 0.0); double r18894 = sqrt(r18893); double r18895 = r18892 / r18894; double r18896 = 2.0; double r18897 = x; double r18898 = fabs(r18897); double r18899 = r18896 * r18898; double r18900 = 3.0; double r18901 = r18896 / r18900; double r18902 = r18898 * r18898; double r18903 = r18902 * r18898; double r18904 = r18901 * r18903; double r18905 = r18899 + r18904; double r18906 = 5.0; double r18907 = r18892 / r18906; double r18908 = r18903 * r18898; double r18909 = r18908 * r18898; double r18910 = r18907 * r18909; double r18911 = r18905 + r18910; double r18912 = 21.0; double r18913 = r18892 / r18912; double r18914 = r18909 * r18898; double r18915 = r18914 * r18898; double r18916 = r18913 * r18915; double r18917 = r18911 + r18916; double r18918 = r18895 * r18917; double r18919 = fabs(r18918); return r18919; } double f_of(float x) { float r18920 = 2.0f; float r18921 = x; float r18922 = fabs(r18921); float r18923 = r18920 * r18922; float r18924 = 3.0f; float r18925 = r18920 / r18924; float r18926 = r18922 * r18925; float r18927 = r18922 * r18922; float r18928 = r18926 * r18927; float r18929 = r18923 + r18928; float r18930 = r18922 * (r18922 * r18922); float r18931 = r18930 * r18930; float r18932 = 21.0f; float r18933 = r18932 / r18922; float r18934 = r18931 / r18933; float r18935 = r18930 * r18927; float r18936 = 5.0f; float r18937 = r18935 / r18936; float r18938 = r18934 + r18937; float r18939 = r18929 + r18938; float r18940 = 1.0f; float r18941 = atan2(1.0, 0.0); float r18942 = sqrt(r18941); float r18943 = r18940 / r18942; float r18944 = r18939 * r18943; float r18945 = fabs(r18944); return r18945; } double f_od(double x) { double r18946 = 2.0; double r18947 = x; double r18948 = fabs(r18947); double r18949 = r18946 * r18948; double r18950 = 3.0; double r18951 = r18946 / r18950; double r18952 = r18948 * r18951; double r18953 = r18948 * r18948; double r18954 = r18952 * r18953; double r18955 = r18949 + r18954; double r18956 = r18948 * (r18948 * r18948); double r18957 = r18956 * r18956; double r18958 = 21.0; double r18959 = r18958 / r18948; double r18960 = r18957 / r18959; double r18961 = r18956 * r18953; double r18962 = 5.0; double r18963 = r18961 / r18962; double r18964 = r18960 + r18963; double r18965 = r18955 + r18964; double r18966 = 1.0; double r18967 = atan2(1.0, 0.0); double r18968 = sqrt(r18967); double r18969 = r18966 / r18968; double r18970 = r18965 * r18969; double r18971 = fabs(r18970); return r18971; } 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 r18972, r18973, r18974, r18975, r18976, r18977, r18978, r18979, r18980, r18981, r18982, r18983, r18984, r18985, r18986, r18987, r18988, r18989, r18990, r18991, r18992, r18993, r18994, r18995, r18996, r18997, r18998, r18999; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18972, "1", 10, MPFR_RNDN); mpfr_init(r18973); mpfr_init(r18974); mpfr_init(r18975); mpfr_init_set_str(r18976, "2", 10, MPFR_RNDN); mpfr_init(r18977); mpfr_init(r18978); mpfr_init(r18979); mpfr_init_set_str(r18980, "3", 10, MPFR_RNDN); mpfr_init(r18981); mpfr_init(r18982); mpfr_init(r18983); mpfr_init(r18984); mpfr_init(r18985); mpfr_init_set_str(r18986, "5", 10, MPFR_RNDN); mpfr_init(r18987); mpfr_init(r18988); mpfr_init(r18989); mpfr_init(r18990); mpfr_init(r18991); mpfr_init_set_str(r18992, "21", 10, MPFR_RNDN); mpfr_init(r18993); mpfr_init(r18994); mpfr_init(r18995); mpfr_init(r18996); mpfr_init(r18997); mpfr_init(r18998); mpfr_init(r18999); } double f_im(double x) { ; mpfr_const_pi(r18973, MPFR_RNDN); mpfr_sqrt(r18974, r18973, MPFR_RNDN); mpfr_div(r18975, r18972, r18974, MPFR_RNDN); ; mpfr_set_d(r18977, x, MPFR_RNDN); mpfr_abs(r18978, r18977, MPFR_RNDN); mpfr_mul(r18979, r18976, r18978, MPFR_RNDN); ; mpfr_div(r18981, r18976, r18980, MPFR_RNDN); mpfr_mul(r18982, r18978, r18978, MPFR_RNDN); mpfr_mul(r18983, r18982, r18978, MPFR_RNDN); mpfr_mul(r18984, r18981, r18983, MPFR_RNDN); mpfr_add(r18985, r18979, r18984, MPFR_RNDN); ; mpfr_div(r18987, r18972, r18986, MPFR_RNDN); mpfr_mul(r18988, r18983, r18978, MPFR_RNDN); mpfr_mul(r18989, r18988, r18978, MPFR_RNDN); mpfr_mul(r18990, r18987, r18989, MPFR_RNDN); mpfr_add(r18991, r18985, r18990, MPFR_RNDN); ; mpfr_div(r18993, r18972, r18992, MPFR_RNDN); mpfr_mul(r18994, r18989, r18978, MPFR_RNDN); mpfr_mul(r18995, r18994, r18978, MPFR_RNDN); mpfr_mul(r18996, r18993, r18995, MPFR_RNDN); mpfr_add(r18997, r18991, r18996, MPFR_RNDN); mpfr_mul(r18998, r18975, r18997, MPFR_RNDN); mpfr_abs(r18999, r18998, MPFR_RNDN); return mpfr_get_d(r18999, MPFR_RNDN); } 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19000, "2", 10, MPFR_RNDN); mpfr_init(r19001); mpfr_init(r19002); mpfr_init(r19003); mpfr_init_set_str(r19004, "3", 10, MPFR_RNDN); mpfr_init(r19005); mpfr_init(r19006); mpfr_init(r19007); mpfr_init(r19008); mpfr_init(r19009); mpfr_init(r19010); mpfr_init(r19011); mpfr_init_set_str(r19012, "21", 10, MPFR_RNDN); mpfr_init(r19013); mpfr_init(r19014); mpfr_init(r19015); mpfr_init_set_str(r19016, "5", 10, MPFR_RNDN); mpfr_init(r19017); mpfr_init(r19018); mpfr_init(r19019); mpfr_init_set_str(r19020, "1", 10, MPFR_RNDN); mpfr_init(r19021); mpfr_init(r19022); mpfr_init(r19023); mpfr_init(r19024); mpfr_init(r19025); } double f_fm(double x) { ; mpfr_set_d(r19001, x, MPFR_RNDN); mpfr_abs(r19002, r19001, MPFR_RNDN); mpfr_mul(r19003, r19000, r19002, MPFR_RNDN); ; mpfr_div(r19005, r19000, r19004, MPFR_RNDN); mpfr_mul(r19006, r19002, r19005, MPFR_RNDN); mpfr_sqr(r19007, r19002, MPFR_RNDN); mpfr_mul(r19008, r19006, r19007, MPFR_RNDN); mpfr_add(r19009, r19003, r19008, MPFR_RNDN); mpfr_mul(r19010, r19002, r19002, MPFR_RNDN); mpfr_mul(r19010, r19010, r19002, MPFR_RNDN); mpfr_sqr(r19011, r19010, MPFR_RNDN); ; mpfr_div(r19013, r19012, r19002, MPFR_RNDN); mpfr_div(r19014, r19011, r19013, MPFR_RNDN); mpfr_mul(r19015, r19010, r19007, MPFR_RNDN); ; mpfr_div(r19017, r19015, r19016, MPFR_RNDN); mpfr_add(r19018, r19014, r19017, MPFR_RNDN); mpfr_add(r19019, r19009, r19018, MPFR_RNDN); ; mpfr_const_pi(r19021, MPFR_RNDN); mpfr_sqrt(r19022, r19021, MPFR_RNDN); mpfr_div(r19023, r19020, r19022, MPFR_RNDN); mpfr_mul(r19024, r19019, r19023, MPFR_RNDN); mpfr_abs(r19025, r19024, MPFR_RNDN); return mpfr_get_d(r19025, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19026, "2", 10, MPFR_RNDN); mpfr_init(r19027); mpfr_init(r19028); mpfr_init(r19029); mpfr_init_set_str(r19030, "3", 10, MPFR_RNDN); mpfr_init(r19031); mpfr_init(r19032); mpfr_init(r19033); mpfr_init(r19034); mpfr_init(r19035); mpfr_init(r19036); mpfr_init(r19037); mpfr_init_set_str(r19038, "21", 10, MPFR_RNDN); mpfr_init(r19039); mpfr_init(r19040); mpfr_init(r19041); mpfr_init_set_str(r19042, "5", 10, MPFR_RNDN); mpfr_init(r19043); mpfr_init(r19044); mpfr_init(r19045); mpfr_init_set_str(r19046, "1", 10, MPFR_RNDN); mpfr_init(r19047); mpfr_init(r19048); mpfr_init(r19049); mpfr_init(r19050); mpfr_init(r19051); } double f_dm(double x) { ; mpfr_set_d(r19027, x, MPFR_RNDN); mpfr_abs(r19028, r19027, MPFR_RNDN); mpfr_mul(r19029, r19026, r19028, MPFR_RNDN); ; mpfr_div(r19031, r19026, r19030, MPFR_RNDN); mpfr_mul(r19032, r19028, r19031, MPFR_RNDN); mpfr_sqr(r19033, r19028, MPFR_RNDN); mpfr_mul(r19034, r19032, r19033, MPFR_RNDN); mpfr_add(r19035, r19029, r19034, MPFR_RNDN); mpfr_mul(r19036, r19028, r19028, MPFR_RNDN); mpfr_mul(r19036, r19036, r19028, MPFR_RNDN); mpfr_sqr(r19037, r19036, MPFR_RNDN); ; mpfr_div(r19039, r19038, r19028, MPFR_RNDN); mpfr_div(r19040, r19037, r19039, MPFR_RNDN); mpfr_mul(r19041, r19036, r19033, MPFR_RNDN); ; mpfr_div(r19043, r19041, r19042, MPFR_RNDN); mpfr_add(r19044, r19040, r19043, MPFR_RNDN); mpfr_add(r19045, r19035, r19044, MPFR_RNDN); ; mpfr_const_pi(r19047, MPFR_RNDN); mpfr_sqrt(r19048, r19047, MPFR_RNDN); mpfr_div(r19049, r19046, r19048, MPFR_RNDN); mpfr_mul(r19050, r19045, r19049, MPFR_RNDN); mpfr_abs(r19051, r19050, MPFR_RNDN); return mpfr_get_d(r19051, MPFR_RNDN); }