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