#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 r26558 = 1; float r26559 = atan2(1.0, 0.0); float r26560 = sqrt(r26559); float r26561 = r26558 / r26560; float r26562 = 2; float r26563 = x; float r26564 = fabs(r26563); float r26565 = r26562 * r26564; float r26566 = 3; float r26567 = r26562 / r26566; float r26568 = r26564 * r26564; float r26569 = r26568 * r26564; float r26570 = r26567 * r26569; float r26571 = r26565 + r26570; float r26572 = 5; float r26573 = r26558 / r26572; float r26574 = r26569 * r26564; float r26575 = r26574 * r26564; float r26576 = r26573 * r26575; float r26577 = r26571 + r26576; float r26578 = 21; float r26579 = r26558 / r26578; float r26580 = r26575 * r26564; float r26581 = r26580 * r26564; float r26582 = r26579 * r26581; float r26583 = r26577 + r26582; float r26584 = r26561 * r26583; float r26585 = fabs(r26584); return r26585; } double f_id(double x) { double r26586 = 1; double r26587 = atan2(1.0, 0.0); double r26588 = sqrt(r26587); double r26589 = r26586 / r26588; double r26590 = 2; double r26591 = x; double r26592 = fabs(r26591); double r26593 = r26590 * r26592; double r26594 = 3; double r26595 = r26590 / r26594; double r26596 = r26592 * r26592; double r26597 = r26596 * r26592; double r26598 = r26595 * r26597; double r26599 = r26593 + r26598; double r26600 = 5; double r26601 = r26586 / r26600; double r26602 = r26597 * r26592; double r26603 = r26602 * r26592; double r26604 = r26601 * r26603; double r26605 = r26599 + r26604; double r26606 = 21; double r26607 = r26586 / r26606; double r26608 = r26603 * r26592; double r26609 = r26608 * r26592; double r26610 = r26607 * r26609; double r26611 = r26605 + r26610; double r26612 = r26589 * r26611; double r26613 = fabs(r26612); return r26613; } double f_of(float x) { float r26614 = 1; float r26615 = atan2(1.0, 0.0); float r26616 = sqrt(r26615); float r26617 = r26614 / r26616; float r26618 = 2; float r26619 = x; float r26620 = fabs(r26619); float r26621 = r26618 * r26620; float r26622 = 3; float r26623 = r26618 / r26622; float r26624 = r26620 * r26620; float r26625 = r26624 * r26620; float r26626 = r26623 * r26625; float r26627 = r26621 + r26626; float r26628 = 5; float r26629 = r26614 / r26628; float r26630 = r26625 * r26620; float r26631 = r26630 * r26620; float r26632 = r26629 * r26631; float r26633 = r26627 + r26632; float r26634 = 21; float r26635 = r26614 / r26634; float r26636 = r26631 * r26620; float r26637 = r26636 * r26620; float r26638 = r26635 * r26637; float r26639 = r26633 + r26638; float r26640 = r26617 * r26639; float r26641 = fabs(r26640); return r26641; } double f_od(double x) { double r26642 = 1; double r26643 = atan2(1.0, 0.0); double r26644 = sqrt(r26643); double r26645 = r26642 / r26644; double r26646 = 2; double r26647 = x; double r26648 = fabs(r26647); double r26649 = r26646 * r26648; double r26650 = 3; double r26651 = r26646 / r26650; double r26652 = r26648 * r26648; double r26653 = r26652 * r26648; double r26654 = r26651 * r26653; double r26655 = r26649 + r26654; double r26656 = 5; double r26657 = r26642 / r26656; double r26658 = r26653 * r26648; double r26659 = r26658 * r26648; double r26660 = r26657 * r26659; double r26661 = r26655 + r26660; double r26662 = 21; double r26663 = r26642 / r26662; double r26664 = r26659 * r26648; double r26665 = r26664 * r26648; double r26666 = r26663 * r26665; double r26667 = r26661 + r26666; double r26668 = r26645 * r26667; double r26669 = fabs(r26668); return r26669; } 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 r26670, r26671, r26672, r26673, r26674, r26675, r26676, r26677, r26678, r26679, r26680, r26681, r26682, r26683, r26684, r26685, r26686, r26687, r26688, r26689, r26690, r26691, r26692, r26693, r26694, r26695, r26696, r26697; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r26670, "1", 10, MPFR_RNDN); mpfr_init(r26671); mpfr_init(r26672); mpfr_init(r26673); mpfr_init_set_str(r26674, "2", 10, MPFR_RNDN); mpfr_init(r26675); mpfr_init(r26676); mpfr_init(r26677); mpfr_init_set_str(r26678, "3", 10, MPFR_RNDN); mpfr_init(r26679); mpfr_init(r26680); mpfr_init(r26681); mpfr_init(r26682); mpfr_init(r26683); mpfr_init_set_str(r26684, "5", 10, MPFR_RNDN); mpfr_init(r26685); mpfr_init(r26686); mpfr_init(r26687); mpfr_init(r26688); mpfr_init(r26689); mpfr_init_set_str(r26690, "21", 10, MPFR_RNDN); mpfr_init(r26691); mpfr_init(r26692); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init(r26696); mpfr_init(r26697); } double f_im(double x) { ; mpfr_const_pi(r26671, MPFR_RNDN); mpfr_sqrt(r26672, r26671, MPFR_RNDN); mpfr_div(r26673, r26670, r26672, MPFR_RNDN); ; mpfr_set_d(r26675, x, MPFR_RNDN); mpfr_abs(r26676, r26675, MPFR_RNDN); mpfr_mul(r26677, r26674, r26676, MPFR_RNDN); ; mpfr_div(r26679, r26674, r26678, MPFR_RNDN); mpfr_mul(r26680, r26676, r26676, MPFR_RNDN); mpfr_mul(r26681, r26680, r26676, MPFR_RNDN); mpfr_mul(r26682, r26679, r26681, MPFR_RNDN); mpfr_add(r26683, r26677, r26682, MPFR_RNDN); ; mpfr_div(r26685, r26670, r26684, MPFR_RNDN); mpfr_mul(r26686, r26681, r26676, MPFR_RNDN); mpfr_mul(r26687, r26686, r26676, MPFR_RNDN); mpfr_mul(r26688, r26685, r26687, MPFR_RNDN); mpfr_add(r26689, r26683, r26688, MPFR_RNDN); ; mpfr_div(r26691, r26670, r26690, MPFR_RNDN); mpfr_mul(r26692, r26687, r26676, MPFR_RNDN); mpfr_mul(r26693, r26692, r26676, MPFR_RNDN); mpfr_mul(r26694, r26691, r26693, MPFR_RNDN); mpfr_add(r26695, r26689, r26694, MPFR_RNDN); mpfr_mul(r26696, r26673, r26695, MPFR_RNDN); mpfr_abs(r26697, r26696, MPFR_RNDN); return mpfr_get_d(r26697, MPFR_RNDN); } static mpfr_t r26698, r26699, r26700, r26701, r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26698, "1", 10, MPFR_RNDN); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); mpfr_init_set_str(r26702, "2", 10, MPFR_RNDN); mpfr_init(r26703); mpfr_init(r26704); mpfr_init(r26705); mpfr_init_set_str(r26706, "3", 10, MPFR_RNDN); mpfr_init(r26707); mpfr_init(r26708); mpfr_init(r26709); mpfr_init(r26710); mpfr_init(r26711); mpfr_init_set_str(r26712, "5", 10, MPFR_RNDN); mpfr_init(r26713); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init_set_str(r26718, "21", 10, MPFR_RNDN); mpfr_init(r26719); mpfr_init(r26720); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); mpfr_init(r26725); } double f_fm(double x) { ; mpfr_const_pi(r26699, MPFR_RNDN); mpfr_sqrt(r26700, r26699, MPFR_RNDN); mpfr_div(r26701, r26698, r26700, MPFR_RNDN); ; mpfr_set_d(r26703, x, MPFR_RNDN); mpfr_abs(r26704, r26703, MPFR_RNDN); mpfr_mul(r26705, r26702, r26704, MPFR_RNDN); ; mpfr_div(r26707, r26702, r26706, MPFR_RNDN); mpfr_mul(r26708, r26704, r26704, MPFR_RNDN); mpfr_mul(r26709, r26708, r26704, MPFR_RNDN); mpfr_mul(r26710, r26707, r26709, MPFR_RNDN); mpfr_add(r26711, r26705, r26710, MPFR_RNDN); ; mpfr_div(r26713, r26698, r26712, MPFR_RNDN); mpfr_mul(r26714, r26709, r26704, MPFR_RNDN); mpfr_mul(r26715, r26714, r26704, MPFR_RNDN); mpfr_mul(r26716, r26713, r26715, MPFR_RNDN); mpfr_add(r26717, r26711, r26716, MPFR_RNDN); ; mpfr_div(r26719, r26698, r26718, MPFR_RNDN); mpfr_mul(r26720, r26715, r26704, MPFR_RNDN); mpfr_mul(r26721, r26720, r26704, MPFR_RNDN); mpfr_mul(r26722, r26719, r26721, MPFR_RNDN); mpfr_add(r26723, r26717, r26722, MPFR_RNDN); mpfr_mul(r26724, r26701, r26723, MPFR_RNDN); mpfr_abs(r26725, r26724, MPFR_RNDN); return mpfr_get_d(r26725, MPFR_RNDN); } static mpfr_t r26726, r26727, r26728, r26729, r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r26726, "1", 10, MPFR_RNDN); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); mpfr_init_set_str(r26730, "2", 10, MPFR_RNDN); mpfr_init(r26731); mpfr_init(r26732); mpfr_init(r26733); mpfr_init_set_str(r26734, "3", 10, MPFR_RNDN); mpfr_init(r26735); mpfr_init(r26736); mpfr_init(r26737); mpfr_init(r26738); mpfr_init(r26739); mpfr_init_set_str(r26740, "5", 10, MPFR_RNDN); mpfr_init(r26741); mpfr_init(r26742); mpfr_init(r26743); mpfr_init(r26744); mpfr_init(r26745); mpfr_init_set_str(r26746, "21", 10, MPFR_RNDN); mpfr_init(r26747); mpfr_init(r26748); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); } double f_dm(double x) { ; mpfr_const_pi(r26727, MPFR_RNDN); mpfr_sqrt(r26728, r26727, MPFR_RNDN); mpfr_div(r26729, r26726, r26728, MPFR_RNDN); ; mpfr_set_d(r26731, x, MPFR_RNDN); mpfr_abs(r26732, r26731, MPFR_RNDN); mpfr_mul(r26733, r26730, r26732, MPFR_RNDN); ; mpfr_div(r26735, r26730, r26734, MPFR_RNDN); mpfr_mul(r26736, r26732, r26732, MPFR_RNDN); mpfr_mul(r26737, r26736, r26732, MPFR_RNDN); mpfr_mul(r26738, r26735, r26737, MPFR_RNDN); mpfr_add(r26739, r26733, r26738, MPFR_RNDN); ; mpfr_div(r26741, r26726, r26740, MPFR_RNDN); mpfr_mul(r26742, r26737, r26732, MPFR_RNDN); mpfr_mul(r26743, r26742, r26732, MPFR_RNDN); mpfr_mul(r26744, r26741, r26743, MPFR_RNDN); mpfr_add(r26745, r26739, r26744, MPFR_RNDN); ; mpfr_div(r26747, r26726, r26746, MPFR_RNDN); mpfr_mul(r26748, r26743, r26732, MPFR_RNDN); mpfr_mul(r26749, r26748, r26732, MPFR_RNDN); mpfr_mul(r26750, r26747, r26749, MPFR_RNDN); mpfr_add(r26751, r26745, r26750, MPFR_RNDN); mpfr_mul(r26752, r26729, r26751, MPFR_RNDN); mpfr_abs(r26753, r26752, MPFR_RNDN); return mpfr_get_d(r26753, MPFR_RNDN); }