#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 r23715 = 1; float r23716 = atan2(1.0, 0.0); float r23717 = sqrt(r23716); float r23718 = r23715 / r23717; float r23719 = 2; float r23720 = x; float r23721 = fabs(r23720); float r23722 = r23719 * r23721; float r23723 = 3; float r23724 = r23719 / r23723; float r23725 = r23721 * r23721; float r23726 = r23725 * r23721; float r23727 = r23724 * r23726; float r23728 = r23722 + r23727; float r23729 = 5; float r23730 = r23715 / r23729; float r23731 = r23726 * r23721; float r23732 = r23731 * r23721; float r23733 = r23730 * r23732; float r23734 = r23728 + r23733; float r23735 = 21; float r23736 = r23715 / r23735; float r23737 = r23732 * r23721; float r23738 = r23737 * r23721; float r23739 = r23736 * r23738; float r23740 = r23734 + r23739; float r23741 = r23718 * r23740; float r23742 = fabs(r23741); return r23742; } double f_id(double x) { double r23743 = 1; double r23744 = atan2(1.0, 0.0); double r23745 = sqrt(r23744); double r23746 = r23743 / r23745; double r23747 = 2; double r23748 = x; double r23749 = fabs(r23748); double r23750 = r23747 * r23749; double r23751 = 3; double r23752 = r23747 / r23751; double r23753 = r23749 * r23749; double r23754 = r23753 * r23749; double r23755 = r23752 * r23754; double r23756 = r23750 + r23755; double r23757 = 5; double r23758 = r23743 / r23757; double r23759 = r23754 * r23749; double r23760 = r23759 * r23749; double r23761 = r23758 * r23760; double r23762 = r23756 + r23761; double r23763 = 21; double r23764 = r23743 / r23763; double r23765 = r23760 * r23749; double r23766 = r23765 * r23749; double r23767 = r23764 * r23766; double r23768 = r23762 + r23767; double r23769 = r23746 * r23768; double r23770 = fabs(r23769); return r23770; } double f_of(float x) { float r23771 = 1; float r23772 = atan2(1.0, 0.0); float r23773 = sqrt(r23772); float r23774 = r23771 / r23773; float r23775 = 2; float r23776 = x; float r23777 = fabs(r23776); float r23778 = r23775 * r23777; float r23779 = 3; float r23780 = r23775 / r23779; float r23781 = r23777 * r23777; float r23782 = r23781 * r23777; float r23783 = r23780 * r23782; float r23784 = r23778 + r23783; float r23785 = 5; float r23786 = r23771 / r23785; float r23787 = r23782 * r23777; float r23788 = r23787 * r23777; float r23789 = r23786 * r23788; float r23790 = r23784 + r23789; float r23791 = 21; float r23792 = r23771 / r23791; float r23793 = r23788 * r23777; float r23794 = r23793 * r23777; float r23795 = r23792 * r23794; float r23796 = r23790 + r23795; float r23797 = r23774 * r23796; float r23798 = fabs(r23797); return r23798; } double f_od(double x) { double r23799 = 1; double r23800 = atan2(1.0, 0.0); double r23801 = sqrt(r23800); double r23802 = r23799 / r23801; double r23803 = 2; double r23804 = x; double r23805 = fabs(r23804); double r23806 = r23803 * r23805; double r23807 = 3; double r23808 = r23803 / r23807; double r23809 = r23805 * r23805; double r23810 = r23809 * r23805; double r23811 = r23808 * r23810; double r23812 = r23806 + r23811; double r23813 = 5; double r23814 = r23799 / r23813; double r23815 = r23810 * r23805; double r23816 = r23815 * r23805; double r23817 = r23814 * r23816; double r23818 = r23812 + r23817; double r23819 = 21; double r23820 = r23799 / r23819; double r23821 = r23816 * r23805; double r23822 = r23821 * r23805; double r23823 = r23820 * r23822; double r23824 = r23818 + r23823; double r23825 = r23802 * r23824; double r23826 = fabs(r23825); return r23826; } 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 r23827, r23828, r23829, r23830, r23831, r23832, r23833, r23834, r23835, r23836, r23837, r23838, r23839, r23840, r23841, r23842, r23843, r23844, r23845, r23846, r23847, r23848, r23849, r23850, r23851, r23852, r23853, r23854; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23827, "1", 10, MPFR_RNDN); mpfr_init(r23828); mpfr_init(r23829); mpfr_init(r23830); mpfr_init_set_str(r23831, "2", 10, MPFR_RNDN); mpfr_init(r23832); mpfr_init(r23833); mpfr_init(r23834); mpfr_init_set_str(r23835, "3", 10, MPFR_RNDN); mpfr_init(r23836); mpfr_init(r23837); mpfr_init(r23838); mpfr_init(r23839); mpfr_init(r23840); mpfr_init_set_str(r23841, "5", 10, MPFR_RNDN); mpfr_init(r23842); mpfr_init(r23843); mpfr_init(r23844); mpfr_init(r23845); mpfr_init(r23846); mpfr_init_set_str(r23847, "21", 10, MPFR_RNDN); mpfr_init(r23848); mpfr_init(r23849); mpfr_init(r23850); mpfr_init(r23851); mpfr_init(r23852); mpfr_init(r23853); mpfr_init(r23854); } double f_im(double x) { ; mpfr_const_pi(r23828, MPFR_RNDN); mpfr_sqrt(r23829, r23828, MPFR_RNDN); mpfr_div(r23830, r23827, r23829, MPFR_RNDN); ; mpfr_set_d(r23832, x, MPFR_RNDN); mpfr_abs(r23833, r23832, MPFR_RNDN); mpfr_mul(r23834, r23831, r23833, MPFR_RNDN); ; mpfr_div(r23836, r23831, r23835, MPFR_RNDN); mpfr_mul(r23837, r23833, r23833, MPFR_RNDN); mpfr_mul(r23838, r23837, r23833, MPFR_RNDN); mpfr_mul(r23839, r23836, r23838, MPFR_RNDN); mpfr_add(r23840, r23834, r23839, MPFR_RNDN); ; mpfr_div(r23842, r23827, r23841, MPFR_RNDN); mpfr_mul(r23843, r23838, r23833, MPFR_RNDN); mpfr_mul(r23844, r23843, r23833, MPFR_RNDN); mpfr_mul(r23845, r23842, r23844, MPFR_RNDN); mpfr_add(r23846, r23840, r23845, MPFR_RNDN); ; mpfr_div(r23848, r23827, r23847, MPFR_RNDN); mpfr_mul(r23849, r23844, r23833, MPFR_RNDN); mpfr_mul(r23850, r23849, r23833, MPFR_RNDN); mpfr_mul(r23851, r23848, r23850, MPFR_RNDN); mpfr_add(r23852, r23846, r23851, MPFR_RNDN); mpfr_mul(r23853, r23830, r23852, MPFR_RNDN); mpfr_abs(r23854, r23853, MPFR_RNDN); return mpfr_get_d(r23854, MPFR_RNDN); } static mpfr_t r23855, r23856, r23857, r23858, r23859, r23860, r23861, r23862, r23863, r23864, r23865, r23866, r23867, r23868, r23869, r23870, r23871, r23872, r23873, r23874, r23875, r23876, r23877, r23878, r23879, r23880, r23881, r23882; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23855, "1", 10, MPFR_RNDN); mpfr_init(r23856); mpfr_init(r23857); mpfr_init(r23858); mpfr_init_set_str(r23859, "2", 10, MPFR_RNDN); mpfr_init(r23860); mpfr_init(r23861); mpfr_init(r23862); mpfr_init_set_str(r23863, "3", 10, MPFR_RNDN); mpfr_init(r23864); mpfr_init(r23865); mpfr_init(r23866); mpfr_init(r23867); mpfr_init(r23868); mpfr_init_set_str(r23869, "5", 10, MPFR_RNDN); mpfr_init(r23870); mpfr_init(r23871); mpfr_init(r23872); mpfr_init(r23873); mpfr_init(r23874); mpfr_init_set_str(r23875, "21", 10, MPFR_RNDN); mpfr_init(r23876); mpfr_init(r23877); mpfr_init(r23878); mpfr_init(r23879); mpfr_init(r23880); mpfr_init(r23881); mpfr_init(r23882); } double f_fm(double x) { ; mpfr_const_pi(r23856, MPFR_RNDN); mpfr_sqrt(r23857, r23856, MPFR_RNDN); mpfr_div(r23858, r23855, r23857, MPFR_RNDN); ; mpfr_set_d(r23860, x, MPFR_RNDN); mpfr_abs(r23861, r23860, MPFR_RNDN); mpfr_mul(r23862, r23859, r23861, MPFR_RNDN); ; mpfr_div(r23864, r23859, r23863, MPFR_RNDN); mpfr_mul(r23865, r23861, r23861, MPFR_RNDN); mpfr_mul(r23866, r23865, r23861, MPFR_RNDN); mpfr_mul(r23867, r23864, r23866, MPFR_RNDN); mpfr_add(r23868, r23862, r23867, MPFR_RNDN); ; mpfr_div(r23870, r23855, r23869, MPFR_RNDN); mpfr_mul(r23871, r23866, r23861, MPFR_RNDN); mpfr_mul(r23872, r23871, r23861, MPFR_RNDN); mpfr_mul(r23873, r23870, r23872, MPFR_RNDN); mpfr_add(r23874, r23868, r23873, MPFR_RNDN); ; mpfr_div(r23876, r23855, r23875, MPFR_RNDN); mpfr_mul(r23877, r23872, r23861, MPFR_RNDN); mpfr_mul(r23878, r23877, r23861, MPFR_RNDN); mpfr_mul(r23879, r23876, r23878, MPFR_RNDN); mpfr_add(r23880, r23874, r23879, MPFR_RNDN); mpfr_mul(r23881, r23858, r23880, MPFR_RNDN); mpfr_abs(r23882, r23881, MPFR_RNDN); return mpfr_get_d(r23882, MPFR_RNDN); } static mpfr_t r23883, r23884, r23885, r23886, r23887, r23888, r23889, r23890, r23891, r23892, r23893, r23894, r23895, r23896, r23897, r23898, r23899, r23900, r23901, r23902, r23903, r23904, r23905, r23906, r23907, r23908, r23909, r23910; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23883, "1", 10, MPFR_RNDN); mpfr_init(r23884); mpfr_init(r23885); mpfr_init(r23886); mpfr_init_set_str(r23887, "2", 10, MPFR_RNDN); mpfr_init(r23888); mpfr_init(r23889); mpfr_init(r23890); mpfr_init_set_str(r23891, "3", 10, MPFR_RNDN); mpfr_init(r23892); mpfr_init(r23893); mpfr_init(r23894); mpfr_init(r23895); mpfr_init(r23896); mpfr_init_set_str(r23897, "5", 10, MPFR_RNDN); mpfr_init(r23898); mpfr_init(r23899); mpfr_init(r23900); mpfr_init(r23901); mpfr_init(r23902); mpfr_init_set_str(r23903, "21", 10, MPFR_RNDN); mpfr_init(r23904); mpfr_init(r23905); mpfr_init(r23906); mpfr_init(r23907); mpfr_init(r23908); mpfr_init(r23909); mpfr_init(r23910); } double f_dm(double x) { ; mpfr_const_pi(r23884, MPFR_RNDN); mpfr_sqrt(r23885, r23884, MPFR_RNDN); mpfr_div(r23886, r23883, r23885, MPFR_RNDN); ; mpfr_set_d(r23888, x, MPFR_RNDN); mpfr_abs(r23889, r23888, MPFR_RNDN); mpfr_mul(r23890, r23887, r23889, MPFR_RNDN); ; mpfr_div(r23892, r23887, r23891, MPFR_RNDN); mpfr_mul(r23893, r23889, r23889, MPFR_RNDN); mpfr_mul(r23894, r23893, r23889, MPFR_RNDN); mpfr_mul(r23895, r23892, r23894, MPFR_RNDN); mpfr_add(r23896, r23890, r23895, MPFR_RNDN); ; mpfr_div(r23898, r23883, r23897, MPFR_RNDN); mpfr_mul(r23899, r23894, r23889, MPFR_RNDN); mpfr_mul(r23900, r23899, r23889, MPFR_RNDN); mpfr_mul(r23901, r23898, r23900, MPFR_RNDN); mpfr_add(r23902, r23896, r23901, MPFR_RNDN); ; mpfr_div(r23904, r23883, r23903, MPFR_RNDN); mpfr_mul(r23905, r23900, r23889, MPFR_RNDN); mpfr_mul(r23906, r23905, r23889, MPFR_RNDN); mpfr_mul(r23907, r23904, r23906, MPFR_RNDN); mpfr_add(r23908, r23902, r23907, MPFR_RNDN); mpfr_mul(r23909, r23886, r23908, MPFR_RNDN); mpfr_abs(r23910, r23909, MPFR_RNDN); return mpfr_get_d(r23910, MPFR_RNDN); }