#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Jmat.Real.gamma, branch z less than 0.5";

double f_if(float z) {
        float r10513 = atan2(1.0, 0.0);
        float r10514 = z;
        float r10515 = r10513 * r10514;
        float r10516 = sin(r10515);
        float r10517 = r10513 / r10516;
        float r10518 = 2;
        float r10519 = r10513 * r10518;
        float r10520 = sqrt(r10519);
        float r10521 = 1;
        float r10522 = r10521 - r10514;
        float r10523 = r10522 - r10521;
        float r10524 = 7;
        float r10525 = r10523 + r10524;
        float r10526 = 0.5;
        float r10527 = r10525 + r10526;
        float r10528 = r10523 + r10526;
        float r10529 = pow(r10527, r10528);
        float r10530 = r10520 * r10529;
        float r10531 = -r10527;
        float r10532 = exp(r10531);
        float r10533 = r10530 * r10532;
        float r10534 = 0.9999999999998099;
        float r10535 = 676.5203681218851;
        float r10536 = r10523 + r10521;
        float r10537 = r10535 / r10536;
        float r10538 = r10534 + r10537;
        float r10539 = -1259.1392167224028;
        float r10540 = r10523 + r10518;
        float r10541 = r10539 / r10540;
        float r10542 = r10538 + r10541;
        float r10543 = 771.3234287776531;
        float r10544 = 3;
        float r10545 = r10523 + r10544;
        float r10546 = r10543 / r10545;
        float r10547 = r10542 + r10546;
        float r10548 = -176.6150291621406;
        float r10549 = 4;
        float r10550 = r10523 + r10549;
        float r10551 = r10548 / r10550;
        float r10552 = r10547 + r10551;
        float r10553 = 12.507343278686905;
        float r10554 = 5;
        float r10555 = r10523 + r10554;
        float r10556 = r10553 / r10555;
        float r10557 = r10552 + r10556;
        float r10558 = -0.13857109526572012;
        float r10559 = 6;
        float r10560 = r10523 + r10559;
        float r10561 = r10558 / r10560;
        float r10562 = r10557 + r10561;
        float r10563 = 9.984369578019572e-06;
        float r10564 = r10563 / r10525;
        float r10565 = r10562 + r10564;
        float r10566 = 1.5056327351493116e-07;
        float r10567 = 8;
        float r10568 = r10523 + r10567;
        float r10569 = r10566 / r10568;
        float r10570 = r10565 + r10569;
        float r10571 = r10533 * r10570;
        float r10572 = r10517 * r10571;
        return r10572;
}

double f_id(double z) {
        double r10573 = atan2(1.0, 0.0);
        double r10574 = z;
        double r10575 = r10573 * r10574;
        double r10576 = sin(r10575);
        double r10577 = r10573 / r10576;
        double r10578 = 2;
        double r10579 = r10573 * r10578;
        double r10580 = sqrt(r10579);
        double r10581 = 1;
        double r10582 = r10581 - r10574;
        double r10583 = r10582 - r10581;
        double r10584 = 7;
        double r10585 = r10583 + r10584;
        double r10586 = 0.5;
        double r10587 = r10585 + r10586;
        double r10588 = r10583 + r10586;
        double r10589 = pow(r10587, r10588);
        double r10590 = r10580 * r10589;
        double r10591 = -r10587;
        double r10592 = exp(r10591);
        double r10593 = r10590 * r10592;
        double r10594 = 0.9999999999998099;
        double r10595 = 676.5203681218851;
        double r10596 = r10583 + r10581;
        double r10597 = r10595 / r10596;
        double r10598 = r10594 + r10597;
        double r10599 = -1259.1392167224028;
        double r10600 = r10583 + r10578;
        double r10601 = r10599 / r10600;
        double r10602 = r10598 + r10601;
        double r10603 = 771.3234287776531;
        double r10604 = 3;
        double r10605 = r10583 + r10604;
        double r10606 = r10603 / r10605;
        double r10607 = r10602 + r10606;
        double r10608 = -176.6150291621406;
        double r10609 = 4;
        double r10610 = r10583 + r10609;
        double r10611 = r10608 / r10610;
        double r10612 = r10607 + r10611;
        double r10613 = 12.507343278686905;
        double r10614 = 5;
        double r10615 = r10583 + r10614;
        double r10616 = r10613 / r10615;
        double r10617 = r10612 + r10616;
        double r10618 = -0.13857109526572012;
        double r10619 = 6;
        double r10620 = r10583 + r10619;
        double r10621 = r10618 / r10620;
        double r10622 = r10617 + r10621;
        double r10623 = 9.984369578019572e-06;
        double r10624 = r10623 / r10585;
        double r10625 = r10622 + r10624;
        double r10626 = 1.5056327351493116e-07;
        double r10627 = 8;
        double r10628 = r10583 + r10627;
        double r10629 = r10626 / r10628;
        double r10630 = r10625 + r10629;
        double r10631 = r10593 * r10630;
        double r10632 = r10577 * r10631;
        return r10632;
}


double f_of(float z) {
        float r10633 = 7.5;
        float r10634 = z;
        float r10635 = r10633 - r10634;
        float r10636 = 0.5;
        float r10637 = r10636 - r10634;
        float r10638 = pow(r10635, r10637);
        float r10639 = exp(r10635);
        float r10640 = r10638 / r10639;
        float r10641 = 2;
        float r10642 = atan2(1.0, 0.0);
        float r10643 = r10641 * r10642;
        float r10644 = sqrt(r10643);
        float r10645 = r10640 * r10644;
        float r10646 = r10642 * r10634;
        float r10647 = sin(r10646);
        float r10648 = r10642 / r10647;
        float r10649 = r10645 * r10648;
        float r10650 = -176.6150291621406;
        float r10651 = 1;
        float r10652 = r10651 - r10634;
        float r10653 = 4;
        float r10654 = r10651 - r10653;
        float r10655 = r10652 - r10654;
        float r10656 = r10650 / r10655;
        float r10657 = 771.3234287776531;
        float r10658 = 3;
        float r10659 = r10651 - r10658;
        float r10660 = r10652 - r10659;
        float r10661 = r10657 / r10660;
        float r10662 = r10656 + r10661;
        float r10663 = 676.5203681218851;
        float r10664 = r10663 / r10652;
        float r10665 = 0.9999999999998099;
        float r10666 = r10664 + r10665;
        float r10667 = -1259.1392167224028;
        float r10668 = r10641 + r10651;
        float r10669 = r10651 + r10634;
        float r10670 = r10668 - r10669;
        float r10671 = r10667 / r10670;
        float r10672 = r10666 + r10671;
        float r10673 = r10662 + r10672;
        float r10674 = 12.507343278686905;
        float r10675 = 5;
        float r10676 = r10651 - r10675;
        float r10677 = r10652 - r10676;
        float r10678 = r10674 / r10677;
        float r10679 = -0.13857109526572012;
        float r10680 = 6;
        float r10681 = r10651 + r10680;
        float r10682 = r10681 - r10669;
        float r10683 = r10679 / r10682;
        float r10684 = r10678 + r10683;
        float r10685 = 1.5056327351493116e-07;
        float r10686 = 8;
        float r10687 = r10651 + r10686;
        float r10688 = r10634 + r10651;
        float r10689 = r10687 - r10688;
        float r10690 = r10685 / r10689;
        float r10691 = 9.984369578019572e-06;
        float r10692 = 7;
        float r10693 = r10692 + r10651;
        float r10694 = r10693 - r10688;
        float r10695 = r10691 / r10694;
        float r10696 = r10690 + r10695;
        float r10697 = r10684 + r10696;
        float r10698 = r10673 + r10697;
        float r10699 = r10649 * r10698;
        return r10699;
}

double f_od(double z) {
        double r10700 = 7.5;
        double r10701 = z;
        double r10702 = r10700 - r10701;
        double r10703 = 0.5;
        double r10704 = r10703 - r10701;
        double r10705 = pow(r10702, r10704);
        double r10706 = exp(r10702);
        double r10707 = r10705 / r10706;
        double r10708 = 2;
        double r10709 = atan2(1.0, 0.0);
        double r10710 = r10708 * r10709;
        double r10711 = sqrt(r10710);
        double r10712 = r10707 * r10711;
        double r10713 = r10709 * r10701;
        double r10714 = sin(r10713);
        double r10715 = r10709 / r10714;
        double r10716 = r10712 * r10715;
        double r10717 = -176.6150291621406;
        double r10718 = 1;
        double r10719 = r10718 - r10701;
        double r10720 = 4;
        double r10721 = r10718 - r10720;
        double r10722 = r10719 - r10721;
        double r10723 = r10717 / r10722;
        double r10724 = 771.3234287776531;
        double r10725 = 3;
        double r10726 = r10718 - r10725;
        double r10727 = r10719 - r10726;
        double r10728 = r10724 / r10727;
        double r10729 = r10723 + r10728;
        double r10730 = 676.5203681218851;
        double r10731 = r10730 / r10719;
        double r10732 = 0.9999999999998099;
        double r10733 = r10731 + r10732;
        double r10734 = -1259.1392167224028;
        double r10735 = r10708 + r10718;
        double r10736 = r10718 + r10701;
        double r10737 = r10735 - r10736;
        double r10738 = r10734 / r10737;
        double r10739 = r10733 + r10738;
        double r10740 = r10729 + r10739;
        double r10741 = 12.507343278686905;
        double r10742 = 5;
        double r10743 = r10718 - r10742;
        double r10744 = r10719 - r10743;
        double r10745 = r10741 / r10744;
        double r10746 = -0.13857109526572012;
        double r10747 = 6;
        double r10748 = r10718 + r10747;
        double r10749 = r10748 - r10736;
        double r10750 = r10746 / r10749;
        double r10751 = r10745 + r10750;
        double r10752 = 1.5056327351493116e-07;
        double r10753 = 8;
        double r10754 = r10718 + r10753;
        double r10755 = r10701 + r10718;
        double r10756 = r10754 - r10755;
        double r10757 = r10752 / r10756;
        double r10758 = 9.984369578019572e-06;
        double r10759 = 7;
        double r10760 = r10759 + r10718;
        double r10761 = r10760 - r10755;
        double r10762 = r10758 / r10761;
        double r10763 = r10757 + r10762;
        double r10764 = r10751 + r10763;
        double r10765 = r10740 + r10764;
        double r10766 = r10716 * r10765;
        return r10766;
}

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 r10767, r10768, r10769, r10770, r10771, 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;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r10767);
        mpfr_init(r10768);
        mpfr_init(r10769);
        mpfr_init(r10770);
        mpfr_init(r10771);
        mpfr_init_set_str(r10772, "2", 10, MPFR_RNDN);
        mpfr_init(r10773);
        mpfr_init(r10774);
        mpfr_init_set_str(r10775, "1", 10, MPFR_RNDN);
        mpfr_init(r10776);
        mpfr_init(r10777);
        mpfr_init_set_str(r10778, "7", 10, MPFR_RNDN);
        mpfr_init(r10779);
        mpfr_init_set_str(r10780, "0.5", 10, MPFR_RNDN);
        mpfr_init(r10781);
        mpfr_init(r10782);
        mpfr_init(r10783);
        mpfr_init(r10784);
        mpfr_init(r10785);
        mpfr_init(r10786);
        mpfr_init(r10787);
        mpfr_init_set_str(r10788, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init_set_str(r10789, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init(r10790);
        mpfr_init(r10791);
        mpfr_init(r10792);
        mpfr_init_set_str(r10793, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init(r10794);
        mpfr_init(r10795);
        mpfr_init(r10796);
        mpfr_init_set_str(r10797, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r10798, "3", 10, MPFR_RNDN);
        mpfr_init(r10799);
        mpfr_init(r10800);
        mpfr_init(r10801);
        mpfr_init_set_str(r10802, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r10803, "4", 10, MPFR_RNDN);
        mpfr_init(r10804);
        mpfr_init(r10805);
        mpfr_init(r10806);
        mpfr_init_set_str(r10807, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r10808, "5", 10, MPFR_RNDN);
        mpfr_init(r10809);
        mpfr_init(r10810);
        mpfr_init(r10811);
        mpfr_init_set_str(r10812, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r10813, "6", 10, MPFR_RNDN);
        mpfr_init(r10814);
        mpfr_init(r10815);
        mpfr_init(r10816);
        mpfr_init_set_str(r10817, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r10818);
        mpfr_init(r10819);
        mpfr_init_set_str(r10820, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r10821, "8", 10, MPFR_RNDN);
        mpfr_init(r10822);
        mpfr_init(r10823);
        mpfr_init(r10824);
        mpfr_init(r10825);
        mpfr_init(r10826);
}

double f_im(double z) {
        mpfr_const_pi(r10767, MPFR_RNDN);
        mpfr_set_d(r10768, z, MPFR_RNDN);
        mpfr_mul(r10769, r10767, r10768, MPFR_RNDN);
        mpfr_sin(r10770, r10769, MPFR_RNDN);
        mpfr_div(r10771, r10767, r10770, MPFR_RNDN);
        ;
        mpfr_mul(r10773, r10767, r10772, MPFR_RNDN);
        mpfr_sqrt(r10774, r10773, MPFR_RNDN);
        ;
        mpfr_sub(r10776, r10775, r10768, MPFR_RNDN);
        mpfr_sub(r10777, r10776, r10775, MPFR_RNDN);
        ;
        mpfr_add(r10779, r10777, r10778, MPFR_RNDN);
        ;
        mpfr_add(r10781, r10779, r10780, MPFR_RNDN);
        mpfr_add(r10782, r10777, r10780, MPFR_RNDN);
        mpfr_pow(r10783, r10781, r10782, MPFR_RNDN);
        mpfr_mul(r10784, r10774, r10783, MPFR_RNDN);
        mpfr_neg(r10785, r10781, MPFR_RNDN);
        mpfr_exp(r10786, r10785, MPFR_RNDN);
        mpfr_mul(r10787, r10784, r10786, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10790, r10777, r10775, MPFR_RNDN);
        mpfr_div(r10791, r10789, r10790, MPFR_RNDN);
        mpfr_add(r10792, r10788, r10791, MPFR_RNDN);
        ;
        mpfr_add(r10794, r10777, r10772, MPFR_RNDN);
        mpfr_div(r10795, r10793, r10794, MPFR_RNDN);
        mpfr_add(r10796, r10792, r10795, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10799, r10777, r10798, MPFR_RNDN);
        mpfr_div(r10800, r10797, r10799, MPFR_RNDN);
        mpfr_add(r10801, r10796, r10800, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10804, r10777, r10803, MPFR_RNDN);
        mpfr_div(r10805, r10802, r10804, MPFR_RNDN);
        mpfr_add(r10806, r10801, r10805, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10809, r10777, r10808, MPFR_RNDN);
        mpfr_div(r10810, r10807, r10809, MPFR_RNDN);
        mpfr_add(r10811, r10806, r10810, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10814, r10777, r10813, MPFR_RNDN);
        mpfr_div(r10815, r10812, r10814, MPFR_RNDN);
        mpfr_add(r10816, r10811, r10815, MPFR_RNDN);
        ;
        mpfr_div(r10818, r10817, r10779, MPFR_RNDN);
        mpfr_add(r10819, r10816, r10818, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10822, r10777, r10821, MPFR_RNDN);
        mpfr_div(r10823, r10820, r10822, MPFR_RNDN);
        mpfr_add(r10824, r10819, r10823, MPFR_RNDN);
        mpfr_mul(r10825, r10787, r10824, MPFR_RNDN);
        mpfr_mul(r10826, r10771, r10825, MPFR_RNDN);
        return mpfr_get_d(r10826, MPFR_RNDN);
}

static mpfr_t r10827, r10828, r10829, r10830, r10831, r10832, r10833, r10834, r10835, r10836, r10837, 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;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10827, "7.5", 10, MPFR_RNDN);
        mpfr_init(r10828);
        mpfr_init(r10829);
        mpfr_init_set_str(r10830, "0.5", 10, MPFR_RNDN);
        mpfr_init(r10831);
        mpfr_init(r10832);
        mpfr_init(r10833);
        mpfr_init(r10834);
        mpfr_init_set_str(r10835, "2", 10, MPFR_RNDN);
        mpfr_init(r10836);
        mpfr_init(r10837);
        mpfr_init(r10838);
        mpfr_init(r10839);
        mpfr_init(r10840);
        mpfr_init(r10841);
        mpfr_init(r10842);
        mpfr_init(r10843);
        mpfr_init_set_str(r10844, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r10845, "1", 10, MPFR_RNDN);
        mpfr_init(r10846);
        mpfr_init_set_str(r10847, "4", 10, MPFR_RNDN);
        mpfr_init(r10848);
        mpfr_init(r10849);
        mpfr_init(r10850);
        mpfr_init_set_str(r10851, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r10852, "3", 10, MPFR_RNDN);
        mpfr_init(r10853);
        mpfr_init(r10854);
        mpfr_init(r10855);
        mpfr_init(r10856);
        mpfr_init_set_str(r10857, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init(r10858);
        mpfr_init_set_str(r10859, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r10860);
        mpfr_init_set_str(r10861, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init(r10862);
        mpfr_init(r10863);
        mpfr_init(r10864);
        mpfr_init(r10865);
        mpfr_init(r10866);
        mpfr_init(r10867);
        mpfr_init_set_str(r10868, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r10869, "5", 10, MPFR_RNDN);
        mpfr_init(r10870);
        mpfr_init(r10871);
        mpfr_init(r10872);
        mpfr_init_set_str(r10873, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r10874, "6", 10, MPFR_RNDN);
        mpfr_init(r10875);
        mpfr_init(r10876);
        mpfr_init(r10877);
        mpfr_init(r10878);
        mpfr_init_set_str(r10879, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r10880, "8", 10, MPFR_RNDN);
        mpfr_init(r10881);
        mpfr_init(r10882);
        mpfr_init(r10883);
        mpfr_init(r10884);
        mpfr_init_set_str(r10885, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init_set_str(r10886, "7", 10, MPFR_RNDN);
        mpfr_init(r10887);
        mpfr_init(r10888);
        mpfr_init(r10889);
        mpfr_init(r10890);
        mpfr_init(r10891);
        mpfr_init(r10892);
        mpfr_init(r10893);
}

double f_fm(double z) {
        ;
        mpfr_set_d(r10828, z, MPFR_RNDN);
        mpfr_sub(r10829, r10827, r10828, MPFR_RNDN);
        ;
        mpfr_sub(r10831, r10830, r10828, MPFR_RNDN);
        mpfr_pow(r10832, r10829, r10831, MPFR_RNDN);
        mpfr_exp(r10833, r10829, MPFR_RNDN);
        mpfr_div(r10834, r10832, r10833, MPFR_RNDN);
        ;
        mpfr_const_pi(r10836, MPFR_RNDN);
        mpfr_mul(r10837, r10835, r10836, MPFR_RNDN);
        mpfr_sqrt(r10838, r10837, MPFR_RNDN);
        mpfr_mul(r10839, r10834, r10838, MPFR_RNDN);
        mpfr_mul(r10840, r10836, r10828, MPFR_RNDN);
        mpfr_sin(r10841, r10840, MPFR_RNDN);
        mpfr_div(r10842, r10836, r10841, MPFR_RNDN);
        mpfr_mul(r10843, r10839, r10842, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10846, r10845, r10828, MPFR_RNDN);
        ;
        mpfr_sub(r10848, r10845, r10847, MPFR_RNDN);
        mpfr_sub(r10849, r10846, r10848, MPFR_RNDN);
        mpfr_div(r10850, r10844, r10849, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10853, r10845, r10852, MPFR_RNDN);
        mpfr_sub(r10854, r10846, r10853, MPFR_RNDN);
        mpfr_div(r10855, r10851, r10854, MPFR_RNDN);
        mpfr_add(r10856, r10850, r10855, MPFR_RNDN);
        ;
        mpfr_div(r10858, r10857, r10846, MPFR_RNDN);
        ;
        mpfr_add(r10860, r10858, r10859, MPFR_RNDN);
        ;
        mpfr_add(r10862, r10835, r10845, MPFR_RNDN);
        mpfr_add(r10863, r10845, r10828, MPFR_RNDN);
        mpfr_sub(r10864, r10862, r10863, MPFR_RNDN);
        mpfr_div(r10865, r10861, r10864, MPFR_RNDN);
        mpfr_add(r10866, r10860, r10865, MPFR_RNDN);
        mpfr_add(r10867, r10856, r10866, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10870, r10845, r10869, MPFR_RNDN);
        mpfr_sub(r10871, r10846, r10870, MPFR_RNDN);
        mpfr_div(r10872, r10868, r10871, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10875, r10845, r10874, MPFR_RNDN);
        mpfr_sub(r10876, r10875, r10863, MPFR_RNDN);
        mpfr_div(r10877, r10873, r10876, MPFR_RNDN);
        mpfr_add(r10878, r10872, r10877, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10881, r10845, r10880, MPFR_RNDN);
        mpfr_add(r10882, r10828, r10845, MPFR_RNDN);
        mpfr_sub(r10883, r10881, r10882, MPFR_RNDN);
        mpfr_div(r10884, r10879, r10883, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10887, r10886, r10845, MPFR_RNDN);
        mpfr_sub(r10888, r10887, r10882, MPFR_RNDN);
        mpfr_div(r10889, r10885, r10888, MPFR_RNDN);
        mpfr_add(r10890, r10884, r10889, MPFR_RNDN);
        mpfr_add(r10891, r10878, r10890, MPFR_RNDN);
        mpfr_add(r10892, r10867, r10891, MPFR_RNDN);
        mpfr_mul(r10893, r10843, r10892, MPFR_RNDN);
        return mpfr_get_d(r10893, MPFR_RNDN);
}

static mpfr_t r10894, r10895, r10896, r10897, r10898, r10899, r10900, r10901, r10902, r10903, r10904, r10905, r10906, r10907, r10908, r10909, r10910, r10911, r10912, r10913, r10914, r10915, r10916, r10917, r10918, r10919, r10920, r10921, r10922, r10923, r10924, r10925, r10926, r10927, r10928, r10929, r10930, r10931, r10932, r10933, r10934, r10935, r10936, r10937, r10938, r10939, r10940, r10941, r10942, r10943, r10944, r10945, r10946, r10947, r10948, r10949, r10950, r10951, r10952, r10953, r10954, r10955, r10956, r10957, r10958, r10959, r10960;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r10894, "7.5", 10, MPFR_RNDN);
        mpfr_init(r10895);
        mpfr_init(r10896);
        mpfr_init_set_str(r10897, "0.5", 10, MPFR_RNDN);
        mpfr_init(r10898);
        mpfr_init(r10899);
        mpfr_init(r10900);
        mpfr_init(r10901);
        mpfr_init_set_str(r10902, "2", 10, MPFR_RNDN);
        mpfr_init(r10903);
        mpfr_init(r10904);
        mpfr_init(r10905);
        mpfr_init(r10906);
        mpfr_init(r10907);
        mpfr_init(r10908);
        mpfr_init(r10909);
        mpfr_init(r10910);
        mpfr_init_set_str(r10911, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r10912, "1", 10, MPFR_RNDN);
        mpfr_init(r10913);
        mpfr_init_set_str(r10914, "4", 10, MPFR_RNDN);
        mpfr_init(r10915);
        mpfr_init(r10916);
        mpfr_init(r10917);
        mpfr_init_set_str(r10918, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r10919, "3", 10, MPFR_RNDN);
        mpfr_init(r10920);
        mpfr_init(r10921);
        mpfr_init(r10922);
        mpfr_init(r10923);
        mpfr_init_set_str(r10924, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init(r10925);
        mpfr_init_set_str(r10926, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r10927);
        mpfr_init_set_str(r10928, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init(r10929);
        mpfr_init(r10930);
        mpfr_init(r10931);
        mpfr_init(r10932);
        mpfr_init(r10933);
        mpfr_init(r10934);
        mpfr_init_set_str(r10935, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r10936, "5", 10, MPFR_RNDN);
        mpfr_init(r10937);
        mpfr_init(r10938);
        mpfr_init(r10939);
        mpfr_init_set_str(r10940, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r10941, "6", 10, MPFR_RNDN);
        mpfr_init(r10942);
        mpfr_init(r10943);
        mpfr_init(r10944);
        mpfr_init(r10945);
        mpfr_init_set_str(r10946, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r10947, "8", 10, MPFR_RNDN);
        mpfr_init(r10948);
        mpfr_init(r10949);
        mpfr_init(r10950);
        mpfr_init(r10951);
        mpfr_init_set_str(r10952, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init_set_str(r10953, "7", 10, MPFR_RNDN);
        mpfr_init(r10954);
        mpfr_init(r10955);
        mpfr_init(r10956);
        mpfr_init(r10957);
        mpfr_init(r10958);
        mpfr_init(r10959);
        mpfr_init(r10960);
}

double f_dm(double z) {
        ;
        mpfr_set_d(r10895, z, MPFR_RNDN);
        mpfr_sub(r10896, r10894, r10895, MPFR_RNDN);
        ;
        mpfr_sub(r10898, r10897, r10895, MPFR_RNDN);
        mpfr_pow(r10899, r10896, r10898, MPFR_RNDN);
        mpfr_exp(r10900, r10896, MPFR_RNDN);
        mpfr_div(r10901, r10899, r10900, MPFR_RNDN);
        ;
        mpfr_const_pi(r10903, MPFR_RNDN);
        mpfr_mul(r10904, r10902, r10903, MPFR_RNDN);
        mpfr_sqrt(r10905, r10904, MPFR_RNDN);
        mpfr_mul(r10906, r10901, r10905, MPFR_RNDN);
        mpfr_mul(r10907, r10903, r10895, MPFR_RNDN);
        mpfr_sin(r10908, r10907, MPFR_RNDN);
        mpfr_div(r10909, r10903, r10908, MPFR_RNDN);
        mpfr_mul(r10910, r10906, r10909, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10913, r10912, r10895, MPFR_RNDN);
        ;
        mpfr_sub(r10915, r10912, r10914, MPFR_RNDN);
        mpfr_sub(r10916, r10913, r10915, MPFR_RNDN);
        mpfr_div(r10917, r10911, r10916, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10920, r10912, r10919, MPFR_RNDN);
        mpfr_sub(r10921, r10913, r10920, MPFR_RNDN);
        mpfr_div(r10922, r10918, r10921, MPFR_RNDN);
        mpfr_add(r10923, r10917, r10922, MPFR_RNDN);
        ;
        mpfr_div(r10925, r10924, r10913, MPFR_RNDN);
        ;
        mpfr_add(r10927, r10925, r10926, MPFR_RNDN);
        ;
        mpfr_add(r10929, r10902, r10912, MPFR_RNDN);
        mpfr_add(r10930, r10912, r10895, MPFR_RNDN);
        mpfr_sub(r10931, r10929, r10930, MPFR_RNDN);
        mpfr_div(r10932, r10928, r10931, MPFR_RNDN);
        mpfr_add(r10933, r10927, r10932, MPFR_RNDN);
        mpfr_add(r10934, r10923, r10933, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r10937, r10912, r10936, MPFR_RNDN);
        mpfr_sub(r10938, r10913, r10937, MPFR_RNDN);
        mpfr_div(r10939, r10935, r10938, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10942, r10912, r10941, MPFR_RNDN);
        mpfr_sub(r10943, r10942, r10930, MPFR_RNDN);
        mpfr_div(r10944, r10940, r10943, MPFR_RNDN);
        mpfr_add(r10945, r10939, r10944, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10948, r10912, r10947, MPFR_RNDN);
        mpfr_add(r10949, r10895, r10912, MPFR_RNDN);
        mpfr_sub(r10950, r10948, r10949, MPFR_RNDN);
        mpfr_div(r10951, r10946, r10950, MPFR_RNDN);
        ;
        ;
        mpfr_add(r10954, r10953, r10912, MPFR_RNDN);
        mpfr_sub(r10955, r10954, r10949, MPFR_RNDN);
        mpfr_div(r10956, r10952, r10955, MPFR_RNDN);
        mpfr_add(r10957, r10951, r10956, MPFR_RNDN);
        mpfr_add(r10958, r10945, r10957, MPFR_RNDN);
        mpfr_add(r10959, r10934, r10958, MPFR_RNDN);
        mpfr_mul(r10960, r10910, r10959, MPFR_RNDN);
        return mpfr_get_d(r10960, MPFR_RNDN);
}

