#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 r17616 = 1.0f; float r17617 = atan2(1.0, 0.0); float r17618 = sqrt(r17617); float r17619 = r17616 / r17618; float r17620 = 2.0f; float r17621 = x; float r17622 = fabs(r17621); float r17623 = r17620 * r17622; float r17624 = 3.0f; float r17625 = r17620 / r17624; float r17626 = r17622 * r17622; float r17627 = r17626 * r17622; float r17628 = r17625 * r17627; float r17629 = r17623 + r17628; float r17630 = 5.0f; float r17631 = r17616 / r17630; float r17632 = r17627 * r17622; float r17633 = r17632 * r17622; float r17634 = r17631 * r17633; float r17635 = r17629 + r17634; float r17636 = 21.0f; float r17637 = r17616 / r17636; float r17638 = r17633 * r17622; float r17639 = r17638 * r17622; float r17640 = r17637 * r17639; float r17641 = r17635 + r17640; float r17642 = r17619 * r17641; float r17643 = fabs(r17642); return r17643; } double f_id(double x) { double r17644 = 1.0; double r17645 = atan2(1.0, 0.0); double r17646 = sqrt(r17645); double r17647 = r17644 / r17646; double r17648 = 2.0; double r17649 = x; double r17650 = fabs(r17649); double r17651 = r17648 * r17650; double r17652 = 3.0; double r17653 = r17648 / r17652; double r17654 = r17650 * r17650; double r17655 = r17654 * r17650; double r17656 = r17653 * r17655; double r17657 = r17651 + r17656; double r17658 = 5.0; double r17659 = r17644 / r17658; double r17660 = r17655 * r17650; double r17661 = r17660 * r17650; double r17662 = r17659 * r17661; double r17663 = r17657 + r17662; double r17664 = 21.0; double r17665 = r17644 / r17664; double r17666 = r17661 * r17650; double r17667 = r17666 * r17650; double r17668 = r17665 * r17667; double r17669 = r17663 + r17668; double r17670 = r17647 * r17669; double r17671 = fabs(r17670); return r17671; } double f_of(float x) { float r17672 = x; float r17673 = fabs(r17672); float r17674 = 5.0f; float r17675 = r17673 / r17674; float r17676 = r17673 * (r17673 * r17673); float r17677 = r17675 * r17676; float r17678 = 2.0f; float r17679 = 3.0f; float r17680 = r17678 / r17679; float r17681 = r17678 * r17673; float r17682 = fma(r17680, r17676, r17681); float r17683 = fma(r17677, r17673, r17682); float r17684 = r17673 * r17673; float r17685 = r17684 * (r17684 * r17684); float r17686 = 21.0f; float r17687 = r17686 / r17673; float r17688 = r17685 / r17687; float r17689 = r17683 + r17688; float r17690 = atan2(1.0, 0.0); float r17691 = sqrt(r17690); float r17692 = r17689 / r17691; float r17693 = fabs(r17692); return r17693; } double f_od(double x) { double r17694 = x; double r17695 = fabs(r17694); double r17696 = 5.0; double r17697 = r17695 / r17696; double r17698 = r17695 * (r17695 * r17695); double r17699 = r17697 * r17698; double r17700 = 2.0; double r17701 = 3.0; double r17702 = r17700 / r17701; double r17703 = r17700 * r17695; double r17704 = fma(r17702, r17698, r17703); double r17705 = fma(r17699, r17695, r17704); double r17706 = r17695 * r17695; double r17707 = r17706 * (r17706 * r17706); double r17708 = 21.0; double r17709 = r17708 / r17695; double r17710 = r17707 / r17709; double r17711 = r17705 + r17710; double r17712 = atan2(1.0, 0.0); double r17713 = sqrt(r17712); double r17714 = r17711 / r17713; double r17715 = fabs(r17714); return r17715; } 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 r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731, r17732, r17733, r17734, r17735, r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17716, "1", 10, MPFR_RNDN); mpfr_init(r17717); mpfr_init(r17718); mpfr_init(r17719); mpfr_init_set_str(r17720, "2", 10, MPFR_RNDN); mpfr_init(r17721); mpfr_init(r17722); mpfr_init(r17723); mpfr_init_set_str(r17724, "3", 10, MPFR_RNDN); mpfr_init(r17725); mpfr_init(r17726); mpfr_init(r17727); mpfr_init(r17728); mpfr_init(r17729); mpfr_init_set_str(r17730, "5", 10, MPFR_RNDN); mpfr_init(r17731); mpfr_init(r17732); mpfr_init(r17733); mpfr_init(r17734); mpfr_init(r17735); mpfr_init_set_str(r17736, "21", 10, MPFR_RNDN); mpfr_init(r17737); mpfr_init(r17738); mpfr_init(r17739); mpfr_init(r17740); mpfr_init(r17741); mpfr_init(r17742); mpfr_init(r17743); } double f_im(double x) { ; mpfr_const_pi(r17717, MPFR_RNDN); mpfr_sqrt(r17718, r17717, MPFR_RNDN); mpfr_div(r17719, r17716, r17718, MPFR_RNDN); ; mpfr_set_d(r17721, x, MPFR_RNDN); mpfr_abs(r17722, r17721, MPFR_RNDN); mpfr_mul(r17723, r17720, r17722, MPFR_RNDN); ; mpfr_div(r17725, r17720, r17724, MPFR_RNDN); mpfr_mul(r17726, r17722, r17722, MPFR_RNDN); mpfr_mul(r17727, r17726, r17722, MPFR_RNDN); mpfr_mul(r17728, r17725, r17727, MPFR_RNDN); mpfr_add(r17729, r17723, r17728, MPFR_RNDN); ; mpfr_div(r17731, r17716, r17730, MPFR_RNDN); mpfr_mul(r17732, r17727, r17722, MPFR_RNDN); mpfr_mul(r17733, r17732, r17722, MPFR_RNDN); mpfr_mul(r17734, r17731, r17733, MPFR_RNDN); mpfr_add(r17735, r17729, r17734, MPFR_RNDN); ; mpfr_div(r17737, r17716, r17736, MPFR_RNDN); mpfr_mul(r17738, r17733, r17722, MPFR_RNDN); mpfr_mul(r17739, r17738, r17722, MPFR_RNDN); mpfr_mul(r17740, r17737, r17739, MPFR_RNDN); mpfr_add(r17741, r17735, r17740, MPFR_RNDN); mpfr_mul(r17742, r17719, r17741, MPFR_RNDN); mpfr_abs(r17743, r17742, MPFR_RNDN); return mpfr_get_d(r17743, MPFR_RNDN); } static mpfr_t r17744, r17745, r17746, r17747, r17748, r17749, r17750, r17751, r17752, r17753, r17754, r17755, r17756, r17757, r17758, r17759, r17760, r17761, r17762, r17763, r17764, r17765; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17744); mpfr_init(r17745); mpfr_init_set_str(r17746, "5", 10, MPFR_RNDN); mpfr_init(r17747); mpfr_init(r17748); mpfr_init(r17749); mpfr_init_set_str(r17750, "2", 10, MPFR_RNDN); mpfr_init_set_str(r17751, "3", 10, MPFR_RNDN); mpfr_init(r17752); mpfr_init(r17753); mpfr_init(r17754); mpfr_init(r17755); mpfr_init(r17756); mpfr_init(r17757); mpfr_init_set_str(r17758, "21", 10, MPFR_RNDN); mpfr_init(r17759); mpfr_init(r17760); mpfr_init(r17761); mpfr_init(r17762); mpfr_init(r17763); mpfr_init(r17764); mpfr_init(r17765); } double f_fm(double x) { mpfr_set_d(r17744, x, MPFR_RNDN); mpfr_abs(r17745, r17744, MPFR_RNDN); ; mpfr_div(r17747, r17745, r17746, MPFR_RNDN); mpfr_mul(r17748, r17745, r17745, MPFR_RNDN); mpfr_mul(r17748, r17748, r17745, MPFR_RNDN); mpfr_mul(r17749, r17747, r17748, MPFR_RNDN); ; ; mpfr_div(r17752, r17750, r17751, MPFR_RNDN); mpfr_mul(r17753, r17750, r17745, MPFR_RNDN); mpfr_fma(r17754, r17752, r17748, r17753, MPFR_RNDN); mpfr_fma(r17755, r17749, r17745, r17754, MPFR_RNDN); mpfr_sqr(r17756, r17745, MPFR_RNDN); mpfr_mul(r17757, r17756, r17756, MPFR_RNDN); mpfr_mul(r17757, r17757, r17756, MPFR_RNDN); ; mpfr_div(r17759, r17758, r17745, MPFR_RNDN); mpfr_div(r17760, r17757, r17759, MPFR_RNDN); mpfr_add(r17761, r17755, r17760, MPFR_RNDN); mpfr_const_pi(r17762, MPFR_RNDN); mpfr_sqrt(r17763, r17762, MPFR_RNDN); mpfr_div(r17764, r17761, r17763, MPFR_RNDN); mpfr_abs(r17765, r17764, MPFR_RNDN); return mpfr_get_d(r17765, MPFR_RNDN); } static mpfr_t r17766, r17767, r17768, r17769, r17770, r17771, r17772, r17773, r17774, r17775, r17776, r17777, r17778, r17779, r17780, r17781, r17782, r17783, r17784, r17785, r17786, r17787; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17766); mpfr_init(r17767); mpfr_init_set_str(r17768, "5", 10, MPFR_RNDN); mpfr_init(r17769); mpfr_init(r17770); mpfr_init(r17771); mpfr_init_set_str(r17772, "2", 10, MPFR_RNDN); mpfr_init_set_str(r17773, "3", 10, MPFR_RNDN); mpfr_init(r17774); mpfr_init(r17775); mpfr_init(r17776); mpfr_init(r17777); mpfr_init(r17778); mpfr_init(r17779); mpfr_init_set_str(r17780, "21", 10, MPFR_RNDN); mpfr_init(r17781); mpfr_init(r17782); mpfr_init(r17783); mpfr_init(r17784); mpfr_init(r17785); mpfr_init(r17786); mpfr_init(r17787); } double f_dm(double x) { mpfr_set_d(r17766, x, MPFR_RNDN); mpfr_abs(r17767, r17766, MPFR_RNDN); ; mpfr_div(r17769, r17767, r17768, MPFR_RNDN); mpfr_mul(r17770, r17767, r17767, MPFR_RNDN); mpfr_mul(r17770, r17770, r17767, MPFR_RNDN); mpfr_mul(r17771, r17769, r17770, MPFR_RNDN); ; ; mpfr_div(r17774, r17772, r17773, MPFR_RNDN); mpfr_mul(r17775, r17772, r17767, MPFR_RNDN); mpfr_fma(r17776, r17774, r17770, r17775, MPFR_RNDN); mpfr_fma(r17777, r17771, r17767, r17776, MPFR_RNDN); mpfr_sqr(r17778, r17767, MPFR_RNDN); mpfr_mul(r17779, r17778, r17778, MPFR_RNDN); mpfr_mul(r17779, r17779, r17778, MPFR_RNDN); ; mpfr_div(r17781, r17780, r17767, MPFR_RNDN); mpfr_div(r17782, r17779, r17781, MPFR_RNDN); mpfr_add(r17783, r17777, r17782, MPFR_RNDN); mpfr_const_pi(r17784, MPFR_RNDN); mpfr_sqrt(r17785, r17784, MPFR_RNDN); mpfr_div(r17786, r17783, r17785, MPFR_RNDN); mpfr_abs(r17787, r17786, MPFR_RNDN); return mpfr_get_d(r17787, MPFR_RNDN); }