#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 r24530 = 1; float r24531 = atan2(1.0, 0.0); float r24532 = sqrt(r24531); float r24533 = r24530 / r24532; float r24534 = 2; float r24535 = x; float r24536 = fabs(r24535); float r24537 = r24534 * r24536; float r24538 = 3; float r24539 = r24534 / r24538; float r24540 = r24536 * r24536; float r24541 = r24540 * r24536; float r24542 = r24539 * r24541; float r24543 = r24537 + r24542; float r24544 = 5; float r24545 = r24530 / r24544; float r24546 = r24541 * r24536; float r24547 = r24546 * r24536; float r24548 = r24545 * r24547; float r24549 = r24543 + r24548; float r24550 = 21; float r24551 = r24530 / r24550; float r24552 = r24547 * r24536; float r24553 = r24552 * r24536; float r24554 = r24551 * r24553; float r24555 = r24549 + r24554; float r24556 = r24533 * r24555; float r24557 = fabs(r24556); return r24557; } double f_id(double x) { double r24558 = 1; double r24559 = atan2(1.0, 0.0); double r24560 = sqrt(r24559); double r24561 = r24558 / r24560; double r24562 = 2; double r24563 = x; double r24564 = fabs(r24563); double r24565 = r24562 * r24564; double r24566 = 3; double r24567 = r24562 / r24566; double r24568 = r24564 * r24564; double r24569 = r24568 * r24564; double r24570 = r24567 * r24569; double r24571 = r24565 + r24570; double r24572 = 5; double r24573 = r24558 / r24572; double r24574 = r24569 * r24564; double r24575 = r24574 * r24564; double r24576 = r24573 * r24575; double r24577 = r24571 + r24576; double r24578 = 21; double r24579 = r24558 / r24578; double r24580 = r24575 * r24564; double r24581 = r24580 * r24564; double r24582 = r24579 * r24581; double r24583 = r24577 + r24582; double r24584 = r24561 * r24583; double r24585 = fabs(r24584); return r24585; } double f_of(float x) { float r24586 = 1; float r24587 = atan2(1.0, 0.0); float r24588 = sqrt(r24587); float r24589 = r24586 / r24588; float r24590 = 2; float r24591 = x; float r24592 = fabs(r24591); float r24593 = r24590 * r24592; float r24594 = 3; float r24595 = r24590 / r24594; float r24596 = r24592 * r24592; float r24597 = r24596 * r24592; float r24598 = r24595 * r24597; float r24599 = r24593 + r24598; float r24600 = 5; float r24601 = r24586 / r24600; float r24602 = r24597 * r24592; float r24603 = r24602 * r24592; float r24604 = r24601 * r24603; float r24605 = r24599 + r24604; float r24606 = r24589 * r24605; float r24607 = r24592 / r24588; float r24608 = 21; float r24609 = r24607 / r24608; float r24610 = pow(r24592, r24594); float r24611 = r24610 * r24610; float r24612 = r24609 * r24611; float r24613 = r24606 + r24612; float r24614 = fabs(r24613); return r24614; } double f_od(double x) { double r24615 = 1; double r24616 = atan2(1.0, 0.0); double r24617 = sqrt(r24616); double r24618 = r24615 / r24617; double r24619 = 2; double r24620 = x; double r24621 = fabs(r24620); double r24622 = r24619 * r24621; double r24623 = 3; double r24624 = r24619 / r24623; double r24625 = r24621 * r24621; double r24626 = r24625 * r24621; double r24627 = r24624 * r24626; double r24628 = r24622 + r24627; double r24629 = 5; double r24630 = r24615 / r24629; double r24631 = r24626 * r24621; double r24632 = r24631 * r24621; double r24633 = r24630 * r24632; double r24634 = r24628 + r24633; double r24635 = r24618 * r24634; double r24636 = r24621 / r24617; double r24637 = 21; double r24638 = r24636 / r24637; double r24639 = pow(r24621, r24623); double r24640 = r24639 * r24639; double r24641 = r24638 * r24640; double r24642 = r24635 + r24641; double r24643 = fabs(r24642); return r24643; } 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 r24644, r24645, r24646, r24647, r24648, r24649, r24650, r24651, r24652, r24653, r24654, r24655, r24656, r24657, r24658, r24659, r24660, r24661, r24662, r24663, r24664, r24665, r24666, r24667, r24668, r24669, r24670, r24671; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24644, "1", 10, MPFR_RNDN); mpfr_init(r24645); mpfr_init(r24646); mpfr_init(r24647); mpfr_init_set_str(r24648, "2", 10, MPFR_RNDN); mpfr_init(r24649); mpfr_init(r24650); mpfr_init(r24651); mpfr_init_set_str(r24652, "3", 10, MPFR_RNDN); mpfr_init(r24653); mpfr_init(r24654); mpfr_init(r24655); mpfr_init(r24656); mpfr_init(r24657); mpfr_init_set_str(r24658, "5", 10, MPFR_RNDN); mpfr_init(r24659); mpfr_init(r24660); mpfr_init(r24661); mpfr_init(r24662); mpfr_init(r24663); mpfr_init_set_str(r24664, "21", 10, MPFR_RNDN); mpfr_init(r24665); mpfr_init(r24666); mpfr_init(r24667); mpfr_init(r24668); mpfr_init(r24669); mpfr_init(r24670); mpfr_init(r24671); } double f_im(double x) { ; mpfr_const_pi(r24645, MPFR_RNDN); mpfr_sqrt(r24646, r24645, MPFR_RNDN); mpfr_div(r24647, r24644, r24646, MPFR_RNDN); ; mpfr_set_d(r24649, x, MPFR_RNDN); mpfr_abs(r24650, r24649, MPFR_RNDN); mpfr_mul(r24651, r24648, r24650, MPFR_RNDN); ; mpfr_div(r24653, r24648, r24652, MPFR_RNDN); mpfr_mul(r24654, r24650, r24650, MPFR_RNDN); mpfr_mul(r24655, r24654, r24650, MPFR_RNDN); mpfr_mul(r24656, r24653, r24655, MPFR_RNDN); mpfr_add(r24657, r24651, r24656, MPFR_RNDN); ; mpfr_div(r24659, r24644, r24658, MPFR_RNDN); mpfr_mul(r24660, r24655, r24650, MPFR_RNDN); mpfr_mul(r24661, r24660, r24650, MPFR_RNDN); mpfr_mul(r24662, r24659, r24661, MPFR_RNDN); mpfr_add(r24663, r24657, r24662, MPFR_RNDN); ; mpfr_div(r24665, r24644, r24664, MPFR_RNDN); mpfr_mul(r24666, r24661, r24650, MPFR_RNDN); mpfr_mul(r24667, r24666, r24650, MPFR_RNDN); mpfr_mul(r24668, r24665, r24667, MPFR_RNDN); mpfr_add(r24669, r24663, r24668, MPFR_RNDN); mpfr_mul(r24670, r24647, r24669, MPFR_RNDN); mpfr_abs(r24671, r24670, MPFR_RNDN); return mpfr_get_d(r24671, MPFR_RNDN); } static mpfr_t r24672, r24673, r24674, r24675, r24676, r24677, r24678, r24679, r24680, r24681, r24682, r24683, r24684, r24685, r24686, r24687, r24688, r24689, r24690, r24691, r24692, r24693, r24694, r24695, r24696, r24697, r24698, r24699, r24700; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24672, "1", 10, MPFR_RNDN); mpfr_init(r24673); mpfr_init(r24674); mpfr_init(r24675); mpfr_init_set_str(r24676, "2", 10, MPFR_RNDN); mpfr_init(r24677); mpfr_init(r24678); mpfr_init(r24679); mpfr_init_set_str(r24680, "3", 10, MPFR_RNDN); mpfr_init(r24681); mpfr_init(r24682); mpfr_init(r24683); mpfr_init(r24684); mpfr_init(r24685); mpfr_init_set_str(r24686, "5", 10, MPFR_RNDN); mpfr_init(r24687); mpfr_init(r24688); mpfr_init(r24689); mpfr_init(r24690); mpfr_init(r24691); mpfr_init(r24692); mpfr_init(r24693); mpfr_init_set_str(r24694, "21", 10, MPFR_RNDN); mpfr_init(r24695); mpfr_init(r24696); mpfr_init(r24697); mpfr_init(r24698); mpfr_init(r24699); mpfr_init(r24700); } double f_fm(double x) { ; mpfr_const_pi(r24673, MPFR_RNDN); mpfr_sqrt(r24674, r24673, MPFR_RNDN); mpfr_div(r24675, r24672, r24674, MPFR_RNDN); ; mpfr_set_d(r24677, x, MPFR_RNDN); mpfr_abs(r24678, r24677, MPFR_RNDN); mpfr_mul(r24679, r24676, r24678, MPFR_RNDN); ; mpfr_div(r24681, r24676, r24680, MPFR_RNDN); mpfr_mul(r24682, r24678, r24678, MPFR_RNDN); mpfr_mul(r24683, r24682, r24678, MPFR_RNDN); mpfr_mul(r24684, r24681, r24683, MPFR_RNDN); mpfr_add(r24685, r24679, r24684, MPFR_RNDN); ; mpfr_div(r24687, r24672, r24686, MPFR_RNDN); mpfr_mul(r24688, r24683, r24678, MPFR_RNDN); mpfr_mul(r24689, r24688, r24678, MPFR_RNDN); mpfr_mul(r24690, r24687, r24689, MPFR_RNDN); mpfr_add(r24691, r24685, r24690, MPFR_RNDN); mpfr_mul(r24692, r24675, r24691, MPFR_RNDN); mpfr_div(r24693, r24678, r24674, MPFR_RNDN); ; mpfr_div(r24695, r24693, r24694, MPFR_RNDN); mpfr_pow(r24696, r24678, r24680, MPFR_RNDN); mpfr_mul(r24697, r24696, r24696, MPFR_RNDN); mpfr_mul(r24698, r24695, r24697, MPFR_RNDN); mpfr_add(r24699, r24692, r24698, MPFR_RNDN); mpfr_abs(r24700, r24699, MPFR_RNDN); return mpfr_get_d(r24700, MPFR_RNDN); } static mpfr_t r24701, r24702, r24703, r24704, r24705, r24706, r24707, r24708, r24709, r24710, r24711, r24712, r24713, r24714, r24715, r24716, r24717, r24718, r24719, r24720, r24721, r24722, r24723, r24724, r24725, r24726, r24727, r24728, r24729; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24701, "1", 10, MPFR_RNDN); mpfr_init(r24702); mpfr_init(r24703); mpfr_init(r24704); mpfr_init_set_str(r24705, "2", 10, MPFR_RNDN); mpfr_init(r24706); mpfr_init(r24707); mpfr_init(r24708); mpfr_init_set_str(r24709, "3", 10, MPFR_RNDN); mpfr_init(r24710); mpfr_init(r24711); mpfr_init(r24712); mpfr_init(r24713); mpfr_init(r24714); mpfr_init_set_str(r24715, "5", 10, MPFR_RNDN); mpfr_init(r24716); mpfr_init(r24717); mpfr_init(r24718); mpfr_init(r24719); mpfr_init(r24720); mpfr_init(r24721); mpfr_init(r24722); mpfr_init_set_str(r24723, "21", 10, MPFR_RNDN); mpfr_init(r24724); mpfr_init(r24725); mpfr_init(r24726); mpfr_init(r24727); mpfr_init(r24728); mpfr_init(r24729); } double f_dm(double x) { ; mpfr_const_pi(r24702, MPFR_RNDN); mpfr_sqrt(r24703, r24702, MPFR_RNDN); mpfr_div(r24704, r24701, r24703, MPFR_RNDN); ; mpfr_set_d(r24706, x, MPFR_RNDN); mpfr_abs(r24707, r24706, MPFR_RNDN); mpfr_mul(r24708, r24705, r24707, MPFR_RNDN); ; mpfr_div(r24710, r24705, r24709, MPFR_RNDN); mpfr_mul(r24711, r24707, r24707, MPFR_RNDN); mpfr_mul(r24712, r24711, r24707, MPFR_RNDN); mpfr_mul(r24713, r24710, r24712, MPFR_RNDN); mpfr_add(r24714, r24708, r24713, MPFR_RNDN); ; mpfr_div(r24716, r24701, r24715, MPFR_RNDN); mpfr_mul(r24717, r24712, r24707, MPFR_RNDN); mpfr_mul(r24718, r24717, r24707, MPFR_RNDN); mpfr_mul(r24719, r24716, r24718, MPFR_RNDN); mpfr_add(r24720, r24714, r24719, MPFR_RNDN); mpfr_mul(r24721, r24704, r24720, MPFR_RNDN); mpfr_div(r24722, r24707, r24703, MPFR_RNDN); ; mpfr_div(r24724, r24722, r24723, MPFR_RNDN); mpfr_pow(r24725, r24707, r24709, MPFR_RNDN); mpfr_mul(r24726, r24725, r24725, MPFR_RNDN); mpfr_mul(r24727, r24724, r24726, MPFR_RNDN); mpfr_add(r24728, r24721, r24727, MPFR_RNDN); mpfr_abs(r24729, r24728, MPFR_RNDN); return mpfr_get_d(r24729, MPFR_RNDN); }