#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 r23965 = 1; float r23966 = atan2(1.0, 0.0); float r23967 = sqrt(r23966); float r23968 = r23965 / r23967; float r23969 = 2; float r23970 = x; float r23971 = fabs(r23970); float r23972 = r23969 * r23971; float r23973 = 3; float r23974 = r23969 / r23973; float r23975 = r23971 * r23971; float r23976 = r23975 * r23971; float r23977 = r23974 * r23976; float r23978 = r23972 + r23977; float r23979 = 5; float r23980 = r23965 / r23979; float r23981 = r23976 * r23971; float r23982 = r23981 * r23971; float r23983 = r23980 * r23982; float r23984 = r23978 + r23983; float r23985 = 21; float r23986 = r23965 / r23985; float r23987 = r23982 * r23971; float r23988 = r23987 * r23971; float r23989 = r23986 * r23988; float r23990 = r23984 + r23989; float r23991 = r23968 * r23990; float r23992 = fabs(r23991); return r23992; } double f_id(double x) { double r23993 = 1; double r23994 = atan2(1.0, 0.0); double r23995 = sqrt(r23994); double r23996 = r23993 / r23995; double r23997 = 2; double r23998 = x; double r23999 = fabs(r23998); double r24000 = r23997 * r23999; double r24001 = 3; double r24002 = r23997 / r24001; double r24003 = r23999 * r23999; double r24004 = r24003 * r23999; double r24005 = r24002 * r24004; double r24006 = r24000 + r24005; double r24007 = 5; double r24008 = r23993 / r24007; double r24009 = r24004 * r23999; double r24010 = r24009 * r23999; double r24011 = r24008 * r24010; double r24012 = r24006 + r24011; double r24013 = 21; double r24014 = r23993 / r24013; double r24015 = r24010 * r23999; double r24016 = r24015 * r23999; double r24017 = r24014 * r24016; double r24018 = r24012 + r24017; double r24019 = r23996 * r24018; double r24020 = fabs(r24019); return r24020; } double f_of(float x) { float r24021 = 1; float r24022 = atan2(1.0, 0.0); float r24023 = sqrt(r24022); float r24024 = r24021 / r24023; float r24025 = 2; float r24026 = x; float r24027 = fabs(r24026); float r24028 = r24025 * r24027; float r24029 = 3; float r24030 = r24025 / r24029; float r24031 = r24027 * r24027; float r24032 = r24031 * r24027; float r24033 = r24030 * r24032; float r24034 = r24028 + r24033; float r24035 = 5; float r24036 = r24021 / r24035; float r24037 = r24032 * r24027; float r24038 = r24037 * r24027; float r24039 = r24036 * r24038; float r24040 = r24034 + r24039; float r24041 = 21; float r24042 = r24021 / r24041; float r24043 = r24038 * r24027; float r24044 = r24043 * r24027; float r24045 = r24042 * r24044; float r24046 = r24040 + r24045; float r24047 = r24024 * r24046; float r24048 = fabs(r24047); return r24048; } double f_od(double x) { double r24049 = 1; double r24050 = atan2(1.0, 0.0); double r24051 = sqrt(r24050); double r24052 = r24049 / r24051; double r24053 = 2; double r24054 = x; double r24055 = fabs(r24054); double r24056 = r24053 * r24055; double r24057 = 3; double r24058 = r24053 / r24057; double r24059 = r24055 * r24055; double r24060 = r24059 * r24055; double r24061 = r24058 * r24060; double r24062 = r24056 + r24061; double r24063 = 5; double r24064 = r24049 / r24063; double r24065 = r24060 * r24055; double r24066 = r24065 * r24055; double r24067 = r24064 * r24066; double r24068 = r24062 + r24067; double r24069 = 21; double r24070 = r24049 / r24069; double r24071 = r24066 * r24055; double r24072 = r24071 * r24055; double r24073 = r24070 * r24072; double r24074 = r24068 + r24073; double r24075 = r24052 * r24074; double r24076 = fabs(r24075); return r24076; } 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 r24077, r24078, r24079, r24080, r24081, r24082, r24083, r24084, r24085, r24086, r24087, r24088, r24089, r24090, r24091, r24092, r24093, r24094, r24095, r24096, r24097, r24098, r24099, r24100, r24101, r24102, r24103, r24104; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r24077, "1", 10, MPFR_RNDN); mpfr_init(r24078); mpfr_init(r24079); mpfr_init(r24080); mpfr_init_set_str(r24081, "2", 10, MPFR_RNDN); mpfr_init(r24082); mpfr_init(r24083); mpfr_init(r24084); mpfr_init_set_str(r24085, "3", 10, MPFR_RNDN); mpfr_init(r24086); mpfr_init(r24087); mpfr_init(r24088); mpfr_init(r24089); mpfr_init(r24090); mpfr_init_set_str(r24091, "5", 10, MPFR_RNDN); mpfr_init(r24092); mpfr_init(r24093); mpfr_init(r24094); mpfr_init(r24095); mpfr_init(r24096); mpfr_init_set_str(r24097, "21", 10, MPFR_RNDN); mpfr_init(r24098); mpfr_init(r24099); mpfr_init(r24100); mpfr_init(r24101); mpfr_init(r24102); mpfr_init(r24103); mpfr_init(r24104); } double f_im(double x) { ; mpfr_const_pi(r24078, MPFR_RNDN); mpfr_sqrt(r24079, r24078, MPFR_RNDN); mpfr_div(r24080, r24077, r24079, MPFR_RNDN); ; mpfr_set_d(r24082, x, MPFR_RNDN); mpfr_abs(r24083, r24082, MPFR_RNDN); mpfr_mul(r24084, r24081, r24083, MPFR_RNDN); ; mpfr_div(r24086, r24081, r24085, MPFR_RNDN); mpfr_mul(r24087, r24083, r24083, MPFR_RNDN); mpfr_mul(r24088, r24087, r24083, MPFR_RNDN); mpfr_mul(r24089, r24086, r24088, MPFR_RNDN); mpfr_add(r24090, r24084, r24089, MPFR_RNDN); ; mpfr_div(r24092, r24077, r24091, MPFR_RNDN); mpfr_mul(r24093, r24088, r24083, MPFR_RNDN); mpfr_mul(r24094, r24093, r24083, MPFR_RNDN); mpfr_mul(r24095, r24092, r24094, MPFR_RNDN); mpfr_add(r24096, r24090, r24095, MPFR_RNDN); ; mpfr_div(r24098, r24077, r24097, MPFR_RNDN); mpfr_mul(r24099, r24094, r24083, MPFR_RNDN); mpfr_mul(r24100, r24099, r24083, MPFR_RNDN); mpfr_mul(r24101, r24098, r24100, MPFR_RNDN); mpfr_add(r24102, r24096, r24101, MPFR_RNDN); mpfr_mul(r24103, r24080, r24102, MPFR_RNDN); mpfr_abs(r24104, r24103, MPFR_RNDN); return mpfr_get_d(r24104, MPFR_RNDN); } static mpfr_t r24105, r24106, r24107, r24108, r24109, r24110, r24111, r24112, r24113, r24114, r24115, r24116, r24117, r24118, r24119, r24120, r24121, r24122, r24123, r24124, r24125, r24126, r24127, r24128, r24129, r24130, r24131, r24132; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24105, "1", 10, MPFR_RNDN); mpfr_init(r24106); mpfr_init(r24107); mpfr_init(r24108); mpfr_init_set_str(r24109, "2", 10, MPFR_RNDN); mpfr_init(r24110); mpfr_init(r24111); mpfr_init(r24112); mpfr_init_set_str(r24113, "3", 10, MPFR_RNDN); mpfr_init(r24114); mpfr_init(r24115); mpfr_init(r24116); mpfr_init(r24117); mpfr_init(r24118); mpfr_init_set_str(r24119, "5", 10, MPFR_RNDN); mpfr_init(r24120); mpfr_init(r24121); mpfr_init(r24122); mpfr_init(r24123); mpfr_init(r24124); mpfr_init_set_str(r24125, "21", 10, MPFR_RNDN); mpfr_init(r24126); mpfr_init(r24127); mpfr_init(r24128); mpfr_init(r24129); mpfr_init(r24130); mpfr_init(r24131); mpfr_init(r24132); } double f_fm(double x) { ; mpfr_const_pi(r24106, MPFR_RNDN); mpfr_sqrt(r24107, r24106, MPFR_RNDN); mpfr_div(r24108, r24105, r24107, MPFR_RNDN); ; mpfr_set_d(r24110, x, MPFR_RNDN); mpfr_abs(r24111, r24110, MPFR_RNDN); mpfr_mul(r24112, r24109, r24111, MPFR_RNDN); ; mpfr_div(r24114, r24109, r24113, MPFR_RNDN); mpfr_mul(r24115, r24111, r24111, MPFR_RNDN); mpfr_mul(r24116, r24115, r24111, MPFR_RNDN); mpfr_mul(r24117, r24114, r24116, MPFR_RNDN); mpfr_add(r24118, r24112, r24117, MPFR_RNDN); ; mpfr_div(r24120, r24105, r24119, MPFR_RNDN); mpfr_mul(r24121, r24116, r24111, MPFR_RNDN); mpfr_mul(r24122, r24121, r24111, MPFR_RNDN); mpfr_mul(r24123, r24120, r24122, MPFR_RNDN); mpfr_add(r24124, r24118, r24123, MPFR_RNDN); ; mpfr_div(r24126, r24105, r24125, MPFR_RNDN); mpfr_mul(r24127, r24122, r24111, MPFR_RNDN); mpfr_mul(r24128, r24127, r24111, MPFR_RNDN); mpfr_mul(r24129, r24126, r24128, MPFR_RNDN); mpfr_add(r24130, r24124, r24129, MPFR_RNDN); mpfr_mul(r24131, r24108, r24130, MPFR_RNDN); mpfr_abs(r24132, r24131, MPFR_RNDN); return mpfr_get_d(r24132, MPFR_RNDN); } static mpfr_t r24133, r24134, r24135, r24136, r24137, r24138, r24139, r24140, r24141, r24142, r24143, r24144, r24145, r24146, r24147, r24148, r24149, r24150, r24151, r24152, r24153, r24154, r24155, r24156, r24157, r24158, r24159, r24160; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r24133, "1", 10, MPFR_RNDN); mpfr_init(r24134); mpfr_init(r24135); mpfr_init(r24136); mpfr_init_set_str(r24137, "2", 10, MPFR_RNDN); mpfr_init(r24138); mpfr_init(r24139); mpfr_init(r24140); mpfr_init_set_str(r24141, "3", 10, MPFR_RNDN); mpfr_init(r24142); mpfr_init(r24143); mpfr_init(r24144); mpfr_init(r24145); mpfr_init(r24146); mpfr_init_set_str(r24147, "5", 10, MPFR_RNDN); mpfr_init(r24148); mpfr_init(r24149); mpfr_init(r24150); mpfr_init(r24151); mpfr_init(r24152); mpfr_init_set_str(r24153, "21", 10, MPFR_RNDN); mpfr_init(r24154); mpfr_init(r24155); mpfr_init(r24156); mpfr_init(r24157); mpfr_init(r24158); mpfr_init(r24159); mpfr_init(r24160); } double f_dm(double x) { ; mpfr_const_pi(r24134, MPFR_RNDN); mpfr_sqrt(r24135, r24134, MPFR_RNDN); mpfr_div(r24136, r24133, r24135, MPFR_RNDN); ; mpfr_set_d(r24138, x, MPFR_RNDN); mpfr_abs(r24139, r24138, MPFR_RNDN); mpfr_mul(r24140, r24137, r24139, MPFR_RNDN); ; mpfr_div(r24142, r24137, r24141, MPFR_RNDN); mpfr_mul(r24143, r24139, r24139, MPFR_RNDN); mpfr_mul(r24144, r24143, r24139, MPFR_RNDN); mpfr_mul(r24145, r24142, r24144, MPFR_RNDN); mpfr_add(r24146, r24140, r24145, MPFR_RNDN); ; mpfr_div(r24148, r24133, r24147, MPFR_RNDN); mpfr_mul(r24149, r24144, r24139, MPFR_RNDN); mpfr_mul(r24150, r24149, r24139, MPFR_RNDN); mpfr_mul(r24151, r24148, r24150, MPFR_RNDN); mpfr_add(r24152, r24146, r24151, MPFR_RNDN); ; mpfr_div(r24154, r24133, r24153, MPFR_RNDN); mpfr_mul(r24155, r24150, r24139, MPFR_RNDN); mpfr_mul(r24156, r24155, r24139, MPFR_RNDN); mpfr_mul(r24157, r24154, r24156, MPFR_RNDN); mpfr_add(r24158, r24152, r24157, MPFR_RNDN); mpfr_mul(r24159, r24136, r24158, MPFR_RNDN); mpfr_abs(r24160, r24159, MPFR_RNDN); return mpfr_get_d(r24160, MPFR_RNDN); }