#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 r18820 = 1.0f; float r18821 = atan2(1.0, 0.0); float r18822 = sqrt(r18821); float r18823 = r18820 / r18822; float r18824 = 2.0f; float r18825 = x; float r18826 = fabs(r18825); float r18827 = r18824 * r18826; float r18828 = 3.0f; float r18829 = r18824 / r18828; float r18830 = r18826 * r18826; float r18831 = r18830 * r18826; float r18832 = r18829 * r18831; float r18833 = r18827 + r18832; float r18834 = 5.0f; float r18835 = r18820 / r18834; float r18836 = r18831 * r18826; float r18837 = r18836 * r18826; float r18838 = r18835 * r18837; float r18839 = r18833 + r18838; float r18840 = 21.0f; float r18841 = r18820 / r18840; float r18842 = r18837 * r18826; float r18843 = r18842 * r18826; float r18844 = r18841 * r18843; float r18845 = r18839 + r18844; float r18846 = r18823 * r18845; float r18847 = fabs(r18846); return r18847; } double f_id(double x) { double r18848 = 1.0; double r18849 = atan2(1.0, 0.0); double r18850 = sqrt(r18849); double r18851 = r18848 / r18850; double r18852 = 2.0; double r18853 = x; double r18854 = fabs(r18853); double r18855 = r18852 * r18854; double r18856 = 3.0; double r18857 = r18852 / r18856; double r18858 = r18854 * r18854; double r18859 = r18858 * r18854; double r18860 = r18857 * r18859; double r18861 = r18855 + r18860; double r18862 = 5.0; double r18863 = r18848 / r18862; double r18864 = r18859 * r18854; double r18865 = r18864 * r18854; double r18866 = r18863 * r18865; double r18867 = r18861 + r18866; double r18868 = 21.0; double r18869 = r18848 / r18868; double r18870 = r18865 * r18854; double r18871 = r18870 * r18854; double r18872 = r18869 * r18871; double r18873 = r18867 + r18872; double r18874 = r18851 * r18873; double r18875 = fabs(r18874); return r18875; } double f_of(float x) { float r18876 = 2.0f; float r18877 = x; float r18878 = fabs(r18877); float r18879 = r18876 * r18878; float r18880 = 3.0f; float r18881 = r18876 / r18880; float r18882 = r18878 * r18881; float r18883 = r18878 * r18878; float r18884 = r18882 * r18883; float r18885 = r18879 + r18884; float r18886 = r18878 * (r18878 * r18878); float r18887 = r18886 * r18886; float r18888 = 1.0f; float r18889 = pow(r18887, r18888); float r18890 = 21.0f; float r18891 = r18890 / r18878; float r18892 = r18889 / r18891; float r18893 = r18886 * r18883; float r18894 = 5.0f; float r18895 = r18893 / r18894; float r18896 = r18892 + r18895; float r18897 = r18885 + r18896; float r18898 = atan2(1.0, 0.0); float r18899 = sqrt(r18898); float r18900 = r18897 / r18899; float r18901 = fabs(r18900); return r18901; } double f_od(double x) { double r18902 = 2.0; double r18903 = x; double r18904 = fabs(r18903); double r18905 = r18902 * r18904; double r18906 = 3.0; double r18907 = r18902 / r18906; double r18908 = r18904 * r18907; double r18909 = r18904 * r18904; double r18910 = r18908 * r18909; double r18911 = r18905 + r18910; double r18912 = r18904 * (r18904 * r18904); double r18913 = r18912 * r18912; double r18914 = 1.0; double r18915 = pow(r18913, r18914); double r18916 = 21.0; double r18917 = r18916 / r18904; double r18918 = r18915 / r18917; double r18919 = r18912 * r18909; double r18920 = 5.0; double r18921 = r18919 / r18920; double r18922 = r18918 + r18921; double r18923 = r18911 + r18922; double r18924 = atan2(1.0, 0.0); double r18925 = sqrt(r18924); double r18926 = r18923 / r18925; double r18927 = fabs(r18926); return r18927; } 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 r18928, r18929, r18930, r18931, r18932, r18933, r18934, r18935, r18936, r18937, r18938, r18939, r18940, r18941, r18942, r18943, r18944, r18945, r18946, r18947, r18948, r18949, r18950, r18951, r18952, r18953, r18954, r18955; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18928, "1", 10, MPFR_RNDN); mpfr_init(r18929); mpfr_init(r18930); mpfr_init(r18931); mpfr_init_set_str(r18932, "2", 10, MPFR_RNDN); mpfr_init(r18933); mpfr_init(r18934); mpfr_init(r18935); mpfr_init_set_str(r18936, "3", 10, MPFR_RNDN); mpfr_init(r18937); mpfr_init(r18938); mpfr_init(r18939); mpfr_init(r18940); mpfr_init(r18941); mpfr_init_set_str(r18942, "5", 10, MPFR_RNDN); mpfr_init(r18943); mpfr_init(r18944); mpfr_init(r18945); mpfr_init(r18946); mpfr_init(r18947); mpfr_init_set_str(r18948, "21", 10, MPFR_RNDN); mpfr_init(r18949); mpfr_init(r18950); mpfr_init(r18951); mpfr_init(r18952); mpfr_init(r18953); mpfr_init(r18954); mpfr_init(r18955); } double f_im(double x) { ; mpfr_const_pi(r18929, MPFR_RNDN); mpfr_sqrt(r18930, r18929, MPFR_RNDN); mpfr_div(r18931, r18928, r18930, MPFR_RNDN); ; mpfr_set_d(r18933, x, MPFR_RNDN); mpfr_abs(r18934, r18933, MPFR_RNDN); mpfr_mul(r18935, r18932, r18934, MPFR_RNDN); ; mpfr_div(r18937, r18932, r18936, MPFR_RNDN); mpfr_mul(r18938, r18934, r18934, MPFR_RNDN); mpfr_mul(r18939, r18938, r18934, MPFR_RNDN); mpfr_mul(r18940, r18937, r18939, MPFR_RNDN); mpfr_add(r18941, r18935, r18940, MPFR_RNDN); ; mpfr_div(r18943, r18928, r18942, MPFR_RNDN); mpfr_mul(r18944, r18939, r18934, MPFR_RNDN); mpfr_mul(r18945, r18944, r18934, MPFR_RNDN); mpfr_mul(r18946, r18943, r18945, MPFR_RNDN); mpfr_add(r18947, r18941, r18946, MPFR_RNDN); ; mpfr_div(r18949, r18928, r18948, MPFR_RNDN); mpfr_mul(r18950, r18945, r18934, MPFR_RNDN); mpfr_mul(r18951, r18950, r18934, MPFR_RNDN); mpfr_mul(r18952, r18949, r18951, MPFR_RNDN); mpfr_add(r18953, r18947, r18952, MPFR_RNDN); mpfr_mul(r18954, r18931, r18953, MPFR_RNDN); mpfr_abs(r18955, r18954, MPFR_RNDN); return mpfr_get_d(r18955, MPFR_RNDN); } static mpfr_t 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18956, "2", 10, MPFR_RNDN); mpfr_init(r18957); mpfr_init(r18958); mpfr_init(r18959); mpfr_init_set_str(r18960, "3", 10, MPFR_RNDN); mpfr_init(r18961); mpfr_init(r18962); mpfr_init(r18963); mpfr_init(r18964); mpfr_init(r18965); mpfr_init(r18966); mpfr_init(r18967); mpfr_init_set_str(r18968, "1", 10, MPFR_RNDN); mpfr_init(r18969); mpfr_init_set_str(r18970, "21", 10, MPFR_RNDN); mpfr_init(r18971); mpfr_init(r18972); mpfr_init(r18973); mpfr_init_set_str(r18974, "5", 10, MPFR_RNDN); mpfr_init(r18975); mpfr_init(r18976); mpfr_init(r18977); mpfr_init(r18978); mpfr_init(r18979); mpfr_init(r18980); mpfr_init(r18981); } double f_fm(double x) { ; mpfr_set_d(r18957, x, MPFR_RNDN); mpfr_abs(r18958, r18957, MPFR_RNDN); mpfr_mul(r18959, r18956, r18958, MPFR_RNDN); ; mpfr_div(r18961, r18956, r18960, MPFR_RNDN); mpfr_mul(r18962, r18958, r18961, MPFR_RNDN); mpfr_sqr(r18963, r18958, MPFR_RNDN); mpfr_mul(r18964, r18962, r18963, MPFR_RNDN); mpfr_add(r18965, r18959, r18964, MPFR_RNDN); mpfr_mul(r18966, r18958, r18958, MPFR_RNDN); mpfr_mul(r18966, r18966, r18958, MPFR_RNDN); mpfr_sqr(r18967, r18966, MPFR_RNDN); ; mpfr_pow(r18969, r18967, r18968, MPFR_RNDN); ; mpfr_div(r18971, r18970, r18958, MPFR_RNDN); mpfr_div(r18972, r18969, r18971, MPFR_RNDN); mpfr_mul(r18973, r18966, r18963, MPFR_RNDN); ; mpfr_div(r18975, r18973, r18974, MPFR_RNDN); mpfr_add(r18976, r18972, r18975, MPFR_RNDN); mpfr_add(r18977, r18965, r18976, MPFR_RNDN); mpfr_const_pi(r18978, MPFR_RNDN); mpfr_sqrt(r18979, r18978, MPFR_RNDN); mpfr_div(r18980, r18977, r18979, MPFR_RNDN); mpfr_abs(r18981, r18980, MPFR_RNDN); return mpfr_get_d(r18981, MPFR_RNDN); } static mpfr_t r18982, r18983, r18984, r18985, r18986, r18987, r18988, r18989, r18990, r18991, r18992, r18993, r18994, r18995, r18996, r18997, r18998, r18999, r19000, r19001, r19002, r19003, r19004, r19005, r19006, r19007; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18982, "2", 10, MPFR_RNDN); mpfr_init(r18983); mpfr_init(r18984); mpfr_init(r18985); mpfr_init_set_str(r18986, "3", 10, MPFR_RNDN); mpfr_init(r18987); mpfr_init(r18988); mpfr_init(r18989); mpfr_init(r18990); mpfr_init(r18991); mpfr_init(r18992); mpfr_init(r18993); mpfr_init_set_str(r18994, "1", 10, MPFR_RNDN); mpfr_init(r18995); mpfr_init_set_str(r18996, "21", 10, MPFR_RNDN); mpfr_init(r18997); mpfr_init(r18998); mpfr_init(r18999); mpfr_init_set_str(r19000, "5", 10, MPFR_RNDN); mpfr_init(r19001); mpfr_init(r19002); mpfr_init(r19003); mpfr_init(r19004); mpfr_init(r19005); mpfr_init(r19006); mpfr_init(r19007); } double f_dm(double x) { ; mpfr_set_d(r18983, x, MPFR_RNDN); mpfr_abs(r18984, r18983, MPFR_RNDN); mpfr_mul(r18985, r18982, r18984, MPFR_RNDN); ; mpfr_div(r18987, r18982, r18986, MPFR_RNDN); mpfr_mul(r18988, r18984, r18987, MPFR_RNDN); mpfr_sqr(r18989, r18984, MPFR_RNDN); mpfr_mul(r18990, r18988, r18989, MPFR_RNDN); mpfr_add(r18991, r18985, r18990, MPFR_RNDN); mpfr_mul(r18992, r18984, r18984, MPFR_RNDN); mpfr_mul(r18992, r18992, r18984, MPFR_RNDN); mpfr_sqr(r18993, r18992, MPFR_RNDN); ; mpfr_pow(r18995, r18993, r18994, MPFR_RNDN); ; mpfr_div(r18997, r18996, r18984, MPFR_RNDN); mpfr_div(r18998, r18995, r18997, MPFR_RNDN); mpfr_mul(r18999, r18992, r18989, MPFR_RNDN); ; mpfr_div(r19001, r18999, r19000, MPFR_RNDN); mpfr_add(r19002, r18998, r19001, MPFR_RNDN); mpfr_add(r19003, r18991, r19002, MPFR_RNDN); mpfr_const_pi(r19004, MPFR_RNDN); mpfr_sqrt(r19005, r19004, MPFR_RNDN); mpfr_div(r19006, r19003, r19005, MPFR_RNDN); mpfr_abs(r19007, r19006, MPFR_RNDN); return mpfr_get_d(r19007, MPFR_RNDN); }