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