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