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