#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 r18548 = 1.0f; float r18549 = atan2(1.0, 0.0); float r18550 = sqrt(r18549); float r18551 = r18548 / r18550; float r18552 = 2.0f; float r18553 = x; float r18554 = fabs(r18553); float r18555 = r18552 * r18554; float r18556 = 3.0f; float r18557 = r18552 / r18556; float r18558 = r18554 * r18554; float r18559 = r18558 * r18554; float r18560 = r18557 * r18559; float r18561 = r18555 + r18560; float r18562 = 5.0f; float r18563 = r18548 / r18562; float r18564 = r18559 * r18554; float r18565 = r18564 * r18554; float r18566 = r18563 * r18565; float r18567 = r18561 + r18566; float r18568 = 21.0f; float r18569 = r18548 / r18568; float r18570 = r18565 * r18554; float r18571 = r18570 * r18554; float r18572 = r18569 * r18571; float r18573 = r18567 + r18572; float r18574 = r18551 * r18573; float r18575 = fabs(r18574); return r18575; } double f_id(double x) { double r18576 = 1.0; double r18577 = atan2(1.0, 0.0); double r18578 = sqrt(r18577); double r18579 = r18576 / r18578; double r18580 = 2.0; double r18581 = x; double r18582 = fabs(r18581); double r18583 = r18580 * r18582; double r18584 = 3.0; double r18585 = r18580 / r18584; double r18586 = r18582 * r18582; double r18587 = r18586 * r18582; double r18588 = r18585 * r18587; double r18589 = r18583 + r18588; double r18590 = 5.0; double r18591 = r18576 / r18590; double r18592 = r18587 * r18582; double r18593 = r18592 * r18582; double r18594 = r18591 * r18593; double r18595 = r18589 + r18594; double r18596 = 21.0; double r18597 = r18576 / r18596; double r18598 = r18593 * r18582; double r18599 = r18598 * r18582; double r18600 = r18597 * r18599; double r18601 = r18595 + r18600; double r18602 = r18579 * r18601; double r18603 = fabs(r18602); return r18603; } double f_of(float x) { float r18604 = 2.0f; float r18605 = x; float r18606 = fabs(r18605); float r18607 = r18604 * r18606; float r18608 = 3.0f; float r18609 = r18604 / r18608; float r18610 = r18606 * r18609; float r18611 = r18606 * r18606; float r18612 = r18610 * r18611; float r18613 = r18607 + r18612; float r18614 = r18606 * (r18606 * r18606); float r18615 = pow(r18614, r18604); float r18616 = 21.0f; float r18617 = r18616 / r18606; float r18618 = r18615 / r18617; float r18619 = r18614 * r18611; float r18620 = 5.0f; float r18621 = r18619 / r18620; float r18622 = r18618 + r18621; float r18623 = r18613 + r18622; float r18624 = atan2(1.0, 0.0); float r18625 = sqrt(r18624); float r18626 = r18623 / r18625; float r18627 = fabs(r18626); return r18627; } double f_od(double x) { double r18628 = 2.0; double r18629 = x; double r18630 = fabs(r18629); double r18631 = r18628 * r18630; double r18632 = 3.0; double r18633 = r18628 / r18632; double r18634 = r18630 * r18633; double r18635 = r18630 * r18630; double r18636 = r18634 * r18635; double r18637 = r18631 + r18636; double r18638 = r18630 * (r18630 * r18630); double r18639 = pow(r18638, r18628); double r18640 = 21.0; double r18641 = r18640 / r18630; double r18642 = r18639 / r18641; double r18643 = r18638 * r18635; double r18644 = 5.0; double r18645 = r18643 / r18644; double r18646 = r18642 + r18645; double r18647 = r18637 + r18646; double r18648 = atan2(1.0, 0.0); double r18649 = sqrt(r18648); double r18650 = r18647 / r18649; double r18651 = fabs(r18650); return r18651; } 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 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_im() { 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, "2", 10, MPFR_RNDN); mpfr_init(r18657); mpfr_init(r18658); mpfr_init(r18659); mpfr_init_set_str(r18660, "3", 10, MPFR_RNDN); mpfr_init(r18661); mpfr_init(r18662); mpfr_init(r18663); mpfr_init(r18664); mpfr_init(r18665); mpfr_init_set_str(r18666, "5", 10, MPFR_RNDN); mpfr_init(r18667); mpfr_init(r18668); mpfr_init(r18669); mpfr_init(r18670); mpfr_init(r18671); mpfr_init_set_str(r18672, "21", 10, MPFR_RNDN); mpfr_init(r18673); mpfr_init(r18674); mpfr_init(r18675); mpfr_init(r18676); mpfr_init(r18677); mpfr_init(r18678); mpfr_init(r18679); } double f_im(double x) { ; mpfr_const_pi(r18653, MPFR_RNDN); mpfr_sqrt(r18654, r18653, MPFR_RNDN); mpfr_div(r18655, r18652, r18654, MPFR_RNDN); ; mpfr_set_d(r18657, x, MPFR_RNDN); mpfr_abs(r18658, r18657, MPFR_RNDN); mpfr_mul(r18659, r18656, r18658, MPFR_RNDN); ; mpfr_div(r18661, r18656, r18660, MPFR_RNDN); mpfr_mul(r18662, r18658, r18658, MPFR_RNDN); mpfr_mul(r18663, r18662, r18658, MPFR_RNDN); mpfr_mul(r18664, r18661, r18663, MPFR_RNDN); mpfr_add(r18665, r18659, r18664, MPFR_RNDN); ; mpfr_div(r18667, r18652, r18666, MPFR_RNDN); mpfr_mul(r18668, r18663, r18658, MPFR_RNDN); mpfr_mul(r18669, r18668, r18658, MPFR_RNDN); mpfr_mul(r18670, r18667, r18669, MPFR_RNDN); mpfr_add(r18671, r18665, r18670, MPFR_RNDN); ; mpfr_div(r18673, r18652, r18672, MPFR_RNDN); mpfr_mul(r18674, r18669, r18658, MPFR_RNDN); mpfr_mul(r18675, r18674, r18658, MPFR_RNDN); mpfr_mul(r18676, r18673, r18675, MPFR_RNDN); mpfr_add(r18677, r18671, r18676, MPFR_RNDN); mpfr_mul(r18678, r18655, r18677, MPFR_RNDN); mpfr_abs(r18679, r18678, MPFR_RNDN); return mpfr_get_d(r18679, MPFR_RNDN); } static mpfr_t r18680, r18681, r18682, r18683, r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18680, "2", 10, MPFR_RNDN); mpfr_init(r18681); mpfr_init(r18682); mpfr_init(r18683); mpfr_init_set_str(r18684, "3", 10, MPFR_RNDN); mpfr_init(r18685); mpfr_init(r18686); mpfr_init(r18687); mpfr_init(r18688); mpfr_init(r18689); mpfr_init(r18690); mpfr_init(r18691); mpfr_init_set_str(r18692, "21", 10, MPFR_RNDN); mpfr_init(r18693); mpfr_init(r18694); mpfr_init(r18695); mpfr_init_set_str(r18696, "5", 10, MPFR_RNDN); mpfr_init(r18697); mpfr_init(r18698); mpfr_init(r18699); mpfr_init(r18700); mpfr_init(r18701); mpfr_init(r18702); mpfr_init(r18703); } double f_fm(double x) { ; mpfr_set_d(r18681, x, MPFR_RNDN); mpfr_abs(r18682, r18681, MPFR_RNDN); mpfr_mul(r18683, r18680, r18682, MPFR_RNDN); ; mpfr_div(r18685, r18680, r18684, MPFR_RNDN); mpfr_mul(r18686, r18682, r18685, MPFR_RNDN); mpfr_sqr(r18687, r18682, MPFR_RNDN); mpfr_mul(r18688, r18686, r18687, MPFR_RNDN); mpfr_add(r18689, r18683, r18688, MPFR_RNDN); mpfr_mul(r18690, r18682, r18682, MPFR_RNDN); mpfr_mul(r18690, r18690, r18682, MPFR_RNDN); mpfr_pow(r18691, r18690, r18680, MPFR_RNDN); ; mpfr_div(r18693, r18692, r18682, MPFR_RNDN); mpfr_div(r18694, r18691, r18693, MPFR_RNDN); mpfr_mul(r18695, r18690, r18687, MPFR_RNDN); ; mpfr_div(r18697, r18695, r18696, MPFR_RNDN); mpfr_add(r18698, r18694, r18697, MPFR_RNDN); mpfr_add(r18699, r18689, r18698, MPFR_RNDN); mpfr_const_pi(r18700, MPFR_RNDN); mpfr_sqrt(r18701, r18700, MPFR_RNDN); mpfr_div(r18702, r18699, r18701, MPFR_RNDN); mpfr_abs(r18703, r18702, MPFR_RNDN); return mpfr_get_d(r18703, MPFR_RNDN); } static mpfr_t r18704, r18705, r18706, r18707, r18708, r18709, r18710, r18711, r18712, r18713, r18714, r18715, r18716, r18717, r18718, r18719, r18720, r18721, r18722, r18723, r18724, r18725, r18726, r18727; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18704, "2", 10, MPFR_RNDN); mpfr_init(r18705); mpfr_init(r18706); mpfr_init(r18707); mpfr_init_set_str(r18708, "3", 10, MPFR_RNDN); mpfr_init(r18709); mpfr_init(r18710); mpfr_init(r18711); mpfr_init(r18712); mpfr_init(r18713); mpfr_init(r18714); mpfr_init(r18715); mpfr_init_set_str(r18716, "21", 10, MPFR_RNDN); mpfr_init(r18717); mpfr_init(r18718); mpfr_init(r18719); mpfr_init_set_str(r18720, "5", 10, MPFR_RNDN); mpfr_init(r18721); mpfr_init(r18722); mpfr_init(r18723); mpfr_init(r18724); mpfr_init(r18725); mpfr_init(r18726); mpfr_init(r18727); } double f_dm(double x) { ; mpfr_set_d(r18705, x, MPFR_RNDN); mpfr_abs(r18706, r18705, MPFR_RNDN); mpfr_mul(r18707, r18704, r18706, MPFR_RNDN); ; mpfr_div(r18709, r18704, r18708, MPFR_RNDN); mpfr_mul(r18710, r18706, r18709, MPFR_RNDN); mpfr_sqr(r18711, r18706, MPFR_RNDN); mpfr_mul(r18712, r18710, r18711, MPFR_RNDN); mpfr_add(r18713, r18707, r18712, MPFR_RNDN); mpfr_mul(r18714, r18706, r18706, MPFR_RNDN); mpfr_mul(r18714, r18714, r18706, MPFR_RNDN); mpfr_pow(r18715, r18714, r18704, MPFR_RNDN); ; mpfr_div(r18717, r18716, r18706, MPFR_RNDN); mpfr_div(r18718, r18715, r18717, MPFR_RNDN); mpfr_mul(r18719, r18714, r18711, MPFR_RNDN); ; mpfr_div(r18721, r18719, r18720, MPFR_RNDN); mpfr_add(r18722, r18718, r18721, MPFR_RNDN); mpfr_add(r18723, r18713, r18722, MPFR_RNDN); mpfr_const_pi(r18724, MPFR_RNDN); mpfr_sqrt(r18725, r18724, MPFR_RNDN); mpfr_div(r18726, r18723, r18725, MPFR_RNDN); mpfr_abs(r18727, r18726, MPFR_RNDN); return mpfr_get_d(r18727, MPFR_RNDN); }