#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 r8652 = 1;
        float r8653 = -r8652;
        float r8654 = 2;
        float r8655 = r8653 / r8654;
        float r8656 = 3;
        float r8657 = sqrt(r8656);
        float r8658 = r8657 / r8654;
        float r8659 = /* ERROR: no complex support in C */;
        float r8660 = r8659 * r8659;
        float r8661 = r8660 * r8659;
        float r8662 = r8661 * r8659;
        float r8663 = -r8654;
        float r8664 = 0;
        float r8665 = /* ERROR: no complex support in C */;
        float r8666 = r8665 * r8659;
        float r8667 = r8666 * r8659;
        float r8668 = r8667 * r8659;
        float r8669 = r8662 + r8668;
        float r8670 = 5;
        float r8671 = /* ERROR: no complex support in C */;
        float r8672 = r8671 * r8659;
        float r8673 = r8672 * r8659;
        float r8674 = r8669 + r8673;
        float r8675 = 4;
        float r8676 = /* ERROR: no complex support in C */;
        float r8677 = r8676 * r8659;
        float r8678 = r8674 + r8677;
        float r8679 = 7;
        float r8680 = /* ERROR: no complex support in C */;
        float r8681 = r8678 + r8680;
        float r8682 = /* ERROR: no complex support in C */;
        return r8682;
}

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


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

double f_od() {
        double r8745 = 1;
        double r8746 = -r8745;
        double r8747 = 2;
        double r8748 = r8746 / r8747;
        double r8749 = 3;
        double r8750 = sqrt(r8749);
        double r8751 = r8750 / r8747;
        double r8752 = /* ERROR: no complex support in C */;
        double r8753 = r8752 * r8752;
        double r8754 = r8753 * r8752;
        double r8755 = r8754 * r8752;
        double r8756 = -r8747;
        double r8757 = 0;
        double r8758 = /* ERROR: no complex support in C */;
        double r8759 = r8758 * r8752;
        double r8760 = r8759 * r8752;
        double r8761 = r8760 * r8752;
        double r8762 = r8755 + r8761;
        double r8763 = 5;
        double r8764 = /* ERROR: no complex support in C */;
        double r8765 = r8764 * r8752;
        double r8766 = r8765 * r8752;
        double r8767 = r8762 + r8766;
        double r8768 = 4;
        double r8769 = /* ERROR: no complex support in C */;
        double r8770 = r8769 * r8752;
        double r8771 = r8767 + r8770;
        double r8772 = 7;
        double r8773 = /* ERROR: no complex support in C */;
        double r8774 = r8771 + r8773;
        double r8775 = /* ERROR: no complex support in C */;
        return r8775;
}

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 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, r8802, r8803, r8804, r8805, r8806;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8776, "1", 10, MPFR_RNDN);
        mpfr_init(r8777);
        mpfr_init_set_str(r8778, "2", 10, MPFR_RNDN);
        mpfr_init(r8779);
        mpfr_init_set_str(r8780, "3", 10, MPFR_RNDN);
        mpfr_init(r8781);
        mpfr_init(r8782);
        mpfr_init(r8783);
        mpfr_init(r8784);
        mpfr_init(r8785);
        mpfr_init(r8786);
        mpfr_init(r8787);
        mpfr_init_set_str(r8788, "0", 10, MPFR_RNDN);
        mpfr_init(r8789);
        mpfr_init(r8790);
        mpfr_init(r8791);
        mpfr_init(r8792);
        mpfr_init(r8793);
        mpfr_init_set_str(r8794, "5", 10, MPFR_RNDN);
        mpfr_init(r8795);
        mpfr_init(r8796);
        mpfr_init(r8797);
        mpfr_init(r8798);
        mpfr_init_set_str(r8799, "4", 10, MPFR_RNDN);
        mpfr_init(r8800);
        mpfr_init(r8801);
        mpfr_init(r8802);
        mpfr_init_set_str(r8803, "7", 10, MPFR_RNDN);
        mpfr_init(r8804);
        mpfr_init(r8805);
        mpfr_init(r8806);
}

double f_im() {
        ;
        mpfr_neg(r8777, r8776, MPFR_RNDN);
        ;
        mpfr_div(r8779, r8777, r8778, MPFR_RNDN);
        ;
        mpfr_sqrt(r8781, r8780, MPFR_RNDN);
        mpfr_div(r8782, r8781, r8778, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8784, r8783, r8783, MPFR_RNDN);
        mpfr_mul(r8785, r8784, r8783, MPFR_RNDN);
        mpfr_mul(r8786, r8785, r8783, MPFR_RNDN);
        mpfr_neg(r8787, r8778, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8790, r8789, r8783, MPFR_RNDN);
        mpfr_mul(r8791, r8790, r8783, MPFR_RNDN);
        mpfr_mul(r8792, r8791, r8783, MPFR_RNDN);
        mpfr_add(r8793, r8786, r8792, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8796, r8795, r8783, MPFR_RNDN);
        mpfr_mul(r8797, r8796, r8783, MPFR_RNDN);
        mpfr_add(r8798, r8793, r8797, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8801, r8800, r8783, MPFR_RNDN);
        mpfr_add(r8802, r8798, r8801, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8805, r8802, r8804, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8806, MPFR_RNDN);
}

static mpfr_t 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, r8833, r8834, r8835, r8836, r8837;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8807, "1", 10, MPFR_RNDN);
        mpfr_init(r8808);
        mpfr_init_set_str(r8809, "2", 10, MPFR_RNDN);
        mpfr_init(r8810);
        mpfr_init_set_str(r8811, "3", 10, MPFR_RNDN);
        mpfr_init(r8812);
        mpfr_init(r8813);
        mpfr_init(r8814);
        mpfr_init(r8815);
        mpfr_init(r8816);
        mpfr_init(r8817);
        mpfr_init(r8818);
        mpfr_init_set_str(r8819, "0", 10, MPFR_RNDN);
        mpfr_init(r8820);
        mpfr_init(r8821);
        mpfr_init(r8822);
        mpfr_init(r8823);
        mpfr_init(r8824);
        mpfr_init_set_str(r8825, "5", 10, MPFR_RNDN);
        mpfr_init(r8826);
        mpfr_init(r8827);
        mpfr_init(r8828);
        mpfr_init(r8829);
        mpfr_init_set_str(r8830, "4", 10, MPFR_RNDN);
        mpfr_init(r8831);
        mpfr_init(r8832);
        mpfr_init(r8833);
        mpfr_init_set_str(r8834, "7", 10, MPFR_RNDN);
        mpfr_init(r8835);
        mpfr_init(r8836);
        mpfr_init(r8837);
}

double f_fm() {
        ;
        mpfr_neg(r8808, r8807, MPFR_RNDN);
        ;
        mpfr_div(r8810, r8808, r8809, MPFR_RNDN);
        ;
        mpfr_sqrt(r8812, r8811, MPFR_RNDN);
        mpfr_div(r8813, r8812, r8809, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8815, r8814, r8814, MPFR_RNDN);
        mpfr_mul(r8816, r8815, r8814, MPFR_RNDN);
        mpfr_mul(r8817, r8816, r8814, MPFR_RNDN);
        mpfr_neg(r8818, r8809, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8821, r8820, r8814, MPFR_RNDN);
        mpfr_mul(r8822, r8821, r8814, MPFR_RNDN);
        mpfr_mul(r8823, r8822, r8814, MPFR_RNDN);
        mpfr_add(r8824, r8817, r8823, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8827, r8826, r8814, MPFR_RNDN);
        mpfr_mul(r8828, r8827, r8814, MPFR_RNDN);
        mpfr_add(r8829, r8824, r8828, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8832, r8831, r8814, MPFR_RNDN);
        mpfr_add(r8833, r8829, r8832, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8836, r8833, r8835, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8837, MPFR_RNDN);
}

static mpfr_t r8838, r8839, r8840, r8841, r8842, r8843, r8844, r8845, r8846, r8847, r8848, r8849, r8850, r8851, r8852, r8853, r8854, r8855, r8856, r8857, r8858, r8859, r8860, r8861, r8862, r8863, r8864, r8865, r8866, r8867, r8868;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r8838, "1", 10, MPFR_RNDN);
        mpfr_init(r8839);
        mpfr_init_set_str(r8840, "2", 10, MPFR_RNDN);
        mpfr_init(r8841);
        mpfr_init_set_str(r8842, "3", 10, MPFR_RNDN);
        mpfr_init(r8843);
        mpfr_init(r8844);
        mpfr_init(r8845);
        mpfr_init(r8846);
        mpfr_init(r8847);
        mpfr_init(r8848);
        mpfr_init(r8849);
        mpfr_init_set_str(r8850, "0", 10, MPFR_RNDN);
        mpfr_init(r8851);
        mpfr_init(r8852);
        mpfr_init(r8853);
        mpfr_init(r8854);
        mpfr_init(r8855);
        mpfr_init_set_str(r8856, "5", 10, MPFR_RNDN);
        mpfr_init(r8857);
        mpfr_init(r8858);
        mpfr_init(r8859);
        mpfr_init(r8860);
        mpfr_init_set_str(r8861, "4", 10, MPFR_RNDN);
        mpfr_init(r8862);
        mpfr_init(r8863);
        mpfr_init(r8864);
        mpfr_init_set_str(r8865, "7", 10, MPFR_RNDN);
        mpfr_init(r8866);
        mpfr_init(r8867);
        mpfr_init(r8868);
}

double f_dm() {
        ;
        mpfr_neg(r8839, r8838, MPFR_RNDN);
        ;
        mpfr_div(r8841, r8839, r8840, MPFR_RNDN);
        ;
        mpfr_sqrt(r8843, r8842, MPFR_RNDN);
        mpfr_div(r8844, r8843, r8840, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        mpfr_mul(r8846, r8845, r8845, MPFR_RNDN);
        mpfr_mul(r8847, r8846, r8845, MPFR_RNDN);
        mpfr_mul(r8848, r8847, r8845, MPFR_RNDN);
        mpfr_neg(r8849, r8840, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8852, r8851, r8845, MPFR_RNDN);
        mpfr_mul(r8853, r8852, r8845, MPFR_RNDN);
        mpfr_mul(r8854, r8853, r8845, MPFR_RNDN);
        mpfr_add(r8855, r8848, r8854, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8858, r8857, r8845, MPFR_RNDN);
        mpfr_mul(r8859, r8858, r8845, MPFR_RNDN);
        mpfr_add(r8860, r8855, r8859, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_mul(r8863, r8862, r8845, MPFR_RNDN);
        mpfr_add(r8864, r8860, r8863, MPFR_RNDN);
        ;
        /* ERROR: no complex support in C */;
        mpfr_add(r8867, r8864, r8866, MPFR_RNDN);
        /* ERROR: no complex support in C */;
        return mpfr_get_d(r8868, MPFR_RNDN);
}

