#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 r27636 = 1; float r27637 = atan2(1.0, 0.0); float r27638 = sqrt(r27637); float r27639 = r27636 / r27638; float r27640 = 2; float r27641 = x; float r27642 = fabs(r27641); float r27643 = r27640 * r27642; float r27644 = 3; float r27645 = r27640 / r27644; float r27646 = r27642 * r27642; float r27647 = r27646 * r27642; float r27648 = r27645 * r27647; float r27649 = r27643 + r27648; float r27650 = 5; float r27651 = r27636 / r27650; float r27652 = r27647 * r27642; float r27653 = r27652 * r27642; float r27654 = r27651 * r27653; float r27655 = r27649 + r27654; float r27656 = 21; float r27657 = r27636 / r27656; float r27658 = r27653 * r27642; float r27659 = r27658 * r27642; float r27660 = r27657 * r27659; float r27661 = r27655 + r27660; float r27662 = r27639 * r27661; float r27663 = fabs(r27662); return r27663; } double f_id(double x) { double r27664 = 1; double r27665 = atan2(1.0, 0.0); double r27666 = sqrt(r27665); double r27667 = r27664 / r27666; double r27668 = 2; double r27669 = x; double r27670 = fabs(r27669); double r27671 = r27668 * r27670; double r27672 = 3; double r27673 = r27668 / r27672; double r27674 = r27670 * r27670; double r27675 = r27674 * r27670; double r27676 = r27673 * r27675; double r27677 = r27671 + r27676; double r27678 = 5; double r27679 = r27664 / r27678; double r27680 = r27675 * r27670; double r27681 = r27680 * r27670; double r27682 = r27679 * r27681; double r27683 = r27677 + r27682; double r27684 = 21; double r27685 = r27664 / r27684; double r27686 = r27681 * r27670; double r27687 = r27686 * r27670; double r27688 = r27685 * r27687; double r27689 = r27683 + r27688; double r27690 = r27667 * r27689; double r27691 = fabs(r27690); return r27691; } double f_of(float x) { float r27692 = 1; float r27693 = atan2(1.0, 0.0); float r27694 = sqrt(r27693); float r27695 = r27692 / r27694; float r27696 = 2; float r27697 = x; float r27698 = fabs(r27697); float r27699 = r27696 * r27698; float r27700 = 3; float r27701 = r27696 / r27700; float r27702 = r27698 * r27698; float r27703 = r27702 * r27698; float r27704 = r27701 * r27703; float r27705 = r27699 + r27704; float r27706 = 5; float r27707 = r27692 / r27706; float r27708 = r27703 * r27698; float r27709 = r27708 * r27698; float r27710 = r27707 * r27709; float r27711 = r27705 + r27710; float r27712 = 21; float r27713 = r27698 / r27712; float r27714 = r27700 + r27700; float r27715 = pow(r27698, r27714); float r27716 = r27713 * r27715; float r27717 = pow(r27716, r27692); float r27718 = r27711 + r27717; float r27719 = r27695 * r27718; float r27720 = fabs(r27719); return r27720; } double f_od(double x) { double r27721 = 1; double r27722 = atan2(1.0, 0.0); double r27723 = sqrt(r27722); double r27724 = r27721 / r27723; double r27725 = 2; double r27726 = x; double r27727 = fabs(r27726); double r27728 = r27725 * r27727; double r27729 = 3; double r27730 = r27725 / r27729; double r27731 = r27727 * r27727; double r27732 = r27731 * r27727; double r27733 = r27730 * r27732; double r27734 = r27728 + r27733; double r27735 = 5; double r27736 = r27721 / r27735; double r27737 = r27732 * r27727; double r27738 = r27737 * r27727; double r27739 = r27736 * r27738; double r27740 = r27734 + r27739; double r27741 = 21; double r27742 = r27727 / r27741; double r27743 = r27729 + r27729; double r27744 = pow(r27727, r27743); double r27745 = r27742 * r27744; double r27746 = pow(r27745, r27721); double r27747 = r27740 + r27746; double r27748 = r27724 * r27747; double r27749 = fabs(r27748); return r27749; } 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 r27750, r27751, r27752, r27753, r27754, r27755, r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r27750, "1", 10, MPFR_RNDN); mpfr_init(r27751); mpfr_init(r27752); mpfr_init(r27753); mpfr_init_set_str(r27754, "2", 10, MPFR_RNDN); mpfr_init(r27755); mpfr_init(r27756); mpfr_init(r27757); mpfr_init_set_str(r27758, "3", 10, MPFR_RNDN); mpfr_init(r27759); mpfr_init(r27760); mpfr_init(r27761); mpfr_init(r27762); mpfr_init(r27763); mpfr_init_set_str(r27764, "5", 10, MPFR_RNDN); mpfr_init(r27765); mpfr_init(r27766); mpfr_init(r27767); mpfr_init(r27768); mpfr_init(r27769); mpfr_init_set_str(r27770, "21", 10, MPFR_RNDN); mpfr_init(r27771); mpfr_init(r27772); mpfr_init(r27773); mpfr_init(r27774); mpfr_init(r27775); mpfr_init(r27776); mpfr_init(r27777); } double f_im(double x) { ; mpfr_const_pi(r27751, MPFR_RNDN); mpfr_sqrt(r27752, r27751, MPFR_RNDN); mpfr_div(r27753, r27750, r27752, MPFR_RNDN); ; mpfr_set_d(r27755, x, MPFR_RNDN); mpfr_abs(r27756, r27755, MPFR_RNDN); mpfr_mul(r27757, r27754, r27756, MPFR_RNDN); ; mpfr_div(r27759, r27754, r27758, MPFR_RNDN); mpfr_mul(r27760, r27756, r27756, MPFR_RNDN); mpfr_mul(r27761, r27760, r27756, MPFR_RNDN); mpfr_mul(r27762, r27759, r27761, MPFR_RNDN); mpfr_add(r27763, r27757, r27762, MPFR_RNDN); ; mpfr_div(r27765, r27750, r27764, MPFR_RNDN); mpfr_mul(r27766, r27761, r27756, MPFR_RNDN); mpfr_mul(r27767, r27766, r27756, MPFR_RNDN); mpfr_mul(r27768, r27765, r27767, MPFR_RNDN); mpfr_add(r27769, r27763, r27768, MPFR_RNDN); ; mpfr_div(r27771, r27750, r27770, MPFR_RNDN); mpfr_mul(r27772, r27767, r27756, MPFR_RNDN); mpfr_mul(r27773, r27772, r27756, MPFR_RNDN); mpfr_mul(r27774, r27771, r27773, MPFR_RNDN); mpfr_add(r27775, r27769, r27774, MPFR_RNDN); mpfr_mul(r27776, r27753, r27775, MPFR_RNDN); mpfr_abs(r27777, r27776, MPFR_RNDN); return mpfr_get_d(r27777, MPFR_RNDN); } static mpfr_t r27778, r27779, r27780, r27781, r27782, r27783, r27784, r27785, r27786, r27787, r27788, r27789, r27790, r27791, r27792, r27793, r27794, r27795, r27796, r27797, r27798, r27799, r27800, r27801, r27802, r27803, r27804, r27805, r27806; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27778, "1", 10, MPFR_RNDN); mpfr_init(r27779); mpfr_init(r27780); mpfr_init(r27781); mpfr_init_set_str(r27782, "2", 10, MPFR_RNDN); mpfr_init(r27783); mpfr_init(r27784); mpfr_init(r27785); mpfr_init_set_str(r27786, "3", 10, MPFR_RNDN); mpfr_init(r27787); mpfr_init(r27788); mpfr_init(r27789); mpfr_init(r27790); mpfr_init(r27791); mpfr_init_set_str(r27792, "5", 10, MPFR_RNDN); mpfr_init(r27793); mpfr_init(r27794); mpfr_init(r27795); mpfr_init(r27796); mpfr_init(r27797); mpfr_init_set_str(r27798, "21", 10, MPFR_RNDN); mpfr_init(r27799); mpfr_init(r27800); mpfr_init(r27801); mpfr_init(r27802); mpfr_init(r27803); mpfr_init(r27804); mpfr_init(r27805); mpfr_init(r27806); } double f_fm(double x) { ; mpfr_const_pi(r27779, MPFR_RNDN); mpfr_sqrt(r27780, r27779, MPFR_RNDN); mpfr_div(r27781, r27778, r27780, MPFR_RNDN); ; mpfr_set_d(r27783, x, MPFR_RNDN); mpfr_abs(r27784, r27783, MPFR_RNDN); mpfr_mul(r27785, r27782, r27784, MPFR_RNDN); ; mpfr_div(r27787, r27782, r27786, MPFR_RNDN); mpfr_mul(r27788, r27784, r27784, MPFR_RNDN); mpfr_mul(r27789, r27788, r27784, MPFR_RNDN); mpfr_mul(r27790, r27787, r27789, MPFR_RNDN); mpfr_add(r27791, r27785, r27790, MPFR_RNDN); ; mpfr_div(r27793, r27778, r27792, MPFR_RNDN); mpfr_mul(r27794, r27789, r27784, MPFR_RNDN); mpfr_mul(r27795, r27794, r27784, MPFR_RNDN); mpfr_mul(r27796, r27793, r27795, MPFR_RNDN); mpfr_add(r27797, r27791, r27796, MPFR_RNDN); ; mpfr_div(r27799, r27784, r27798, MPFR_RNDN); mpfr_add(r27800, r27786, r27786, MPFR_RNDN); mpfr_pow(r27801, r27784, r27800, MPFR_RNDN); mpfr_mul(r27802, r27799, r27801, MPFR_RNDN); mpfr_pow(r27803, r27802, r27778, MPFR_RNDN); mpfr_add(r27804, r27797, r27803, MPFR_RNDN); mpfr_mul(r27805, r27781, r27804, MPFR_RNDN); mpfr_abs(r27806, r27805, MPFR_RNDN); return mpfr_get_d(r27806, MPFR_RNDN); } static mpfr_t r27807, r27808, r27809, r27810, r27811, r27812, r27813, r27814, r27815, r27816, r27817, r27818, r27819, r27820, r27821, r27822, r27823, r27824, r27825, r27826, r27827, r27828, r27829, r27830, r27831, r27832, r27833, r27834, r27835; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r27807, "1", 10, MPFR_RNDN); mpfr_init(r27808); mpfr_init(r27809); mpfr_init(r27810); mpfr_init_set_str(r27811, "2", 10, MPFR_RNDN); mpfr_init(r27812); mpfr_init(r27813); mpfr_init(r27814); mpfr_init_set_str(r27815, "3", 10, MPFR_RNDN); mpfr_init(r27816); mpfr_init(r27817); mpfr_init(r27818); mpfr_init(r27819); mpfr_init(r27820); mpfr_init_set_str(r27821, "5", 10, MPFR_RNDN); mpfr_init(r27822); mpfr_init(r27823); mpfr_init(r27824); mpfr_init(r27825); mpfr_init(r27826); mpfr_init_set_str(r27827, "21", 10, MPFR_RNDN); mpfr_init(r27828); mpfr_init(r27829); mpfr_init(r27830); mpfr_init(r27831); mpfr_init(r27832); mpfr_init(r27833); mpfr_init(r27834); mpfr_init(r27835); } double f_dm(double x) { ; mpfr_const_pi(r27808, MPFR_RNDN); mpfr_sqrt(r27809, r27808, MPFR_RNDN); mpfr_div(r27810, r27807, r27809, MPFR_RNDN); ; mpfr_set_d(r27812, x, MPFR_RNDN); mpfr_abs(r27813, r27812, MPFR_RNDN); mpfr_mul(r27814, r27811, r27813, MPFR_RNDN); ; mpfr_div(r27816, r27811, r27815, MPFR_RNDN); mpfr_mul(r27817, r27813, r27813, MPFR_RNDN); mpfr_mul(r27818, r27817, r27813, MPFR_RNDN); mpfr_mul(r27819, r27816, r27818, MPFR_RNDN); mpfr_add(r27820, r27814, r27819, MPFR_RNDN); ; mpfr_div(r27822, r27807, r27821, MPFR_RNDN); mpfr_mul(r27823, r27818, r27813, MPFR_RNDN); mpfr_mul(r27824, r27823, r27813, MPFR_RNDN); mpfr_mul(r27825, r27822, r27824, MPFR_RNDN); mpfr_add(r27826, r27820, r27825, MPFR_RNDN); ; mpfr_div(r27828, r27813, r27827, MPFR_RNDN); mpfr_add(r27829, r27815, r27815, MPFR_RNDN); mpfr_pow(r27830, r27813, r27829, MPFR_RNDN); mpfr_mul(r27831, r27828, r27830, MPFR_RNDN); mpfr_pow(r27832, r27831, r27807, MPFR_RNDN); mpfr_add(r27833, r27826, r27832, MPFR_RNDN); mpfr_mul(r27834, r27810, r27833, MPFR_RNDN); mpfr_abs(r27835, r27834, MPFR_RNDN); return mpfr_get_d(r27835, MPFR_RNDN); }