#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 r18253 = 1.0f; float r18254 = atan2(1.0, 0.0); float r18255 = sqrt(r18254); float r18256 = r18253 / r18255; float r18257 = 2.0f; float r18258 = x; float r18259 = fabs(r18258); float r18260 = r18257 * r18259; float r18261 = 3.0f; float r18262 = r18257 / r18261; float r18263 = r18259 * r18259; float r18264 = r18263 * r18259; float r18265 = r18262 * r18264; float r18266 = r18260 + r18265; float r18267 = 5.0f; float r18268 = r18253 / r18267; float r18269 = r18264 * r18259; float r18270 = r18269 * r18259; float r18271 = r18268 * r18270; float r18272 = r18266 + r18271; float r18273 = 21.0f; float r18274 = r18253 / r18273; float r18275 = r18270 * r18259; float r18276 = r18275 * r18259; float r18277 = r18274 * r18276; float r18278 = r18272 + r18277; float r18279 = r18256 * r18278; float r18280 = fabs(r18279); return r18280; } double f_id(double x) { double r18281 = 1.0; double r18282 = atan2(1.0, 0.0); double r18283 = sqrt(r18282); double r18284 = r18281 / r18283; double r18285 = 2.0; double r18286 = x; double r18287 = fabs(r18286); double r18288 = r18285 * r18287; double r18289 = 3.0; double r18290 = r18285 / r18289; double r18291 = r18287 * r18287; double r18292 = r18291 * r18287; double r18293 = r18290 * r18292; double r18294 = r18288 + r18293; double r18295 = 5.0; double r18296 = r18281 / r18295; double r18297 = r18292 * r18287; double r18298 = r18297 * r18287; double r18299 = r18296 * r18298; double r18300 = r18294 + r18299; double r18301 = 21.0; double r18302 = r18281 / r18301; double r18303 = r18298 * r18287; double r18304 = r18303 * r18287; double r18305 = r18302 * r18304; double r18306 = r18300 + r18305; double r18307 = r18284 * r18306; double r18308 = fabs(r18307); return r18308; } double f_of(float x) { float r18309 = 2.0f; float r18310 = x; float r18311 = fabs(r18310); float r18312 = r18309 * r18311; float r18313 = cbrt(r18311); float r18314 = 3.0f; float r18315 = r18309 / r18314; float r18316 = cbrt(r18315); float r18317 = r18313 * r18316; float r18318 = r18317 * (r18317 * r18317); float r18319 = r18311 * r18311; float r18320 = r18318 * r18319; float r18321 = r18312 + r18320; float r18322 = r18311 * (r18311 * r18311); float r18323 = r18322 * r18322; float r18324 = 21.0f; float r18325 = r18324 / r18311; float r18326 = r18323 / r18325; float r18327 = r18322 * r18319; float r18328 = 5.0f; float r18329 = r18327 / r18328; float r18330 = r18326 + r18329; float r18331 = r18321 + r18330; float r18332 = atan2(1.0, 0.0); float r18333 = sqrt(r18332); float r18334 = r18331 / r18333; float r18335 = fabs(r18334); return r18335; } double f_od(double x) { double r18336 = 2.0; double r18337 = x; double r18338 = fabs(r18337); double r18339 = r18336 * r18338; double r18340 = cbrt(r18338); double r18341 = 3.0; double r18342 = r18336 / r18341; double r18343 = cbrt(r18342); double r18344 = r18340 * r18343; double r18345 = r18344 * (r18344 * r18344); double r18346 = r18338 * r18338; double r18347 = r18345 * r18346; double r18348 = r18339 + r18347; double r18349 = r18338 * (r18338 * r18338); double r18350 = r18349 * r18349; double r18351 = 21.0; double r18352 = r18351 / r18338; double r18353 = r18350 / r18352; double r18354 = r18349 * r18346; double r18355 = 5.0; double r18356 = r18354 / r18355; double r18357 = r18353 + r18356; double r18358 = r18348 + r18357; double r18359 = atan2(1.0, 0.0); double r18360 = sqrt(r18359); double r18361 = r18358 / r18360; double r18362 = fabs(r18361); return r18362; } 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 r18363, r18364, r18365, r18366, r18367, r18368, r18369, r18370, r18371, r18372, r18373, r18374, r18375, r18376, r18377, r18378, r18379, r18380, r18381, r18382, r18383, r18384, r18385, r18386, r18387, r18388, r18389, r18390; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18363, "1", 10, MPFR_RNDN); mpfr_init(r18364); mpfr_init(r18365); mpfr_init(r18366); mpfr_init_set_str(r18367, "2", 10, MPFR_RNDN); mpfr_init(r18368); mpfr_init(r18369); mpfr_init(r18370); mpfr_init_set_str(r18371, "3", 10, MPFR_RNDN); mpfr_init(r18372); mpfr_init(r18373); mpfr_init(r18374); mpfr_init(r18375); mpfr_init(r18376); mpfr_init_set_str(r18377, "5", 10, MPFR_RNDN); mpfr_init(r18378); mpfr_init(r18379); mpfr_init(r18380); mpfr_init(r18381); mpfr_init(r18382); mpfr_init_set_str(r18383, "21", 10, MPFR_RNDN); mpfr_init(r18384); mpfr_init(r18385); mpfr_init(r18386); mpfr_init(r18387); mpfr_init(r18388); mpfr_init(r18389); mpfr_init(r18390); } double f_im(double x) { ; mpfr_const_pi(r18364, MPFR_RNDN); mpfr_sqrt(r18365, r18364, MPFR_RNDN); mpfr_div(r18366, r18363, r18365, MPFR_RNDN); ; mpfr_set_d(r18368, x, MPFR_RNDN); mpfr_abs(r18369, r18368, MPFR_RNDN); mpfr_mul(r18370, r18367, r18369, MPFR_RNDN); ; mpfr_div(r18372, r18367, r18371, MPFR_RNDN); mpfr_mul(r18373, r18369, r18369, MPFR_RNDN); mpfr_mul(r18374, r18373, r18369, MPFR_RNDN); mpfr_mul(r18375, r18372, r18374, MPFR_RNDN); mpfr_add(r18376, r18370, r18375, MPFR_RNDN); ; mpfr_div(r18378, r18363, r18377, MPFR_RNDN); mpfr_mul(r18379, r18374, r18369, MPFR_RNDN); mpfr_mul(r18380, r18379, r18369, MPFR_RNDN); mpfr_mul(r18381, r18378, r18380, MPFR_RNDN); mpfr_add(r18382, r18376, r18381, MPFR_RNDN); ; mpfr_div(r18384, r18363, r18383, MPFR_RNDN); mpfr_mul(r18385, r18380, r18369, MPFR_RNDN); mpfr_mul(r18386, r18385, r18369, MPFR_RNDN); mpfr_mul(r18387, r18384, r18386, MPFR_RNDN); mpfr_add(r18388, r18382, r18387, MPFR_RNDN); mpfr_mul(r18389, r18366, r18388, MPFR_RNDN); mpfr_abs(r18390, r18389, MPFR_RNDN); return mpfr_get_d(r18390, MPFR_RNDN); } static mpfr_t r18391, r18392, r18393, r18394, r18395, r18396, r18397, r18398, r18399, r18400, r18401, r18402, r18403, r18404, r18405, r18406, r18407, r18408, r18409, r18410, r18411, r18412, r18413, r18414, r18415, r18416, r18417; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18391, "2", 10, MPFR_RNDN); mpfr_init(r18392); mpfr_init(r18393); mpfr_init(r18394); mpfr_init(r18395); mpfr_init_set_str(r18396, "3", 10, MPFR_RNDN); mpfr_init(r18397); mpfr_init(r18398); mpfr_init(r18399); mpfr_init(r18400); mpfr_init(r18401); mpfr_init(r18402); mpfr_init(r18403); mpfr_init(r18404); mpfr_init(r18405); mpfr_init_set_str(r18406, "21", 10, MPFR_RNDN); mpfr_init(r18407); mpfr_init(r18408); mpfr_init(r18409); mpfr_init_set_str(r18410, "5", 10, MPFR_RNDN); mpfr_init(r18411); mpfr_init(r18412); mpfr_init(r18413); mpfr_init(r18414); mpfr_init(r18415); mpfr_init(r18416); mpfr_init(r18417); } double f_fm(double x) { ; mpfr_set_d(r18392, x, MPFR_RNDN); mpfr_abs(r18393, r18392, MPFR_RNDN); mpfr_mul(r18394, r18391, r18393, MPFR_RNDN); mpfr_cbrt(r18395, r18393, MPFR_RNDN); ; mpfr_div(r18397, r18391, r18396, MPFR_RNDN); mpfr_cbrt(r18398, r18397, MPFR_RNDN); mpfr_mul(r18399, r18395, r18398, MPFR_RNDN); mpfr_mul(r18400, r18399, r18399, MPFR_RNDN); mpfr_mul(r18400, r18400, r18399, MPFR_RNDN); mpfr_sqr(r18401, r18393, MPFR_RNDN); mpfr_mul(r18402, r18400, r18401, MPFR_RNDN); mpfr_add(r18403, r18394, r18402, MPFR_RNDN); mpfr_mul(r18404, r18393, r18393, MPFR_RNDN); mpfr_mul(r18404, r18404, r18393, MPFR_RNDN); mpfr_sqr(r18405, r18404, MPFR_RNDN); ; mpfr_div(r18407, r18406, r18393, MPFR_RNDN); mpfr_div(r18408, r18405, r18407, MPFR_RNDN); mpfr_mul(r18409, r18404, r18401, MPFR_RNDN); ; mpfr_div(r18411, r18409, r18410, MPFR_RNDN); mpfr_add(r18412, r18408, r18411, MPFR_RNDN); mpfr_add(r18413, r18403, r18412, MPFR_RNDN); mpfr_const_pi(r18414, MPFR_RNDN); mpfr_sqrt(r18415, r18414, MPFR_RNDN); mpfr_div(r18416, r18413, r18415, MPFR_RNDN); mpfr_abs(r18417, r18416, MPFR_RNDN); return mpfr_get_d(r18417, MPFR_RNDN); } static mpfr_t r18418, r18419, r18420, r18421, r18422, r18423, r18424, r18425, r18426, r18427, r18428, r18429, r18430, r18431, r18432, r18433, r18434, r18435, r18436, r18437, r18438, r18439, r18440, r18441, r18442, r18443, r18444; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18418, "2", 10, MPFR_RNDN); mpfr_init(r18419); mpfr_init(r18420); mpfr_init(r18421); mpfr_init(r18422); mpfr_init_set_str(r18423, "3", 10, MPFR_RNDN); mpfr_init(r18424); mpfr_init(r18425); mpfr_init(r18426); mpfr_init(r18427); mpfr_init(r18428); mpfr_init(r18429); mpfr_init(r18430); mpfr_init(r18431); mpfr_init(r18432); mpfr_init_set_str(r18433, "21", 10, MPFR_RNDN); mpfr_init(r18434); mpfr_init(r18435); mpfr_init(r18436); mpfr_init_set_str(r18437, "5", 10, MPFR_RNDN); mpfr_init(r18438); mpfr_init(r18439); mpfr_init(r18440); mpfr_init(r18441); mpfr_init(r18442); mpfr_init(r18443); mpfr_init(r18444); } double f_dm(double x) { ; mpfr_set_d(r18419, x, MPFR_RNDN); mpfr_abs(r18420, r18419, MPFR_RNDN); mpfr_mul(r18421, r18418, r18420, MPFR_RNDN); mpfr_cbrt(r18422, r18420, MPFR_RNDN); ; mpfr_div(r18424, r18418, r18423, MPFR_RNDN); mpfr_cbrt(r18425, r18424, MPFR_RNDN); mpfr_mul(r18426, r18422, r18425, MPFR_RNDN); mpfr_mul(r18427, r18426, r18426, MPFR_RNDN); mpfr_mul(r18427, r18427, r18426, MPFR_RNDN); mpfr_sqr(r18428, r18420, MPFR_RNDN); mpfr_mul(r18429, r18427, r18428, MPFR_RNDN); mpfr_add(r18430, r18421, r18429, MPFR_RNDN); mpfr_mul(r18431, r18420, r18420, MPFR_RNDN); mpfr_mul(r18431, r18431, r18420, MPFR_RNDN); mpfr_sqr(r18432, r18431, MPFR_RNDN); ; mpfr_div(r18434, r18433, r18420, MPFR_RNDN); mpfr_div(r18435, r18432, r18434, MPFR_RNDN); mpfr_mul(r18436, r18431, r18428, MPFR_RNDN); ; mpfr_div(r18438, r18436, r18437, MPFR_RNDN); mpfr_add(r18439, r18435, r18438, MPFR_RNDN); mpfr_add(r18440, r18430, r18439, MPFR_RNDN); mpfr_const_pi(r18441, MPFR_RNDN); mpfr_sqrt(r18442, r18441, MPFR_RNDN); mpfr_div(r18443, r18440, r18442, MPFR_RNDN); mpfr_abs(r18444, r18443, MPFR_RNDN); return mpfr_get_d(r18444, MPFR_RNDN); }