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