#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 r18380 = 1.0f; float r18381 = atan2(1.0, 0.0); float r18382 = sqrt(r18381); float r18383 = r18380 / r18382; float r18384 = 2.0f; float r18385 = x; float r18386 = fabs(r18385); float r18387 = r18384 * r18386; float r18388 = 3.0f; float r18389 = r18384 / r18388; float r18390 = r18386 * r18386; float r18391 = r18390 * r18386; float r18392 = r18389 * r18391; float r18393 = r18387 + r18392; float r18394 = 5.0f; float r18395 = r18380 / r18394; float r18396 = r18391 * r18386; float r18397 = r18396 * r18386; float r18398 = r18395 * r18397; float r18399 = r18393 + r18398; float r18400 = 21.0f; float r18401 = r18380 / r18400; float r18402 = r18397 * r18386; float r18403 = r18402 * r18386; float r18404 = r18401 * r18403; float r18405 = r18399 + r18404; float r18406 = r18383 * r18405; float r18407 = fabs(r18406); return r18407; } double f_id(double x) { double r18408 = 1.0; double r18409 = atan2(1.0, 0.0); double r18410 = sqrt(r18409); double r18411 = r18408 / r18410; double r18412 = 2.0; double r18413 = x; double r18414 = fabs(r18413); double r18415 = r18412 * r18414; double r18416 = 3.0; double r18417 = r18412 / r18416; double r18418 = r18414 * r18414; double r18419 = r18418 * r18414; double r18420 = r18417 * r18419; double r18421 = r18415 + r18420; double r18422 = 5.0; double r18423 = r18408 / r18422; double r18424 = r18419 * r18414; double r18425 = r18424 * r18414; double r18426 = r18423 * r18425; double r18427 = r18421 + r18426; double r18428 = 21.0; double r18429 = r18408 / r18428; double r18430 = r18425 * r18414; double r18431 = r18430 * r18414; double r18432 = r18429 * r18431; double r18433 = r18427 + r18432; double r18434 = r18411 * r18433; double r18435 = fabs(r18434); return r18435; } double f_of(float x) { float r18436 = 1.0f; float r18437 = atan2(1.0, 0.0); float r18438 = r18436 / r18437; float r18439 = sqrt(r18438); float r18440 = 0.047619047619047616f; float r18441 = x; float r18442 = fabs(r18441); float r18443 = r18442 * (r18442 * r18442); float r18444 = r18443 * r18443; float r18445 = r18442 * r18444; float r18446 = r18440 * r18445; float r18447 = 0.6666666666666666f; float r18448 = 3.0f; float r18449 = pow(r18442, r18448); float r18450 = r18447 * r18449; float r18451 = 2.0f; float r18452 = r18451 * r18442; float r18453 = 0.2f; float r18454 = r18442 * r18442; float r18455 = r18454 * r18443; float r18456 = r18453 * r18455; float r18457 = r18452 + r18456; float r18458 = r18450 + r18457; float r18459 = r18446 + r18458; float r18460 = r18439 * r18459; float r18461 = fabs(r18460); return r18461; } double f_od(double x) { double r18462 = 1.0; double r18463 = atan2(1.0, 0.0); double r18464 = r18462 / r18463; double r18465 = sqrt(r18464); double r18466 = 0.047619047619047616; double r18467 = x; double r18468 = fabs(r18467); double r18469 = r18468 * (r18468 * r18468); double r18470 = r18469 * r18469; double r18471 = r18468 * r18470; double r18472 = r18466 * r18471; double r18473 = 0.6666666666666666; double r18474 = 3.0; double r18475 = pow(r18468, r18474); double r18476 = r18473 * r18475; double r18477 = 2.0; double r18478 = r18477 * r18468; double r18479 = 0.2; double r18480 = r18468 * r18468; double r18481 = r18480 * r18469; double r18482 = r18479 * r18481; double r18483 = r18478 + r18482; double r18484 = r18476 + r18483; double r18485 = r18472 + r18484; double r18486 = r18465 * r18485; double r18487 = fabs(r18486); return r18487; } 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 r18488, r18489, r18490, r18491, r18492, r18493, r18494, r18495, r18496, r18497, r18498, r18499, r18500, r18501, r18502, r18503, r18504, r18505, r18506, r18507, r18508, r18509, r18510, r18511, r18512, r18513, r18514, r18515; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18488, "1", 10, MPFR_RNDN); mpfr_init(r18489); mpfr_init(r18490); mpfr_init(r18491); mpfr_init_set_str(r18492, "2", 10, MPFR_RNDN); mpfr_init(r18493); mpfr_init(r18494); mpfr_init(r18495); mpfr_init_set_str(r18496, "3", 10, MPFR_RNDN); mpfr_init(r18497); mpfr_init(r18498); mpfr_init(r18499); mpfr_init(r18500); mpfr_init(r18501); mpfr_init_set_str(r18502, "5", 10, MPFR_RNDN); mpfr_init(r18503); mpfr_init(r18504); mpfr_init(r18505); mpfr_init(r18506); mpfr_init(r18507); mpfr_init_set_str(r18508, "21", 10, MPFR_RNDN); mpfr_init(r18509); mpfr_init(r18510); mpfr_init(r18511); mpfr_init(r18512); mpfr_init(r18513); mpfr_init(r18514); mpfr_init(r18515); } double f_im(double x) { ; mpfr_const_pi(r18489, MPFR_RNDN); mpfr_sqrt(r18490, r18489, MPFR_RNDN); mpfr_div(r18491, r18488, r18490, MPFR_RNDN); ; mpfr_set_d(r18493, x, MPFR_RNDN); mpfr_abs(r18494, r18493, MPFR_RNDN); mpfr_mul(r18495, r18492, r18494, MPFR_RNDN); ; mpfr_div(r18497, r18492, r18496, MPFR_RNDN); mpfr_mul(r18498, r18494, r18494, MPFR_RNDN); mpfr_mul(r18499, r18498, r18494, MPFR_RNDN); mpfr_mul(r18500, r18497, r18499, MPFR_RNDN); mpfr_add(r18501, r18495, r18500, MPFR_RNDN); ; mpfr_div(r18503, r18488, r18502, MPFR_RNDN); mpfr_mul(r18504, r18499, r18494, MPFR_RNDN); mpfr_mul(r18505, r18504, r18494, MPFR_RNDN); mpfr_mul(r18506, r18503, r18505, MPFR_RNDN); mpfr_add(r18507, r18501, r18506, MPFR_RNDN); ; mpfr_div(r18509, r18488, r18508, MPFR_RNDN); mpfr_mul(r18510, r18505, r18494, MPFR_RNDN); mpfr_mul(r18511, r18510, r18494, MPFR_RNDN); mpfr_mul(r18512, r18509, r18511, MPFR_RNDN); mpfr_add(r18513, r18507, r18512, MPFR_RNDN); mpfr_mul(r18514, r18491, r18513, MPFR_RNDN); mpfr_abs(r18515, r18514, MPFR_RNDN); return mpfr_get_d(r18515, MPFR_RNDN); } static mpfr_t r18516, r18517, r18518, r18519, r18520, r18521, r18522, r18523, r18524, r18525, r18526, r18527, r18528, r18529, r18530, r18531, r18532, r18533, r18534, r18535, r18536, r18537, r18538, r18539, r18540, r18541; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18516, "1", 10, MPFR_RNDN); mpfr_init(r18517); mpfr_init(r18518); mpfr_init(r18519); mpfr_init_set_str(r18520, "1/21", 10, MPFR_RNDN); mpfr_init(r18521); mpfr_init(r18522); mpfr_init(r18523); mpfr_init(r18524); mpfr_init(r18525); mpfr_init(r18526); mpfr_init_set_str(r18527, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18528, "3", 10, MPFR_RNDN); mpfr_init(r18529); mpfr_init(r18530); mpfr_init_set_str(r18531, "2", 10, MPFR_RNDN); mpfr_init(r18532); mpfr_init_set_str(r18533, "1/5", 10, MPFR_RNDN); mpfr_init(r18534); mpfr_init(r18535); mpfr_init(r18536); mpfr_init(r18537); mpfr_init(r18538); mpfr_init(r18539); mpfr_init(r18540); mpfr_init(r18541); } double f_fm(double x) { ; mpfr_const_pi(r18517, MPFR_RNDN); mpfr_div(r18518, r18516, r18517, MPFR_RNDN); mpfr_sqrt(r18519, r18518, MPFR_RNDN); ; mpfr_set_d(r18521, x, MPFR_RNDN); mpfr_abs(r18522, r18521, MPFR_RNDN); mpfr_mul(r18523, r18522, r18522, MPFR_RNDN); mpfr_mul(r18523, r18523, r18522, MPFR_RNDN); mpfr_sqr(r18524, r18523, MPFR_RNDN); mpfr_mul(r18525, r18522, r18524, MPFR_RNDN); mpfr_mul(r18526, r18520, r18525, MPFR_RNDN); ; ; mpfr_pow(r18529, r18522, r18528, MPFR_RNDN); mpfr_mul(r18530, r18527, r18529, MPFR_RNDN); ; mpfr_mul(r18532, r18531, r18522, MPFR_RNDN); ; mpfr_sqr(r18534, r18522, MPFR_RNDN); mpfr_mul(r18535, r18534, r18523, MPFR_RNDN); mpfr_mul(r18536, r18533, r18535, MPFR_RNDN); mpfr_add(r18537, r18532, r18536, MPFR_RNDN); mpfr_add(r18538, r18530, r18537, MPFR_RNDN); mpfr_add(r18539, r18526, r18538, MPFR_RNDN); mpfr_mul(r18540, r18519, r18539, MPFR_RNDN); mpfr_abs(r18541, r18540, MPFR_RNDN); return mpfr_get_d(r18541, MPFR_RNDN); } static mpfr_t r18542, r18543, r18544, r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558, r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18542, "1", 10, MPFR_RNDN); mpfr_init(r18543); mpfr_init(r18544); mpfr_init(r18545); mpfr_init_set_str(r18546, "1/21", 10, MPFR_RNDN); mpfr_init(r18547); mpfr_init(r18548); mpfr_init(r18549); mpfr_init(r18550); mpfr_init(r18551); mpfr_init(r18552); mpfr_init_set_str(r18553, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18554, "3", 10, MPFR_RNDN); mpfr_init(r18555); mpfr_init(r18556); mpfr_init_set_str(r18557, "2", 10, MPFR_RNDN); mpfr_init(r18558); mpfr_init_set_str(r18559, "1/5", 10, MPFR_RNDN); mpfr_init(r18560); mpfr_init(r18561); mpfr_init(r18562); mpfr_init(r18563); mpfr_init(r18564); mpfr_init(r18565); mpfr_init(r18566); mpfr_init(r18567); } double f_dm(double x) { ; mpfr_const_pi(r18543, MPFR_RNDN); mpfr_div(r18544, r18542, r18543, MPFR_RNDN); mpfr_sqrt(r18545, r18544, MPFR_RNDN); ; mpfr_set_d(r18547, x, MPFR_RNDN); mpfr_abs(r18548, r18547, MPFR_RNDN); mpfr_mul(r18549, r18548, r18548, MPFR_RNDN); mpfr_mul(r18549, r18549, r18548, MPFR_RNDN); mpfr_sqr(r18550, r18549, MPFR_RNDN); mpfr_mul(r18551, r18548, r18550, MPFR_RNDN); mpfr_mul(r18552, r18546, r18551, MPFR_RNDN); ; ; mpfr_pow(r18555, r18548, r18554, MPFR_RNDN); mpfr_mul(r18556, r18553, r18555, MPFR_RNDN); ; mpfr_mul(r18558, r18557, r18548, MPFR_RNDN); ; mpfr_sqr(r18560, r18548, MPFR_RNDN); mpfr_mul(r18561, r18560, r18549, MPFR_RNDN); mpfr_mul(r18562, r18559, r18561, MPFR_RNDN); mpfr_add(r18563, r18558, r18562, MPFR_RNDN); mpfr_add(r18564, r18556, r18563, MPFR_RNDN); mpfr_add(r18565, r18552, r18564, MPFR_RNDN); mpfr_mul(r18566, r18545, r18565, MPFR_RNDN); mpfr_abs(r18567, r18566, MPFR_RNDN); return mpfr_get_d(r18567, MPFR_RNDN); }