#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 = 1.0f; float r18745 = atan2(1.0, 0.0); float r18746 = r18744 / r18745; float r18747 = sqrt(r18746); float r18748 = 0.047619047619047616f; float r18749 = x; float r18750 = fabs(r18749); float r18751 = r18750 * (r18750 * r18750); float r18752 = r18751 * r18751; float r18753 = r18750 * r18752; float r18754 = r18748 * r18753; float r18755 = 0.6666666666666666f; float r18756 = 3.0f; float r18757 = pow(r18750, r18756); float r18758 = r18755 * r18757; float r18759 = 2.0f; float r18760 = r18759 * r18750; float r18761 = 0.2f; float r18762 = r18750 * r18750; float r18763 = r18762 * r18751; float r18764 = r18761 * r18763; float r18765 = r18760 + r18764; float r18766 = r18758 + r18765; float r18767 = r18754 + r18766; float r18768 = r18747 * r18767; float r18769 = fabs(r18768); return r18769; } double f_od(double x) { double r18770 = 1.0; double r18771 = atan2(1.0, 0.0); double r18772 = r18770 / r18771; double r18773 = sqrt(r18772); double r18774 = 0.047619047619047616; double r18775 = x; double r18776 = fabs(r18775); double r18777 = r18776 * (r18776 * r18776); double r18778 = r18777 * r18777; double r18779 = r18776 * r18778; double r18780 = r18774 * r18779; double r18781 = 0.6666666666666666; double r18782 = 3.0; double r18783 = pow(r18776, r18782); double r18784 = r18781 * r18783; double r18785 = 2.0; double r18786 = r18785 * r18776; double r18787 = 0.2; double r18788 = r18776 * r18776; double r18789 = r18788 * r18777; double r18790 = r18787 * r18789; double r18791 = r18786 + r18790; double r18792 = r18784 + r18791; double r18793 = r18780 + r18792; double r18794 = r18773 * 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, "1", 10, MPFR_RNDN); mpfr_init(r18825); mpfr_init(r18826); mpfr_init(r18827); mpfr_init_set_str(r18828, "1/21", 10, MPFR_RNDN); mpfr_init(r18829); mpfr_init(r18830); mpfr_init(r18831); mpfr_init(r18832); mpfr_init(r18833); mpfr_init(r18834); mpfr_init_set_str(r18835, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18836, "3", 10, MPFR_RNDN); mpfr_init(r18837); mpfr_init(r18838); mpfr_init_set_str(r18839, "2", 10, MPFR_RNDN); mpfr_init(r18840); mpfr_init_set_str(r18841, "1/5", 10, MPFR_RNDN); 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_const_pi(r18825, MPFR_RNDN); mpfr_div(r18826, r18824, r18825, MPFR_RNDN); mpfr_sqrt(r18827, r18826, MPFR_RNDN); ; mpfr_set_d(r18829, x, MPFR_RNDN); mpfr_abs(r18830, r18829, MPFR_RNDN); mpfr_mul(r18831, r18830, r18830, MPFR_RNDN); mpfr_mul(r18831, r18831, r18830, MPFR_RNDN); mpfr_sqr(r18832, r18831, MPFR_RNDN); mpfr_mul(r18833, r18830, r18832, MPFR_RNDN); mpfr_mul(r18834, r18828, r18833, MPFR_RNDN); ; ; mpfr_pow(r18837, r18830, r18836, MPFR_RNDN); mpfr_mul(r18838, r18835, r18837, MPFR_RNDN); ; mpfr_mul(r18840, r18839, r18830, MPFR_RNDN); ; mpfr_sqr(r18842, r18830, MPFR_RNDN); mpfr_mul(r18843, r18842, r18831, MPFR_RNDN); mpfr_mul(r18844, r18841, r18843, MPFR_RNDN); mpfr_add(r18845, r18840, r18844, MPFR_RNDN); mpfr_add(r18846, r18838, r18845, MPFR_RNDN); mpfr_add(r18847, r18834, r18846, MPFR_RNDN); mpfr_mul(r18848, r18827, 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, "1", 10, MPFR_RNDN); mpfr_init(r18851); mpfr_init(r18852); mpfr_init(r18853); mpfr_init_set_str(r18854, "1/21", 10, MPFR_RNDN); mpfr_init(r18855); mpfr_init(r18856); mpfr_init(r18857); mpfr_init(r18858); mpfr_init(r18859); mpfr_init(r18860); mpfr_init_set_str(r18861, "2/3", 10, MPFR_RNDN); mpfr_init_set_str(r18862, "3", 10, MPFR_RNDN); mpfr_init(r18863); mpfr_init(r18864); mpfr_init_set_str(r18865, "2", 10, MPFR_RNDN); mpfr_init(r18866); mpfr_init_set_str(r18867, "1/5", 10, MPFR_RNDN); 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_const_pi(r18851, MPFR_RNDN); mpfr_div(r18852, r18850, r18851, MPFR_RNDN); mpfr_sqrt(r18853, r18852, MPFR_RNDN); ; mpfr_set_d(r18855, x, MPFR_RNDN); mpfr_abs(r18856, r18855, MPFR_RNDN); mpfr_mul(r18857, r18856, r18856, MPFR_RNDN); mpfr_mul(r18857, r18857, r18856, MPFR_RNDN); mpfr_sqr(r18858, r18857, MPFR_RNDN); mpfr_mul(r18859, r18856, r18858, MPFR_RNDN); mpfr_mul(r18860, r18854, r18859, MPFR_RNDN); ; ; mpfr_pow(r18863, r18856, r18862, MPFR_RNDN); mpfr_mul(r18864, r18861, r18863, MPFR_RNDN); ; mpfr_mul(r18866, r18865, r18856, MPFR_RNDN); ; mpfr_sqr(r18868, r18856, MPFR_RNDN); mpfr_mul(r18869, r18868, r18857, MPFR_RNDN); mpfr_mul(r18870, r18867, r18869, MPFR_RNDN); mpfr_add(r18871, r18866, r18870, MPFR_RNDN); mpfr_add(r18872, r18864, r18871, MPFR_RNDN); mpfr_add(r18873, r18860, r18872, MPFR_RNDN); mpfr_mul(r18874, r18853, r18873, MPFR_RNDN); mpfr_abs(r18875, r18874, MPFR_RNDN); return mpfr_get_d(r18875, MPFR_RNDN); }