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