#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 r23642 = 1; float r23643 = atan2(1.0, 0.0); float r23644 = sqrt(r23643); float r23645 = r23642 / r23644; float r23646 = 2; float r23647 = x; float r23648 = fabs(r23647); float r23649 = r23646 * r23648; float r23650 = 3; float r23651 = r23646 / r23650; float r23652 = r23648 * r23648; float r23653 = r23652 * r23648; float r23654 = r23651 * r23653; float r23655 = r23649 + r23654; float r23656 = 5; float r23657 = r23642 / r23656; float r23658 = r23653 * r23648; float r23659 = r23658 * r23648; float r23660 = r23657 * r23659; float r23661 = r23655 + r23660; float r23662 = 21; float r23663 = r23642 / r23662; float r23664 = r23659 * r23648; float r23665 = r23664 * r23648; float r23666 = r23663 * r23665; float r23667 = r23661 + r23666; float r23668 = r23645 * r23667; float r23669 = fabs(r23668); return r23669; } double f_id(double x) { double r23670 = 1; double r23671 = atan2(1.0, 0.0); double r23672 = sqrt(r23671); double r23673 = r23670 / r23672; double r23674 = 2; double r23675 = x; double r23676 = fabs(r23675); double r23677 = r23674 * r23676; double r23678 = 3; double r23679 = r23674 / r23678; double r23680 = r23676 * r23676; double r23681 = r23680 * r23676; double r23682 = r23679 * r23681; double r23683 = r23677 + r23682; double r23684 = 5; double r23685 = r23670 / r23684; double r23686 = r23681 * r23676; double r23687 = r23686 * r23676; double r23688 = r23685 * r23687; double r23689 = r23683 + r23688; double r23690 = 21; double r23691 = r23670 / r23690; double r23692 = r23687 * r23676; double r23693 = r23692 * r23676; double r23694 = r23691 * r23693; double r23695 = r23689 + r23694; double r23696 = r23673 * r23695; double r23697 = fabs(r23696); return r23697; } double f_of(float x) { float r23698 = 1; float r23699 = 5; float r23700 = r23698 / r23699; float r23701 = x; float r23702 = fabs(r23701); float r23703 = r23700 * r23702; float r23704 = r23702 * r23702; float r23705 = r23703 * r23704; float r23706 = 2; float r23707 = 3; float r23708 = r23706 / r23707; float r23709 = r23708 * r23702; float r23710 = r23702 + r23702; float r23711 = fma(r23709, r23704, r23710); float r23712 = fma(r23705, r23704, r23711); float r23713 = atan2(1.0, 0.0); float r23714 = sqrt(r23713); float r23715 = r23698 / r23714; float r23716 = pow(r23702, r23707); float r23717 = r23716 * r23716; float r23718 = 21; float r23719 = r23702 / r23718; float r23720 = sqrt(r23719); float r23721 = r23720 * r23720; float r23722 = r23721 * r23715; float r23723 = r23717 * r23722; float r23724 = fma(r23712, r23715, r23723); float r23725 = fabs(r23724); return r23725; } double f_od(double x) { double r23726 = 1; double r23727 = 5; double r23728 = r23726 / r23727; double r23729 = x; double r23730 = fabs(r23729); double r23731 = r23728 * r23730; double r23732 = r23730 * r23730; double r23733 = r23731 * r23732; double r23734 = 2; double r23735 = 3; double r23736 = r23734 / r23735; double r23737 = r23736 * r23730; double r23738 = r23730 + r23730; double r23739 = fma(r23737, r23732, r23738); double r23740 = fma(r23733, r23732, r23739); double r23741 = atan2(1.0, 0.0); double r23742 = sqrt(r23741); double r23743 = r23726 / r23742; double r23744 = pow(r23730, r23735); double r23745 = r23744 * r23744; double r23746 = 21; double r23747 = r23730 / r23746; double r23748 = sqrt(r23747); double r23749 = r23748 * r23748; double r23750 = r23749 * r23743; double r23751 = r23745 * r23750; double r23752 = fma(r23740, r23743, r23751); double r23753 = fabs(r23752); return r23753; } 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 r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762, r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778, r23779, r23780, r23781; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r23754, "1", 10, MPFR_RNDN); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init_set_str(r23758, "2", 10, MPFR_RNDN); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init_set_str(r23762, "3", 10, MPFR_RNDN); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init(r23767); mpfr_init_set_str(r23768, "5", 10, MPFR_RNDN); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init_set_str(r23774, "21", 10, MPFR_RNDN); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); mpfr_init(r23778); mpfr_init(r23779); mpfr_init(r23780); mpfr_init(r23781); } double f_im(double x) { ; mpfr_const_pi(r23755, MPFR_RNDN); mpfr_sqrt(r23756, r23755, MPFR_RNDN); mpfr_div(r23757, r23754, r23756, MPFR_RNDN); ; mpfr_set_d(r23759, x, MPFR_RNDN); mpfr_abs(r23760, r23759, MPFR_RNDN); mpfr_mul(r23761, r23758, r23760, MPFR_RNDN); ; mpfr_div(r23763, r23758, r23762, MPFR_RNDN); mpfr_mul(r23764, r23760, r23760, MPFR_RNDN); mpfr_mul(r23765, r23764, r23760, MPFR_RNDN); mpfr_mul(r23766, r23763, r23765, MPFR_RNDN); mpfr_add(r23767, r23761, r23766, MPFR_RNDN); ; mpfr_div(r23769, r23754, r23768, MPFR_RNDN); mpfr_mul(r23770, r23765, r23760, MPFR_RNDN); mpfr_mul(r23771, r23770, r23760, MPFR_RNDN); mpfr_mul(r23772, r23769, r23771, MPFR_RNDN); mpfr_add(r23773, r23767, r23772, MPFR_RNDN); ; mpfr_div(r23775, r23754, r23774, MPFR_RNDN); mpfr_mul(r23776, r23771, r23760, MPFR_RNDN); mpfr_mul(r23777, r23776, r23760, MPFR_RNDN); mpfr_mul(r23778, r23775, r23777, MPFR_RNDN); mpfr_add(r23779, r23773, r23778, MPFR_RNDN); mpfr_mul(r23780, r23757, r23779, MPFR_RNDN); mpfr_abs(r23781, r23780, MPFR_RNDN); return mpfr_get_d(r23781, MPFR_RNDN); } static mpfr_t r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807, r23808, r23809; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23782, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23783, "5", 10, MPFR_RNDN); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init_set_str(r23790, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23791, "3", 10, MPFR_RNDN); mpfr_init(r23792); mpfr_init(r23793); mpfr_init(r23794); mpfr_init(r23795); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init(r23801); mpfr_init_set_str(r23802, "21", 10, MPFR_RNDN); mpfr_init(r23803); mpfr_init(r23804); mpfr_init(r23805); mpfr_init(r23806); mpfr_init(r23807); mpfr_init(r23808); mpfr_init(r23809); } double f_fm(double x) { ; ; mpfr_div(r23784, r23782, r23783, MPFR_RNDN); mpfr_set_d(r23785, x, MPFR_RNDN); mpfr_abs(r23786, r23785, MPFR_RNDN); mpfr_mul(r23787, r23784, r23786, MPFR_RNDN); mpfr_mul(r23788, r23786, r23786, MPFR_RNDN); mpfr_mul(r23789, r23787, r23788, MPFR_RNDN); ; ; mpfr_div(r23792, r23790, r23791, MPFR_RNDN); mpfr_mul(r23793, r23792, r23786, MPFR_RNDN); mpfr_add(r23794, r23786, r23786, MPFR_RNDN); mpfr_fma(r23795, r23793, r23788, r23794, MPFR_RNDN); mpfr_fma(r23796, r23789, r23788, r23795, MPFR_RNDN); mpfr_const_pi(r23797, MPFR_RNDN); mpfr_sqrt(r23798, r23797, MPFR_RNDN); mpfr_div(r23799, r23782, r23798, MPFR_RNDN); mpfr_pow(r23800, r23786, r23791, MPFR_RNDN); mpfr_mul(r23801, r23800, r23800, MPFR_RNDN); ; mpfr_div(r23803, r23786, r23802, MPFR_RNDN); mpfr_sqrt(r23804, r23803, MPFR_RNDN); mpfr_mul(r23805, r23804, r23804, MPFR_RNDN); mpfr_mul(r23806, r23805, r23799, MPFR_RNDN); mpfr_mul(r23807, r23801, r23806, MPFR_RNDN); mpfr_fma(r23808, r23796, r23799, r23807, MPFR_RNDN); mpfr_abs(r23809, r23808, MPFR_RNDN); return mpfr_get_d(r23809, MPFR_RNDN); } static mpfr_t r23810, r23811, r23812, r23813, r23814, r23815, r23816, r23817, r23818, r23819, r23820, r23821, r23822, r23823, r23824, r23825, r23826, r23827, r23828, r23829, r23830, r23831, r23832, r23833, r23834, r23835, r23836, r23837; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r23810, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23811, "5", 10, MPFR_RNDN); mpfr_init(r23812); mpfr_init(r23813); mpfr_init(r23814); mpfr_init(r23815); mpfr_init(r23816); mpfr_init(r23817); mpfr_init_set_str(r23818, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23819, "3", 10, MPFR_RNDN); mpfr_init(r23820); mpfr_init(r23821); mpfr_init(r23822); mpfr_init(r23823); mpfr_init(r23824); mpfr_init(r23825); mpfr_init(r23826); mpfr_init(r23827); mpfr_init(r23828); mpfr_init(r23829); mpfr_init_set_str(r23830, "21", 10, MPFR_RNDN); mpfr_init(r23831); mpfr_init(r23832); mpfr_init(r23833); mpfr_init(r23834); mpfr_init(r23835); mpfr_init(r23836); mpfr_init(r23837); } double f_dm(double x) { ; ; mpfr_div(r23812, r23810, r23811, MPFR_RNDN); mpfr_set_d(r23813, x, MPFR_RNDN); mpfr_abs(r23814, r23813, MPFR_RNDN); mpfr_mul(r23815, r23812, r23814, MPFR_RNDN); mpfr_mul(r23816, r23814, r23814, MPFR_RNDN); mpfr_mul(r23817, r23815, r23816, MPFR_RNDN); ; ; mpfr_div(r23820, r23818, r23819, MPFR_RNDN); mpfr_mul(r23821, r23820, r23814, MPFR_RNDN); mpfr_add(r23822, r23814, r23814, MPFR_RNDN); mpfr_fma(r23823, r23821, r23816, r23822, MPFR_RNDN); mpfr_fma(r23824, r23817, r23816, r23823, MPFR_RNDN); mpfr_const_pi(r23825, MPFR_RNDN); mpfr_sqrt(r23826, r23825, MPFR_RNDN); mpfr_div(r23827, r23810, r23826, MPFR_RNDN); mpfr_pow(r23828, r23814, r23819, MPFR_RNDN); mpfr_mul(r23829, r23828, r23828, MPFR_RNDN); ; mpfr_div(r23831, r23814, r23830, MPFR_RNDN); mpfr_sqrt(r23832, r23831, MPFR_RNDN); mpfr_mul(r23833, r23832, r23832, MPFR_RNDN); mpfr_mul(r23834, r23833, r23827, MPFR_RNDN); mpfr_mul(r23835, r23829, r23834, MPFR_RNDN); mpfr_fma(r23836, r23824, r23827, r23835, MPFR_RNDN); mpfr_abs(r23837, r23836, MPFR_RNDN); return mpfr_get_d(r23837, MPFR_RNDN); }