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