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