#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 = 1.0f; float r18605 = atan2(1.0, 0.0); float r18606 = r18604 / r18605; float r18607 = sqrt(r18606); float r18608 = 0.047619047619047616f; float r18609 = x; float r18610 = fabs(r18609); float r18611 = r18610 * (r18610 * r18610); float r18612 = r18611 * r18611; float r18613 = r18610 * r18612; float r18614 = r18608 * r18613; float r18615 = 0.6666666666666666f; float r18616 = 3.0f; float r18617 = pow(r18610, r18616); float r18618 = r18615 * r18617; float r18619 = 2.0f; float r18620 = r18619 * r18610; float r18621 = 0.2f; float r18622 = r18610 * r18610; float r18623 = r18622 * r18611; float r18624 = r18621 * r18623; float r18625 = r18620 + r18624; float r18626 = r18618 + r18625; float r18627 = r18614 + r18626; float r18628 = r18607 * r18627; float r18629 = fabs(r18628); return r18629; } double f_od(double x) { double r18630 = 1.0; double r18631 = atan2(1.0, 0.0); double r18632 = r18630 / r18631; double r18633 = sqrt(r18632); double r18634 = 0.047619047619047616; double r18635 = x; double r18636 = fabs(r18635); double r18637 = r18636 * (r18636 * r18636); double r18638 = r18637 * r18637; double r18639 = r18636 * r18638; double r18640 = r18634 * r18639; double r18641 = 0.6666666666666666; double r18642 = 3.0; double r18643 = pow(r18636, r18642); double r18644 = r18641 * r18643; double r18645 = 2.0; double r18646 = r18645 * r18636; double r18647 = 0.2; double r18648 = r18636 * r18636; double r18649 = r18648 * r18637; double r18650 = r18647 * r18649; double r18651 = r18646 + r18650; double r18652 = r18644 + r18651; double r18653 = r18640 + r18652; double r18654 = r18633 * r18653; double r18655 = fabs(r18654); return r18655; } 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 r18656, r18657, r18658, r18659, r18660, r18661, r18662, r18663, r18664, r18665, r18666, r18667, r18668, r18669, r18670, r18671, r18672, r18673, r18674, r18675, r18676, r18677, r18678, r18679, r18680, r18681, r18682, r18683; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18656, "1", 10, MPFR_RNDN); mpfr_init(r18657); mpfr_init(r18658); mpfr_init(r18659); mpfr_init_set_str(r18660, "2", 10, MPFR_RNDN); mpfr_init(r18661); mpfr_init(r18662); mpfr_init(r18663); mpfr_init_set_str(r18664, "3", 10, MPFR_RNDN); mpfr_init(r18665); mpfr_init(r18666); mpfr_init(r18667); mpfr_init(r18668); mpfr_init(r18669); mpfr_init_set_str(r18670, "5", 10, MPFR_RNDN); mpfr_init(r18671); mpfr_init(r18672); mpfr_init(r18673); mpfr_init(r18674); mpfr_init(r18675); mpfr_init_set_str(r18676, "21", 10, MPFR_RNDN); mpfr_init(r18677); mpfr_init(r18678); mpfr_init(r18679); mpfr_init(r18680); mpfr_init(r18681); mpfr_init(r18682); mpfr_init(r18683); } double f_im(double x) { ; mpfr_const_pi(r18657, MPFR_RNDN); mpfr_sqrt(r18658, r18657, MPFR_RNDN); mpfr_div(r18659, r18656, r18658, MPFR_RNDN); ; mpfr_set_d(r18661, x, MPFR_RNDN); mpfr_abs(r18662, r18661, MPFR_RNDN); mpfr_mul(r18663, r18660, r18662, MPFR_RNDN); ; mpfr_div(r18665, r18660, r18664, MPFR_RNDN); mpfr_mul(r18666, r18662, r18662, MPFR_RNDN); mpfr_mul(r18667, r18666, r18662, MPFR_RNDN); mpfr_mul(r18668, r18665, r18667, MPFR_RNDN); mpfr_add(r18669, r18663, r18668, MPFR_RNDN); ; mpfr_div(r18671, r18656, r18670, MPFR_RNDN); mpfr_mul(r18672, r18667, r18662, MPFR_RNDN); mpfr_mul(r18673, r18672, r18662, MPFR_RNDN); mpfr_mul(r18674, r18671, r18673, MPFR_RNDN); mpfr_add(r18675, r18669, r18674, MPFR_RNDN); ; mpfr_div(r18677, r18656, r18676, MPFR_RNDN); mpfr_mul(r18678, r18673, r18662, MPFR_RNDN); mpfr_mul(r18679, r18678, r18662, MPFR_RNDN); mpfr_mul(r18680, r18677, r18679, MPFR_RNDN); mpfr_add(r18681, r18675, r18680, MPFR_RNDN); mpfr_mul(r18682, r18659, r18681, MPFR_RNDN); mpfr_abs(r18683, r18682, MPFR_RNDN); return mpfr_get_d(r18683, MPFR_RNDN); } static mpfr_t r18684, r18685, r18686, r18687, r18688, r18689, r18690, r18691, r18692, r18693, r18694, r18695, r18696, r18697, r18698, r18699, r18700, r18701, r18702, r18703, r18704, r18705, r18706, r18707, r18708, r18709; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18684, "1", 10, MPFR_RNDN); mpfr_init(r18685); mpfr_init(r18686); mpfr_init(r18687); mpfr_init_set_str(r18688, "1/21", 10, MPFR_RNDN); mpfr_init(r18689); mpfr_init(r18690); mpfr_init(r18691); mpfr_init(r18692); mpfr_init(r18693); mpfr_init(r18694); mpfr_init_set_str(r18695, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18696, "3", 10, MPFR_RNDN); mpfr_init(r18697); mpfr_init(r18698); mpfr_init_set_str(r18699, "2", 10, MPFR_RNDN); mpfr_init(r18700); mpfr_init_set_str(r18701, "1/5", 10, MPFR_RNDN); mpfr_init(r18702); mpfr_init(r18703); mpfr_init(r18704); mpfr_init(r18705); mpfr_init(r18706); mpfr_init(r18707); mpfr_init(r18708); mpfr_init(r18709); } double f_fm(double x) { ; mpfr_const_pi(r18685, MPFR_RNDN); mpfr_div(r18686, r18684, r18685, MPFR_RNDN); mpfr_sqrt(r18687, r18686, MPFR_RNDN); ; mpfr_set_d(r18689, x, MPFR_RNDN); mpfr_abs(r18690, r18689, MPFR_RNDN); mpfr_mul(r18691, r18690, r18690, MPFR_RNDN); mpfr_mul(r18691, r18691, r18690, MPFR_RNDN); mpfr_sqr(r18692, r18691, MPFR_RNDN); mpfr_mul(r18693, r18690, r18692, MPFR_RNDN); mpfr_mul(r18694, r18688, r18693, MPFR_RNDN); ; ; mpfr_pow(r18697, r18690, r18696, MPFR_RNDN); mpfr_mul(r18698, r18695, r18697, MPFR_RNDN); ; mpfr_mul(r18700, r18699, r18690, MPFR_RNDN); ; mpfr_sqr(r18702, r18690, MPFR_RNDN); mpfr_mul(r18703, r18702, r18691, MPFR_RNDN); mpfr_mul(r18704, r18701, r18703, MPFR_RNDN); mpfr_add(r18705, r18700, r18704, MPFR_RNDN); mpfr_add(r18706, r18698, r18705, MPFR_RNDN); mpfr_add(r18707, r18694, r18706, MPFR_RNDN); mpfr_mul(r18708, r18687, r18707, MPFR_RNDN); mpfr_abs(r18709, r18708, MPFR_RNDN); return mpfr_get_d(r18709, MPFR_RNDN); } static mpfr_t r18710, r18711, r18712, r18713, r18714, r18715, r18716, r18717, r18718, r18719, r18720, r18721, r18722, r18723, r18724, r18725, r18726, r18727, r18728, r18729, r18730, r18731, r18732, r18733, r18734, r18735; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18710, "1", 10, MPFR_RNDN); mpfr_init(r18711); mpfr_init(r18712); mpfr_init(r18713); mpfr_init_set_str(r18714, "1/21", 10, MPFR_RNDN); mpfr_init(r18715); mpfr_init(r18716); mpfr_init(r18717); mpfr_init(r18718); mpfr_init(r18719); mpfr_init(r18720); mpfr_init_set_str(r18721, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18722, "3", 10, MPFR_RNDN); mpfr_init(r18723); mpfr_init(r18724); mpfr_init_set_str(r18725, "2", 10, MPFR_RNDN); mpfr_init(r18726); mpfr_init_set_str(r18727, "1/5", 10, MPFR_RNDN); mpfr_init(r18728); mpfr_init(r18729); mpfr_init(r18730); mpfr_init(r18731); mpfr_init(r18732); mpfr_init(r18733); mpfr_init(r18734); mpfr_init(r18735); } double f_dm(double x) { ; mpfr_const_pi(r18711, MPFR_RNDN); mpfr_div(r18712, r18710, r18711, MPFR_RNDN); mpfr_sqrt(r18713, r18712, MPFR_RNDN); ; mpfr_set_d(r18715, x, MPFR_RNDN); mpfr_abs(r18716, r18715, MPFR_RNDN); mpfr_mul(r18717, r18716, r18716, MPFR_RNDN); mpfr_mul(r18717, r18717, r18716, MPFR_RNDN); mpfr_sqr(r18718, r18717, MPFR_RNDN); mpfr_mul(r18719, r18716, r18718, MPFR_RNDN); mpfr_mul(r18720, r18714, r18719, MPFR_RNDN); ; ; mpfr_pow(r18723, r18716, r18722, MPFR_RNDN); mpfr_mul(r18724, r18721, r18723, MPFR_RNDN); ; mpfr_mul(r18726, r18725, r18716, MPFR_RNDN); ; mpfr_sqr(r18728, r18716, MPFR_RNDN); mpfr_mul(r18729, r18728, r18717, MPFR_RNDN); mpfr_mul(r18730, r18727, r18729, MPFR_RNDN); mpfr_add(r18731, r18726, r18730, MPFR_RNDN); mpfr_add(r18732, r18724, r18731, MPFR_RNDN); mpfr_add(r18733, r18720, r18732, MPFR_RNDN); mpfr_mul(r18734, r18713, r18733, MPFR_RNDN); mpfr_abs(r18735, r18734, MPFR_RNDN); return mpfr_get_d(r18735, MPFR_RNDN); }