#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 r17576 = 1.0f; float r17577 = atan2(1.0, 0.0); float r17578 = sqrt(r17577); float r17579 = r17576 / r17578; float r17580 = 2.0f; float r17581 = x; float r17582 = fabs(r17581); float r17583 = r17580 * r17582; float r17584 = 3.0f; float r17585 = r17580 / r17584; float r17586 = r17582 * r17582; float r17587 = r17586 * r17582; float r17588 = r17585 * r17587; float r17589 = r17583 + r17588; float r17590 = 5.0f; float r17591 = r17576 / r17590; float r17592 = r17587 * r17582; float r17593 = r17592 * r17582; float r17594 = r17591 * r17593; float r17595 = r17589 + r17594; float r17596 = 21.0f; float r17597 = r17576 / r17596; float r17598 = r17593 * r17582; float r17599 = r17598 * r17582; float r17600 = r17597 * r17599; float r17601 = r17595 + r17600; float r17602 = r17579 * r17601; float r17603 = fabs(r17602); return r17603; } double f_id(double x) { double r17604 = 1.0; double r17605 = atan2(1.0, 0.0); double r17606 = sqrt(r17605); double r17607 = r17604 / r17606; double r17608 = 2.0; double r17609 = x; double r17610 = fabs(r17609); double r17611 = r17608 * r17610; double r17612 = 3.0; double r17613 = r17608 / r17612; double r17614 = r17610 * r17610; double r17615 = r17614 * r17610; double r17616 = r17613 * r17615; double r17617 = r17611 + r17616; double r17618 = 5.0; double r17619 = r17604 / r17618; double r17620 = r17615 * r17610; double r17621 = r17620 * r17610; double r17622 = r17619 * r17621; double r17623 = r17617 + r17622; double r17624 = 21.0; double r17625 = r17604 / r17624; double r17626 = r17621 * r17610; double r17627 = r17626 * r17610; double r17628 = r17625 * r17627; double r17629 = r17623 + r17628; double r17630 = r17607 * r17629; double r17631 = fabs(r17630); return r17631; } double f_of(float x) { float r17632 = 0.20000000298023224f; float r17633 = x; float r17634 = fabs(r17633); float r17635 = r17634 * (r17634 * r17634); float r17636 = r17634 * r17635; float r17637 = r17632 * r17636; float r17638 = 2.0f; float r17639 = 3.0f; float r17640 = r17638 / r17639; float r17641 = r17638 * r17634; float r17642 = fma(r17640, r17635, r17641); float r17643 = fma(r17637, r17634, r17642); float r17644 = r17634 * r17634; float r17645 = r17644 * (r17644 * r17644); float r17646 = 1.0f; float r17647 = pow(r17645, r17646); float r17648 = 21.0f; float r17649 = r17648 / r17634; float r17650 = r17647 / r17649; float r17651 = r17643 + r17650; float r17652 = atan2(1.0, 0.0); float r17653 = sqrt(r17652); float r17654 = r17651 / r17653; float r17655 = fabs(r17654); return r17655; } double f_od(double x) { double r17656 = 0.20000000298023224; double r17657 = x; double r17658 = fabs(r17657); double r17659 = r17658 * (r17658 * r17658); double r17660 = r17658 * r17659; double r17661 = r17656 * r17660; double r17662 = 2.0; double r17663 = 3.0; double r17664 = r17662 / r17663; double r17665 = r17662 * r17658; double r17666 = fma(r17664, r17659, r17665); double r17667 = fma(r17661, r17658, r17666); double r17668 = r17658 * r17658; double r17669 = r17668 * (r17668 * r17668); double r17670 = 1.0; double r17671 = pow(r17669, r17670); double r17672 = 21.0; double r17673 = r17672 / r17658; double r17674 = r17671 / r17673; double r17675 = r17667 + r17674; double r17676 = atan2(1.0, 0.0); double r17677 = sqrt(r17676); double r17678 = r17675 / r17677; double r17679 = fabs(r17678); return r17679; } 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 r17680, r17681, r17682, r17683, r17684, r17685, r17686, r17687, r17688, r17689, r17690, r17691, r17692, r17693, r17694, r17695, r17696, r17697, r17698, r17699, r17700, r17701, r17702, r17703, r17704, r17705, r17706, r17707; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17680, "1", 10, MPFR_RNDN); mpfr_init(r17681); mpfr_init(r17682); mpfr_init(r17683); mpfr_init_set_str(r17684, "2", 10, MPFR_RNDN); mpfr_init(r17685); mpfr_init(r17686); mpfr_init(r17687); mpfr_init_set_str(r17688, "3", 10, MPFR_RNDN); mpfr_init(r17689); mpfr_init(r17690); mpfr_init(r17691); mpfr_init(r17692); mpfr_init(r17693); mpfr_init_set_str(r17694, "5", 10, MPFR_RNDN); mpfr_init(r17695); mpfr_init(r17696); mpfr_init(r17697); mpfr_init(r17698); mpfr_init(r17699); mpfr_init_set_str(r17700, "21", 10, MPFR_RNDN); mpfr_init(r17701); mpfr_init(r17702); mpfr_init(r17703); mpfr_init(r17704); mpfr_init(r17705); mpfr_init(r17706); mpfr_init(r17707); } double f_im(double x) { ; mpfr_const_pi(r17681, MPFR_RNDN); mpfr_sqrt(r17682, r17681, MPFR_RNDN); mpfr_div(r17683, r17680, r17682, MPFR_RNDN); ; mpfr_set_d(r17685, x, MPFR_RNDN); mpfr_abs(r17686, r17685, MPFR_RNDN); mpfr_mul(r17687, r17684, r17686, MPFR_RNDN); ; mpfr_div(r17689, r17684, r17688, MPFR_RNDN); mpfr_mul(r17690, r17686, r17686, MPFR_RNDN); mpfr_mul(r17691, r17690, r17686, MPFR_RNDN); mpfr_mul(r17692, r17689, r17691, MPFR_RNDN); mpfr_add(r17693, r17687, r17692, MPFR_RNDN); ; mpfr_div(r17695, r17680, r17694, MPFR_RNDN); mpfr_mul(r17696, r17691, r17686, MPFR_RNDN); mpfr_mul(r17697, r17696, r17686, MPFR_RNDN); mpfr_mul(r17698, r17695, r17697, MPFR_RNDN); mpfr_add(r17699, r17693, r17698, MPFR_RNDN); ; mpfr_div(r17701, r17680, r17700, MPFR_RNDN); mpfr_mul(r17702, r17697, r17686, MPFR_RNDN); mpfr_mul(r17703, r17702, r17686, MPFR_RNDN); mpfr_mul(r17704, r17701, r17703, MPFR_RNDN); mpfr_add(r17705, r17699, r17704, MPFR_RNDN); mpfr_mul(r17706, r17683, r17705, MPFR_RNDN); mpfr_abs(r17707, r17706, MPFR_RNDN); return mpfr_get_d(r17707, MPFR_RNDN); } static mpfr_t r17708, r17709, r17710, r17711, r17712, r17713, r17714, r17715, r17716, r17717, r17718, r17719, r17720, r17721, r17722, r17723, r17724, r17725, r17726, r17727, r17728, r17729, r17730, r17731; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17708, "1/5", 10, MPFR_RNDN); mpfr_init(r17709); mpfr_init(r17710); mpfr_init(r17711); mpfr_init(r17712); mpfr_init(r17713); mpfr_init_set_str(r17714, "2", 10, MPFR_RNDN); mpfr_init_set_str(r17715, "3", 10, MPFR_RNDN); mpfr_init(r17716); mpfr_init(r17717); mpfr_init(r17718); mpfr_init(r17719); mpfr_init(r17720); mpfr_init(r17721); mpfr_init_set_str(r17722, "1", 10, MPFR_RNDN); mpfr_init(r17723); mpfr_init_set_str(r17724, "21", 10, MPFR_RNDN); mpfr_init(r17725); mpfr_init(r17726); mpfr_init(r17727); mpfr_init(r17728); mpfr_init(r17729); mpfr_init(r17730); mpfr_init(r17731); } double f_fm(double x) { ; mpfr_set_d(r17709, x, MPFR_RNDN); mpfr_abs(r17710, r17709, MPFR_RNDN); mpfr_mul(r17711, r17710, r17710, MPFR_RNDN); mpfr_mul(r17711, r17711, r17710, MPFR_RNDN); mpfr_mul(r17712, r17710, r17711, MPFR_RNDN); mpfr_mul(r17713, r17708, r17712, MPFR_RNDN); ; ; mpfr_div(r17716, r17714, r17715, MPFR_RNDN); mpfr_mul(r17717, r17714, r17710, MPFR_RNDN); mpfr_fma(r17718, r17716, r17711, r17717, MPFR_RNDN); mpfr_fma(r17719, r17713, r17710, r17718, MPFR_RNDN); mpfr_sqr(r17720, r17710, MPFR_RNDN); mpfr_mul(r17721, r17720, r17720, MPFR_RNDN); mpfr_mul(r17721, r17721, r17720, MPFR_RNDN); ; mpfr_pow(r17723, r17721, r17722, MPFR_RNDN); ; mpfr_div(r17725, r17724, r17710, MPFR_RNDN); mpfr_div(r17726, r17723, r17725, MPFR_RNDN); mpfr_add(r17727, r17719, r17726, MPFR_RNDN); mpfr_const_pi(r17728, MPFR_RNDN); mpfr_sqrt(r17729, r17728, MPFR_RNDN); mpfr_div(r17730, r17727, r17729, MPFR_RNDN); mpfr_abs(r17731, r17730, MPFR_RNDN); return mpfr_get_d(r17731, MPFR_RNDN); } static mpfr_t r17732, r17733, r17734, r17735, r17736, r17737, r17738, r17739, r17740, r17741, r17742, r17743, r17744, r17745, r17746, r17747, r17748, r17749, r17750, r17751, r17752, r17753, r17754, r17755; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r17732, "1/5", 10, MPFR_RNDN); mpfr_init(r17733); mpfr_init(r17734); mpfr_init(r17735); mpfr_init(r17736); mpfr_init(r17737); mpfr_init_set_str(r17738, "2", 10, MPFR_RNDN); mpfr_init_set_str(r17739, "3", 10, MPFR_RNDN); mpfr_init(r17740); mpfr_init(r17741); mpfr_init(r17742); mpfr_init(r17743); mpfr_init(r17744); mpfr_init(r17745); mpfr_init_set_str(r17746, "1", 10, MPFR_RNDN); mpfr_init(r17747); mpfr_init_set_str(r17748, "21", 10, MPFR_RNDN); mpfr_init(r17749); mpfr_init(r17750); mpfr_init(r17751); mpfr_init(r17752); mpfr_init(r17753); mpfr_init(r17754); mpfr_init(r17755); } double f_dm(double x) { ; mpfr_set_d(r17733, x, MPFR_RNDN); mpfr_abs(r17734, r17733, MPFR_RNDN); mpfr_mul(r17735, r17734, r17734, MPFR_RNDN); mpfr_mul(r17735, r17735, r17734, MPFR_RNDN); mpfr_mul(r17736, r17734, r17735, MPFR_RNDN); mpfr_mul(r17737, r17732, r17736, MPFR_RNDN); ; ; mpfr_div(r17740, r17738, r17739, MPFR_RNDN); mpfr_mul(r17741, r17738, r17734, MPFR_RNDN); mpfr_fma(r17742, r17740, r17735, r17741, MPFR_RNDN); mpfr_fma(r17743, r17737, r17734, r17742, MPFR_RNDN); mpfr_sqr(r17744, r17734, MPFR_RNDN); mpfr_mul(r17745, r17744, r17744, MPFR_RNDN); mpfr_mul(r17745, r17745, r17744, MPFR_RNDN); ; mpfr_pow(r17747, r17745, r17746, MPFR_RNDN); ; mpfr_div(r17749, r17748, r17734, MPFR_RNDN); mpfr_div(r17750, r17747, r17749, MPFR_RNDN); mpfr_add(r17751, r17743, r17750, MPFR_RNDN); mpfr_const_pi(r17752, MPFR_RNDN); mpfr_sqrt(r17753, r17752, MPFR_RNDN); mpfr_div(r17754, r17751, r17753, MPFR_RNDN); mpfr_abs(r17755, r17754, MPFR_RNDN); return mpfr_get_d(r17755, MPFR_RNDN); }