#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 r18172 = 1.0f; float r18173 = atan2(1.0, 0.0); float r18174 = sqrt(r18173); float r18175 = r18172 / r18174; float r18176 = 2.0f; float r18177 = x; float r18178 = fabs(r18177); float r18179 = r18176 * r18178; float r18180 = 3.0f; float r18181 = r18176 / r18180; float r18182 = r18178 * r18178; float r18183 = r18182 * r18178; float r18184 = r18181 * r18183; float r18185 = r18179 + r18184; float r18186 = 5.0f; float r18187 = r18172 / r18186; float r18188 = r18183 * r18178; float r18189 = r18188 * r18178; float r18190 = r18187 * r18189; float r18191 = r18185 + r18190; float r18192 = 21.0f; float r18193 = r18172 / r18192; float r18194 = r18189 * r18178; float r18195 = r18194 * r18178; float r18196 = r18193 * r18195; float r18197 = r18191 + r18196; float r18198 = r18175 * r18197; float r18199 = fabs(r18198); return r18199; } double f_id(double x) { double r18200 = 1.0; double r18201 = atan2(1.0, 0.0); double r18202 = sqrt(r18201); double r18203 = r18200 / r18202; double r18204 = 2.0; double r18205 = x; double r18206 = fabs(r18205); double r18207 = r18204 * r18206; double r18208 = 3.0; double r18209 = r18204 / r18208; double r18210 = r18206 * r18206; double r18211 = r18210 * r18206; double r18212 = r18209 * r18211; double r18213 = r18207 + r18212; double r18214 = 5.0; double r18215 = r18200 / r18214; double r18216 = r18211 * r18206; double r18217 = r18216 * r18206; double r18218 = r18215 * r18217; double r18219 = r18213 + r18218; double r18220 = 21.0; double r18221 = r18200 / r18220; double r18222 = r18217 * r18206; double r18223 = r18222 * r18206; double r18224 = r18221 * r18223; double r18225 = r18219 + r18224; double r18226 = r18203 * r18225; double r18227 = fabs(r18226); return r18227; } double f_of(float x) { float r18228 = 2.0f; float r18229 = x; float r18230 = fabs(r18229); float r18231 = r18228 * r18230; float r18232 = cbrt(r18230); float r18233 = 3.0f; float r18234 = r18228 / r18233; float r18235 = cbrt(r18234); float r18236 = r18232 * r18235; float r18237 = r18236 * (r18236 * r18236); float r18238 = r18230 * r18230; float r18239 = r18237 * r18238; float r18240 = r18231 + r18239; float r18241 = r18230 * (r18230 * r18230); float r18242 = r18241 * r18241; float r18243 = 21.0f; float r18244 = r18243 / r18230; float r18245 = r18242 / r18244; float r18246 = r18241 * r18238; float r18247 = 5.0f; float r18248 = r18246 / r18247; float r18249 = r18245 + r18248; float r18250 = r18240 + r18249; float r18251 = atan2(1.0, 0.0); float r18252 = sqrt(r18251); float r18253 = r18250 / r18252; float r18254 = fabs(r18253); return r18254; } double f_od(double x) { double r18255 = 2.0; double r18256 = x; double r18257 = fabs(r18256); double r18258 = r18255 * r18257; double r18259 = cbrt(r18257); double r18260 = 3.0; double r18261 = r18255 / r18260; double r18262 = cbrt(r18261); double r18263 = r18259 * r18262; double r18264 = r18263 * (r18263 * r18263); double r18265 = r18257 * r18257; double r18266 = r18264 * r18265; double r18267 = r18258 + r18266; double r18268 = r18257 * (r18257 * r18257); double r18269 = r18268 * r18268; double r18270 = 21.0; double r18271 = r18270 / r18257; double r18272 = r18269 / r18271; double r18273 = r18268 * r18265; double r18274 = 5.0; double r18275 = r18273 / r18274; double r18276 = r18272 + r18275; double r18277 = r18267 + r18276; double r18278 = atan2(1.0, 0.0); double r18279 = sqrt(r18278); double r18280 = r18277 / r18279; double r18281 = fabs(r18280); return r18281; } 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 r18282, r18283, r18284, r18285, r18286, r18287, r18288, r18289, r18290, r18291, r18292, r18293, r18294, r18295, r18296, r18297, r18298, r18299, r18300, r18301, r18302, r18303, r18304, r18305, r18306, r18307, r18308, r18309; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18282, "1", 10, MPFR_RNDN); mpfr_init(r18283); mpfr_init(r18284); mpfr_init(r18285); mpfr_init_set_str(r18286, "2", 10, MPFR_RNDN); mpfr_init(r18287); mpfr_init(r18288); mpfr_init(r18289); mpfr_init_set_str(r18290, "3", 10, MPFR_RNDN); mpfr_init(r18291); mpfr_init(r18292); mpfr_init(r18293); mpfr_init(r18294); mpfr_init(r18295); mpfr_init_set_str(r18296, "5", 10, MPFR_RNDN); mpfr_init(r18297); mpfr_init(r18298); mpfr_init(r18299); mpfr_init(r18300); mpfr_init(r18301); mpfr_init_set_str(r18302, "21", 10, MPFR_RNDN); mpfr_init(r18303); mpfr_init(r18304); mpfr_init(r18305); mpfr_init(r18306); mpfr_init(r18307); mpfr_init(r18308); mpfr_init(r18309); } double f_im(double x) { ; mpfr_const_pi(r18283, MPFR_RNDN); mpfr_sqrt(r18284, r18283, MPFR_RNDN); mpfr_div(r18285, r18282, r18284, MPFR_RNDN); ; mpfr_set_d(r18287, x, MPFR_RNDN); mpfr_abs(r18288, r18287, MPFR_RNDN); mpfr_mul(r18289, r18286, r18288, MPFR_RNDN); ; mpfr_div(r18291, r18286, r18290, MPFR_RNDN); mpfr_mul(r18292, r18288, r18288, MPFR_RNDN); mpfr_mul(r18293, r18292, r18288, MPFR_RNDN); mpfr_mul(r18294, r18291, r18293, MPFR_RNDN); mpfr_add(r18295, r18289, r18294, MPFR_RNDN); ; mpfr_div(r18297, r18282, r18296, MPFR_RNDN); mpfr_mul(r18298, r18293, r18288, MPFR_RNDN); mpfr_mul(r18299, r18298, r18288, MPFR_RNDN); mpfr_mul(r18300, r18297, r18299, MPFR_RNDN); mpfr_add(r18301, r18295, r18300, MPFR_RNDN); ; mpfr_div(r18303, r18282, r18302, MPFR_RNDN); mpfr_mul(r18304, r18299, r18288, MPFR_RNDN); mpfr_mul(r18305, r18304, r18288, MPFR_RNDN); mpfr_mul(r18306, r18303, r18305, MPFR_RNDN); mpfr_add(r18307, r18301, r18306, MPFR_RNDN); mpfr_mul(r18308, r18285, r18307, MPFR_RNDN); mpfr_abs(r18309, r18308, MPFR_RNDN); return mpfr_get_d(r18309, MPFR_RNDN); } static mpfr_t r18310, r18311, r18312, r18313, r18314, r18315, r18316, r18317, r18318, r18319, r18320, r18321, r18322, r18323, r18324, r18325, r18326, r18327, r18328, r18329, r18330, r18331, r18332, r18333, r18334, r18335, r18336; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18310, "2", 10, MPFR_RNDN); mpfr_init(r18311); mpfr_init(r18312); mpfr_init(r18313); mpfr_init(r18314); mpfr_init_set_str(r18315, "3", 10, MPFR_RNDN); mpfr_init(r18316); mpfr_init(r18317); mpfr_init(r18318); mpfr_init(r18319); mpfr_init(r18320); mpfr_init(r18321); mpfr_init(r18322); mpfr_init(r18323); mpfr_init(r18324); mpfr_init_set_str(r18325, "21", 10, MPFR_RNDN); mpfr_init(r18326); mpfr_init(r18327); mpfr_init(r18328); mpfr_init_set_str(r18329, "5", 10, MPFR_RNDN); mpfr_init(r18330); mpfr_init(r18331); mpfr_init(r18332); mpfr_init(r18333); mpfr_init(r18334); mpfr_init(r18335); mpfr_init(r18336); } double f_fm(double x) { ; mpfr_set_d(r18311, x, MPFR_RNDN); mpfr_abs(r18312, r18311, MPFR_RNDN); mpfr_mul(r18313, r18310, r18312, MPFR_RNDN); mpfr_cbrt(r18314, r18312, MPFR_RNDN); ; mpfr_div(r18316, r18310, r18315, MPFR_RNDN); mpfr_cbrt(r18317, r18316, MPFR_RNDN); mpfr_mul(r18318, r18314, r18317, MPFR_RNDN); mpfr_mul(r18319, r18318, r18318, MPFR_RNDN); mpfr_mul(r18319, r18319, r18318, MPFR_RNDN); mpfr_sqr(r18320, r18312, MPFR_RNDN); mpfr_mul(r18321, r18319, r18320, MPFR_RNDN); mpfr_add(r18322, r18313, r18321, MPFR_RNDN); mpfr_mul(r18323, r18312, r18312, MPFR_RNDN); mpfr_mul(r18323, r18323, r18312, MPFR_RNDN); mpfr_sqr(r18324, r18323, MPFR_RNDN); ; mpfr_div(r18326, r18325, r18312, MPFR_RNDN); mpfr_div(r18327, r18324, r18326, MPFR_RNDN); mpfr_mul(r18328, r18323, r18320, MPFR_RNDN); ; mpfr_div(r18330, r18328, r18329, MPFR_RNDN); mpfr_add(r18331, r18327, r18330, MPFR_RNDN); mpfr_add(r18332, r18322, r18331, MPFR_RNDN); mpfr_const_pi(r18333, MPFR_RNDN); mpfr_sqrt(r18334, r18333, MPFR_RNDN); mpfr_div(r18335, r18332, r18334, MPFR_RNDN); mpfr_abs(r18336, r18335, MPFR_RNDN); return mpfr_get_d(r18336, MPFR_RNDN); } static mpfr_t r18337, r18338, r18339, r18340, r18341, r18342, r18343, r18344, r18345, r18346, r18347, r18348, r18349, r18350, r18351, r18352, r18353, r18354, r18355, r18356, r18357, r18358, r18359, r18360, r18361, r18362, r18363; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18337, "2", 10, MPFR_RNDN); mpfr_init(r18338); mpfr_init(r18339); mpfr_init(r18340); mpfr_init(r18341); mpfr_init_set_str(r18342, "3", 10, MPFR_RNDN); mpfr_init(r18343); mpfr_init(r18344); mpfr_init(r18345); mpfr_init(r18346); mpfr_init(r18347); mpfr_init(r18348); mpfr_init(r18349); mpfr_init(r18350); mpfr_init(r18351); mpfr_init_set_str(r18352, "21", 10, MPFR_RNDN); mpfr_init(r18353); mpfr_init(r18354); mpfr_init(r18355); mpfr_init_set_str(r18356, "5", 10, MPFR_RNDN); mpfr_init(r18357); mpfr_init(r18358); mpfr_init(r18359); mpfr_init(r18360); mpfr_init(r18361); mpfr_init(r18362); mpfr_init(r18363); } double f_dm(double x) { ; mpfr_set_d(r18338, x, MPFR_RNDN); mpfr_abs(r18339, r18338, MPFR_RNDN); mpfr_mul(r18340, r18337, r18339, MPFR_RNDN); mpfr_cbrt(r18341, r18339, MPFR_RNDN); ; mpfr_div(r18343, r18337, r18342, MPFR_RNDN); mpfr_cbrt(r18344, r18343, MPFR_RNDN); mpfr_mul(r18345, r18341, r18344, MPFR_RNDN); mpfr_mul(r18346, r18345, r18345, MPFR_RNDN); mpfr_mul(r18346, r18346, r18345, MPFR_RNDN); mpfr_sqr(r18347, r18339, MPFR_RNDN); mpfr_mul(r18348, r18346, r18347, MPFR_RNDN); mpfr_add(r18349, r18340, r18348, MPFR_RNDN); mpfr_mul(r18350, r18339, r18339, MPFR_RNDN); mpfr_mul(r18350, r18350, r18339, MPFR_RNDN); mpfr_sqr(r18351, r18350, MPFR_RNDN); ; mpfr_div(r18353, r18352, r18339, MPFR_RNDN); mpfr_div(r18354, r18351, r18353, MPFR_RNDN); mpfr_mul(r18355, r18350, r18347, MPFR_RNDN); ; mpfr_div(r18357, r18355, r18356, MPFR_RNDN); mpfr_add(r18358, r18354, r18357, MPFR_RNDN); mpfr_add(r18359, r18349, r18358, MPFR_RNDN); mpfr_const_pi(r18360, MPFR_RNDN); mpfr_sqrt(r18361, r18360, MPFR_RNDN); mpfr_div(r18362, r18359, r18361, MPFR_RNDN); mpfr_abs(r18363, r18362, MPFR_RNDN); return mpfr_get_d(r18363, MPFR_RNDN); }