#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 r22662 = 1; float r22663 = atan2(1.0, 0.0); float r22664 = sqrt(r22663); float r22665 = r22662 / r22664; float r22666 = 2; float r22667 = x; float r22668 = fabs(r22667); float r22669 = r22666 * r22668; float r22670 = 3; float r22671 = r22666 / r22670; float r22672 = r22668 * r22668; float r22673 = r22672 * r22668; float r22674 = r22671 * r22673; float r22675 = r22669 + r22674; float r22676 = 5; float r22677 = r22662 / r22676; float r22678 = r22673 * r22668; float r22679 = r22678 * r22668; float r22680 = r22677 * r22679; float r22681 = r22675 + r22680; float r22682 = 21; float r22683 = r22662 / r22682; float r22684 = r22679 * r22668; float r22685 = r22684 * r22668; float r22686 = r22683 * r22685; float r22687 = r22681 + r22686; float r22688 = r22665 * r22687; float r22689 = fabs(r22688); return r22689; } double f_id(double x) { double r22690 = 1; double r22691 = atan2(1.0, 0.0); double r22692 = sqrt(r22691); double r22693 = r22690 / r22692; double r22694 = 2; double r22695 = x; double r22696 = fabs(r22695); double r22697 = r22694 * r22696; double r22698 = 3; double r22699 = r22694 / r22698; double r22700 = r22696 * r22696; double r22701 = r22700 * r22696; double r22702 = r22699 * r22701; double r22703 = r22697 + r22702; double r22704 = 5; double r22705 = r22690 / r22704; double r22706 = r22701 * r22696; double r22707 = r22706 * r22696; double r22708 = r22705 * r22707; double r22709 = r22703 + r22708; double r22710 = 21; double r22711 = r22690 / r22710; double r22712 = r22707 * r22696; double r22713 = r22712 * r22696; double r22714 = r22711 * r22713; double r22715 = r22709 + r22714; double r22716 = r22693 * r22715; double r22717 = fabs(r22716); return r22717; } double f_of(float x) { float r22718 = x; float r22719 = fabs(r22718); float r22720 = 2; float r22721 = 3; float r22722 = r22720 / r22721; float r22723 = r22719 * r22722; float r22724 = r22719 * r22719; float r22725 = r22723 * r22724; float r22726 = r22719 + r22719; float r22727 = r22725 + r22726; float r22728 = pow(r22719, r22721); float r22729 = r22728 * r22724; float r22730 = 5; float r22731 = r22729 / r22730; float r22732 = r22727 + r22731; float r22733 = r22728 * r22728; float r22734 = r22733 * r22719; float r22735 = 21; float r22736 = r22734 / r22735; float r22737 = r22732 + r22736; float r22738 = atan2(1.0, 0.0); float r22739 = sqrt(r22738); float r22740 = sqrt(r22739); float r22741 = r22740 * r22740; float r22742 = r22737 / r22741; float r22743 = fabs(r22742); return r22743; } double f_od(double x) { double r22744 = x; double r22745 = fabs(r22744); double r22746 = 2; double r22747 = 3; double r22748 = r22746 / r22747; double r22749 = r22745 * r22748; double r22750 = r22745 * r22745; double r22751 = r22749 * r22750; double r22752 = r22745 + r22745; double r22753 = r22751 + r22752; double r22754 = pow(r22745, r22747); double r22755 = r22754 * r22750; double r22756 = 5; double r22757 = r22755 / r22756; double r22758 = r22753 + r22757; double r22759 = r22754 * r22754; double r22760 = r22759 * r22745; double r22761 = 21; double r22762 = r22760 / r22761; double r22763 = r22758 + r22762; double r22764 = atan2(1.0, 0.0); double r22765 = sqrt(r22764); double r22766 = sqrt(r22765); double r22767 = r22766 * r22766; double r22768 = r22763 / r22767; double r22769 = fabs(r22768); return r22769; } 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 r22770, r22771, r22772, r22773, r22774, r22775, r22776, r22777, r22778, r22779, r22780, r22781, r22782, r22783, r22784, r22785, r22786, r22787, r22788, r22789, r22790, r22791, r22792, r22793, r22794, r22795, r22796, r22797; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22770, "1", 10, MPFR_RNDN); mpfr_init(r22771); mpfr_init(r22772); mpfr_init(r22773); mpfr_init_set_str(r22774, "2", 10, MPFR_RNDN); mpfr_init(r22775); mpfr_init(r22776); mpfr_init(r22777); mpfr_init_set_str(r22778, "3", 10, MPFR_RNDN); mpfr_init(r22779); mpfr_init(r22780); mpfr_init(r22781); mpfr_init(r22782); mpfr_init(r22783); mpfr_init_set_str(r22784, "5", 10, MPFR_RNDN); mpfr_init(r22785); mpfr_init(r22786); mpfr_init(r22787); mpfr_init(r22788); mpfr_init(r22789); mpfr_init_set_str(r22790, "21", 10, MPFR_RNDN); mpfr_init(r22791); mpfr_init(r22792); mpfr_init(r22793); mpfr_init(r22794); mpfr_init(r22795); mpfr_init(r22796); mpfr_init(r22797); } double f_im(double x) { ; mpfr_const_pi(r22771, MPFR_RNDN); mpfr_sqrt(r22772, r22771, MPFR_RNDN); mpfr_div(r22773, r22770, r22772, MPFR_RNDN); ; mpfr_set_d(r22775, x, MPFR_RNDN); mpfr_abs(r22776, r22775, MPFR_RNDN); mpfr_mul(r22777, r22774, r22776, MPFR_RNDN); ; mpfr_div(r22779, r22774, r22778, MPFR_RNDN); mpfr_mul(r22780, r22776, r22776, MPFR_RNDN); mpfr_mul(r22781, r22780, r22776, MPFR_RNDN); mpfr_mul(r22782, r22779, r22781, MPFR_RNDN); mpfr_add(r22783, r22777, r22782, MPFR_RNDN); ; mpfr_div(r22785, r22770, r22784, MPFR_RNDN); mpfr_mul(r22786, r22781, r22776, MPFR_RNDN); mpfr_mul(r22787, r22786, r22776, MPFR_RNDN); mpfr_mul(r22788, r22785, r22787, MPFR_RNDN); mpfr_add(r22789, r22783, r22788, MPFR_RNDN); ; mpfr_div(r22791, r22770, r22790, MPFR_RNDN); mpfr_mul(r22792, r22787, r22776, MPFR_RNDN); mpfr_mul(r22793, r22792, r22776, MPFR_RNDN); mpfr_mul(r22794, r22791, r22793, MPFR_RNDN); mpfr_add(r22795, r22789, r22794, MPFR_RNDN); mpfr_mul(r22796, r22773, r22795, MPFR_RNDN); mpfr_abs(r22797, r22796, MPFR_RNDN); return mpfr_get_d(r22797, MPFR_RNDN); } static mpfr_t r22798, r22799, r22800, r22801, r22802, r22803, r22804, r22805, r22806, r22807, r22808, r22809, r22810, r22811, r22812, r22813, r22814, r22815, r22816, r22817, r22818, r22819, r22820, r22821, r22822, r22823; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22798); mpfr_init(r22799); mpfr_init_set_str(r22800, "2", 10, MPFR_RNDN); mpfr_init_set_str(r22801, "3", 10, MPFR_RNDN); mpfr_init(r22802); mpfr_init(r22803); mpfr_init(r22804); mpfr_init(r22805); mpfr_init(r22806); mpfr_init(r22807); mpfr_init(r22808); mpfr_init(r22809); mpfr_init_set_str(r22810, "5", 10, MPFR_RNDN); mpfr_init(r22811); mpfr_init(r22812); mpfr_init(r22813); mpfr_init(r22814); mpfr_init_set_str(r22815, "21", 10, MPFR_RNDN); mpfr_init(r22816); mpfr_init(r22817); mpfr_init(r22818); mpfr_init(r22819); mpfr_init(r22820); mpfr_init(r22821); mpfr_init(r22822); mpfr_init(r22823); } double f_fm(double x) { mpfr_set_d(r22798, x, MPFR_RNDN); mpfr_abs(r22799, r22798, MPFR_RNDN); ; ; mpfr_div(r22802, r22800, r22801, MPFR_RNDN); mpfr_mul(r22803, r22799, r22802, MPFR_RNDN); mpfr_mul(r22804, r22799, r22799, MPFR_RNDN); mpfr_mul(r22805, r22803, r22804, MPFR_RNDN); mpfr_add(r22806, r22799, r22799, MPFR_RNDN); mpfr_add(r22807, r22805, r22806, MPFR_RNDN); mpfr_pow(r22808, r22799, r22801, MPFR_RNDN); mpfr_mul(r22809, r22808, r22804, MPFR_RNDN); ; mpfr_div(r22811, r22809, r22810, MPFR_RNDN); mpfr_add(r22812, r22807, r22811, MPFR_RNDN); mpfr_mul(r22813, r22808, r22808, MPFR_RNDN); mpfr_mul(r22814, r22813, r22799, MPFR_RNDN); ; mpfr_div(r22816, r22814, r22815, MPFR_RNDN); mpfr_add(r22817, r22812, r22816, MPFR_RNDN); mpfr_const_pi(r22818, MPFR_RNDN); mpfr_sqrt(r22819, r22818, MPFR_RNDN); mpfr_sqrt(r22820, r22819, MPFR_RNDN); mpfr_mul(r22821, r22820, r22820, MPFR_RNDN); mpfr_div(r22822, r22817, r22821, MPFR_RNDN); mpfr_abs(r22823, r22822, MPFR_RNDN); return mpfr_get_d(r22823, MPFR_RNDN); } static mpfr_t r22824, r22825, r22826, r22827, r22828, r22829, r22830, r22831, r22832, r22833, r22834, r22835, r22836, r22837, r22838, r22839, r22840, r22841, r22842, r22843, r22844, r22845, r22846, r22847, r22848, r22849; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22824); mpfr_init(r22825); mpfr_init_set_str(r22826, "2", 10, MPFR_RNDN); mpfr_init_set_str(r22827, "3", 10, MPFR_RNDN); mpfr_init(r22828); mpfr_init(r22829); mpfr_init(r22830); mpfr_init(r22831); mpfr_init(r22832); mpfr_init(r22833); mpfr_init(r22834); mpfr_init(r22835); mpfr_init_set_str(r22836, "5", 10, MPFR_RNDN); mpfr_init(r22837); mpfr_init(r22838); mpfr_init(r22839); mpfr_init(r22840); mpfr_init_set_str(r22841, "21", 10, MPFR_RNDN); mpfr_init(r22842); mpfr_init(r22843); mpfr_init(r22844); mpfr_init(r22845); mpfr_init(r22846); mpfr_init(r22847); mpfr_init(r22848); mpfr_init(r22849); } double f_dm(double x) { mpfr_set_d(r22824, x, MPFR_RNDN); mpfr_abs(r22825, r22824, MPFR_RNDN); ; ; mpfr_div(r22828, r22826, r22827, MPFR_RNDN); mpfr_mul(r22829, r22825, r22828, MPFR_RNDN); mpfr_mul(r22830, r22825, r22825, MPFR_RNDN); mpfr_mul(r22831, r22829, r22830, MPFR_RNDN); mpfr_add(r22832, r22825, r22825, MPFR_RNDN); mpfr_add(r22833, r22831, r22832, MPFR_RNDN); mpfr_pow(r22834, r22825, r22827, MPFR_RNDN); mpfr_mul(r22835, r22834, r22830, MPFR_RNDN); ; mpfr_div(r22837, r22835, r22836, MPFR_RNDN); mpfr_add(r22838, r22833, r22837, MPFR_RNDN); mpfr_mul(r22839, r22834, r22834, MPFR_RNDN); mpfr_mul(r22840, r22839, r22825, MPFR_RNDN); ; mpfr_div(r22842, r22840, r22841, MPFR_RNDN); mpfr_add(r22843, r22838, r22842, MPFR_RNDN); mpfr_const_pi(r22844, MPFR_RNDN); mpfr_sqrt(r22845, r22844, MPFR_RNDN); mpfr_sqrt(r22846, r22845, MPFR_RNDN); mpfr_mul(r22847, r22846, r22846, MPFR_RNDN); mpfr_div(r22848, r22843, r22847, MPFR_RNDN); mpfr_abs(r22849, r22848, MPFR_RNDN); return mpfr_get_d(r22849, MPFR_RNDN); }