#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 r18476 = 1.0f; float r18477 = atan2(1.0, 0.0); float r18478 = sqrt(r18477); float r18479 = r18476 / r18478; float r18480 = 2.0f; float r18481 = x; float r18482 = fabs(r18481); float r18483 = r18480 * r18482; float r18484 = 3.0f; float r18485 = r18480 / r18484; float r18486 = r18482 * r18482; float r18487 = r18486 * r18482; float r18488 = r18485 * r18487; float r18489 = r18483 + r18488; float r18490 = 5.0f; float r18491 = r18476 / r18490; float r18492 = r18487 * r18482; float r18493 = r18492 * r18482; float r18494 = r18491 * r18493; float r18495 = r18489 + r18494; float r18496 = 21.0f; float r18497 = r18476 / r18496; float r18498 = r18493 * r18482; float r18499 = r18498 * r18482; float r18500 = r18497 * r18499; float r18501 = r18495 + r18500; float r18502 = r18479 * r18501; float r18503 = fabs(r18502); return r18503; } double f_id(double x) { double r18504 = 1.0; double r18505 = atan2(1.0, 0.0); double r18506 = sqrt(r18505); double r18507 = r18504 / r18506; double r18508 = 2.0; double r18509 = x; double r18510 = fabs(r18509); double r18511 = r18508 * r18510; double r18512 = 3.0; double r18513 = r18508 / r18512; double r18514 = r18510 * r18510; double r18515 = r18514 * r18510; double r18516 = r18513 * r18515; double r18517 = r18511 + r18516; double r18518 = 5.0; double r18519 = r18504 / r18518; double r18520 = r18515 * r18510; double r18521 = r18520 * r18510; double r18522 = r18519 * r18521; double r18523 = r18517 + r18522; double r18524 = 21.0; double r18525 = r18504 / r18524; double r18526 = r18521 * r18510; double r18527 = r18526 * r18510; double r18528 = r18525 * r18527; double r18529 = r18523 + r18528; double r18530 = r18507 * r18529; double r18531 = fabs(r18530); return r18531; } double f_of(float x) { float r18532 = 1.0f; float r18533 = atan2(1.0, 0.0); float r18534 = r18532 / r18533; float r18535 = sqrt(r18534); float r18536 = 0.047619047619047616f; float r18537 = x; float r18538 = fabs(r18537); float r18539 = r18538 * (r18538 * r18538); float r18540 = r18539 * r18539; float r18541 = r18538 * r18540; float r18542 = r18536 * r18541; float r18543 = 0.6666666666666666f; float r18544 = 3.0f; float r18545 = pow(r18538, r18544); float r18546 = r18543 * r18545; float r18547 = 2.0f; float r18548 = r18547 * r18538; float r18549 = 0.2f; float r18550 = r18538 * r18538; float r18551 = r18550 * r18539; float r18552 = r18549 * r18551; float r18553 = r18548 + r18552; float r18554 = r18546 + r18553; float r18555 = r18542 + r18554; float r18556 = r18535 * r18555; float r18557 = fabs(r18556); return r18557; } double f_od(double x) { double r18558 = 1.0; double r18559 = atan2(1.0, 0.0); double r18560 = r18558 / r18559; double r18561 = sqrt(r18560); double r18562 = 0.047619047619047616; double r18563 = x; double r18564 = fabs(r18563); double r18565 = r18564 * (r18564 * r18564); double r18566 = r18565 * r18565; double r18567 = r18564 * r18566; double r18568 = r18562 * r18567; double r18569 = 0.6666666666666666; double r18570 = 3.0; double r18571 = pow(r18564, r18570); double r18572 = r18569 * r18571; double r18573 = 2.0; double r18574 = r18573 * r18564; double r18575 = 0.2; double r18576 = r18564 * r18564; double r18577 = r18576 * r18565; double r18578 = r18575 * r18577; double r18579 = r18574 + r18578; double r18580 = r18572 + r18579; double r18581 = r18568 + r18580; double r18582 = r18561 * r18581; double r18583 = fabs(r18582); return r18583; } 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 r18584, r18585, r18586, r18587, r18588, r18589, r18590, r18591, r18592, r18593, r18594, r18595, r18596, r18597, r18598, r18599, r18600, r18601, r18602, r18603, r18604, r18605, r18606, r18607, r18608, r18609, r18610, r18611; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18584, "1", 10, MPFR_RNDN); mpfr_init(r18585); mpfr_init(r18586); mpfr_init(r18587); mpfr_init_set_str(r18588, "2", 10, MPFR_RNDN); mpfr_init(r18589); mpfr_init(r18590); mpfr_init(r18591); mpfr_init_set_str(r18592, "3", 10, MPFR_RNDN); mpfr_init(r18593); mpfr_init(r18594); mpfr_init(r18595); mpfr_init(r18596); mpfr_init(r18597); mpfr_init_set_str(r18598, "5", 10, MPFR_RNDN); mpfr_init(r18599); mpfr_init(r18600); mpfr_init(r18601); mpfr_init(r18602); mpfr_init(r18603); mpfr_init_set_str(r18604, "21", 10, MPFR_RNDN); mpfr_init(r18605); mpfr_init(r18606); mpfr_init(r18607); mpfr_init(r18608); mpfr_init(r18609); mpfr_init(r18610); mpfr_init(r18611); } double f_im(double x) { ; mpfr_const_pi(r18585, MPFR_RNDN); mpfr_sqrt(r18586, r18585, MPFR_RNDN); mpfr_div(r18587, r18584, r18586, MPFR_RNDN); ; mpfr_set_d(r18589, x, MPFR_RNDN); mpfr_abs(r18590, r18589, MPFR_RNDN); mpfr_mul(r18591, r18588, r18590, MPFR_RNDN); ; mpfr_div(r18593, r18588, r18592, MPFR_RNDN); mpfr_mul(r18594, r18590, r18590, MPFR_RNDN); mpfr_mul(r18595, r18594, r18590, MPFR_RNDN); mpfr_mul(r18596, r18593, r18595, MPFR_RNDN); mpfr_add(r18597, r18591, r18596, MPFR_RNDN); ; mpfr_div(r18599, r18584, r18598, MPFR_RNDN); mpfr_mul(r18600, r18595, r18590, MPFR_RNDN); mpfr_mul(r18601, r18600, r18590, MPFR_RNDN); mpfr_mul(r18602, r18599, r18601, MPFR_RNDN); mpfr_add(r18603, r18597, r18602, MPFR_RNDN); ; mpfr_div(r18605, r18584, r18604, MPFR_RNDN); mpfr_mul(r18606, r18601, r18590, MPFR_RNDN); mpfr_mul(r18607, r18606, r18590, MPFR_RNDN); mpfr_mul(r18608, r18605, r18607, MPFR_RNDN); mpfr_add(r18609, r18603, r18608, MPFR_RNDN); mpfr_mul(r18610, r18587, r18609, MPFR_RNDN); mpfr_abs(r18611, r18610, MPFR_RNDN); return mpfr_get_d(r18611, MPFR_RNDN); } static mpfr_t r18612, r18613, r18614, r18615, r18616, r18617, r18618, r18619, r18620, r18621, r18622, r18623, r18624, r18625, r18626, r18627, r18628, r18629, r18630, r18631, r18632, r18633, r18634, r18635, r18636, r18637; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18612, "1", 10, MPFR_RNDN); mpfr_init(r18613); mpfr_init(r18614); mpfr_init(r18615); mpfr_init_set_str(r18616, "1/21", 10, MPFR_RNDN); mpfr_init(r18617); mpfr_init(r18618); mpfr_init(r18619); mpfr_init(r18620); mpfr_init(r18621); mpfr_init(r18622); mpfr_init_set_str(r18623, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18624, "3", 10, MPFR_RNDN); mpfr_init(r18625); mpfr_init(r18626); mpfr_init_set_str(r18627, "2", 10, MPFR_RNDN); mpfr_init(r18628); mpfr_init_set_str(r18629, "1/5", 10, MPFR_RNDN); mpfr_init(r18630); mpfr_init(r18631); mpfr_init(r18632); mpfr_init(r18633); mpfr_init(r18634); mpfr_init(r18635); mpfr_init(r18636); mpfr_init(r18637); } double f_fm(double x) { ; mpfr_const_pi(r18613, MPFR_RNDN); mpfr_div(r18614, r18612, r18613, MPFR_RNDN); mpfr_sqrt(r18615, r18614, MPFR_RNDN); ; mpfr_set_d(r18617, x, MPFR_RNDN); mpfr_abs(r18618, r18617, MPFR_RNDN); mpfr_mul(r18619, r18618, r18618, MPFR_RNDN); mpfr_mul(r18619, r18619, r18618, MPFR_RNDN); mpfr_sqr(r18620, r18619, MPFR_RNDN); mpfr_mul(r18621, r18618, r18620, MPFR_RNDN); mpfr_mul(r18622, r18616, r18621, MPFR_RNDN); ; ; mpfr_pow(r18625, r18618, r18624, MPFR_RNDN); mpfr_mul(r18626, r18623, r18625, MPFR_RNDN); ; mpfr_mul(r18628, r18627, r18618, MPFR_RNDN); ; mpfr_sqr(r18630, r18618, MPFR_RNDN); mpfr_mul(r18631, r18630, r18619, MPFR_RNDN); mpfr_mul(r18632, r18629, r18631, MPFR_RNDN); mpfr_add(r18633, r18628, r18632, MPFR_RNDN); mpfr_add(r18634, r18626, r18633, MPFR_RNDN); mpfr_add(r18635, r18622, r18634, MPFR_RNDN); mpfr_mul(r18636, r18615, r18635, MPFR_RNDN); mpfr_abs(r18637, r18636, MPFR_RNDN); return mpfr_get_d(r18637, MPFR_RNDN); } static mpfr_t r18638, r18639, r18640, r18641, r18642, r18643, r18644, r18645, r18646, r18647, r18648, r18649, r18650, r18651, r18652, r18653, r18654, r18655, r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18638, "1", 10, MPFR_RNDN); mpfr_init(r18639); mpfr_init(r18640); mpfr_init(r18641); mpfr_init_set_str(r18642, "1/21", 10, MPFR_RNDN); mpfr_init(r18643); mpfr_init(r18644); mpfr_init(r18645); mpfr_init(r18646); mpfr_init(r18647); mpfr_init(r18648); mpfr_init_set_str(r18649, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18650, "3", 10, MPFR_RNDN); mpfr_init(r18651); mpfr_init(r18652); mpfr_init_set_str(r18653, "2", 10, MPFR_RNDN); mpfr_init(r18654); mpfr_init_set_str(r18655, "1/5", 10, MPFR_RNDN); mpfr_init(r18656); mpfr_init(r18657); mpfr_init(r18658); mpfr_init(r18659); mpfr_init(r18660); mpfr_init(r18661); mpfr_init(r18662); mpfr_init(r18663); } double f_dm(double x) { ; mpfr_const_pi(r18639, MPFR_RNDN); mpfr_div(r18640, r18638, r18639, MPFR_RNDN); mpfr_sqrt(r18641, r18640, MPFR_RNDN); ; mpfr_set_d(r18643, x, MPFR_RNDN); mpfr_abs(r18644, r18643, MPFR_RNDN); mpfr_mul(r18645, r18644, r18644, MPFR_RNDN); mpfr_mul(r18645, r18645, r18644, MPFR_RNDN); mpfr_sqr(r18646, r18645, MPFR_RNDN); mpfr_mul(r18647, r18644, r18646, MPFR_RNDN); mpfr_mul(r18648, r18642, r18647, MPFR_RNDN); ; ; mpfr_pow(r18651, r18644, r18650, MPFR_RNDN); mpfr_mul(r18652, r18649, r18651, MPFR_RNDN); ; mpfr_mul(r18654, r18653, r18644, MPFR_RNDN); ; mpfr_sqr(r18656, r18644, MPFR_RNDN); mpfr_mul(r18657, r18656, r18645, MPFR_RNDN); mpfr_mul(r18658, r18655, r18657, MPFR_RNDN); mpfr_add(r18659, r18654, r18658, MPFR_RNDN); mpfr_add(r18660, r18652, r18659, MPFR_RNDN); mpfr_add(r18661, r18648, r18660, MPFR_RNDN); mpfr_mul(r18662, r18641, r18661, MPFR_RNDN); mpfr_abs(r18663, r18662, MPFR_RNDN); return mpfr_get_d(r18663, MPFR_RNDN); }