#include #include #include #include #include char *name = "Jmat.Real.erf"; double f_if(float x) { float r18941 = 1.0f; float r18942 = 0.32759109139442444f; float r18943 = x; float r18944 = fabs(r18943); float r18945 = r18942 * r18944; float r18946 = r18941 + r18945; float r18947 = r18941 / r18946; float r18948 = 0.2548295855522156f; float r18949 = -0.2844967246055603f; float r18950 = 1.421413779258728f; float r18951 = -1.453152060508728f; float r18952 = 1.0614054203033447f; float r18953 = r18947 * r18952; float r18954 = r18951 + r18953; float r18955 = r18947 * r18954; float r18956 = r18950 + r18955; float r18957 = r18947 * r18956; float r18958 = r18949 + r18957; float r18959 = r18947 * r18958; float r18960 = r18948 + r18959; float r18961 = r18947 * r18960; float r18962 = r18944 * r18944; float r18963 = -r18962; float r18964 = exp(r18963); float r18965 = r18961 * r18964; float r18966 = r18941 - r18965; return r18966; } double f_id(double x) { double r18967 = 1.0; double r18968 = 0.32759109139442444; double r18969 = x; double r18970 = fabs(r18969); double r18971 = r18968 * r18970; double r18972 = r18967 + r18971; double r18973 = r18967 / r18972; double r18974 = 0.2548295855522156; double r18975 = -0.2844967246055603; double r18976 = 1.421413779258728; double r18977 = -1.453152060508728; double r18978 = 1.0614054203033447; double r18979 = r18973 * r18978; double r18980 = r18977 + r18979; double r18981 = r18973 * r18980; double r18982 = r18976 + r18981; double r18983 = r18973 * r18982; double r18984 = r18975 + r18983; double r18985 = r18973 * r18984; double r18986 = r18974 + r18985; double r18987 = r18973 * r18986; double r18988 = r18970 * r18970; double r18989 = -r18988; double r18990 = exp(r18989); double r18991 = r18987 * r18990; double r18992 = r18967 - r18991; return r18992; } double f_of(float x) { float r18993 = 1.0f; float r18994 = 0.32759109139442444f; float r18995 = x; float r18996 = fabs(r18995); float r18997 = r18994 * r18996; float r18998 = r18993 + r18997; float r18999 = r18993 / r18998; float r19000 = 0.2548295855522156f; float r19001 = -0.2844967246055603f; float r19002 = 1.421413779258728f; float r19003 = -1.453152060508728f; float r19004 = r19003 * r19003; float r19005 = 1.0614054203033447f; float r19006 = r18999 * r19005; float r19007 = r19006 * r19006; float r19008 = r19004 - r19007; float r19009 = pow(r19008, r18993); float r19010 = r18996 * r18994; float r19011 = r18993 + r19010; float r19012 = r19005 / r19011; float r19013 = r19003 - r19012; float r19014 = r19009 / r19013; float r19015 = r18999 * r19014; float r19016 = r19002 + r19015; float r19017 = r18999 * r19016; float r19018 = r19001 + r19017; float r19019 = r18999 * r19018; float r19020 = r19000 + r19019; float r19021 = r18999 * r19020; float r19022 = r18996 * r18996; float r19023 = -r19022; float r19024 = exp(r19023); float r19025 = r19021 * r19024; float r19026 = r18993 - r19025; return r19026; } double f_od(double x) { double r19027 = 1.0; double r19028 = 0.32759109139442444; double r19029 = x; double r19030 = fabs(r19029); double r19031 = r19028 * r19030; double r19032 = r19027 + r19031; double r19033 = r19027 / r19032; double r19034 = 0.2548295855522156; double r19035 = -0.2844967246055603; double r19036 = 1.421413779258728; double r19037 = -1.453152060508728; double r19038 = r19037 * r19037; double r19039 = 1.0614054203033447; double r19040 = r19033 * r19039; double r19041 = r19040 * r19040; double r19042 = r19038 - r19041; double r19043 = pow(r19042, r19027); double r19044 = r19030 * r19028; double r19045 = r19027 + r19044; double r19046 = r19039 / r19045; double r19047 = r19037 - r19046; double r19048 = r19043 / r19047; double r19049 = r19033 * r19048; double r19050 = r19036 + r19049; double r19051 = r19033 * r19050; double r19052 = r19035 + r19051; double r19053 = r19033 * r19052; double r19054 = r19034 + r19053; double r19055 = r19033 * r19054; double r19056 = r19030 * r19030; double r19057 = -r19056; double r19058 = exp(r19057); double r19059 = r19055 * r19058; double r19060 = r19027 - r19059; return r19060; } 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 r19061, r19062, r19063, r19064, r19065, r19066, r19067, r19068, r19069, r19070, r19071, r19072, r19073, r19074, r19075, r19076, r19077, r19078, r19079, r19080, r19081, r19082, r19083, r19084, r19085, r19086; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19061, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19062, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19063); mpfr_init(r19064); mpfr_init(r19065); mpfr_init(r19066); mpfr_init(r19067); mpfr_init_set_str(r19068, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r19069, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19070, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r19071, "-1.453152027", 10, MPFR_RNDN); mpfr_init_set_str(r19072, "1.061405429", 10, MPFR_RNDN); mpfr_init(r19073); mpfr_init(r19074); mpfr_init(r19075); mpfr_init(r19076); mpfr_init(r19077); mpfr_init(r19078); mpfr_init(r19079); mpfr_init(r19080); mpfr_init(r19081); mpfr_init(r19082); mpfr_init(r19083); mpfr_init(r19084); mpfr_init(r19085); mpfr_init(r19086); } double f_im(double x) { ; ; mpfr_set_d(r19063, x, MPFR_RNDN); mpfr_abs(r19064, r19063, MPFR_RNDN); mpfr_mul(r19065, r19062, r19064, MPFR_RNDN); mpfr_add(r19066, r19061, r19065, MPFR_RNDN); mpfr_div(r19067, r19061, r19066, MPFR_RNDN); ; ; ; ; ; mpfr_mul(r19073, r19067, r19072, MPFR_RNDN); mpfr_add(r19074, r19071, r19073, MPFR_RNDN); mpfr_mul(r19075, r19067, r19074, MPFR_RNDN); mpfr_add(r19076, r19070, r19075, MPFR_RNDN); mpfr_mul(r19077, r19067, r19076, MPFR_RNDN); mpfr_add(r19078, r19069, r19077, MPFR_RNDN); mpfr_mul(r19079, r19067, r19078, MPFR_RNDN); mpfr_add(r19080, r19068, r19079, MPFR_RNDN); mpfr_mul(r19081, r19067, r19080, MPFR_RNDN); mpfr_mul(r19082, r19064, r19064, MPFR_RNDN); mpfr_neg(r19083, r19082, MPFR_RNDN); mpfr_exp(r19084, r19083, MPFR_RNDN); mpfr_mul(r19085, r19081, r19084, MPFR_RNDN); mpfr_sub(r19086, r19061, r19085, MPFR_RNDN); return mpfr_get_d(r19086, MPFR_RNDN); } static mpfr_t r19087, r19088, r19089, r19090, r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098, r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116, r19117, r19118, r19119, r19120; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19087, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19088, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19089); mpfr_init(r19090); mpfr_init(r19091); mpfr_init(r19092); mpfr_init(r19093); mpfr_init_set_str(r19094, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r19095, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19096, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r19097, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19098); mpfr_init_set_str(r19099, "1.061405429", 10, MPFR_RNDN); mpfr_init(r19100); mpfr_init(r19101); mpfr_init(r19102); mpfr_init(r19103); mpfr_init(r19104); mpfr_init(r19105); mpfr_init(r19106); mpfr_init(r19107); mpfr_init(r19108); mpfr_init(r19109); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init(r19113); mpfr_init(r19114); mpfr_init(r19115); mpfr_init(r19116); mpfr_init(r19117); mpfr_init(r19118); mpfr_init(r19119); mpfr_init(r19120); } double f_fm(double x) { ; ; mpfr_set_d(r19089, x, MPFR_RNDN); mpfr_abs(r19090, r19089, MPFR_RNDN); mpfr_mul(r19091, r19088, r19090, MPFR_RNDN); mpfr_add(r19092, r19087, r19091, MPFR_RNDN); mpfr_div(r19093, r19087, r19092, MPFR_RNDN); ; ; ; ; mpfr_sqr(r19098, r19097, MPFR_RNDN); ; mpfr_mul(r19100, r19093, r19099, MPFR_RNDN); mpfr_sqr(r19101, r19100, MPFR_RNDN); mpfr_sub(r19102, r19098, r19101, MPFR_RNDN); mpfr_pow(r19103, r19102, r19087, MPFR_RNDN); mpfr_mul(r19104, r19090, r19088, MPFR_RNDN); mpfr_add(r19105, r19087, r19104, MPFR_RNDN); mpfr_div(r19106, r19099, r19105, MPFR_RNDN); mpfr_sub(r19107, r19097, r19106, MPFR_RNDN); mpfr_div(r19108, r19103, r19107, MPFR_RNDN); mpfr_mul(r19109, r19093, r19108, MPFR_RNDN); mpfr_add(r19110, r19096, r19109, MPFR_RNDN); mpfr_mul(r19111, r19093, r19110, MPFR_RNDN); mpfr_add(r19112, r19095, r19111, MPFR_RNDN); mpfr_mul(r19113, r19093, r19112, MPFR_RNDN); mpfr_add(r19114, r19094, r19113, MPFR_RNDN); mpfr_mul(r19115, r19093, r19114, MPFR_RNDN); mpfr_mul(r19116, r19090, r19090, MPFR_RNDN); mpfr_neg(r19117, r19116, MPFR_RNDN); mpfr_exp(r19118, r19117, MPFR_RNDN); mpfr_mul(r19119, r19115, r19118, MPFR_RNDN); mpfr_sub(r19120, r19087, r19119, MPFR_RNDN); return mpfr_get_d(r19120, MPFR_RNDN); } static mpfr_t r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134, r19135, r19136, r19137, r19138, r19139, r19140, r19141, r19142, r19143, r19144, r19145, r19146, r19147, r19148, r19149, r19150, r19151, r19152, r19153, r19154; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19121, "1", 10, MPFR_RNDN); mpfr_init_set_str(r19122, "0.3275911", 10, MPFR_RNDN); mpfr_init(r19123); mpfr_init(r19124); mpfr_init(r19125); mpfr_init(r19126); mpfr_init(r19127); mpfr_init_set_str(r19128, "0.254829592", 10, MPFR_RNDN); mpfr_init_set_str(r19129, "-0.284496736", 10, MPFR_RNDN); mpfr_init_set_str(r19130, "1.421413741", 10, MPFR_RNDN); mpfr_init_set_str(r19131, "-1.453152027", 10, MPFR_RNDN); mpfr_init(r19132); mpfr_init_set_str(r19133, "1.061405429", 10, MPFR_RNDN); mpfr_init(r19134); mpfr_init(r19135); mpfr_init(r19136); mpfr_init(r19137); mpfr_init(r19138); mpfr_init(r19139); mpfr_init(r19140); mpfr_init(r19141); mpfr_init(r19142); mpfr_init(r19143); mpfr_init(r19144); mpfr_init(r19145); mpfr_init(r19146); mpfr_init(r19147); mpfr_init(r19148); mpfr_init(r19149); mpfr_init(r19150); mpfr_init(r19151); mpfr_init(r19152); mpfr_init(r19153); mpfr_init(r19154); } double f_dm(double x) { ; ; mpfr_set_d(r19123, x, MPFR_RNDN); mpfr_abs(r19124, r19123, MPFR_RNDN); mpfr_mul(r19125, r19122, r19124, MPFR_RNDN); mpfr_add(r19126, r19121, r19125, MPFR_RNDN); mpfr_div(r19127, r19121, r19126, MPFR_RNDN); ; ; ; ; mpfr_sqr(r19132, r19131, MPFR_RNDN); ; mpfr_mul(r19134, r19127, r19133, MPFR_RNDN); mpfr_sqr(r19135, r19134, MPFR_RNDN); mpfr_sub(r19136, r19132, r19135, MPFR_RNDN); mpfr_pow(r19137, r19136, r19121, MPFR_RNDN); mpfr_mul(r19138, r19124, r19122, MPFR_RNDN); mpfr_add(r19139, r19121, r19138, MPFR_RNDN); mpfr_div(r19140, r19133, r19139, MPFR_RNDN); mpfr_sub(r19141, r19131, r19140, MPFR_RNDN); mpfr_div(r19142, r19137, r19141, MPFR_RNDN); mpfr_mul(r19143, r19127, r19142, MPFR_RNDN); mpfr_add(r19144, r19130, r19143, MPFR_RNDN); mpfr_mul(r19145, r19127, r19144, MPFR_RNDN); mpfr_add(r19146, r19129, r19145, MPFR_RNDN); mpfr_mul(r19147, r19127, r19146, MPFR_RNDN); mpfr_add(r19148, r19128, r19147, MPFR_RNDN); mpfr_mul(r19149, r19127, r19148, MPFR_RNDN); mpfr_mul(r19150, r19124, r19124, MPFR_RNDN); mpfr_neg(r19151, r19150, MPFR_RNDN); mpfr_exp(r19152, r19151, MPFR_RNDN); mpfr_mul(r19153, r19149, r19152, MPFR_RNDN); mpfr_sub(r19154, r19121, r19153, MPFR_RNDN); return mpfr_get_d(r19154, MPFR_RNDN); }