#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 r18664 = 1.0f; float r18665 = atan2(1.0, 0.0); float r18666 = sqrt(r18665); float r18667 = r18664 / r18666; float r18668 = 2.0f; float r18669 = x; float r18670 = fabs(r18669); float r18671 = r18668 * r18670; float r18672 = 3.0f; float r18673 = r18668 / r18672; float r18674 = r18670 * r18670; float r18675 = r18674 * r18670; float r18676 = r18673 * r18675; float r18677 = r18671 + r18676; float r18678 = 5.0f; float r18679 = r18664 / r18678; float r18680 = r18675 * r18670; float r18681 = r18680 * r18670; float r18682 = r18679 * r18681; float r18683 = r18677 + r18682; float r18684 = 21.0f; float r18685 = r18664 / r18684; float r18686 = r18681 * r18670; float r18687 = r18686 * r18670; float r18688 = r18685 * r18687; float r18689 = r18683 + r18688; float r18690 = r18667 * r18689; float r18691 = fabs(r18690); return r18691; } double f_id(double x) { double r18692 = 1.0; double r18693 = atan2(1.0, 0.0); double r18694 = sqrt(r18693); double r18695 = r18692 / r18694; double r18696 = 2.0; double r18697 = x; double r18698 = fabs(r18697); double r18699 = r18696 * r18698; double r18700 = 3.0; double r18701 = r18696 / r18700; double r18702 = r18698 * r18698; double r18703 = r18702 * r18698; double r18704 = r18701 * r18703; double r18705 = r18699 + r18704; double r18706 = 5.0; double r18707 = r18692 / r18706; double r18708 = r18703 * r18698; double r18709 = r18708 * r18698; double r18710 = r18707 * r18709; double r18711 = r18705 + r18710; double r18712 = 21.0; double r18713 = r18692 / r18712; double r18714 = r18709 * r18698; double r18715 = r18714 * r18698; double r18716 = r18713 * r18715; double r18717 = r18711 + r18716; double r18718 = r18695 * r18717; double r18719 = fabs(r18718); return r18719; } double f_of(float x) { float r18720 = 2.0f; float r18721 = x; float r18722 = fabs(r18721); float r18723 = r18720 * r18722; float r18724 = 3.0f; float r18725 = r18720 / r18724; float r18726 = r18722 * r18725; float r18727 = cbrt(r18726); float r18728 = r18727 * (r18727 * r18727); float r18729 = r18722 * r18722; float r18730 = r18728 * r18729; float r18731 = r18723 + r18730; float r18732 = r18722 * (r18722 * r18722); float r18733 = r18732 * r18732; float r18734 = 21.0f; float r18735 = r18734 / r18722; float r18736 = r18733 / r18735; float r18737 = r18732 * r18729; float r18738 = 5.0f; float r18739 = r18737 / r18738; float r18740 = r18736 + r18739; float r18741 = r18731 + r18740; float r18742 = atan2(1.0, 0.0); float r18743 = sqrt(r18742); float r18744 = r18741 / r18743; float r18745 = fabs(r18744); return r18745; } double f_od(double x) { double r18746 = 2.0; double r18747 = x; double r18748 = fabs(r18747); double r18749 = r18746 * r18748; double r18750 = 3.0; double r18751 = r18746 / r18750; double r18752 = r18748 * r18751; double r18753 = cbrt(r18752); double r18754 = r18753 * (r18753 * r18753); double r18755 = r18748 * r18748; double r18756 = r18754 * r18755; double r18757 = r18749 + r18756; double r18758 = r18748 * (r18748 * r18748); double r18759 = r18758 * r18758; double r18760 = 21.0; double r18761 = r18760 / r18748; double r18762 = r18759 / r18761; double r18763 = r18758 * r18755; double r18764 = 5.0; double r18765 = r18763 / r18764; double r18766 = r18762 + r18765; double r18767 = r18757 + r18766; double r18768 = atan2(1.0, 0.0); double r18769 = sqrt(r18768); double r18770 = r18767 / r18769; double r18771 = fabs(r18770); return r18771; } 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 r18772, r18773, r18774, r18775, r18776, r18777, r18778, r18779, r18780, r18781, r18782, r18783, r18784, r18785, r18786, r18787, r18788, r18789, r18790, r18791, r18792, r18793, r18794, r18795, r18796, r18797, r18798, r18799; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18772, "1", 10, MPFR_RNDN); mpfr_init(r18773); mpfr_init(r18774); mpfr_init(r18775); mpfr_init_set_str(r18776, "2", 10, MPFR_RNDN); mpfr_init(r18777); mpfr_init(r18778); mpfr_init(r18779); mpfr_init_set_str(r18780, "3", 10, MPFR_RNDN); mpfr_init(r18781); mpfr_init(r18782); mpfr_init(r18783); mpfr_init(r18784); mpfr_init(r18785); mpfr_init_set_str(r18786, "5", 10, MPFR_RNDN); mpfr_init(r18787); mpfr_init(r18788); mpfr_init(r18789); mpfr_init(r18790); mpfr_init(r18791); mpfr_init_set_str(r18792, "21", 10, MPFR_RNDN); mpfr_init(r18793); mpfr_init(r18794); mpfr_init(r18795); mpfr_init(r18796); mpfr_init(r18797); mpfr_init(r18798); mpfr_init(r18799); } double f_im(double x) { ; mpfr_const_pi(r18773, MPFR_RNDN); mpfr_sqrt(r18774, r18773, MPFR_RNDN); mpfr_div(r18775, r18772, r18774, MPFR_RNDN); ; mpfr_set_d(r18777, x, MPFR_RNDN); mpfr_abs(r18778, r18777, MPFR_RNDN); mpfr_mul(r18779, r18776, r18778, MPFR_RNDN); ; mpfr_div(r18781, r18776, r18780, MPFR_RNDN); mpfr_mul(r18782, r18778, r18778, MPFR_RNDN); mpfr_mul(r18783, r18782, r18778, MPFR_RNDN); mpfr_mul(r18784, r18781, r18783, MPFR_RNDN); mpfr_add(r18785, r18779, r18784, MPFR_RNDN); ; mpfr_div(r18787, r18772, r18786, MPFR_RNDN); mpfr_mul(r18788, r18783, r18778, MPFR_RNDN); mpfr_mul(r18789, r18788, r18778, MPFR_RNDN); mpfr_mul(r18790, r18787, r18789, MPFR_RNDN); mpfr_add(r18791, r18785, r18790, MPFR_RNDN); ; mpfr_div(r18793, r18772, r18792, MPFR_RNDN); mpfr_mul(r18794, r18789, r18778, MPFR_RNDN); mpfr_mul(r18795, r18794, r18778, MPFR_RNDN); mpfr_mul(r18796, r18793, r18795, MPFR_RNDN); mpfr_add(r18797, r18791, r18796, MPFR_RNDN); mpfr_mul(r18798, r18775, r18797, MPFR_RNDN); mpfr_abs(r18799, r18798, MPFR_RNDN); return mpfr_get_d(r18799, MPFR_RNDN); } static mpfr_t r18800, r18801, r18802, r18803, r18804, r18805, r18806, r18807, r18808, r18809, r18810, r18811, r18812, r18813, r18814, r18815, r18816, r18817, r18818, r18819, r18820, r18821, r18822, r18823, r18824, r18825; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18800, "2", 10, MPFR_RNDN); mpfr_init(r18801); mpfr_init(r18802); mpfr_init(r18803); mpfr_init_set_str(r18804, "3", 10, MPFR_RNDN); mpfr_init(r18805); mpfr_init(r18806); mpfr_init(r18807); mpfr_init(r18808); mpfr_init(r18809); mpfr_init(r18810); mpfr_init(r18811); mpfr_init(r18812); mpfr_init(r18813); mpfr_init_set_str(r18814, "21", 10, MPFR_RNDN); mpfr_init(r18815); mpfr_init(r18816); mpfr_init(r18817); mpfr_init_set_str(r18818, "5", 10, MPFR_RNDN); mpfr_init(r18819); mpfr_init(r18820); mpfr_init(r18821); mpfr_init(r18822); mpfr_init(r18823); mpfr_init(r18824); mpfr_init(r18825); } double f_fm(double x) { ; mpfr_set_d(r18801, x, MPFR_RNDN); mpfr_abs(r18802, r18801, MPFR_RNDN); mpfr_mul(r18803, r18800, r18802, MPFR_RNDN); ; mpfr_div(r18805, r18800, r18804, MPFR_RNDN); mpfr_mul(r18806, r18802, r18805, MPFR_RNDN); mpfr_cbrt(r18807, r18806, MPFR_RNDN); mpfr_mul(r18808, r18807, r18807, MPFR_RNDN); mpfr_mul(r18808, r18808, r18807, MPFR_RNDN); mpfr_sqr(r18809, r18802, MPFR_RNDN); mpfr_mul(r18810, r18808, r18809, MPFR_RNDN); mpfr_add(r18811, r18803, r18810, MPFR_RNDN); mpfr_mul(r18812, r18802, r18802, MPFR_RNDN); mpfr_mul(r18812, r18812, r18802, MPFR_RNDN); mpfr_sqr(r18813, r18812, MPFR_RNDN); ; mpfr_div(r18815, r18814, r18802, MPFR_RNDN); mpfr_div(r18816, r18813, r18815, MPFR_RNDN); mpfr_mul(r18817, r18812, r18809, MPFR_RNDN); ; mpfr_div(r18819, r18817, r18818, MPFR_RNDN); mpfr_add(r18820, r18816, r18819, MPFR_RNDN); mpfr_add(r18821, r18811, r18820, MPFR_RNDN); mpfr_const_pi(r18822, MPFR_RNDN); mpfr_sqrt(r18823, r18822, MPFR_RNDN); mpfr_div(r18824, r18821, r18823, MPFR_RNDN); mpfr_abs(r18825, r18824, MPFR_RNDN); return mpfr_get_d(r18825, MPFR_RNDN); } static mpfr_t r18826, r18827, r18828, r18829, r18830, r18831, r18832, r18833, r18834, r18835, r18836, r18837, r18838, r18839, r18840, r18841, r18842, r18843, r18844, r18845, r18846, r18847, r18848, r18849, r18850, r18851; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r18826, "2", 10, MPFR_RNDN); mpfr_init(r18827); mpfr_init(r18828); mpfr_init(r18829); mpfr_init_set_str(r18830, "3", 10, MPFR_RNDN); mpfr_init(r18831); mpfr_init(r18832); mpfr_init(r18833); mpfr_init(r18834); mpfr_init(r18835); mpfr_init(r18836); mpfr_init(r18837); mpfr_init(r18838); mpfr_init(r18839); mpfr_init_set_str(r18840, "21", 10, MPFR_RNDN); mpfr_init(r18841); mpfr_init(r18842); mpfr_init(r18843); mpfr_init_set_str(r18844, "5", 10, MPFR_RNDN); mpfr_init(r18845); mpfr_init(r18846); mpfr_init(r18847); mpfr_init(r18848); mpfr_init(r18849); mpfr_init(r18850); mpfr_init(r18851); } double f_dm(double x) { ; mpfr_set_d(r18827, x, MPFR_RNDN); mpfr_abs(r18828, r18827, MPFR_RNDN); mpfr_mul(r18829, r18826, r18828, MPFR_RNDN); ; mpfr_div(r18831, r18826, r18830, MPFR_RNDN); mpfr_mul(r18832, r18828, r18831, MPFR_RNDN); mpfr_cbrt(r18833, r18832, MPFR_RNDN); mpfr_mul(r18834, r18833, r18833, MPFR_RNDN); mpfr_mul(r18834, r18834, r18833, MPFR_RNDN); mpfr_sqr(r18835, r18828, MPFR_RNDN); mpfr_mul(r18836, r18834, r18835, MPFR_RNDN); mpfr_add(r18837, r18829, r18836, MPFR_RNDN); mpfr_mul(r18838, r18828, r18828, MPFR_RNDN); mpfr_mul(r18838, r18838, r18828, MPFR_RNDN); mpfr_sqr(r18839, r18838, MPFR_RNDN); ; mpfr_div(r18841, r18840, r18828, MPFR_RNDN); mpfr_div(r18842, r18839, r18841, MPFR_RNDN); mpfr_mul(r18843, r18838, r18835, MPFR_RNDN); ; mpfr_div(r18845, r18843, r18844, MPFR_RNDN); mpfr_add(r18846, r18842, r18845, MPFR_RNDN); mpfr_add(r18847, r18837, r18846, MPFR_RNDN); mpfr_const_pi(r18848, MPFR_RNDN); mpfr_sqrt(r18849, r18848, MPFR_RNDN); mpfr_div(r18850, r18847, r18849, MPFR_RNDN); mpfr_abs(r18851, r18850, MPFR_RNDN); return mpfr_get_d(r18851, MPFR_RNDN); }