#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 r18012 = 1.0f; float r18013 = atan2(1.0, 0.0); float r18014 = sqrt(r18013); float r18015 = r18012 / r18014; float r18016 = 2.0f; float r18017 = x; float r18018 = fabs(r18017); float r18019 = r18016 * r18018; float r18020 = 3.0f; float r18021 = r18016 / r18020; float r18022 = r18018 * r18018; float r18023 = r18022 * r18018; float r18024 = r18021 * r18023; float r18025 = r18019 + r18024; float r18026 = 5.0f; float r18027 = r18012 / r18026; float r18028 = r18023 * r18018; float r18029 = r18028 * r18018; float r18030 = r18027 * r18029; float r18031 = r18025 + r18030; float r18032 = 21.0f; float r18033 = r18012 / r18032; float r18034 = r18029 * r18018; float r18035 = r18034 * r18018; float r18036 = r18033 * r18035; float r18037 = r18031 + r18036; float r18038 = r18015 * r18037; float r18039 = fabs(r18038); return r18039; } double f_id(double x) { double r18040 = 1.0; double r18041 = atan2(1.0, 0.0); double r18042 = sqrt(r18041); double r18043 = r18040 / r18042; double r18044 = 2.0; double r18045 = x; double r18046 = fabs(r18045); double r18047 = r18044 * r18046; double r18048 = 3.0; double r18049 = r18044 / r18048; double r18050 = r18046 * r18046; double r18051 = r18050 * r18046; double r18052 = r18049 * r18051; double r18053 = r18047 + r18052; double r18054 = 5.0; double r18055 = r18040 / r18054; double r18056 = r18051 * r18046; double r18057 = r18056 * r18046; double r18058 = r18055 * r18057; double r18059 = r18053 + r18058; double r18060 = 21.0; double r18061 = r18040 / r18060; double r18062 = r18057 * r18046; double r18063 = r18062 * r18046; double r18064 = r18061 * r18063; double r18065 = r18059 + r18064; double r18066 = r18043 * r18065; double r18067 = fabs(r18066); return r18067; } double f_of(float x) { float r18068 = x; float r18069 = fabs(r18068); float r18070 = 5.0f; float r18071 = r18069 / r18070; float r18072 = r18069 * (r18069 * r18069); float r18073 = r18071 * r18072; float r18074 = 2.0f; float r18075 = 3.0f; float r18076 = r18074 / r18075; float r18077 = r18074 * r18069; float r18078 = fma(r18076, r18072, r18077); float r18079 = fma(r18073, r18069, r18078); float r18080 = r18072 * r18072; float r18081 = 21.0f; float r18082 = r18081 / r18069; float r18083 = r18080 / r18082; float r18084 = r18079 + r18083; float r18085 = atan2(1.0, 0.0); float r18086 = sqrt(r18085); float r18087 = sqrt(r18086); float r18088 = r18087 * r18087; float r18089 = r18084 / r18088; float r18090 = fabs(r18089); return r18090; } double f_od(double x) { double r18091 = x; double r18092 = fabs(r18091); double r18093 = 5.0; double r18094 = r18092 / r18093; double r18095 = r18092 * (r18092 * r18092); double r18096 = r18094 * r18095; double r18097 = 2.0; double r18098 = 3.0; double r18099 = r18097 / r18098; double r18100 = r18097 * r18092; double r18101 = fma(r18099, r18095, r18100); double r18102 = fma(r18096, r18092, r18101); double r18103 = r18095 * r18095; double r18104 = 21.0; double r18105 = r18104 / r18092; double r18106 = r18103 / r18105; double r18107 = r18102 + r18106; double r18108 = atan2(1.0, 0.0); double r18109 = sqrt(r18108); double r18110 = sqrt(r18109); double r18111 = r18110 * r18110; double r18112 = r18107 / r18111; double r18113 = fabs(r18112); return r18113; } 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 r18114, r18115, r18116, r18117, r18118, r18119, r18120, r18121, r18122, r18123, r18124, r18125, r18126, r18127, r18128, r18129, r18130, r18131, r18132, r18133, r18134, r18135, r18136, r18137, r18138, r18139, r18140, r18141; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18114, "1", 10, MPFR_RNDN); mpfr_init(r18115); mpfr_init(r18116); mpfr_init(r18117); mpfr_init_set_str(r18118, "2", 10, MPFR_RNDN); mpfr_init(r18119); mpfr_init(r18120); mpfr_init(r18121); mpfr_init_set_str(r18122, "3", 10, MPFR_RNDN); mpfr_init(r18123); mpfr_init(r18124); mpfr_init(r18125); mpfr_init(r18126); mpfr_init(r18127); mpfr_init_set_str(r18128, "5", 10, MPFR_RNDN); mpfr_init(r18129); mpfr_init(r18130); mpfr_init(r18131); mpfr_init(r18132); mpfr_init(r18133); mpfr_init_set_str(r18134, "21", 10, MPFR_RNDN); mpfr_init(r18135); mpfr_init(r18136); mpfr_init(r18137); mpfr_init(r18138); mpfr_init(r18139); mpfr_init(r18140); mpfr_init(r18141); } double f_im(double x) { ; mpfr_const_pi(r18115, MPFR_RNDN); mpfr_sqrt(r18116, r18115, MPFR_RNDN); mpfr_div(r18117, r18114, r18116, MPFR_RNDN); ; mpfr_set_d(r18119, x, MPFR_RNDN); mpfr_abs(r18120, r18119, MPFR_RNDN); mpfr_mul(r18121, r18118, r18120, MPFR_RNDN); ; mpfr_div(r18123, r18118, r18122, MPFR_RNDN); mpfr_mul(r18124, r18120, r18120, MPFR_RNDN); mpfr_mul(r18125, r18124, r18120, MPFR_RNDN); mpfr_mul(r18126, r18123, r18125, MPFR_RNDN); mpfr_add(r18127, r18121, r18126, MPFR_RNDN); ; mpfr_div(r18129, r18114, r18128, MPFR_RNDN); mpfr_mul(r18130, r18125, r18120, MPFR_RNDN); mpfr_mul(r18131, r18130, r18120, MPFR_RNDN); mpfr_mul(r18132, r18129, r18131, MPFR_RNDN); mpfr_add(r18133, r18127, r18132, MPFR_RNDN); ; mpfr_div(r18135, r18114, r18134, MPFR_RNDN); mpfr_mul(r18136, r18131, r18120, MPFR_RNDN); mpfr_mul(r18137, r18136, r18120, MPFR_RNDN); mpfr_mul(r18138, r18135, r18137, MPFR_RNDN); mpfr_add(r18139, r18133, r18138, MPFR_RNDN); mpfr_mul(r18140, r18117, r18139, MPFR_RNDN); mpfr_abs(r18141, r18140, MPFR_RNDN); return mpfr_get_d(r18141, MPFR_RNDN); } static mpfr_t r18142, r18143, r18144, r18145, r18146, r18147, r18148, r18149, r18150, r18151, r18152, r18153, r18154, r18155, r18156, r18157, r18158, r18159, r18160, r18161, r18162, r18163, r18164; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18142); mpfr_init(r18143); mpfr_init_set_str(r18144, "5", 10, MPFR_RNDN); mpfr_init(r18145); mpfr_init(r18146); mpfr_init(r18147); mpfr_init_set_str(r18148, "2", 10, MPFR_RNDN); mpfr_init_set_str(r18149, "3", 10, MPFR_RNDN); mpfr_init(r18150); mpfr_init(r18151); mpfr_init(r18152); mpfr_init(r18153); mpfr_init(r18154); mpfr_init_set_str(r18155, "21", 10, MPFR_RNDN); mpfr_init(r18156); mpfr_init(r18157); mpfr_init(r18158); mpfr_init(r18159); mpfr_init(r18160); mpfr_init(r18161); mpfr_init(r18162); mpfr_init(r18163); mpfr_init(r18164); } double f_fm(double x) { mpfr_set_d(r18142, x, MPFR_RNDN); mpfr_abs(r18143, r18142, MPFR_RNDN); ; mpfr_div(r18145, r18143, r18144, MPFR_RNDN); mpfr_mul(r18146, r18143, r18143, MPFR_RNDN); mpfr_mul(r18146, r18146, r18143, MPFR_RNDN); mpfr_mul(r18147, r18145, r18146, MPFR_RNDN); ; ; mpfr_div(r18150, r18148, r18149, MPFR_RNDN); mpfr_mul(r18151, r18148, r18143, MPFR_RNDN); mpfr_fma(r18152, r18150, r18146, r18151, MPFR_RNDN); mpfr_fma(r18153, r18147, r18143, r18152, MPFR_RNDN); mpfr_mul(r18154, r18146, r18146, MPFR_RNDN); ; mpfr_div(r18156, r18155, r18143, MPFR_RNDN); mpfr_div(r18157, r18154, r18156, MPFR_RNDN); mpfr_add(r18158, r18153, r18157, MPFR_RNDN); mpfr_const_pi(r18159, MPFR_RNDN); mpfr_sqrt(r18160, r18159, MPFR_RNDN); mpfr_sqrt(r18161, r18160, MPFR_RNDN); mpfr_sqr(r18162, r18161, MPFR_RNDN); mpfr_div(r18163, r18158, r18162, MPFR_RNDN); mpfr_abs(r18164, r18163, MPFR_RNDN); return mpfr_get_d(r18164, MPFR_RNDN); } static mpfr_t r18165, r18166, r18167, r18168, r18169, r18170, r18171, r18172, r18173, r18174, r18175, r18176, r18177, r18178, r18179, r18180, r18181, r18182, r18183, r18184, r18185, r18186, r18187; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18165); mpfr_init(r18166); mpfr_init_set_str(r18167, "5", 10, MPFR_RNDN); mpfr_init(r18168); mpfr_init(r18169); mpfr_init(r18170); mpfr_init_set_str(r18171, "2", 10, MPFR_RNDN); mpfr_init_set_str(r18172, "3", 10, MPFR_RNDN); mpfr_init(r18173); mpfr_init(r18174); mpfr_init(r18175); mpfr_init(r18176); mpfr_init(r18177); mpfr_init_set_str(r18178, "21", 10, MPFR_RNDN); mpfr_init(r18179); mpfr_init(r18180); mpfr_init(r18181); mpfr_init(r18182); mpfr_init(r18183); mpfr_init(r18184); mpfr_init(r18185); mpfr_init(r18186); mpfr_init(r18187); } double f_dm(double x) { mpfr_set_d(r18165, x, MPFR_RNDN); mpfr_abs(r18166, r18165, MPFR_RNDN); ; mpfr_div(r18168, r18166, r18167, MPFR_RNDN); mpfr_mul(r18169, r18166, r18166, MPFR_RNDN); mpfr_mul(r18169, r18169, r18166, MPFR_RNDN); mpfr_mul(r18170, r18168, r18169, MPFR_RNDN); ; ; mpfr_div(r18173, r18171, r18172, MPFR_RNDN); mpfr_mul(r18174, r18171, r18166, MPFR_RNDN); mpfr_fma(r18175, r18173, r18169, r18174, MPFR_RNDN); mpfr_fma(r18176, r18170, r18166, r18175, MPFR_RNDN); mpfr_mul(r18177, r18169, r18169, MPFR_RNDN); ; mpfr_div(r18179, r18178, r18166, MPFR_RNDN); mpfr_div(r18180, r18177, r18179, MPFR_RNDN); mpfr_add(r18181, r18176, r18180, MPFR_RNDN); mpfr_const_pi(r18182, MPFR_RNDN); mpfr_sqrt(r18183, r18182, MPFR_RNDN); mpfr_sqrt(r18184, r18183, MPFR_RNDN); mpfr_sqr(r18185, r18184, MPFR_RNDN); mpfr_div(r18186, r18181, r18185, MPFR_RNDN); mpfr_abs(r18187, r18186, MPFR_RNDN); return mpfr_get_d(r18187, MPFR_RNDN); }