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