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