#include #include #include #include #include char *name = "Jmat.Real.erf"; double f_if(float x) { float r27797 = 1; float r27798 = 0.3275911; float r27799 = x; float r27800 = fabs(r27799); float r27801 = r27798 * r27800; float r27802 = r27797 + r27801; float r27803 = r27797 / r27802; float r27804 = 0.254829592; float r27805 = -0.284496736; float r27806 = 1.421413741; float r27807 = -1.453152027; float r27808 = 1.061405429; float r27809 = r27803 * r27808; float r27810 = r27807 + r27809; float r27811 = r27803 * r27810; float r27812 = r27806 + r27811; float r27813 = r27803 * r27812; float r27814 = r27805 + r27813; float r27815 = r27803 * r27814; float r27816 = r27804 + r27815; float r27817 = r27803 * r27816; float r27818 = r27800 * r27800; float r27819 = -r27818; float r27820 = exp(r27819); float r27821 = r27817 * r27820; float r27822 = r27797 - r27821; return r27822; } double f_id(double x) { double r27823 = 1; double r27824 = 0.3275911; double r27825 = x; double r27826 = fabs(r27825); double r27827 = r27824 * r27826; double r27828 = r27823 + r27827; double r27829 = r27823 / r27828; double r27830 = 0.254829592; double r27831 = -0.284496736; double r27832 = 1.421413741; double r27833 = -1.453152027; double r27834 = 1.061405429; double r27835 = r27829 * r27834; double r27836 = r27833 + r27835; double r27837 = r27829 * r27836; double r27838 = r27832 + r27837; double r27839 = r27829 * r27838; double r27840 = r27831 + r27839; double r27841 = r27829 * r27840; double r27842 = r27830 + r27841; double r27843 = r27829 * r27842; double r27844 = r27826 * r27826; double r27845 = -r27844; double r27846 = exp(r27845); double r27847 = r27843 * r27846; double r27848 = r27823 - r27847; return r27848; } double f_of(float x) { float r27849 = x; float r27850 = fabs(r27849); float r27851 = exp(r27850); float r27852 = -r27850; float r27853 = pow(r27851, r27852); float r27854 = 1.453152027; float r27855 = 0.3275911; float r27856 = 1; float r27857 = fma(r27855, r27850, r27856); float r27858 = 4; float r27859 = pow(r27857, r27858); float r27860 = r27854 / r27859; float r27861 = 1.061405429; float r27862 = 5; float r27863 = pow(r27857, r27862); float r27864 = r27861 / r27863; float r27865 = r27860 - r27864; float r27866 = 0.284496736; float r27867 = r27857 * r27857; float r27868 = r27853 / r27867; float r27869 = fma(r27866, r27868, r27856); float r27870 = fma(r27853, r27865, r27869); float r27871 = 1.421413741; float r27872 = 2; float r27873 = pow(r27850, r27872); float r27874 = exp(r27873); float r27875 = r27855 * r27850; float r27876 = r27875 + r27856; float r27877 = 3; float r27878 = pow(r27876, r27877); float r27879 = r27874 * r27878; float r27880 = r27856 / r27879; float r27881 = r27871 * r27880; float r27882 = 0.254829592; float r27883 = r27874 * r27876; float r27884 = r27856 / r27883; float r27885 = r27882 * r27884; float r27886 = r27881 + r27885; float r27887 = r27870 - r27886; float r27888 = pow(r27887, r27872); float r27889 = 1/3; float r27890 = pow(r27888, r27889); float r27891 = -r27856; float r27892 = r27850 * r27850; float r27893 = exp(r27892); float r27894 = r27891 / r27893; float r27895 = fma(r27850, r27855, r27856); float r27896 = r27861 / r27895; float r27897 = -1.453152027; float r27898 = r27896 + r27897; float r27899 = r27856 / r27895; float r27900 = fma(r27898, r27899, r27871); float r27901 = -0.284496736; float r27902 = fma(r27900, r27899, r27901); float r27903 = r27899 * r27899; float r27904 = r27882 / r27895; float r27905 = fma(r27902, r27903, r27904); float r27906 = fma(r27894, r27905, r27856); float r27907 = cbrt(r27906); float r27908 = r27890 * r27907; return r27908; } double f_od(double x) { double r27909 = x; double r27910 = fabs(r27909); double r27911 = exp(r27910); double r27912 = -r27910; double r27913 = pow(r27911, r27912); double r27914 = 1.453152027; double r27915 = 0.3275911; double r27916 = 1; double r27917 = fma(r27915, r27910, r27916); double r27918 = 4; double r27919 = pow(r27917, r27918); double r27920 = r27914 / r27919; double r27921 = 1.061405429; double r27922 = 5; double r27923 = pow(r27917, r27922); double r27924 = r27921 / r27923; double r27925 = r27920 - r27924; double r27926 = 0.284496736; double r27927 = r27917 * r27917; double r27928 = r27913 / r27927; double r27929 = fma(r27926, r27928, r27916); double r27930 = fma(r27913, r27925, r27929); double r27931 = 1.421413741; double r27932 = 2; double r27933 = pow(r27910, r27932); double r27934 = exp(r27933); double r27935 = r27915 * r27910; double r27936 = r27935 + r27916; double r27937 = 3; double r27938 = pow(r27936, r27937); double r27939 = r27934 * r27938; double r27940 = r27916 / r27939; double r27941 = r27931 * r27940; double r27942 = 0.254829592; double r27943 = r27934 * r27936; double r27944 = r27916 / r27943; double r27945 = r27942 * r27944; double r27946 = r27941 + r27945; double r27947 = r27930 - r27946; double r27948 = pow(r27947, r27932); double r27949 = 1/3; double r27950 = pow(r27948, r27949); double r27951 = -r27916; double r27952 = r27910 * r27910; double r27953 = exp(r27952); double r27954 = r27951 / r27953; double r27955 = fma(r27910, r27915, r27916); double r27956 = r27921 / r27955; double r27957 = -1.453152027; double r27958 = r27956 + r27957; double r27959 = r27916 / r27955; double r27960 = fma(r27958, r27959, r27931); double r27961 = -0.284496736; double r27962 = fma(r27960, r27959, r27961); double r27963 = r27959 * r27959; double r27964 = r27942 / r27955; double r27965 = fma(r27962, r27963, r27964); double r27966 = fma(r27954, r27965, r27916); double r27967 = cbrt(r27966); double r27968 = r27950 * r27967; return r27968; } 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 r27969, r27970, r27971, r27972, r27973, r27974, r27975, r27976, r27977, r27978, r27979, r27980, r27981, r27982, r27983, r27984, r27985, r27986, r27987, r27988, r27989, r27990, r27991, r27992, r27993, r27994; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27969, "1", 10, MPFR_RNDN); mpfr_init_set_str(r27970, "0.3275911", 10, MPFR_RNDN); mpfr_init(r27971); mpfr_init(r27972); mpfr_init(r27973); mpfr_init(r27974); mpfr_init(r27975); mpfr_init_set_str(r27976, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r27977, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r27978, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r27979, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r27980, "1.061405429", 10, MPFR_RNDN); mpfr_init(r27981); mpfr_init(r27982); mpfr_init(r27983); mpfr_init(r27984); mpfr_init(r27985); mpfr_init(r27986); mpfr_init(r27987); mpfr_init(r27988); mpfr_init(r27989); mpfr_init(r27990); mpfr_init(r27991); mpfr_init(r27992); mpfr_init(r27993); mpfr_init(r27994); } double f_im(double x) { ; ; mpfr_set_d(r27971, x, MPFR_RNDN); mpfr_abs(r27972, r27971, MPFR_RNDN); mpfr_mul(r27973, r27970, r27972, MPFR_RNDN); mpfr_add(r27974, r27969, r27973, MPFR_RNDN); mpfr_div(r27975, r27969, r27974, MPFR_RNDN); ; ; ; ; ; mpfr_mul(r27981, r27975, r27980, MPFR_RNDN); mpfr_add(r27982, r27979, r27981, MPFR_RNDN); mpfr_mul(r27983, r27975, r27982, MPFR_RNDN); mpfr_add(r27984, r27978, r27983, MPFR_RNDN); mpfr_mul(r27985, r27975, r27984, MPFR_RNDN); mpfr_add(r27986, r27977, r27985, MPFR_RNDN); mpfr_mul(r27987, r27975, r27986, MPFR_RNDN); mpfr_add(r27988, r27976, r27987, MPFR_RNDN); mpfr_mul(r27989, r27975, r27988, MPFR_RNDN); mpfr_mul(r27990, r27972, r27972, MPFR_RNDN); mpfr_neg(r27991, r27990, MPFR_RNDN); mpfr_exp(r27992, r27991, MPFR_RNDN); mpfr_mul(r27993, r27989, r27992, MPFR_RNDN); mpfr_sub(r27994, r27969, r27993, MPFR_RNDN); return mpfr_get_d(r27994, MPFR_RNDN); } static mpfr_t r27995, r27996, r27997, r27998, r27999, r28000, r28001, r28002, r28003, r28004, r28005, r28006, r28007, r28008, r28009, r28010, r28011, r28012, r28013, r28014, r28015, r28016, r28017, r28018, r28019, r28020, r28021, r28022, r28023, r28024, r28025, r28026, r28027, r28028, r28029, r28030, r28031, r28032, r28033, r28034, r28035, r28036, r28037, r28038, r28039, r28040, r28041, r28042, r28043, r28044, r28045, r28046, r28047, r28048, r28049, r28050, r28051, r28052, r28053, r28054; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r27995); mpfr_init(r27996); mpfr_init(r27997); mpfr_init(r27998); mpfr_init(r27999); mpfr_init_set_str(r28000, "1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r28001, "0.3275911", 10, MPFR_RNDN); mpfr_init_set_str(r28002, "1", 10, MPFR_RNDN); mpfr_init(r28003); mpfr_init_set_str(r28004, "4", 10, MPFR_RNDN); mpfr_init(r28005); mpfr_init(r28006); mpfr_init_set_str(r28007, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r28008, "5", 10, MPFR_RNDN); mpfr_init(r28009); mpfr_init(r28010); mpfr_init(r28011); mpfr_init_set_str(r28012, "0.284496736", 10, MPFR_RNDN); mpfr_init(r28013); mpfr_init(r28014); mpfr_init(r28015); mpfr_init(r28016); mpfr_init_set_str(r28017, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r28018, "2", 10, MPFR_RNDN); mpfr_init(r28019); mpfr_init(r28020); mpfr_init(r28021); mpfr_init(r28022); mpfr_init_set_str(r28023, "3", 10, MPFR_RNDN); mpfr_init(r28024); mpfr_init(r28025); mpfr_init(r28026); mpfr_init(r28027); mpfr_init_set_str(r28028, "0.254829592", 10, MPFR_RNDN); mpfr_init(r28029); mpfr_init(r28030); mpfr_init(r28031); mpfr_init(r28032); mpfr_init(r28033); mpfr_init(r28034); mpfr_init_set_str(r28035, "1/3", 10, MPFR_RNDN); mpfr_init(r28036); mpfr_init(r28037); mpfr_init(r28038); mpfr_init(r28039); mpfr_init(r28040); mpfr_init(r28041); mpfr_init(r28042); mpfr_init_set_str(r28043, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r28044); mpfr_init(r28045); mpfr_init(r28046); mpfr_init_set_str(r28047, "-0.284496736", 10, MPFR_RNDN); mpfr_init(r28048); mpfr_init(r28049); mpfr_init(r28050); mpfr_init(r28051); mpfr_init(r28052); mpfr_init(r28053); mpfr_init(r28054); } double f_fm(double x) { mpfr_set_d(r27995, x, MPFR_RNDN); mpfr_abs(r27996, r27995, MPFR_RNDN); mpfr_exp(r27997, r27996, MPFR_RNDN); mpfr_neg(r27998, r27996, MPFR_RNDN); mpfr_pow(r27999, r27997, r27998, MPFR_RNDN); ; ; ; mpfr_fma(r28003, r28001, r27996, r28002, MPFR_RNDN); ; mpfr_pow(r28005, r28003, r28004, MPFR_RNDN); mpfr_div(r28006, r28000, r28005, MPFR_RNDN); ; ; mpfr_pow(r28009, r28003, r28008, MPFR_RNDN); mpfr_div(r28010, r28007, r28009, MPFR_RNDN); mpfr_sub(r28011, r28006, r28010, MPFR_RNDN); ; mpfr_mul(r28013, r28003, r28003, MPFR_RNDN); mpfr_div(r28014, r27999, r28013, MPFR_RNDN); mpfr_fma(r28015, r28012, r28014, r28002, MPFR_RNDN); mpfr_fma(r28016, r27999, r28011, r28015, MPFR_RNDN); ; ; mpfr_pow(r28019, r27996, r28018, MPFR_RNDN); mpfr_exp(r28020, r28019, MPFR_RNDN); mpfr_mul(r28021, r28001, r27996, MPFR_RNDN); mpfr_add(r28022, r28021, r28002, MPFR_RNDN); ; mpfr_pow(r28024, r28022, r28023, MPFR_RNDN); mpfr_mul(r28025, r28020, r28024, MPFR_RNDN); mpfr_div(r28026, r28002, r28025, MPFR_RNDN); mpfr_mul(r28027, r28017, r28026, MPFR_RNDN); ; mpfr_mul(r28029, r28020, r28022, MPFR_RNDN); mpfr_div(r28030, r28002, r28029, MPFR_RNDN); mpfr_mul(r28031, r28028, r28030, MPFR_RNDN); mpfr_add(r28032, r28027, r28031, MPFR_RNDN); mpfr_sub(r28033, r28016, r28032, MPFR_RNDN); mpfr_pow(r28034, r28033, r28018, MPFR_RNDN); ; mpfr_pow(r28036, r28034, r28035, MPFR_RNDN); mpfr_neg(r28037, r28002, MPFR_RNDN); mpfr_mul(r28038, r27996, r27996, MPFR_RNDN); mpfr_exp(r28039, r28038, MPFR_RNDN); mpfr_div(r28040, r28037, r28039, MPFR_RNDN); mpfr_fma(r28041, r27996, r28001, r28002, MPFR_RNDN); mpfr_div(r28042, r28007, r28041, MPFR_RNDN); ; mpfr_add(r28044, r28042, r28043, MPFR_RNDN); mpfr_div(r28045, r28002, r28041, MPFR_RNDN); mpfr_fma(r28046, r28044, r28045, r28017, MPFR_RNDN); ; mpfr_fma(r28048, r28046, r28045, r28047, MPFR_RNDN); mpfr_mul(r28049, r28045, r28045, MPFR_RNDN); mpfr_div(r28050, r28028, r28041, MPFR_RNDN); mpfr_fma(r28051, r28048, r28049, r28050, MPFR_RNDN); mpfr_fma(r28052, r28040, r28051, r28002, MPFR_RNDN); mpfr_cbrt(r28053, r28052, MPFR_RNDN); mpfr_mul(r28054, r28036, r28053, MPFR_RNDN); return mpfr_get_d(r28054, MPFR_RNDN); } static mpfr_t r28055, r28056, r28057, r28058, r28059, r28060, r28061, r28062, r28063, r28064, r28065, r28066, r28067, r28068, r28069, r28070, r28071, r28072, r28073, r28074, r28075, r28076, r28077, r28078, r28079, r28080, r28081, r28082, r28083, r28084, r28085, r28086, r28087, r28088, r28089, r28090, r28091, r28092, r28093, r28094, r28095, r28096, r28097, r28098, r28099, r28100, r28101, r28102, r28103, r28104, r28105, r28106, r28107, r28108, r28109, r28110, r28111, r28112, r28113, r28114; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r28055); mpfr_init(r28056); mpfr_init(r28057); mpfr_init(r28058); mpfr_init(r28059); mpfr_init_set_str(r28060, "1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r28061, "0.3275911", 10, MPFR_RNDN); mpfr_init_set_str(r28062, "1", 10, MPFR_RNDN); mpfr_init(r28063); mpfr_init_set_str(r28064, "4", 10, MPFR_RNDN); mpfr_init(r28065); mpfr_init(r28066); mpfr_init_set_str(r28067, "1.061405429", 10, MPFR_RNDN); mpfr_init_set_str(r28068, "5", 10, MPFR_RNDN); mpfr_init(r28069); mpfr_init(r28070); mpfr_init(r28071); mpfr_init_set_str(r28072, "0.284496736", 10, MPFR_RNDN); mpfr_init(r28073); mpfr_init(r28074); mpfr_init(r28075); mpfr_init(r28076); mpfr_init_set_str(r28077, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r28078, "2", 10, MPFR_RNDN); mpfr_init(r28079); mpfr_init(r28080); mpfr_init(r28081); mpfr_init(r28082); mpfr_init_set_str(r28083, "3", 10, MPFR_RNDN); mpfr_init(r28084); mpfr_init(r28085); mpfr_init(r28086); mpfr_init(r28087); mpfr_init_set_str(r28088, "0.254829592", 10, MPFR_RNDN); mpfr_init(r28089); mpfr_init(r28090); mpfr_init(r28091); mpfr_init(r28092); mpfr_init(r28093); mpfr_init(r28094); mpfr_init_set_str(r28095, "1/3", 10, MPFR_RNDN); mpfr_init(r28096); mpfr_init(r28097); mpfr_init(r28098); mpfr_init(r28099); mpfr_init(r28100); mpfr_init(r28101); mpfr_init(r28102); mpfr_init_set_str(r28103, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r28104); mpfr_init(r28105); mpfr_init(r28106); mpfr_init_set_str(r28107, "-0.284496736", 10, MPFR_RNDN); mpfr_init(r28108); mpfr_init(r28109); mpfr_init(r28110); mpfr_init(r28111); mpfr_init(r28112); mpfr_init(r28113); mpfr_init(r28114); } double f_dm(double x) { mpfr_set_d(r28055, x, MPFR_RNDN); mpfr_abs(r28056, r28055, MPFR_RNDN); mpfr_exp(r28057, r28056, MPFR_RNDN); mpfr_neg(r28058, r28056, MPFR_RNDN); mpfr_pow(r28059, r28057, r28058, MPFR_RNDN); ; ; ; mpfr_fma(r28063, r28061, r28056, r28062, MPFR_RNDN); ; mpfr_pow(r28065, r28063, r28064, MPFR_RNDN); mpfr_div(r28066, r28060, r28065, MPFR_RNDN); ; ; mpfr_pow(r28069, r28063, r28068, MPFR_RNDN); mpfr_div(r28070, r28067, r28069, MPFR_RNDN); mpfr_sub(r28071, r28066, r28070, MPFR_RNDN); ; mpfr_mul(r28073, r28063, r28063, MPFR_RNDN); mpfr_div(r28074, r28059, r28073, MPFR_RNDN); mpfr_fma(r28075, r28072, r28074, r28062, MPFR_RNDN); mpfr_fma(r28076, r28059, r28071, r28075, MPFR_RNDN); ; ; mpfr_pow(r28079, r28056, r28078, MPFR_RNDN); mpfr_exp(r28080, r28079, MPFR_RNDN); mpfr_mul(r28081, r28061, r28056, MPFR_RNDN); mpfr_add(r28082, r28081, r28062, MPFR_RNDN); ; mpfr_pow(r28084, r28082, r28083, MPFR_RNDN); mpfr_mul(r28085, r28080, r28084, MPFR_RNDN); mpfr_div(r28086, r28062, r28085, MPFR_RNDN); mpfr_mul(r28087, r28077, r28086, MPFR_RNDN); ; mpfr_mul(r28089, r28080, r28082, MPFR_RNDN); mpfr_div(r28090, r28062, r28089, MPFR_RNDN); mpfr_mul(r28091, r28088, r28090, MPFR_RNDN); mpfr_add(r28092, r28087, r28091, MPFR_RNDN); mpfr_sub(r28093, r28076, r28092, MPFR_RNDN); mpfr_pow(r28094, r28093, r28078, MPFR_RNDN); ; mpfr_pow(r28096, r28094, r28095, MPFR_RNDN); mpfr_neg(r28097, r28062, MPFR_RNDN); mpfr_mul(r28098, r28056, r28056, MPFR_RNDN); mpfr_exp(r28099, r28098, MPFR_RNDN); mpfr_div(r28100, r28097, r28099, MPFR_RNDN); mpfr_fma(r28101, r28056, r28061, r28062, MPFR_RNDN); mpfr_div(r28102, r28067, r28101, MPFR_RNDN); ; mpfr_add(r28104, r28102, r28103, MPFR_RNDN); mpfr_div(r28105, r28062, r28101, MPFR_RNDN); mpfr_fma(r28106, r28104, r28105, r28077, MPFR_RNDN); ; mpfr_fma(r28108, r28106, r28105, r28107, MPFR_RNDN); mpfr_mul(r28109, r28105, r28105, MPFR_RNDN); mpfr_div(r28110, r28088, r28101, MPFR_RNDN); mpfr_fma(r28111, r28108, r28109, r28110, MPFR_RNDN); mpfr_fma(r28112, r28100, r28111, r28062, MPFR_RNDN); mpfr_cbrt(r28113, r28112, MPFR_RNDN); mpfr_mul(r28114, r28096, r28113, MPFR_RNDN); return mpfr_get_d(r28114, MPFR_RNDN); }