#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Jmat.Real.erf";

double f_if(float x) {
        float r22850 = 1;
        float r22851 = 0.3275911;
        float r22852 = x;
        float r22853 = fabs(r22852);
        float r22854 = r22851 * r22853;
        float r22855 = r22850 + r22854;
        float r22856 = r22850 / r22855;
        float r22857 = 0.254829592;
        float r22858 = -0.284496736;
        float r22859 = 1.421413741;
        float r22860 = -1.453152027;
        float r22861 = 1.061405429;
        float r22862 = r22856 * r22861;
        float r22863 = r22860 + r22862;
        float r22864 = r22856 * r22863;
        float r22865 = r22859 + r22864;
        float r22866 = r22856 * r22865;
        float r22867 = r22858 + r22866;
        float r22868 = r22856 * r22867;
        float r22869 = r22857 + r22868;
        float r22870 = r22856 * r22869;
        float r22871 = r22853 * r22853;
        float r22872 = -r22871;
        float r22873 = exp(r22872);
        float r22874 = r22870 * r22873;
        float r22875 = r22850 - r22874;
        return r22875;
}

double f_id(double x) {
        double r22876 = 1;
        double r22877 = 0.3275911;
        double r22878 = x;
        double r22879 = fabs(r22878);
        double r22880 = r22877 * r22879;
        double r22881 = r22876 + r22880;
        double r22882 = r22876 / r22881;
        double r22883 = 0.254829592;
        double r22884 = -0.284496736;
        double r22885 = 1.421413741;
        double r22886 = -1.453152027;
        double r22887 = 1.061405429;
        double r22888 = r22882 * r22887;
        double r22889 = r22886 + r22888;
        double r22890 = r22882 * r22889;
        double r22891 = r22885 + r22890;
        double r22892 = r22882 * r22891;
        double r22893 = r22884 + r22892;
        double r22894 = r22882 * r22893;
        double r22895 = r22883 + r22894;
        double r22896 = r22882 * r22895;
        double r22897 = r22879 * r22879;
        double r22898 = -r22897;
        double r22899 = exp(r22898);
        double r22900 = r22896 * r22899;
        double r22901 = r22876 - r22900;
        return r22901;
}


double f_of(float x) {
        float r22902 = 1;
        float r22903 = 0.3275911;
        float r22904 = x;
        float r22905 = fabs(r22904);
        float r22906 = r22903 * r22905;
        float r22907 = r22902 + r22906;
        float r22908 = r22902 / r22907;
        float r22909 = 0.254829592;
        float r22910 = -0.284496736;
        float r22911 = 1.421413741;
        float r22912 = 1.061405429;
        float r22913 = r22905 * r22903;
        float r22914 = r22913 + r22902;
        float r22915 = r22912 / r22914;
        float r22916 = -1.453152027;
        float r22917 = r22915 + r22916;
        float r22918 = r22917 / r22914;
        float r22919 = pow(r22918, r22902);
        float r22920 = r22911 + r22919;
        float r22921 = r22908 * r22920;
        float r22922 = r22910 + r22921;
        float r22923 = r22908 * r22922;
        float r22924 = r22909 + r22923;
        float r22925 = r22908 * r22924;
        float r22926 = r22905 * r22905;
        float r22927 = -r22926;
        float r22928 = exp(r22927);
        float r22929 = r22925 * r22928;
        float r22930 = r22902 - r22929;
        return r22930;
}

double f_od(double x) {
        double r22931 = 1;
        double r22932 = 0.3275911;
        double r22933 = x;
        double r22934 = fabs(r22933);
        double r22935 = r22932 * r22934;
        double r22936 = r22931 + r22935;
        double r22937 = r22931 / r22936;
        double r22938 = 0.254829592;
        double r22939 = -0.284496736;
        double r22940 = 1.421413741;
        double r22941 = 1.061405429;
        double r22942 = r22934 * r22932;
        double r22943 = r22942 + r22931;
        double r22944 = r22941 / r22943;
        double r22945 = -1.453152027;
        double r22946 = r22944 + r22945;
        double r22947 = r22946 / r22943;
        double r22948 = pow(r22947, r22931);
        double r22949 = r22940 + r22948;
        double r22950 = r22937 * r22949;
        double r22951 = r22939 + r22950;
        double r22952 = r22937 * r22951;
        double r22953 = r22938 + r22952;
        double r22954 = r22937 * r22953;
        double r22955 = r22934 * r22934;
        double r22956 = -r22955;
        double r22957 = exp(r22956);
        double r22958 = r22954 * r22957;
        double r22959 = r22931 - r22958;
        return r22959;
}

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 r22960, r22961, r22962, r22963, r22964, r22965, r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980, r22981, r22982, r22983, r22984, r22985;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22960, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22961, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r22962);
        mpfr_init(r22963);
        mpfr_init(r22964);
        mpfr_init(r22965);
        mpfr_init(r22966);
        mpfr_init_set_str(r22967, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r22968, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r22969, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r22970, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init_set_str(r22971, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r22972);
        mpfr_init(r22973);
        mpfr_init(r22974);
        mpfr_init(r22975);
        mpfr_init(r22976);
        mpfr_init(r22977);
        mpfr_init(r22978);
        mpfr_init(r22979);
        mpfr_init(r22980);
        mpfr_init(r22981);
        mpfr_init(r22982);
        mpfr_init(r22983);
        mpfr_init(r22984);
        mpfr_init(r22985);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r22962, x, MPFR_RNDN);
        mpfr_abs(r22963, r22962, MPFR_RNDN);
        mpfr_mul(r22964, r22961, r22963, MPFR_RNDN);
        mpfr_add(r22965, r22960, r22964, MPFR_RNDN);
        mpfr_div(r22966, r22960, r22965, MPFR_RNDN);
        ;
        ;
        ;
        ;
        ;
        mpfr_mul(r22972, r22966, r22971, MPFR_RNDN);
        mpfr_add(r22973, r22970, r22972, MPFR_RNDN);
        mpfr_mul(r22974, r22966, r22973, MPFR_RNDN);
        mpfr_add(r22975, r22969, r22974, MPFR_RNDN);
        mpfr_mul(r22976, r22966, r22975, MPFR_RNDN);
        mpfr_add(r22977, r22968, r22976, MPFR_RNDN);
        mpfr_mul(r22978, r22966, r22977, MPFR_RNDN);
        mpfr_add(r22979, r22967, r22978, MPFR_RNDN);
        mpfr_mul(r22980, r22966, r22979, MPFR_RNDN);
        mpfr_mul(r22981, r22963, r22963, MPFR_RNDN);
        mpfr_neg(r22982, r22981, MPFR_RNDN);
        mpfr_exp(r22983, r22982, MPFR_RNDN);
        mpfr_mul(r22984, r22980, r22983, MPFR_RNDN);
        mpfr_sub(r22985, r22960, r22984, MPFR_RNDN);
        return mpfr_get_d(r22985, MPFR_RNDN);
}

static mpfr_t r22986, r22987, r22988, r22989, r22990, r22991, r22992, r22993, r22994, r22995, r22996, r22997, r22998, r22999, r23000, r23001, r23002, r23003, r23004, r23005, r23006, r23007, r23008, r23009, r23010, r23011, r23012, r23013, r23014;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r22986, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r22987, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r22988);
        mpfr_init(r22989);
        mpfr_init(r22990);
        mpfr_init(r22991);
        mpfr_init(r22992);
        mpfr_init_set_str(r22993, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r22994, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r22995, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r22996, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r22997);
        mpfr_init(r22998);
        mpfr_init(r22999);
        mpfr_init_set_str(r23000, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r23001);
        mpfr_init(r23002);
        mpfr_init(r23003);
        mpfr_init(r23004);
        mpfr_init(r23005);
        mpfr_init(r23006);
        mpfr_init(r23007);
        mpfr_init(r23008);
        mpfr_init(r23009);
        mpfr_init(r23010);
        mpfr_init(r23011);
        mpfr_init(r23012);
        mpfr_init(r23013);
        mpfr_init(r23014);
}

double f_fm(double x) {
        ;
        ;
        mpfr_set_d(r22988, x, MPFR_RNDN);
        mpfr_abs(r22989, r22988, MPFR_RNDN);
        mpfr_mul(r22990, r22987, r22989, MPFR_RNDN);
        mpfr_add(r22991, r22986, r22990, MPFR_RNDN);
        mpfr_div(r22992, r22986, r22991, MPFR_RNDN);
        ;
        ;
        ;
        ;
        mpfr_mul(r22997, r22989, r22987, MPFR_RNDN);
        mpfr_add(r22998, r22997, r22986, MPFR_RNDN);
        mpfr_div(r22999, r22996, r22998, MPFR_RNDN);
        ;
        mpfr_add(r23001, r22999, r23000, MPFR_RNDN);
        mpfr_div(r23002, r23001, r22998, MPFR_RNDN);
        mpfr_pow(r23003, r23002, r22986, MPFR_RNDN);
        mpfr_add(r23004, r22995, r23003, MPFR_RNDN);
        mpfr_mul(r23005, r22992, r23004, MPFR_RNDN);
        mpfr_add(r23006, r22994, r23005, MPFR_RNDN);
        mpfr_mul(r23007, r22992, r23006, MPFR_RNDN);
        mpfr_add(r23008, r22993, r23007, MPFR_RNDN);
        mpfr_mul(r23009, r22992, r23008, MPFR_RNDN);
        mpfr_mul(r23010, r22989, r22989, MPFR_RNDN);
        mpfr_neg(r23011, r23010, MPFR_RNDN);
        mpfr_exp(r23012, r23011, MPFR_RNDN);
        mpfr_mul(r23013, r23009, r23012, MPFR_RNDN);
        mpfr_sub(r23014, r22986, r23013, MPFR_RNDN);
        return mpfr_get_d(r23014, MPFR_RNDN);
}

static mpfr_t r23015, r23016, r23017, r23018, r23019, r23020, r23021, r23022, r23023, r23024, r23025, r23026, r23027, r23028, r23029, r23030, r23031, r23032, r23033, r23034, r23035, r23036, r23037, r23038, r23039, r23040, r23041, r23042, r23043;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23015, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r23016, "0.3275911", 10, MPFR_RNDN);
        mpfr_init(r23017);
        mpfr_init(r23018);
        mpfr_init(r23019);
        mpfr_init(r23020);
        mpfr_init(r23021);
        mpfr_init_set_str(r23022, "0.254829592", 10, MPFR_RNDN);
        mpfr_init_set_str(r23023, "-0.284496736", 10, MPFR_RNDN);
        mpfr_init_set_str(r23024, "1.421413741", 10, MPFR_RNDN);
        mpfr_init_set_str(r23025, "1.061405429", 10, MPFR_RNDN);
        mpfr_init(r23026);
        mpfr_init(r23027);
        mpfr_init(r23028);
        mpfr_init_set_str(r23029, "-1.453152027", 10, MPFR_RNDN);
        mpfr_init(r23030);
        mpfr_init(r23031);
        mpfr_init(r23032);
        mpfr_init(r23033);
        mpfr_init(r23034);
        mpfr_init(r23035);
        mpfr_init(r23036);
        mpfr_init(r23037);
        mpfr_init(r23038);
        mpfr_init(r23039);
        mpfr_init(r23040);
        mpfr_init(r23041);
        mpfr_init(r23042);
        mpfr_init(r23043);
}

double f_dm(double x) {
        ;
        ;
        mpfr_set_d(r23017, x, MPFR_RNDN);
        mpfr_abs(r23018, r23017, MPFR_RNDN);
        mpfr_mul(r23019, r23016, r23018, MPFR_RNDN);
        mpfr_add(r23020, r23015, r23019, MPFR_RNDN);
        mpfr_div(r23021, r23015, r23020, MPFR_RNDN);
        ;
        ;
        ;
        ;
        mpfr_mul(r23026, r23018, r23016, MPFR_RNDN);
        mpfr_add(r23027, r23026, r23015, MPFR_RNDN);
        mpfr_div(r23028, r23025, r23027, MPFR_RNDN);
        ;
        mpfr_add(r23030, r23028, r23029, MPFR_RNDN);
        mpfr_div(r23031, r23030, r23027, MPFR_RNDN);
        mpfr_pow(r23032, r23031, r23015, MPFR_RNDN);
        mpfr_add(r23033, r23024, r23032, MPFR_RNDN);
        mpfr_mul(r23034, r23021, r23033, MPFR_RNDN);
        mpfr_add(r23035, r23023, r23034, MPFR_RNDN);
        mpfr_mul(r23036, r23021, r23035, MPFR_RNDN);
        mpfr_add(r23037, r23022, r23036, MPFR_RNDN);
        mpfr_mul(r23038, r23021, r23037, MPFR_RNDN);
        mpfr_mul(r23039, r23018, r23018, MPFR_RNDN);
        mpfr_neg(r23040, r23039, MPFR_RNDN);
        mpfr_exp(r23041, r23040, MPFR_RNDN);
        mpfr_mul(r23042, r23038, r23041, MPFR_RNDN);
        mpfr_sub(r23043, r23015, r23042, MPFR_RNDN);
        return mpfr_get_d(r23043, MPFR_RNDN);
}