#include #include #include #include #include char *name = "Jmat.Real.gamma, branch z less than 0.5"; double f_if(float z) { float r10460 = atan2(1.0, 0.0); float r10461 = z; float r10462 = r10460 * r10461; float r10463 = sin(r10462); float r10464 = r10460 / r10463; float r10465 = 2.0f; float r10466 = r10460 * r10465; float r10467 = sqrt(r10466); float r10468 = 1.0f; float r10469 = r10468 - r10461; float r10470 = r10469 - r10468; float r10471 = 7.0f; float r10472 = r10470 + r10471; float r10473 = 0.5f; float r10474 = r10472 + r10473; float r10475 = r10470 + r10473; float r10476 = pow(r10474, r10475); float r10477 = r10467 * r10476; float r10478 = -r10474; float r10479 = exp(r10478); float r10480 = r10477 * r10479; float r10481 = 0.9999999999998099f; float r10482 = 676.5203681218851f; float r10483 = r10470 + r10468; float r10484 = r10482 / r10483; float r10485 = r10481 + r10484; float r10486 = -1259.1392167224028f; float r10487 = r10470 + r10465; float r10488 = r10486 / r10487; float r10489 = r10485 + r10488; float r10490 = 771.3234287776531f; float r10491 = 3.0f; float r10492 = r10470 + r10491; float r10493 = r10490 / r10492; float r10494 = r10489 + r10493; float r10495 = -176.6150291621406f; float r10496 = 4.0f; float r10497 = r10470 + r10496; float r10498 = r10495 / r10497; float r10499 = r10494 + r10498; float r10500 = 12.507343278686905f; float r10501 = 5.0f; float r10502 = r10470 + r10501; float r10503 = r10500 / r10502; float r10504 = r10499 + r10503; float r10505 = -0.13857109526572012f; float r10506 = 6.0f; float r10507 = r10470 + r10506; float r10508 = r10505 / r10507; float r10509 = r10504 + r10508; float r10510 = 9.984369578019572e-06f; float r10511 = r10510 / r10472; float r10512 = r10509 + r10511; float r10513 = 1.5056327351493116e-07f; float r10514 = 8.0f; float r10515 = r10470 + r10514; float r10516 = r10513 / r10515; float r10517 = r10512 + r10516; float r10518 = r10480 * r10517; float r10519 = r10464 * r10518; return r10519; } double f_id(double z) { double r10520 = atan2(1.0, 0.0); double r10521 = z; double r10522 = r10520 * r10521; double r10523 = sin(r10522); double r10524 = r10520 / r10523; double r10525 = 2.0; double r10526 = r10520 * r10525; double r10527 = sqrt(r10526); double r10528 = 1.0; double r10529 = r10528 - r10521; double r10530 = r10529 - r10528; double r10531 = 7.0; double r10532 = r10530 + r10531; double r10533 = 0.5; double r10534 = r10532 + r10533; double r10535 = r10530 + r10533; double r10536 = pow(r10534, r10535); double r10537 = r10527 * r10536; double r10538 = -r10534; double r10539 = exp(r10538); double r10540 = r10537 * r10539; double r10541 = 0.9999999999998099; double r10542 = 676.5203681218851; double r10543 = r10530 + r10528; double r10544 = r10542 / r10543; double r10545 = r10541 + r10544; double r10546 = -1259.1392167224028; double r10547 = r10530 + r10525; double r10548 = r10546 / r10547; double r10549 = r10545 + r10548; double r10550 = 771.3234287776531; double r10551 = 3.0; double r10552 = r10530 + r10551; double r10553 = r10550 / r10552; double r10554 = r10549 + r10553; double r10555 = -176.6150291621406; double r10556 = 4.0; double r10557 = r10530 + r10556; double r10558 = r10555 / r10557; double r10559 = r10554 + r10558; double r10560 = 12.507343278686905; double r10561 = 5.0; double r10562 = r10530 + r10561; double r10563 = r10560 / r10562; double r10564 = r10559 + r10563; double r10565 = -0.13857109526572012; double r10566 = 6.0; double r10567 = r10530 + r10566; double r10568 = r10565 / r10567; double r10569 = r10564 + r10568; double r10570 = 9.984369578019572e-06; double r10571 = r10570 / r10532; double r10572 = r10569 + r10571; double r10573 = 1.5056327351493116e-07; double r10574 = 8.0; double r10575 = r10530 + r10574; double r10576 = r10573 / r10575; double r10577 = r10572 + r10576; double r10578 = r10540 * r10577; double r10579 = r10524 * r10578; return r10579; } double f_of(float z) { float r10580 = 12.507343278686905f; float r10581 = 1.0f; float r10582 = z; float r10583 = r10581 - r10582; float r10584 = 5.0f; float r10585 = r10581 - r10584; float r10586 = r10583 - r10585; float r10587 = r10580 / r10586; float r10588 = -0.13857109526572012f; float r10589 = 6.0f; float r10590 = r10589 + r10581; float r10591 = r10581 + r10582; float r10592 = r10590 - r10591; float r10593 = r10588 / r10592; float r10594 = r10587 + r10593; float r10595 = 771.3234287776531f; float r10596 = 3.0f; float r10597 = r10581 - r10596; float r10598 = r10583 - r10597; float r10599 = r10595 / r10598; float r10600 = -176.6150291621406f; float r10601 = 4.0f; float r10602 = r10581 + r10601; float r10603 = r10602 - r10591; float r10604 = r10600 / r10603; float r10605 = r10599 + r10604; float r10606 = r10594 + r10605; float r10607 = 1.5056327351493116e-07f; float r10608 = 8.0f; float r10609 = r10581 + r10608; float r10610 = r10609 - r10591; float r10611 = r10607 / r10610; float r10612 = 9.984369578019572e-06f; float r10613 = 7.0f; float r10614 = r10581 + r10613; float r10615 = r10614 - r10591; float r10616 = r10612 / r10615; float r10617 = r10611 + r10616; float r10618 = 0.9999999999998099f; float r10619 = 676.5203681218851f; float r10620 = r10619 / r10583; float r10621 = r10618 + r10620; float r10622 = -1259.1392167224028f; float r10623 = 2.0f; float r10624 = r10623 + r10581; float r10625 = r10624 - r10591; float r10626 = r10622 / r10625; float r10627 = r10621 + r10626; float r10628 = r10617 + r10627; float r10629 = r10606 + r10628; float r10630 = atan2(1.0, 0.0); float r10631 = r10623 * r10630; float r10632 = sqrt(r10631); float r10633 = r10632 * r10630; float r10634 = r10630 * r10582; float r10635 = sin(r10634); float r10636 = r10633 / r10635; float r10637 = 0.5f; float r10638 = r10615 + r10637; float r10639 = r10581 - r10637; float r10640 = r10583 - r10639; float r10641 = pow(r10638, r10640); float r10642 = r10636 * r10641; float r10643 = exp(r10638); float r10644 = r10642 / r10643; float r10645 = r10629 * r10644; return r10645; } double f_od(double z) { double r10646 = 12.507343278686905; double r10647 = 1.0; double r10648 = z; double r10649 = r10647 - r10648; double r10650 = 5.0; double r10651 = r10647 - r10650; double r10652 = r10649 - r10651; double r10653 = r10646 / r10652; double r10654 = -0.13857109526572012; double r10655 = 6.0; double r10656 = r10655 + r10647; double r10657 = r10647 + r10648; double r10658 = r10656 - r10657; double r10659 = r10654 / r10658; double r10660 = r10653 + r10659; double r10661 = 771.3234287776531; double r10662 = 3.0; double r10663 = r10647 - r10662; double r10664 = r10649 - r10663; double r10665 = r10661 / r10664; double r10666 = -176.6150291621406; double r10667 = 4.0; double r10668 = r10647 + r10667; double r10669 = r10668 - r10657; double r10670 = r10666 / r10669; double r10671 = r10665 + r10670; double r10672 = r10660 + r10671; double r10673 = 1.5056327351493116e-07; double r10674 = 8.0; double r10675 = r10647 + r10674; double r10676 = r10675 - r10657; double r10677 = r10673 / r10676; double r10678 = 9.984369578019572e-06; double r10679 = 7.0; double r10680 = r10647 + r10679; double r10681 = r10680 - r10657; double r10682 = r10678 / r10681; double r10683 = r10677 + r10682; double r10684 = 0.9999999999998099; double r10685 = 676.5203681218851; double r10686 = r10685 / r10649; double r10687 = r10684 + r10686; double r10688 = -1259.1392167224028; double r10689 = 2.0; double r10690 = r10689 + r10647; double r10691 = r10690 - r10657; double r10692 = r10688 / r10691; double r10693 = r10687 + r10692; double r10694 = r10683 + r10693; double r10695 = r10672 + r10694; double r10696 = atan2(1.0, 0.0); double r10697 = r10689 * r10696; double r10698 = sqrt(r10697); double r10699 = r10698 * r10696; double r10700 = r10696 * r10648; double r10701 = sin(r10700); double r10702 = r10699 / r10701; double r10703 = 0.5; double r10704 = r10681 + r10703; double r10705 = r10647 - r10703; double r10706 = r10649 - r10705; double r10707 = pow(r10704, r10706); double r10708 = r10702 * r10707; double r10709 = exp(r10704); double r10710 = r10708 / r10709; double r10711 = r10695 * r10710; return r10711; } 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 r10712, r10713, r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721, r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729, r10730, r10731, r10732, r10733, r10734, r10735, r10736, r10737, r10738, r10739, r10740, r10741, r10742, r10743, r10744, r10745, r10746, r10747, r10748, r10749, r10750, r10751, r10752, r10753, r10754, r10755, r10756, r10757, r10758, r10759, r10760, r10761, r10762, r10763, r10764, r10765, r10766, r10767, r10768, r10769, r10770, r10771; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r10712); mpfr_init(r10713); mpfr_init(r10714); mpfr_init(r10715); mpfr_init(r10716); mpfr_init_set_str(r10717, "2", 10, MPFR_RNDN); mpfr_init(r10718); mpfr_init(r10719); mpfr_init_set_str(r10720, "1", 10, MPFR_RNDN); mpfr_init(r10721); mpfr_init(r10722); mpfr_init_set_str(r10723, "7", 10, MPFR_RNDN); mpfr_init(r10724); mpfr_init_set_str(r10725, "0.5", 10, MPFR_RNDN); mpfr_init(r10726); mpfr_init(r10727); mpfr_init(r10728); mpfr_init(r10729); mpfr_init(r10730); mpfr_init(r10731); mpfr_init(r10732); mpfr_init_set_str(r10733, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r10734, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r10735); mpfr_init(r10736); mpfr_init(r10737); mpfr_init_set_str(r10738, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init(r10739); mpfr_init(r10740); mpfr_init(r10741); mpfr_init_set_str(r10742, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r10743, "3", 10, MPFR_RNDN); mpfr_init(r10744); mpfr_init(r10745); mpfr_init(r10746); mpfr_init_set_str(r10747, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r10748, "4", 10, MPFR_RNDN); mpfr_init(r10749); mpfr_init(r10750); mpfr_init(r10751); mpfr_init_set_str(r10752, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r10753, "5", 10, MPFR_RNDN); mpfr_init(r10754); mpfr_init(r10755); mpfr_init(r10756); mpfr_init_set_str(r10757, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r10758, "6", 10, MPFR_RNDN); mpfr_init(r10759); mpfr_init(r10760); mpfr_init(r10761); mpfr_init_set_str(r10762, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init(r10763); mpfr_init(r10764); mpfr_init_set_str(r10765, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10766, "8", 10, MPFR_RNDN); mpfr_init(r10767); mpfr_init(r10768); mpfr_init(r10769); mpfr_init(r10770); mpfr_init(r10771); } double f_im(double z) { mpfr_const_pi(r10712, MPFR_RNDN); mpfr_set_d(r10713, z, MPFR_RNDN); mpfr_mul(r10714, r10712, r10713, MPFR_RNDN); mpfr_sin(r10715, r10714, MPFR_RNDN); mpfr_div(r10716, r10712, r10715, MPFR_RNDN); ; mpfr_mul(r10718, r10712, r10717, MPFR_RNDN); mpfr_sqrt(r10719, r10718, MPFR_RNDN); ; mpfr_sub(r10721, r10720, r10713, MPFR_RNDN); mpfr_sub(r10722, r10721, r10720, MPFR_RNDN); ; mpfr_add(r10724, r10722, r10723, MPFR_RNDN); ; mpfr_add(r10726, r10724, r10725, MPFR_RNDN); mpfr_add(r10727, r10722, r10725, MPFR_RNDN); mpfr_pow(r10728, r10726, r10727, MPFR_RNDN); mpfr_mul(r10729, r10719, r10728, MPFR_RNDN); mpfr_neg(r10730, r10726, MPFR_RNDN); mpfr_exp(r10731, r10730, MPFR_RNDN); mpfr_mul(r10732, r10729, r10731, MPFR_RNDN); ; ; mpfr_add(r10735, r10722, r10720, MPFR_RNDN); mpfr_div(r10736, r10734, r10735, MPFR_RNDN); mpfr_add(r10737, r10733, r10736, MPFR_RNDN); ; mpfr_add(r10739, r10722, r10717, MPFR_RNDN); mpfr_div(r10740, r10738, r10739, MPFR_RNDN); mpfr_add(r10741, r10737, r10740, MPFR_RNDN); ; ; mpfr_add(r10744, r10722, r10743, MPFR_RNDN); mpfr_div(r10745, r10742, r10744, MPFR_RNDN); mpfr_add(r10746, r10741, r10745, MPFR_RNDN); ; ; mpfr_add(r10749, r10722, r10748, MPFR_RNDN); mpfr_div(r10750, r10747, r10749, MPFR_RNDN); mpfr_add(r10751, r10746, r10750, MPFR_RNDN); ; ; mpfr_add(r10754, r10722, r10753, MPFR_RNDN); mpfr_div(r10755, r10752, r10754, MPFR_RNDN); mpfr_add(r10756, r10751, r10755, MPFR_RNDN); ; ; mpfr_add(r10759, r10722, r10758, MPFR_RNDN); mpfr_div(r10760, r10757, r10759, MPFR_RNDN); mpfr_add(r10761, r10756, r10760, MPFR_RNDN); ; mpfr_div(r10763, r10762, r10724, MPFR_RNDN); mpfr_add(r10764, r10761, r10763, MPFR_RNDN); ; ; mpfr_add(r10767, r10722, r10766, MPFR_RNDN); mpfr_div(r10768, r10765, r10767, MPFR_RNDN); mpfr_add(r10769, r10764, r10768, MPFR_RNDN); mpfr_mul(r10770, r10732, r10769, MPFR_RNDN); mpfr_mul(r10771, r10716, r10770, MPFR_RNDN); return mpfr_get_d(r10771, MPFR_RNDN); } static mpfr_t r10772, r10773, r10774, r10775, r10776, r10777, r10778, r10779, r10780, r10781, r10782, r10783, r10784, r10785, r10786, r10787, r10788, r10789, r10790, r10791, r10792, r10793, r10794, r10795, r10796, r10797, r10798, r10799, r10800, r10801, r10802, r10803, r10804, r10805, r10806, r10807, r10808, r10809, r10810, r10811, r10812, r10813, r10814, r10815, r10816, r10817, r10818, r10819, r10820, r10821, r10822, r10823, r10824, r10825, r10826, r10827, r10828, r10829, r10830, r10831, r10832, r10833, r10834, r10835, r10836, r10837; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10772, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r10773, "1", 10, MPFR_RNDN); mpfr_init(r10774); mpfr_init(r10775); mpfr_init_set_str(r10776, "5", 10, MPFR_RNDN); mpfr_init(r10777); mpfr_init(r10778); mpfr_init(r10779); mpfr_init_set_str(r10780, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r10781, "6", 10, MPFR_RNDN); mpfr_init(r10782); mpfr_init(r10783); mpfr_init(r10784); mpfr_init(r10785); mpfr_init(r10786); mpfr_init_set_str(r10787, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r10788, "3", 10, MPFR_RNDN); mpfr_init(r10789); mpfr_init(r10790); mpfr_init(r10791); mpfr_init_set_str(r10792, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r10793, "4", 10, MPFR_RNDN); mpfr_init(r10794); mpfr_init(r10795); mpfr_init(r10796); mpfr_init(r10797); mpfr_init(r10798); mpfr_init_set_str(r10799, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10800, "8", 10, MPFR_RNDN); mpfr_init(r10801); mpfr_init(r10802); mpfr_init(r10803); mpfr_init_set_str(r10804, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r10805, "7", 10, MPFR_RNDN); mpfr_init(r10806); mpfr_init(r10807); mpfr_init(r10808); mpfr_init(r10809); mpfr_init_set_str(r10810, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r10811, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r10812); mpfr_init(r10813); mpfr_init_set_str(r10814, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r10815, "2", 10, MPFR_RNDN); mpfr_init(r10816); mpfr_init(r10817); mpfr_init(r10818); mpfr_init(r10819); mpfr_init(r10820); mpfr_init(r10821); mpfr_init(r10822); mpfr_init(r10823); mpfr_init(r10824); mpfr_init(r10825); mpfr_init(r10826); mpfr_init(r10827); mpfr_init(r10828); mpfr_init_set_str(r10829, "0.5", 10, MPFR_RNDN); mpfr_init(r10830); mpfr_init(r10831); mpfr_init(r10832); mpfr_init(r10833); mpfr_init(r10834); mpfr_init(r10835); mpfr_init(r10836); mpfr_init(r10837); } double f_fm(double z) { ; ; mpfr_set_d(r10774, z, MPFR_RNDN); mpfr_sub(r10775, r10773, r10774, MPFR_RNDN); ; mpfr_sub(r10777, r10773, r10776, MPFR_RNDN); mpfr_sub(r10778, r10775, r10777, MPFR_RNDN); mpfr_div(r10779, r10772, r10778, MPFR_RNDN); ; ; mpfr_add(r10782, r10781, r10773, MPFR_RNDN); mpfr_add(r10783, r10773, r10774, MPFR_RNDN); mpfr_sub(r10784, r10782, r10783, MPFR_RNDN); mpfr_div(r10785, r10780, r10784, MPFR_RNDN); mpfr_add(r10786, r10779, r10785, MPFR_RNDN); ; ; mpfr_sub(r10789, r10773, r10788, MPFR_RNDN); mpfr_sub(r10790, r10775, r10789, MPFR_RNDN); mpfr_div(r10791, r10787, r10790, MPFR_RNDN); ; ; mpfr_add(r10794, r10773, r10793, MPFR_RNDN); mpfr_sub(r10795, r10794, r10783, MPFR_RNDN); mpfr_div(r10796, r10792, r10795, MPFR_RNDN); mpfr_add(r10797, r10791, r10796, MPFR_RNDN); mpfr_add(r10798, r10786, r10797, MPFR_RNDN); ; ; mpfr_add(r10801, r10773, r10800, MPFR_RNDN); mpfr_sub(r10802, r10801, r10783, MPFR_RNDN); mpfr_div(r10803, r10799, r10802, MPFR_RNDN); ; ; mpfr_add(r10806, r10773, r10805, MPFR_RNDN); mpfr_sub(r10807, r10806, r10783, MPFR_RNDN); mpfr_div(r10808, r10804, r10807, MPFR_RNDN); mpfr_add(r10809, r10803, r10808, MPFR_RNDN); ; ; mpfr_div(r10812, r10811, r10775, MPFR_RNDN); mpfr_add(r10813, r10810, r10812, MPFR_RNDN); ; ; mpfr_add(r10816, r10815, r10773, MPFR_RNDN); mpfr_sub(r10817, r10816, r10783, MPFR_RNDN); mpfr_div(r10818, r10814, r10817, MPFR_RNDN); mpfr_add(r10819, r10813, r10818, MPFR_RNDN); mpfr_add(r10820, r10809, r10819, MPFR_RNDN); mpfr_add(r10821, r10798, r10820, MPFR_RNDN); mpfr_const_pi(r10822, MPFR_RNDN); mpfr_mul(r10823, r10815, r10822, MPFR_RNDN); mpfr_sqrt(r10824, r10823, MPFR_RNDN); mpfr_mul(r10825, r10824, r10822, MPFR_RNDN); mpfr_mul(r10826, r10822, r10774, MPFR_RNDN); mpfr_sin(r10827, r10826, MPFR_RNDN); mpfr_div(r10828, r10825, r10827, MPFR_RNDN); ; mpfr_add(r10830, r10807, r10829, MPFR_RNDN); mpfr_sub(r10831, r10773, r10829, MPFR_RNDN); mpfr_sub(r10832, r10775, r10831, MPFR_RNDN); mpfr_pow(r10833, r10830, r10832, MPFR_RNDN); mpfr_mul(r10834, r10828, r10833, MPFR_RNDN); mpfr_exp(r10835, r10830, MPFR_RNDN); mpfr_div(r10836, r10834, r10835, MPFR_RNDN); mpfr_mul(r10837, r10821, r10836, MPFR_RNDN); return mpfr_get_d(r10837, MPFR_RNDN); } static mpfr_t r10838, r10839, r10840, r10841, r10842, r10843, r10844, r10845, r10846, r10847, r10848, r10849, r10850, r10851, r10852, r10853, r10854, r10855, r10856, r10857, r10858, r10859, r10860, r10861, r10862, r10863, r10864, r10865, r10866, r10867, r10868, r10869, r10870, r10871, r10872, r10873, r10874, r10875, r10876, r10877, r10878, r10879, r10880, r10881, r10882, r10883, r10884, r10885, r10886, r10887, r10888, r10889, r10890, r10891, r10892, r10893, r10894, r10895, r10896, r10897, r10898, r10899, r10900, r10901, r10902, r10903; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r10838, "12.507343278686905", 10, MPFR_RNDN); mpfr_init_set_str(r10839, "1", 10, MPFR_RNDN); mpfr_init(r10840); mpfr_init(r10841); mpfr_init_set_str(r10842, "5", 10, MPFR_RNDN); mpfr_init(r10843); mpfr_init(r10844); mpfr_init(r10845); mpfr_init_set_str(r10846, "-0.13857109526572012", 10, MPFR_RNDN); mpfr_init_set_str(r10847, "6", 10, MPFR_RNDN); mpfr_init(r10848); mpfr_init(r10849); mpfr_init(r10850); mpfr_init(r10851); mpfr_init(r10852); mpfr_init_set_str(r10853, "771.3234287776531", 10, MPFR_RNDN); mpfr_init_set_str(r10854, "3", 10, MPFR_RNDN); mpfr_init(r10855); mpfr_init(r10856); mpfr_init(r10857); mpfr_init_set_str(r10858, "-176.6150291621406", 10, MPFR_RNDN); mpfr_init_set_str(r10859, "4", 10, MPFR_RNDN); mpfr_init(r10860); mpfr_init(r10861); mpfr_init(r10862); mpfr_init(r10863); mpfr_init(r10864); mpfr_init_set_str(r10865, "1.5056327351493116e-07", 10, MPFR_RNDN); mpfr_init_set_str(r10866, "8", 10, MPFR_RNDN); mpfr_init(r10867); mpfr_init(r10868); mpfr_init(r10869); mpfr_init_set_str(r10870, "9.984369578019572e-06", 10, MPFR_RNDN); mpfr_init_set_str(r10871, "7", 10, MPFR_RNDN); mpfr_init(r10872); mpfr_init(r10873); mpfr_init(r10874); mpfr_init(r10875); mpfr_init_set_str(r10876, "0.9999999999998099", 10, MPFR_RNDN); mpfr_init_set_str(r10877, "676.5203681218851", 10, MPFR_RNDN); mpfr_init(r10878); mpfr_init(r10879); mpfr_init_set_str(r10880, "-1259.1392167224028", 10, MPFR_RNDN); mpfr_init_set_str(r10881, "2", 10, MPFR_RNDN); mpfr_init(r10882); mpfr_init(r10883); mpfr_init(r10884); mpfr_init(r10885); mpfr_init(r10886); mpfr_init(r10887); mpfr_init(r10888); mpfr_init(r10889); mpfr_init(r10890); mpfr_init(r10891); mpfr_init(r10892); mpfr_init(r10893); mpfr_init(r10894); mpfr_init_set_str(r10895, "0.5", 10, MPFR_RNDN); mpfr_init(r10896); mpfr_init(r10897); mpfr_init(r10898); mpfr_init(r10899); mpfr_init(r10900); mpfr_init(r10901); mpfr_init(r10902); mpfr_init(r10903); } double f_dm(double z) { ; ; mpfr_set_d(r10840, z, MPFR_RNDN); mpfr_sub(r10841, r10839, r10840, MPFR_RNDN); ; mpfr_sub(r10843, r10839, r10842, MPFR_RNDN); mpfr_sub(r10844, r10841, r10843, MPFR_RNDN); mpfr_div(r10845, r10838, r10844, MPFR_RNDN); ; ; mpfr_add(r10848, r10847, r10839, MPFR_RNDN); mpfr_add(r10849, r10839, r10840, MPFR_RNDN); mpfr_sub(r10850, r10848, r10849, MPFR_RNDN); mpfr_div(r10851, r10846, r10850, MPFR_RNDN); mpfr_add(r10852, r10845, r10851, MPFR_RNDN); ; ; mpfr_sub(r10855, r10839, r10854, MPFR_RNDN); mpfr_sub(r10856, r10841, r10855, MPFR_RNDN); mpfr_div(r10857, r10853, r10856, MPFR_RNDN); ; ; mpfr_add(r10860, r10839, r10859, MPFR_RNDN); mpfr_sub(r10861, r10860, r10849, MPFR_RNDN); mpfr_div(r10862, r10858, r10861, MPFR_RNDN); mpfr_add(r10863, r10857, r10862, MPFR_RNDN); mpfr_add(r10864, r10852, r10863, MPFR_RNDN); ; ; mpfr_add(r10867, r10839, r10866, MPFR_RNDN); mpfr_sub(r10868, r10867, r10849, MPFR_RNDN); mpfr_div(r10869, r10865, r10868, MPFR_RNDN); ; ; mpfr_add(r10872, r10839, r10871, MPFR_RNDN); mpfr_sub(r10873, r10872, r10849, MPFR_RNDN); mpfr_div(r10874, r10870, r10873, MPFR_RNDN); mpfr_add(r10875, r10869, r10874, MPFR_RNDN); ; ; mpfr_div(r10878, r10877, r10841, MPFR_RNDN); mpfr_add(r10879, r10876, r10878, MPFR_RNDN); ; ; mpfr_add(r10882, r10881, r10839, MPFR_RNDN); mpfr_sub(r10883, r10882, r10849, MPFR_RNDN); mpfr_div(r10884, r10880, r10883, MPFR_RNDN); mpfr_add(r10885, r10879, r10884, MPFR_RNDN); mpfr_add(r10886, r10875, r10885, MPFR_RNDN); mpfr_add(r10887, r10864, r10886, MPFR_RNDN); mpfr_const_pi(r10888, MPFR_RNDN); mpfr_mul(r10889, r10881, r10888, MPFR_RNDN); mpfr_sqrt(r10890, r10889, MPFR_RNDN); mpfr_mul(r10891, r10890, r10888, MPFR_RNDN); mpfr_mul(r10892, r10888, r10840, MPFR_RNDN); mpfr_sin(r10893, r10892, MPFR_RNDN); mpfr_div(r10894, r10891, r10893, MPFR_RNDN); ; mpfr_add(r10896, r10873, r10895, MPFR_RNDN); mpfr_sub(r10897, r10839, r10895, MPFR_RNDN); mpfr_sub(r10898, r10841, r10897, MPFR_RNDN); mpfr_pow(r10899, r10896, r10898, MPFR_RNDN); mpfr_mul(r10900, r10894, r10899, MPFR_RNDN); mpfr_exp(r10901, r10896, MPFR_RNDN); mpfr_div(r10902, r10900, r10901, MPFR_RNDN); mpfr_mul(r10903, r10887, r10902, MPFR_RNDN); return mpfr_get_d(r10903, MPFR_RNDN); }