#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 r18348 = 1.0f; float r18349 = atan2(1.0, 0.0); float r18350 = sqrt(r18349); float r18351 = r18348 / r18350; float r18352 = 2.0f; float r18353 = x; float r18354 = fabs(r18353); float r18355 = r18352 * r18354; float r18356 = 3.0f; float r18357 = r18352 / r18356; float r18358 = r18354 * r18354; float r18359 = r18358 * r18354; float r18360 = r18357 * r18359; float r18361 = r18355 + r18360; float r18362 = 5.0f; float r18363 = r18348 / r18362; float r18364 = r18359 * r18354; float r18365 = r18364 * r18354; float r18366 = r18363 * r18365; float r18367 = r18361 + r18366; float r18368 = 21.0f; float r18369 = r18348 / r18368; float r18370 = r18365 * r18354; float r18371 = r18370 * r18354; float r18372 = r18369 * r18371; float r18373 = r18367 + r18372; float r18374 = r18351 * r18373; float r18375 = fabs(r18374); return r18375; } double f_id(double x) { double r18376 = 1.0; double r18377 = atan2(1.0, 0.0); double r18378 = sqrt(r18377); double r18379 = r18376 / r18378; double r18380 = 2.0; double r18381 = x; double r18382 = fabs(r18381); double r18383 = r18380 * r18382; double r18384 = 3.0; double r18385 = r18380 / r18384; double r18386 = r18382 * r18382; double r18387 = r18386 * r18382; double r18388 = r18385 * r18387; double r18389 = r18383 + r18388; double r18390 = 5.0; double r18391 = r18376 / r18390; double r18392 = r18387 * r18382; double r18393 = r18392 * r18382; double r18394 = r18391 * r18393; double r18395 = r18389 + r18394; double r18396 = 21.0; double r18397 = r18376 / r18396; double r18398 = r18393 * r18382; double r18399 = r18398 * r18382; double r18400 = r18397 * r18399; double r18401 = r18395 + r18400; double r18402 = r18379 * r18401; double r18403 = fabs(r18402); return r18403; } double f_of(float x) { float r18404 = 2.0f; float r18405 = x; float r18406 = fabs(r18405); float r18407 = r18404 * r18406; float r18408 = 3.0f; float r18409 = r18404 / r18408; float r18410 = r18406 * r18409; float r18411 = r18406 * r18406; float r18412 = r18410 * r18411; float r18413 = r18407 + r18412; float r18414 = r18406 * (r18406 * r18406); float r18415 = pow(r18414, r18404); float r18416 = 21.0f; float r18417 = r18416 / r18406; float r18418 = r18415 / r18417; float r18419 = r18414 * r18411; float r18420 = 5.0f; float r18421 = r18419 / r18420; float r18422 = r18418 + r18421; float r18423 = r18413 + r18422; float r18424 = atan2(1.0, 0.0); float r18425 = sqrt(r18424); float r18426 = r18423 / r18425; float r18427 = fabs(r18426); return r18427; } double f_od(double x) { double r18428 = 2.0; double r18429 = x; double r18430 = fabs(r18429); double r18431 = r18428 * r18430; double r18432 = 3.0; double r18433 = r18428 / r18432; double r18434 = r18430 * r18433; double r18435 = r18430 * r18430; double r18436 = r18434 * r18435; double r18437 = r18431 + r18436; double r18438 = r18430 * (r18430 * r18430); double r18439 = pow(r18438, r18428); double r18440 = 21.0; double r18441 = r18440 / r18430; double r18442 = r18439 / r18441; double r18443 = r18438 * r18435; double r18444 = 5.0; double r18445 = r18443 / r18444; double r18446 = r18442 + r18445; double r18447 = r18437 + r18446; double r18448 = atan2(1.0, 0.0); double r18449 = sqrt(r18448); double r18450 = r18447 / r18449; double r18451 = fabs(r18450); return r18451; } 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 r18452, r18453, r18454, r18455, r18456, r18457, r18458, r18459, r18460, r18461, r18462, r18463, r18464, r18465, r18466, r18467, r18468, r18469, r18470, r18471, r18472, r18473, r18474, r18475, r18476, r18477, r18478, r18479; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18452, "1", 10, MPFR_RNDN); mpfr_init(r18453); mpfr_init(r18454); mpfr_init(r18455); mpfr_init_set_str(r18456, "2", 10, MPFR_RNDN); mpfr_init(r18457); mpfr_init(r18458); mpfr_init(r18459); mpfr_init_set_str(r18460, "3", 10, MPFR_RNDN); mpfr_init(r18461); mpfr_init(r18462); mpfr_init(r18463); mpfr_init(r18464); mpfr_init(r18465); mpfr_init_set_str(r18466, "5", 10, MPFR_RNDN); mpfr_init(r18467); mpfr_init(r18468); mpfr_init(r18469); mpfr_init(r18470); mpfr_init(r18471); mpfr_init_set_str(r18472, "21", 10, MPFR_RNDN); mpfr_init(r18473); mpfr_init(r18474); mpfr_init(r18475); mpfr_init(r18476); mpfr_init(r18477); mpfr_init(r18478); mpfr_init(r18479); } double f_im(double x) { ; mpfr_const_pi(r18453, MPFR_RNDN); mpfr_sqrt(r18454, r18453, MPFR_RNDN); mpfr_div(r18455, r18452, r18454, MPFR_RNDN); ; mpfr_set_d(r18457, x, MPFR_RNDN); mpfr_abs(r18458, r18457, MPFR_RNDN); mpfr_mul(r18459, r18456, r18458, MPFR_RNDN); ; mpfr_div(r18461, r18456, r18460, MPFR_RNDN); mpfr_mul(r18462, r18458, r18458, MPFR_RNDN); mpfr_mul(r18463, r18462, r18458, MPFR_RNDN); mpfr_mul(r18464, r18461, r18463, MPFR_RNDN); mpfr_add(r18465, r18459, r18464, MPFR_RNDN); ; mpfr_div(r18467, r18452, r18466, MPFR_RNDN); mpfr_mul(r18468, r18463, r18458, MPFR_RNDN); mpfr_mul(r18469, r18468, r18458, MPFR_RNDN); mpfr_mul(r18470, r18467, r18469, MPFR_RNDN); mpfr_add(r18471, r18465, r18470, MPFR_RNDN); ; mpfr_div(r18473, r18452, r18472, MPFR_RNDN); mpfr_mul(r18474, r18469, r18458, MPFR_RNDN); mpfr_mul(r18475, r18474, r18458, MPFR_RNDN); mpfr_mul(r18476, r18473, r18475, MPFR_RNDN); mpfr_add(r18477, r18471, r18476, MPFR_RNDN); mpfr_mul(r18478, r18455, r18477, MPFR_RNDN); mpfr_abs(r18479, r18478, MPFR_RNDN); return mpfr_get_d(r18479, MPFR_RNDN); } static mpfr_t r18480, r18481, r18482, r18483, r18484, r18485, r18486, r18487, r18488, r18489, r18490, r18491, r18492, r18493, r18494, r18495, r18496, r18497, r18498, r18499, r18500, r18501, r18502, r18503; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18480, "2", 10, MPFR_RNDN); mpfr_init(r18481); mpfr_init(r18482); mpfr_init(r18483); mpfr_init_set_str(r18484, "3", 10, MPFR_RNDN); mpfr_init(r18485); mpfr_init(r18486); mpfr_init(r18487); mpfr_init(r18488); mpfr_init(r18489); mpfr_init(r18490); mpfr_init(r18491); mpfr_init_set_str(r18492, "21", 10, MPFR_RNDN); mpfr_init(r18493); mpfr_init(r18494); mpfr_init(r18495); mpfr_init_set_str(r18496, "5", 10, MPFR_RNDN); mpfr_init(r18497); mpfr_init(r18498); mpfr_init(r18499); mpfr_init(r18500); mpfr_init(r18501); mpfr_init(r18502); mpfr_init(r18503); } double f_fm(double x) { ; mpfr_set_d(r18481, x, MPFR_RNDN); mpfr_abs(r18482, r18481, MPFR_RNDN); mpfr_mul(r18483, r18480, r18482, MPFR_RNDN); ; mpfr_div(r18485, r18480, r18484, MPFR_RNDN); mpfr_mul(r18486, r18482, r18485, MPFR_RNDN); mpfr_sqr(r18487, r18482, MPFR_RNDN); mpfr_mul(r18488, r18486, r18487, MPFR_RNDN); mpfr_add(r18489, r18483, r18488, MPFR_RNDN); mpfr_mul(r18490, r18482, r18482, MPFR_RNDN); mpfr_mul(r18490, r18490, r18482, MPFR_RNDN); mpfr_pow(r18491, r18490, r18480, MPFR_RNDN); ; mpfr_div(r18493, r18492, r18482, MPFR_RNDN); mpfr_div(r18494, r18491, r18493, MPFR_RNDN); mpfr_mul(r18495, r18490, r18487, MPFR_RNDN); ; mpfr_div(r18497, r18495, r18496, MPFR_RNDN); mpfr_add(r18498, r18494, r18497, MPFR_RNDN); mpfr_add(r18499, r18489, r18498, MPFR_RNDN); mpfr_const_pi(r18500, MPFR_RNDN); mpfr_sqrt(r18501, r18500, MPFR_RNDN); mpfr_div(r18502, r18499, r18501, MPFR_RNDN); mpfr_abs(r18503, r18502, MPFR_RNDN); return mpfr_get_d(r18503, MPFR_RNDN); } static mpfr_t r18504, r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513, r18514, r18515, r18516, r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525, r18526, r18527; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18504, "2", 10, MPFR_RNDN); mpfr_init(r18505); mpfr_init(r18506); mpfr_init(r18507); mpfr_init_set_str(r18508, "3", 10, MPFR_RNDN); mpfr_init(r18509); mpfr_init(r18510); mpfr_init(r18511); mpfr_init(r18512); mpfr_init(r18513); mpfr_init(r18514); mpfr_init(r18515); mpfr_init_set_str(r18516, "21", 10, MPFR_RNDN); mpfr_init(r18517); mpfr_init(r18518); mpfr_init(r18519); mpfr_init_set_str(r18520, "5", 10, MPFR_RNDN); mpfr_init(r18521); mpfr_init(r18522); mpfr_init(r18523); mpfr_init(r18524); mpfr_init(r18525); mpfr_init(r18526); mpfr_init(r18527); } double f_dm(double x) { ; mpfr_set_d(r18505, x, MPFR_RNDN); mpfr_abs(r18506, r18505, MPFR_RNDN); mpfr_mul(r18507, r18504, r18506, MPFR_RNDN); ; mpfr_div(r18509, r18504, r18508, MPFR_RNDN); mpfr_mul(r18510, r18506, r18509, MPFR_RNDN); mpfr_sqr(r18511, r18506, MPFR_RNDN); mpfr_mul(r18512, r18510, r18511, MPFR_RNDN); mpfr_add(r18513, r18507, r18512, MPFR_RNDN); mpfr_mul(r18514, r18506, r18506, MPFR_RNDN); mpfr_mul(r18514, r18514, r18506, MPFR_RNDN); mpfr_pow(r18515, r18514, r18504, MPFR_RNDN); ; mpfr_div(r18517, r18516, r18506, MPFR_RNDN); mpfr_div(r18518, r18515, r18517, MPFR_RNDN); mpfr_mul(r18519, r18514, r18511, MPFR_RNDN); ; mpfr_div(r18521, r18519, r18520, MPFR_RNDN); mpfr_add(r18522, r18518, r18521, MPFR_RNDN); mpfr_add(r18523, r18513, r18522, MPFR_RNDN); mpfr_const_pi(r18524, MPFR_RNDN); mpfr_sqrt(r18525, r18524, MPFR_RNDN); mpfr_div(r18526, r18523, r18525, MPFR_RNDN); mpfr_abs(r18527, r18526, MPFR_RNDN); return mpfr_get_d(r18527, MPFR_RNDN); }