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

char *name = "3.9.2 real part (p56)";

double f_if() {
        float r8616 = 1;
        float r8617 = -r8616;
        float r8618 = 2;
        float r8619 = r8617 / r8618;
        float r8620 = 3;
        float r8621 = sqrt(r8620);
        float r8622 = r8621 / r8618;
        float r8623 = /* ERROR: no complex support in C */;
        float r8624 = r8623 * r8623;
        float r8625 = r8624 * r8623;
        float r8626 = r8625 * r8623;
        float r8627 = -r8618;
        float r8628 = 0;
        float r8629 = /* ERROR: no complex support in C */;
        float r8630 = r8629 * r8623;
        float r8631 = r8630 * r8623;
        float r8632 = r8631 * r8623;
        float r8633 = r8626 + r8632;
        float r8634 = 5;
        float r8635 = /* ERROR: no complex support in C */;
        float r8636 = r8635 * r8623;
        float r8637 = r8636 * r8623;
        float r8638 = r8633 + r8637;
        float r8639 = 4;
        float r8640 = /* ERROR: no complex support in C */;
        float r8641 = r8640 * r8623;
        float r8642 = r8638 + r8641;
        float r8643 = 7;
        float r8644 = /* ERROR: no complex support in C */;
        float r8645 = r8642 + r8644;
        float r8646 = /* ERROR: no complex support in C */;
        return r8646;
}

double f_id() {
        double r8647 = 1;
        double r8648 = -r8647;
        double r8649 = 2;
        double r8650 = r8648 / r8649;
        double r8651 = 3;
        double r8652 = sqrt(r8651);
        double r8653 = r8652 / r8649;
        double r8654 = /* ERROR: no complex support in C */;
        double r8655 = r8654 * r8654;
        double r8656 = r8655 * r8654;
        double r8657 = r8656 * r8654;
        double r8658 = -r8649;
        double r8659 = 0;
        double r8660 = /* ERROR: no complex support in C */;
        double r8661 = r8660 * r8654;
        double r8662 = r8661 * r8654;
        double r8663 = r8662 * r8654;
        double r8664 = r8657 + r8663;
        double r8665 = 5;
        double r8666 = /* ERROR: no complex support in C */;
        double r8667 = r8666 * r8654;
        double r8668 = r8667 * r8654;
        double r8669 = r8664 + r8668;
        double r8670 = 4;
        double r8671 = /* ERROR: no complex support in C */;
        double r8672 = r8671 * r8654;
        double r8673 = r8669 + r8672;
        double r8674 = 7;
        double r8675 = /* ERROR: no complex support in C */;
        double r8676 = r8673 + r8675;
        double r8677 = /* ERROR: no complex support in C */;
        return r8677;
}


double f_of() {
        float r8678 = 1;
        float r8679 = -r8678;
        float r8680 = 2;
        float r8681 = r8679 / r8680;
        float r8682 = 3;
        float r8683 = sqrt(r8682);
        float r8684 = r8683 / r8680;
        float r8685 = /* ERROR: no complex support in C */;
        float r8686 = r8685 * r8685;
        float r8687 = r8686 * r8685;
        float r8688 = r8687 * r8685;
        float r8689 = -r8680;
        float r8690 = 0;
        float r8691 = /* ERROR: no complex support in C */;
        float r8692 = r8691 * r8685;
        float r8693 = r8692 * r8685;
        float r8694 = r8693 * r8685;
        float r8695 = r8688 + r8694;
        float r8696 = 5;
        float r8697 = /* ERROR: no complex support in C */;
        float r8698 = r8697 * r8685;
        float r8699 = r8698 * r8685;
        float r8700 = r8695 + r8699;
        float r8701 = 4;
        float r8702 = /* ERROR: no complex support in C */;
        float r8703 = r8702 * r8685;
        float r8704 = r8700 + r8703;
        float r8705 = 7;
        float r8706 = /* ERROR: no complex support in C */;
        float r8707 = r8704 + r8706;
        float r8708 = /* ERROR: no complex support in C */;
        return r8708;
}

double f_od() {
        double r8709 = 1;
        double r8710 = -r8709;
        double r8711 = 2;
        double r8712 = r8710 / r8711;
        double r8713 = 3;
        double r8714 = sqrt(r8713);
        double r8715 = r8714 / r8711;
        double r8716 = /* ERROR: no complex support in C */;
        double r8717 = r8716 * r8716;
        double r8718 = r8717 * r8716;
        double r8719 = r8718 * r8716;
        double r8720 = -r8711;
        double r8721 = 0;
        double r8722 = /* ERROR: no complex support in C */;
        double r8723 = r8722 * r8716;
        double r8724 = r8723 * r8716;
        double r8725 = r8724 * r8716;
        double r8726 = r8719 + r8725;
        double r8727 = 5;
        double r8728 = /* ERROR: no complex support in C */;
        double r8729 = r8728 * r8716;
        double r8730 = r8729 * r8716;
        double r8731 = r8726 + r8730;
        double r8732 = 4;
        double r8733 = /* ERROR: no complex support in C */;
        double r8734 = r8733 * r8716;
        double r8735 = r8731 + r8734;
        double r8736 = 7;
        double r8737 = /* ERROR: no complex support in C */;
        double r8738 = r8735 + r8737;
        double r8739 = /* ERROR: no complex support in C */;
        return r8739;
}

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 r8740, r8741, r8742, r8743, r8744, r8745, r8746, r8747, r8748, r8749, r8750, r8751, r8752, r8753, r8754, r8755, r8756, r8757, r8758, r8759, r8760, r8761, r8762, r8763, r8764, r8765, r8766, r8767, r8768, r8769, r8770;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8740, "1", 10, MPFR_RNDN);
        mpfr_init(r8741);
        mpfr_init_set_str(r8742, "2", 10, MPFR_RNDN);
        mpfr_init(r8743);
        mpfr_init_set_str(r8744, "3", 10, MPFR_RNDN);
        mpfr_init(r8745);
        mpfr_init(r8746);
        mpfr_init(r8747);
        mpfr_init(r8748);
        mpfr_init(r8749);
        mpfr_init(r8750);
        mpfr_init(r8751);
        mpfr_init_set_str(r8752, "0", 10, MPFR_RNDN);
        mpfr_init(r8753);
        mpfr_init(r8754);
        mpfr_init(r8755);
        mpfr_init(r8756);
        mpfr_init(r8757);
        mpfr_init_set_str(r8758, "5", 10, MPFR_RNDN);
        mpfr_init(r8759);
        mpfr_init(r8760);
        mpfr_init(r8761);
        mpfr_init(r8762);
        mpfr_init_set_str(r8763, "4", 10, MPFR_RNDN);
        mpfr_init(r8764);
        mpfr_init(r8765);
        mpfr_init(r8766);
        mpfr_init_set_str(r8767, "7", 10, MPFR_RNDN);
        mpfr_init(r8768);
        mpfr_init(r8769);
        mpfr_init(r8770);
}

double f_im() {
        ;
        mpfr_neg(r8741, r8740, MPFR_RNDN);
        ;
        mpfr_div(r8743, r8741, r8742, MPFR_RNDN);
        ;
        mpfr_sqrt(r8745, r8744, MPFR_RNDN);
        mpfr_div(r8746, r8745, r8742, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8748, r8747, r8747, MPFR_RNDN);
        mpfr_mul(r8749, r8748, r8747, MPFR_RNDN);
        mpfr_mul(r8750, r8749, r8747, MPFR_RNDN);
        mpfr_neg(r8751, r8742, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8754, r8753, r8747, MPFR_RNDN);
        mpfr_mul(r8755, r8754, r8747, MPFR_RNDN);
        mpfr_mul(r8756, r8755, r8747, MPFR_RNDN);
        mpfr_add(r8757, r8750, r8756, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8760, r8759, r8747, MPFR_RNDN);
        mpfr_mul(r8761, r8760, r8747, MPFR_RNDN);
        mpfr_add(r8762, r8757, r8761, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8765, r8764, r8747, MPFR_RNDN);
        mpfr_add(r8766, r8762, r8765, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8769, r8766, r8768, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8770, MPFR_RNDN);
}

static mpfr_t r8771, r8772, r8773, r8774, r8775, r8776, r8777, r8778, r8779, r8780, r8781, r8782, r8783, r8784, r8785, r8786, r8787, r8788, r8789, r8790, r8791, r8792, r8793, r8794, r8795, r8796, r8797, r8798, r8799, r8800, r8801;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8771, "1", 10, MPFR_RNDN);
        mpfr_init(r8772);
        mpfr_init_set_str(r8773, "2", 10, MPFR_RNDN);
        mpfr_init(r8774);
        mpfr_init_set_str(r8775, "3", 10, MPFR_RNDN);
        mpfr_init(r8776);
        mpfr_init(r8777);
        mpfr_init(r8778);
        mpfr_init(r8779);
        mpfr_init(r8780);
        mpfr_init(r8781);
        mpfr_init(r8782);
        mpfr_init_set_str(r8783, "0", 10, MPFR_RNDN);
        mpfr_init(r8784);
        mpfr_init(r8785);
        mpfr_init(r8786);
        mpfr_init(r8787);
        mpfr_init(r8788);
        mpfr_init_set_str(r8789, "5", 10, MPFR_RNDN);
        mpfr_init(r8790);
        mpfr_init(r8791);
        mpfr_init(r8792);
        mpfr_init(r8793);
        mpfr_init_set_str(r8794, "4", 10, MPFR_RNDN);
        mpfr_init(r8795);
        mpfr_init(r8796);
        mpfr_init(r8797);
        mpfr_init_set_str(r8798, "7", 10, MPFR_RNDN);
        mpfr_init(r8799);
        mpfr_init(r8800);
        mpfr_init(r8801);
}

double f_fm() {
        ;
        mpfr_neg(r8772, r8771, MPFR_RNDN);
        ;
        mpfr_div(r8774, r8772, r8773, MPFR_RNDN);
        ;
        mpfr_sqrt(r8776, r8775, MPFR_RNDN);
        mpfr_div(r8777, r8776, r8773, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8779, r8778, r8778, MPFR_RNDN);
        mpfr_mul(r8780, r8779, r8778, MPFR_RNDN);
        mpfr_mul(r8781, r8780, r8778, MPFR_RNDN);
        mpfr_neg(r8782, r8773, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8785, r8784, r8778, MPFR_RNDN);
        mpfr_mul(r8786, r8785, r8778, MPFR_RNDN);
        mpfr_mul(r8787, r8786, r8778, MPFR_RNDN);
        mpfr_add(r8788, r8781, r8787, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8791, r8790, r8778, MPFR_RNDN);
        mpfr_mul(r8792, r8791, r8778, MPFR_RNDN);
        mpfr_add(r8793, r8788, r8792, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8796, r8795, r8778, MPFR_RNDN);
        mpfr_add(r8797, r8793, r8796, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8800, r8797, r8799, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8801, MPFR_RNDN);
}

static mpfr_t r8802, r8803, r8804, r8805, r8806, r8807, r8808, r8809, r8810, r8811, r8812, r8813, r8814, r8815, r8816, r8817, r8818, r8819, r8820, r8821, r8822, r8823, r8824, r8825, r8826, r8827, r8828, r8829, r8830, r8831, r8832;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8802, "1", 10, MPFR_RNDN);
        mpfr_init(r8803);
        mpfr_init_set_str(r8804, "2", 10, MPFR_RNDN);
        mpfr_init(r8805);
        mpfr_init_set_str(r8806, "3", 10, MPFR_RNDN);
        mpfr_init(r8807);
        mpfr_init(r8808);
        mpfr_init(r8809);
        mpfr_init(r8810);
        mpfr_init(r8811);
        mpfr_init(r8812);
        mpfr_init(r8813);
        mpfr_init_set_str(r8814, "0", 10, MPFR_RNDN);
        mpfr_init(r8815);
        mpfr_init(r8816);
        mpfr_init(r8817);
        mpfr_init(r8818);
        mpfr_init(r8819);
        mpfr_init_set_str(r8820, "5", 10, MPFR_RNDN);
        mpfr_init(r8821);
        mpfr_init(r8822);
        mpfr_init(r8823);
        mpfr_init(r8824);
        mpfr_init_set_str(r8825, "4", 10, MPFR_RNDN);
        mpfr_init(r8826);
        mpfr_init(r8827);
        mpfr_init(r8828);
        mpfr_init_set_str(r8829, "7", 10, MPFR_RNDN);
        mpfr_init(r8830);
        mpfr_init(r8831);
        mpfr_init(r8832);
}

double f_dm() {
        ;
        mpfr_neg(r8803, r8802, MPFR_RNDN);
        ;
        mpfr_div(r8805, r8803, r8804, MPFR_RNDN);
        ;
        mpfr_sqrt(r8807, r8806, MPFR_RNDN);
        mpfr_div(r8808, r8807, r8804, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8810, r8809, r8809, MPFR_RNDN);
        mpfr_mul(r8811, r8810, r8809, MPFR_RNDN);
        mpfr_mul(r8812, r8811, r8809, MPFR_RNDN);
        mpfr_neg(r8813, r8804, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8816, r8815, r8809, MPFR_RNDN);
        mpfr_mul(r8817, r8816, r8809, MPFR_RNDN);
        mpfr_mul(r8818, r8817, r8809, MPFR_RNDN);
        mpfr_add(r8819, r8812, r8818, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8822, r8821, r8809, MPFR_RNDN);
        mpfr_mul(r8823, r8822, r8809, MPFR_RNDN);
        mpfr_add(r8824, r8819, r8823, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8827, r8826, r8809, MPFR_RNDN);
        mpfr_add(r8828, r8824, r8827, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8831, r8828, r8830, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8832, MPFR_RNDN);
}

