#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 r18592 = 1.0f; float r18593 = atan2(1.0, 0.0); float r18594 = sqrt(r18593); float r18595 = r18592 / r18594; float r18596 = 2.0f; float r18597 = x; float r18598 = fabs(r18597); float r18599 = r18596 * r18598; float r18600 = 3.0f; float r18601 = r18596 / r18600; float r18602 = r18598 * r18598; float r18603 = r18602 * r18598; float r18604 = r18601 * r18603; float r18605 = r18599 + r18604; float r18606 = 5.0f; float r18607 = r18592 / r18606; float r18608 = r18603 * r18598; float r18609 = r18608 * r18598; float r18610 = r18607 * r18609; float r18611 = r18605 + r18610; float r18612 = 21.0f; float r18613 = r18592 / r18612; float r18614 = r18609 * r18598; float r18615 = r18614 * r18598; float r18616 = r18613 * r18615; float r18617 = r18611 + r18616; float r18618 = r18595 * r18617; float r18619 = fabs(r18618); return r18619; } double f_id(double x) { double r18620 = 1.0; double r18621 = atan2(1.0, 0.0); double r18622 = sqrt(r18621); double r18623 = r18620 / r18622; double r18624 = 2.0; double r18625 = x; double r18626 = fabs(r18625); double r18627 = r18624 * r18626; double r18628 = 3.0; double r18629 = r18624 / r18628; double r18630 = r18626 * r18626; double r18631 = r18630 * r18626; double r18632 = r18629 * r18631; double r18633 = r18627 + r18632; double r18634 = 5.0; double r18635 = r18620 / r18634; double r18636 = r18631 * r18626; double r18637 = r18636 * r18626; double r18638 = r18635 * r18637; double r18639 = r18633 + r18638; double r18640 = 21.0; double r18641 = r18620 / r18640; double r18642 = r18637 * r18626; double r18643 = r18642 * r18626; double r18644 = r18641 * r18643; double r18645 = r18639 + r18644; double r18646 = r18623 * r18645; double r18647 = fabs(r18646); return r18647; } double f_of(float x) { float r18648 = 2.0f; float r18649 = x; float r18650 = fabs(r18649); float r18651 = r18648 * r18650; float r18652 = 3.0f; float r18653 = r18648 / r18652; float r18654 = r18650 * r18653; float r18655 = r18650 * r18650; float r18656 = r18654 * r18655; float r18657 = r18651 + r18656; float r18658 = r18650 * (r18650 * r18650); float r18659 = r18658 * r18658; float r18660 = 21.0f; float r18661 = r18660 / r18650; float r18662 = r18659 / r18661; float r18663 = r18658 * r18655; float r18664 = 5.0f; float r18665 = r18663 / r18664; float r18666 = r18662 + r18665; float r18667 = r18657 + r18666; float r18668 = atan2(1.0, 0.0); float r18669 = sqrt(r18668); float r18670 = sqrt(r18669); float r18671 = r18670 * r18670; float r18672 = r18667 / r18671; float r18673 = fabs(r18672); return r18673; } double f_od(double x) { double r18674 = 2.0; double r18675 = x; double r18676 = fabs(r18675); double r18677 = r18674 * r18676; double r18678 = 3.0; double r18679 = r18674 / r18678; double r18680 = r18676 * r18679; double r18681 = r18676 * r18676; double r18682 = r18680 * r18681; double r18683 = r18677 + r18682; double r18684 = r18676 * (r18676 * r18676); double r18685 = r18684 * r18684; double r18686 = 21.0; double r18687 = r18686 / r18676; double r18688 = r18685 / r18687; double r18689 = r18684 * r18681; double r18690 = 5.0; double r18691 = r18689 / r18690; double r18692 = r18688 + r18691; double r18693 = r18683 + r18692; double r18694 = atan2(1.0, 0.0); double r18695 = sqrt(r18694); double r18696 = sqrt(r18695); double r18697 = r18696 * r18696; double r18698 = r18693 / r18697; double r18699 = fabs(r18698); return r18699; } 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 r18700, r18701, r18702, r18703, 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_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18700, "1", 10, MPFR_RNDN); mpfr_init(r18701); mpfr_init(r18702); mpfr_init(r18703); 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_set_str(r18714, "5", 10, MPFR_RNDN); mpfr_init(r18715); mpfr_init(r18716); mpfr_init(r18717); mpfr_init(r18718); mpfr_init(r18719); mpfr_init_set_str(r18720, "21", 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_im(double x) { ; mpfr_const_pi(r18701, MPFR_RNDN); mpfr_sqrt(r18702, r18701, MPFR_RNDN); mpfr_div(r18703, r18700, r18702, MPFR_RNDN); ; 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, r18706, MPFR_RNDN); mpfr_mul(r18711, r18710, r18706, MPFR_RNDN); mpfr_mul(r18712, r18709, r18711, MPFR_RNDN); mpfr_add(r18713, r18707, r18712, MPFR_RNDN); ; mpfr_div(r18715, r18700, r18714, MPFR_RNDN); mpfr_mul(r18716, r18711, r18706, MPFR_RNDN); mpfr_mul(r18717, r18716, r18706, MPFR_RNDN); mpfr_mul(r18718, r18715, r18717, MPFR_RNDN); mpfr_add(r18719, r18713, r18718, MPFR_RNDN); ; mpfr_div(r18721, r18700, r18720, MPFR_RNDN); mpfr_mul(r18722, r18717, r18706, MPFR_RNDN); mpfr_mul(r18723, r18722, r18706, MPFR_RNDN); mpfr_mul(r18724, r18721, r18723, MPFR_RNDN); mpfr_add(r18725, r18719, r18724, MPFR_RNDN); mpfr_mul(r18726, r18703, r18725, MPFR_RNDN); mpfr_abs(r18727, r18726, MPFR_RNDN); return mpfr_get_d(r18727, MPFR_RNDN); } static mpfr_t r18728, r18729, r18730, r18731, r18732, r18733, r18734, r18735, r18736, r18737, r18738, r18739, r18740, r18741, r18742, r18743, r18744, r18745, r18746, r18747, r18748, r18749, r18750, r18751, r18752, r18753; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18728, "2", 10, MPFR_RNDN); mpfr_init(r18729); mpfr_init(r18730); mpfr_init(r18731); mpfr_init_set_str(r18732, "3", 10, MPFR_RNDN); mpfr_init(r18733); mpfr_init(r18734); mpfr_init(r18735); mpfr_init(r18736); mpfr_init(r18737); mpfr_init(r18738); mpfr_init(r18739); mpfr_init_set_str(r18740, "21", 10, MPFR_RNDN); mpfr_init(r18741); mpfr_init(r18742); mpfr_init(r18743); mpfr_init_set_str(r18744, "5", 10, MPFR_RNDN); mpfr_init(r18745); mpfr_init(r18746); mpfr_init(r18747); mpfr_init(r18748); mpfr_init(r18749); mpfr_init(r18750); mpfr_init(r18751); mpfr_init(r18752); mpfr_init(r18753); } double f_fm(double x) { ; mpfr_set_d(r18729, x, MPFR_RNDN); mpfr_abs(r18730, r18729, MPFR_RNDN); mpfr_mul(r18731, r18728, r18730, MPFR_RNDN); ; mpfr_div(r18733, r18728, r18732, MPFR_RNDN); mpfr_mul(r18734, r18730, r18733, MPFR_RNDN); mpfr_sqr(r18735, r18730, MPFR_RNDN); mpfr_mul(r18736, r18734, r18735, MPFR_RNDN); mpfr_add(r18737, r18731, r18736, MPFR_RNDN); mpfr_mul(r18738, r18730, r18730, MPFR_RNDN); mpfr_mul(r18738, r18738, r18730, MPFR_RNDN); mpfr_sqr(r18739, r18738, MPFR_RNDN); ; mpfr_div(r18741, r18740, r18730, MPFR_RNDN); mpfr_div(r18742, r18739, r18741, MPFR_RNDN); mpfr_mul(r18743, r18738, r18735, MPFR_RNDN); ; mpfr_div(r18745, r18743, r18744, MPFR_RNDN); mpfr_add(r18746, r18742, r18745, MPFR_RNDN); mpfr_add(r18747, r18737, r18746, MPFR_RNDN); mpfr_const_pi(r18748, MPFR_RNDN); mpfr_sqrt(r18749, r18748, MPFR_RNDN); mpfr_sqrt(r18750, r18749, MPFR_RNDN); mpfr_sqr(r18751, r18750, MPFR_RNDN); mpfr_div(r18752, r18747, r18751, MPFR_RNDN); mpfr_abs(r18753, r18752, MPFR_RNDN); return mpfr_get_d(r18753, MPFR_RNDN); } static mpfr_t r18754, r18755, r18756, r18757, r18758, r18759, r18760, r18761, r18762, r18763, r18764, r18765, r18766, r18767, r18768, r18769, r18770, r18771, r18772, r18773, r18774, r18775, r18776, r18777, r18778, r18779; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18754, "2", 10, MPFR_RNDN); mpfr_init(r18755); mpfr_init(r18756); mpfr_init(r18757); mpfr_init_set_str(r18758, "3", 10, MPFR_RNDN); mpfr_init(r18759); mpfr_init(r18760); mpfr_init(r18761); mpfr_init(r18762); mpfr_init(r18763); mpfr_init(r18764); mpfr_init(r18765); mpfr_init_set_str(r18766, "21", 10, MPFR_RNDN); mpfr_init(r18767); mpfr_init(r18768); mpfr_init(r18769); mpfr_init_set_str(r18770, "5", 10, MPFR_RNDN); mpfr_init(r18771); mpfr_init(r18772); mpfr_init(r18773); mpfr_init(r18774); mpfr_init(r18775); mpfr_init(r18776); mpfr_init(r18777); mpfr_init(r18778); mpfr_init(r18779); } double f_dm(double x) { ; mpfr_set_d(r18755, x, MPFR_RNDN); mpfr_abs(r18756, r18755, MPFR_RNDN); mpfr_mul(r18757, r18754, r18756, MPFR_RNDN); ; mpfr_div(r18759, r18754, r18758, MPFR_RNDN); mpfr_mul(r18760, r18756, r18759, MPFR_RNDN); mpfr_sqr(r18761, r18756, MPFR_RNDN); mpfr_mul(r18762, r18760, r18761, MPFR_RNDN); mpfr_add(r18763, r18757, r18762, MPFR_RNDN); mpfr_mul(r18764, r18756, r18756, MPFR_RNDN); mpfr_mul(r18764, r18764, r18756, MPFR_RNDN); mpfr_sqr(r18765, r18764, MPFR_RNDN); ; mpfr_div(r18767, r18766, r18756, MPFR_RNDN); mpfr_div(r18768, r18765, r18767, MPFR_RNDN); mpfr_mul(r18769, r18764, r18761, MPFR_RNDN); ; mpfr_div(r18771, r18769, r18770, MPFR_RNDN); mpfr_add(r18772, r18768, r18771, MPFR_RNDN); mpfr_add(r18773, r18763, r18772, MPFR_RNDN); mpfr_const_pi(r18774, MPFR_RNDN); mpfr_sqrt(r18775, r18774, MPFR_RNDN); mpfr_sqrt(r18776, r18775, MPFR_RNDN); mpfr_sqr(r18777, r18776, MPFR_RNDN); mpfr_div(r18778, r18773, r18777, MPFR_RNDN); mpfr_abs(r18779, r18778, MPFR_RNDN); return mpfr_get_d(r18779, MPFR_RNDN); }