#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 r24571 = 1; float r24572 = atan2(1.0, 0.0); float r24573 = sqrt(r24572); float r24574 = r24571 / r24573; float r24575 = 2; float r24576 = x; float r24577 = fabs(r24576); float r24578 = r24575 * r24577; float r24579 = 3; float r24580 = r24575 / r24579; float r24581 = r24577 * r24577; float r24582 = r24581 * r24577; float r24583 = r24580 * r24582; float r24584 = r24578 + r24583; float r24585 = 5; float r24586 = r24571 / r24585; float r24587 = r24582 * r24577; float r24588 = r24587 * r24577; float r24589 = r24586 * r24588; float r24590 = r24584 + r24589; float r24591 = 21; float r24592 = r24571 / r24591; float r24593 = r24588 * r24577; float r24594 = r24593 * r24577; float r24595 = r24592 * r24594; float r24596 = r24590 + r24595; float r24597 = r24574 * r24596; float r24598 = fabs(r24597); return r24598; } double f_id(double x) { double r24599 = 1; double r24600 = atan2(1.0, 0.0); double r24601 = sqrt(r24600); double r24602 = r24599 / r24601; double r24603 = 2; double r24604 = x; double r24605 = fabs(r24604); double r24606 = r24603 * r24605; double r24607 = 3; double r24608 = r24603 / r24607; double r24609 = r24605 * r24605; double r24610 = r24609 * r24605; double r24611 = r24608 * r24610; double r24612 = r24606 + r24611; double r24613 = 5; double r24614 = r24599 / r24613; double r24615 = r24610 * r24605; double r24616 = r24615 * r24605; double r24617 = r24614 * r24616; double r24618 = r24612 + r24617; double r24619 = 21; double r24620 = r24599 / r24619; double r24621 = r24616 * r24605; double r24622 = r24621 * r24605; double r24623 = r24620 * r24622; double r24624 = r24618 + r24623; double r24625 = r24602 * r24624; double r24626 = fabs(r24625); return r24626; } double f_of(float x) { float r24627 = 1; float r24628 = atan2(1.0, 0.0); float r24629 = sqrt(r24628); float r24630 = r24627 / r24629; float r24631 = 2; float r24632 = x; float r24633 = fabs(r24632); float r24634 = r24631 * r24633; float r24635 = 3; float r24636 = r24631 / r24635; float r24637 = r24633 * r24633; float r24638 = r24637 * r24633; float r24639 = r24636 * r24638; float r24640 = r24634 + r24639; float r24641 = 5; float r24642 = r24627 / r24641; float r24643 = r24638 * r24633; float r24644 = r24643 * r24633; float r24645 = r24642 * r24644; float r24646 = r24640 + r24645; float r24647 = 21; float r24648 = r24627 / r24647; float r24649 = r24644 * r24633; float r24650 = r24649 * r24633; float r24651 = r24648 * r24650; float r24652 = r24646 + r24651; float r24653 = r24630 * r24652; float r24654 = fabs(r24653); return r24654; } double f_od(double x) { double r24655 = 1; double r24656 = atan2(1.0, 0.0); double r24657 = sqrt(r24656); double r24658 = r24655 / r24657; double r24659 = 2; double r24660 = x; double r24661 = fabs(r24660); double r24662 = r24659 * r24661; double r24663 = 3; double r24664 = r24659 / r24663; double r24665 = r24661 * r24661; double r24666 = r24665 * r24661; double r24667 = r24664 * r24666; double r24668 = r24662 + r24667; double r24669 = 5; double r24670 = r24655 / r24669; double r24671 = r24666 * r24661; double r24672 = r24671 * r24661; double r24673 = r24670 * r24672; double r24674 = r24668 + r24673; double r24675 = 21; double r24676 = r24655 / r24675; double r24677 = r24672 * r24661; double r24678 = r24677 * r24661; double r24679 = r24676 * r24678; double r24680 = r24674 + r24679; double r24681 = r24658 * r24680; double r24682 = fabs(r24681); return r24682; } 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 r24683, r24684, r24685, r24686, r24687, r24688, r24689, r24690, r24691, r24692, r24693, r24694, r24695, r24696, r24697, r24698, r24699, r24700, r24701, r24702, r24703, r24704, r24705, r24706, r24707, r24708, r24709, r24710; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24683, "1", 10, MPFR_RNDN); mpfr_init(r24684); mpfr_init(r24685); mpfr_init(r24686); mpfr_init_set_str(r24687, "2", 10, MPFR_RNDN); mpfr_init(r24688); mpfr_init(r24689); mpfr_init(r24690); mpfr_init_set_str(r24691, "3", 10, MPFR_RNDN); mpfr_init(r24692); mpfr_init(r24693); mpfr_init(r24694); mpfr_init(r24695); mpfr_init(r24696); mpfr_init_set_str(r24697, "5", 10, MPFR_RNDN); mpfr_init(r24698); mpfr_init(r24699); mpfr_init(r24700); mpfr_init(r24701); mpfr_init(r24702); mpfr_init_set_str(r24703, "21", 10, MPFR_RNDN); mpfr_init(r24704); mpfr_init(r24705); mpfr_init(r24706); mpfr_init(r24707); mpfr_init(r24708); mpfr_init(r24709); mpfr_init(r24710); } double f_im(double x) { ; mpfr_const_pi(r24684, MPFR_RNDN); mpfr_sqrt(r24685, r24684, MPFR_RNDN); mpfr_div(r24686, r24683, r24685, MPFR_RNDN); ; mpfr_set_d(r24688, x, MPFR_RNDN); mpfr_abs(r24689, r24688, MPFR_RNDN); mpfr_mul(r24690, r24687, r24689, MPFR_RNDN); ; mpfr_div(r24692, r24687, r24691, MPFR_RNDN); mpfr_mul(r24693, r24689, r24689, MPFR_RNDN); mpfr_mul(r24694, r24693, r24689, MPFR_RNDN); mpfr_mul(r24695, r24692, r24694, MPFR_RNDN); mpfr_add(r24696, r24690, r24695, MPFR_RNDN); ; mpfr_div(r24698, r24683, r24697, MPFR_RNDN); mpfr_mul(r24699, r24694, r24689, MPFR_RNDN); mpfr_mul(r24700, r24699, r24689, MPFR_RNDN); mpfr_mul(r24701, r24698, r24700, MPFR_RNDN); mpfr_add(r24702, r24696, r24701, MPFR_RNDN); ; mpfr_div(r24704, r24683, r24703, MPFR_RNDN); mpfr_mul(r24705, r24700, r24689, MPFR_RNDN); mpfr_mul(r24706, r24705, r24689, MPFR_RNDN); mpfr_mul(r24707, r24704, r24706, MPFR_RNDN); mpfr_add(r24708, r24702, r24707, MPFR_RNDN); mpfr_mul(r24709, r24686, r24708, MPFR_RNDN); mpfr_abs(r24710, r24709, MPFR_RNDN); return mpfr_get_d(r24710, MPFR_RNDN); } static mpfr_t r24711, r24712, r24713, r24714, r24715, r24716, r24717, r24718, r24719, r24720, r24721, r24722, r24723, r24724, r24725, r24726, r24727, r24728, r24729, r24730, r24731, r24732, r24733, r24734, r24735, r24736, r24737, r24738; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24711, "1", 10, MPFR_RNDN); mpfr_init(r24712); mpfr_init(r24713); mpfr_init(r24714); mpfr_init_set_str(r24715, "2", 10, MPFR_RNDN); mpfr_init(r24716); mpfr_init(r24717); mpfr_init(r24718); mpfr_init_set_str(r24719, "3", 10, MPFR_RNDN); mpfr_init(r24720); mpfr_init(r24721); mpfr_init(r24722); mpfr_init(r24723); mpfr_init(r24724); mpfr_init_set_str(r24725, "5", 10, MPFR_RNDN); mpfr_init(r24726); mpfr_init(r24727); mpfr_init(r24728); mpfr_init(r24729); mpfr_init(r24730); mpfr_init_set_str(r24731, "21", 10, MPFR_RNDN); mpfr_init(r24732); mpfr_init(r24733); mpfr_init(r24734); mpfr_init(r24735); mpfr_init(r24736); mpfr_init(r24737); mpfr_init(r24738); } double f_fm(double x) { ; mpfr_const_pi(r24712, MPFR_RNDN); mpfr_sqrt(r24713, r24712, MPFR_RNDN); mpfr_div(r24714, r24711, r24713, MPFR_RNDN); ; mpfr_set_d(r24716, x, MPFR_RNDN); mpfr_abs(r24717, r24716, MPFR_RNDN); mpfr_mul(r24718, r24715, r24717, MPFR_RNDN); ; mpfr_div(r24720, r24715, r24719, MPFR_RNDN); mpfr_mul(r24721, r24717, r24717, MPFR_RNDN); mpfr_mul(r24722, r24721, r24717, MPFR_RNDN); mpfr_mul(r24723, r24720, r24722, MPFR_RNDN); mpfr_add(r24724, r24718, r24723, MPFR_RNDN); ; mpfr_div(r24726, r24711, r24725, MPFR_RNDN); mpfr_mul(r24727, r24722, r24717, MPFR_RNDN); mpfr_mul(r24728, r24727, r24717, MPFR_RNDN); mpfr_mul(r24729, r24726, r24728, MPFR_RNDN); mpfr_add(r24730, r24724, r24729, MPFR_RNDN); ; mpfr_div(r24732, r24711, r24731, MPFR_RNDN); mpfr_mul(r24733, r24728, r24717, MPFR_RNDN); mpfr_mul(r24734, r24733, r24717, MPFR_RNDN); mpfr_mul(r24735, r24732, r24734, MPFR_RNDN); mpfr_add(r24736, r24730, r24735, MPFR_RNDN); mpfr_mul(r24737, r24714, r24736, MPFR_RNDN); mpfr_abs(r24738, r24737, MPFR_RNDN); return mpfr_get_d(r24738, MPFR_RNDN); } static mpfr_t r24739, r24740, r24741, r24742, r24743, r24744, r24745, r24746, r24747, r24748, r24749, r24750, r24751, r24752, r24753, r24754, r24755, r24756, r24757, r24758, r24759, r24760, r24761, r24762, r24763, r24764, r24765, r24766; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24739, "1", 10, MPFR_RNDN); mpfr_init(r24740); mpfr_init(r24741); mpfr_init(r24742); mpfr_init_set_str(r24743, "2", 10, MPFR_RNDN); mpfr_init(r24744); mpfr_init(r24745); mpfr_init(r24746); mpfr_init_set_str(r24747, "3", 10, MPFR_RNDN); mpfr_init(r24748); mpfr_init(r24749); mpfr_init(r24750); mpfr_init(r24751); mpfr_init(r24752); mpfr_init_set_str(r24753, "5", 10, MPFR_RNDN); mpfr_init(r24754); mpfr_init(r24755); mpfr_init(r24756); mpfr_init(r24757); mpfr_init(r24758); mpfr_init_set_str(r24759, "21", 10, MPFR_RNDN); mpfr_init(r24760); mpfr_init(r24761); mpfr_init(r24762); mpfr_init(r24763); mpfr_init(r24764); mpfr_init(r24765); mpfr_init(r24766); } double f_dm(double x) { ; mpfr_const_pi(r24740, MPFR_RNDN); mpfr_sqrt(r24741, r24740, MPFR_RNDN); mpfr_div(r24742, r24739, r24741, MPFR_RNDN); ; mpfr_set_d(r24744, x, MPFR_RNDN); mpfr_abs(r24745, r24744, MPFR_RNDN); mpfr_mul(r24746, r24743, r24745, MPFR_RNDN); ; mpfr_div(r24748, r24743, r24747, MPFR_RNDN); mpfr_mul(r24749, r24745, r24745, MPFR_RNDN); mpfr_mul(r24750, r24749, r24745, MPFR_RNDN); mpfr_mul(r24751, r24748, r24750, MPFR_RNDN); mpfr_add(r24752, r24746, r24751, MPFR_RNDN); ; mpfr_div(r24754, r24739, r24753, MPFR_RNDN); mpfr_mul(r24755, r24750, r24745, MPFR_RNDN); mpfr_mul(r24756, r24755, r24745, MPFR_RNDN); mpfr_mul(r24757, r24754, r24756, MPFR_RNDN); mpfr_add(r24758, r24752, r24757, MPFR_RNDN); ; mpfr_div(r24760, r24739, r24759, MPFR_RNDN); mpfr_mul(r24761, r24756, r24745, MPFR_RNDN); mpfr_mul(r24762, r24761, r24745, MPFR_RNDN); mpfr_mul(r24763, r24760, r24762, MPFR_RNDN); mpfr_add(r24764, r24758, r24763, MPFR_RNDN); mpfr_mul(r24765, r24742, r24764, MPFR_RNDN); mpfr_abs(r24766, r24765, MPFR_RNDN); return mpfr_get_d(r24766, MPFR_RNDN); }