#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 r18096 = 1.0f; float r18097 = atan2(1.0, 0.0); float r18098 = sqrt(r18097); float r18099 = r18096 / r18098; float r18100 = 2.0f; float r18101 = x; float r18102 = fabs(r18101); float r18103 = r18100 * r18102; float r18104 = 3.0f; float r18105 = r18100 / r18104; float r18106 = r18102 * r18102; float r18107 = r18106 * r18102; float r18108 = r18105 * r18107; float r18109 = r18103 + r18108; float r18110 = 5.0f; float r18111 = r18096 / r18110; float r18112 = r18107 * r18102; float r18113 = r18112 * r18102; float r18114 = r18111 * r18113; float r18115 = r18109 + r18114; float r18116 = 21.0f; float r18117 = r18096 / r18116; float r18118 = r18113 * r18102; float r18119 = r18118 * r18102; float r18120 = r18117 * r18119; float r18121 = r18115 + r18120; float r18122 = r18099 * r18121; float r18123 = fabs(r18122); return r18123; } double f_id(double x) { double r18124 = 1.0; double r18125 = atan2(1.0, 0.0); double r18126 = sqrt(r18125); double r18127 = r18124 / r18126; double r18128 = 2.0; double r18129 = x; double r18130 = fabs(r18129); double r18131 = r18128 * r18130; double r18132 = 3.0; double r18133 = r18128 / r18132; double r18134 = r18130 * r18130; double r18135 = r18134 * r18130; double r18136 = r18133 * r18135; double r18137 = r18131 + r18136; double r18138 = 5.0; double r18139 = r18124 / r18138; double r18140 = r18135 * r18130; double r18141 = r18140 * r18130; double r18142 = r18139 * r18141; double r18143 = r18137 + r18142; double r18144 = 21.0; double r18145 = r18124 / r18144; double r18146 = r18141 * r18130; double r18147 = r18146 * r18130; double r18148 = r18145 * r18147; double r18149 = r18143 + r18148; double r18150 = r18127 * r18149; double r18151 = fabs(r18150); return r18151; } double f_of(float x) { float r18152 = 0.047619047619047616f; float r18153 = x; float r18154 = fabs(r18153); float r18155 = r18154 * (r18154 * r18154); float r18156 = r18152 * r18155; float r18157 = r18156 * r18155; float r18158 = 0.2f; float r18159 = r18154 * r18158; float r18160 = r18155 * r18159; float r18161 = 0.6666666666666666f; float r18162 = 2.0f; float r18163 = r18162 * r18154; float r18164 = fma(r18161, r18155, r18163); float r18165 = fma(r18160, r18154, r18164); float r18166 = fma(r18157, r18154, r18165); float r18167 = 1.0f; float r18168 = atan2(1.0, 0.0); float r18169 = r18167 / r18168; float r18170 = sqrt(r18169); float r18171 = r18166 * r18170; float r18172 = fabs(r18171); return r18172; } double f_od(double x) { double r18173 = 0.047619047619047616; double r18174 = x; double r18175 = fabs(r18174); double r18176 = r18175 * (r18175 * r18175); double r18177 = r18173 * r18176; double r18178 = r18177 * r18176; double r18179 = 0.2; double r18180 = r18175 * r18179; double r18181 = r18176 * r18180; double r18182 = 0.6666666666666666; double r18183 = 2.0; double r18184 = r18183 * r18175; double r18185 = fma(r18182, r18176, r18184); double r18186 = fma(r18181, r18175, r18185); double r18187 = fma(r18178, r18175, r18186); double r18188 = 1.0; double r18189 = atan2(1.0, 0.0); double r18190 = r18188 / r18189; double r18191 = sqrt(r18190); double r18192 = r18187 * r18191; double r18193 = fabs(r18192); return r18193; } 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 r18194, r18195, r18196, r18197, r18198, r18199, r18200, r18201, r18202, r18203, r18204, r18205, r18206, r18207, r18208, r18209, r18210, r18211, r18212, r18213, r18214, r18215, r18216, r18217, r18218, r18219, r18220, r18221; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18194, "1", 10, MPFR_RNDN); mpfr_init(r18195); mpfr_init(r18196); mpfr_init(r18197); mpfr_init_set_str(r18198, "2", 10, MPFR_RNDN); mpfr_init(r18199); mpfr_init(r18200); mpfr_init(r18201); mpfr_init_set_str(r18202, "3", 10, MPFR_RNDN); mpfr_init(r18203); mpfr_init(r18204); mpfr_init(r18205); mpfr_init(r18206); mpfr_init(r18207); mpfr_init_set_str(r18208, "5", 10, MPFR_RNDN); mpfr_init(r18209); mpfr_init(r18210); mpfr_init(r18211); mpfr_init(r18212); mpfr_init(r18213); mpfr_init_set_str(r18214, "21", 10, MPFR_RNDN); mpfr_init(r18215); mpfr_init(r18216); mpfr_init(r18217); mpfr_init(r18218); mpfr_init(r18219); mpfr_init(r18220); mpfr_init(r18221); } double f_im(double x) { ; mpfr_const_pi(r18195, MPFR_RNDN); mpfr_sqrt(r18196, r18195, MPFR_RNDN); mpfr_div(r18197, r18194, r18196, MPFR_RNDN); ; mpfr_set_d(r18199, x, MPFR_RNDN); mpfr_abs(r18200, r18199, MPFR_RNDN); mpfr_mul(r18201, r18198, r18200, MPFR_RNDN); ; mpfr_div(r18203, r18198, r18202, MPFR_RNDN); mpfr_mul(r18204, r18200, r18200, MPFR_RNDN); mpfr_mul(r18205, r18204, r18200, MPFR_RNDN); mpfr_mul(r18206, r18203, r18205, MPFR_RNDN); mpfr_add(r18207, r18201, r18206, MPFR_RNDN); ; mpfr_div(r18209, r18194, r18208, MPFR_RNDN); mpfr_mul(r18210, r18205, r18200, MPFR_RNDN); mpfr_mul(r18211, r18210, r18200, MPFR_RNDN); mpfr_mul(r18212, r18209, r18211, MPFR_RNDN); mpfr_add(r18213, r18207, r18212, MPFR_RNDN); ; mpfr_div(r18215, r18194, r18214, MPFR_RNDN); mpfr_mul(r18216, r18211, r18200, MPFR_RNDN); mpfr_mul(r18217, r18216, r18200, MPFR_RNDN); mpfr_mul(r18218, r18215, r18217, MPFR_RNDN); mpfr_add(r18219, r18213, r18218, MPFR_RNDN); mpfr_mul(r18220, r18197, r18219, MPFR_RNDN); mpfr_abs(r18221, r18220, MPFR_RNDN); return mpfr_get_d(r18221, MPFR_RNDN); } static mpfr_t r18222, r18223, r18224, r18225, r18226, r18227, r18228, r18229, r18230, r18231, r18232, r18233, r18234, r18235, r18236, r18237, r18238, r18239, r18240, r18241, r18242; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18222, "1/21", 10, MPFR_RNDN); mpfr_init(r18223); mpfr_init(r18224); mpfr_init(r18225); mpfr_init(r18226); mpfr_init(r18227); mpfr_init_set_str(r18228, "1/5", 10, MPFR_RNDN); mpfr_init(r18229); mpfr_init(r18230); mpfr_init_set_str(r18231, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18232, "2", 10, MPFR_RNDN); mpfr_init(r18233); mpfr_init(r18234); mpfr_init(r18235); mpfr_init(r18236); mpfr_init_set_str(r18237, "1", 10, MPFR_RNDN); mpfr_init(r18238); mpfr_init(r18239); mpfr_init(r18240); mpfr_init(r18241); mpfr_init(r18242); } double f_fm(double x) { ; mpfr_set_d(r18223, x, MPFR_RNDN); mpfr_abs(r18224, r18223, MPFR_RNDN); mpfr_mul(r18225, r18224, r18224, MPFR_RNDN); mpfr_mul(r18225, r18225, r18224, MPFR_RNDN); mpfr_mul(r18226, r18222, r18225, MPFR_RNDN); mpfr_mul(r18227, r18226, r18225, MPFR_RNDN); ; mpfr_mul(r18229, r18224, r18228, MPFR_RNDN); mpfr_mul(r18230, r18225, r18229, MPFR_RNDN); ; ; mpfr_mul(r18233, r18232, r18224, MPFR_RNDN); mpfr_fma(r18234, r18231, r18225, r18233, MPFR_RNDN); mpfr_fma(r18235, r18230, r18224, r18234, MPFR_RNDN); mpfr_fma(r18236, r18227, r18224, r18235, MPFR_RNDN); ; mpfr_const_pi(r18238, MPFR_RNDN); mpfr_div(r18239, r18237, r18238, MPFR_RNDN); mpfr_sqrt(r18240, r18239, MPFR_RNDN); mpfr_mul(r18241, r18236, r18240, MPFR_RNDN); mpfr_abs(r18242, r18241, MPFR_RNDN); return mpfr_get_d(r18242, MPFR_RNDN); } static mpfr_t r18243, r18244, r18245, r18246, r18247, r18248, r18249, r18250, r18251, r18252, r18253, r18254, r18255, r18256, r18257, r18258, r18259, r18260, r18261, r18262, r18263; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18243, "1/21", 10, MPFR_RNDN); mpfr_init(r18244); mpfr_init(r18245); mpfr_init(r18246); mpfr_init(r18247); mpfr_init(r18248); mpfr_init_set_str(r18249, "1/5", 10, MPFR_RNDN); mpfr_init(r18250); mpfr_init(r18251); mpfr_init_set_str(r18252, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18253, "2", 10, MPFR_RNDN); mpfr_init(r18254); mpfr_init(r18255); mpfr_init(r18256); mpfr_init(r18257); mpfr_init_set_str(r18258, "1", 10, MPFR_RNDN); mpfr_init(r18259); mpfr_init(r18260); mpfr_init(r18261); mpfr_init(r18262); mpfr_init(r18263); } double f_dm(double x) { ; mpfr_set_d(r18244, x, MPFR_RNDN); mpfr_abs(r18245, r18244, MPFR_RNDN); mpfr_mul(r18246, r18245, r18245, MPFR_RNDN); mpfr_mul(r18246, r18246, r18245, MPFR_RNDN); mpfr_mul(r18247, r18243, r18246, MPFR_RNDN); mpfr_mul(r18248, r18247, r18246, MPFR_RNDN); ; mpfr_mul(r18250, r18245, r18249, MPFR_RNDN); mpfr_mul(r18251, r18246, r18250, MPFR_RNDN); ; ; mpfr_mul(r18254, r18253, r18245, MPFR_RNDN); mpfr_fma(r18255, r18252, r18246, r18254, MPFR_RNDN); mpfr_fma(r18256, r18251, r18245, r18255, MPFR_RNDN); mpfr_fma(r18257, r18248, r18245, r18256, MPFR_RNDN); ; mpfr_const_pi(r18259, MPFR_RNDN); mpfr_div(r18260, r18258, r18259, MPFR_RNDN); mpfr_sqrt(r18261, r18260, MPFR_RNDN); mpfr_mul(r18262, r18257, r18261, MPFR_RNDN); mpfr_abs(r18263, r18262, MPFR_RNDN); return mpfr_get_d(r18263, MPFR_RNDN); }