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

char *name = "3frac (problem 3.3.3)";

double f_if(float x) {
        float r24728 = 1;
        float r24729 = x;
        float r24730 = r24729 + r24728;
        float r24731 = r24728 / r24730;
        float r24732 = 2;
        float r24733 = r24732 / r24729;
        float r24734 = r24731 - r24733;
        float r24735 = r24729 - r24728;
        float r24736 = r24728 / r24735;
        float r24737 = r24734 + r24736;
        return r24737;
}

double f_id(double x) {
        double r24738 = 1;
        double r24739 = x;
        double r24740 = r24739 + r24738;
        double r24741 = r24738 / r24740;
        double r24742 = 2;
        double r24743 = r24742 / r24739;
        double r24744 = r24741 - r24743;
        double r24745 = r24739 - r24738;
        double r24746 = r24738 / r24745;
        double r24747 = r24744 + r24746;
        return r24747;
}


double f_of(float x) {
        float r24748 = 2;
        float r24749 = x;
        float r24750 = 1;
        float r24751 = r24749 + r24750;
        float r24752 = r24751 * r24749;
        float r24753 = r24748 / r24752;
        float r24754 = r24749 - r24750;
        float r24755 = r24753 / r24754;
        return r24755;
}

double f_od(double x) {
        double r24756 = 2;
        double r24757 = x;
        double r24758 = 1;
        double r24759 = r24757 + r24758;
        double r24760 = r24759 * r24757;
        double r24761 = r24756 / r24760;
        double r24762 = r24757 - r24758;
        double r24763 = r24761 / r24762;
        return r24763;
}

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 r24764, r24765, r24766, r24767, r24768, r24769, r24770, r24771, r24772, r24773;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r24764, "1", 10, MPFR_RNDN);
        mpfr_init(r24765);
        mpfr_init(r24766);
        mpfr_init(r24767);
        mpfr_init_set_str(r24768, "2", 10, MPFR_RNDN);
        mpfr_init(r24769);
        mpfr_init(r24770);
        mpfr_init(r24771);
        mpfr_init(r24772);
        mpfr_init(r24773);
}

double f_im(double x) {
        ;
        mpfr_set_d(r24765, x, MPFR_RNDN);
        mpfr_add(r24766, r24765, r24764, MPFR_RNDN);
        mpfr_div(r24767, r24764, r24766, MPFR_RNDN);
        ;
        mpfr_div(r24769, r24768, r24765, MPFR_RNDN);
        mpfr_sub(r24770, r24767, r24769, MPFR_RNDN);
        mpfr_sub(r24771, r24765, r24764, MPFR_RNDN);
        mpfr_div(r24772, r24764, r24771, MPFR_RNDN);
        mpfr_add(r24773, r24770, r24772, MPFR_RNDN);
        return mpfr_get_d(r24773, MPFR_RNDN);
}

static mpfr_t r24774, r24775, r24776, r24777, r24778, r24779, r24780, r24781;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r24774, "2", 10, MPFR_RNDN);
        mpfr_init(r24775);
        mpfr_init_set_str(r24776, "1", 10, MPFR_RNDN);
        mpfr_init(r24777);
        mpfr_init(r24778);
        mpfr_init(r24779);
        mpfr_init(r24780);
        mpfr_init(r24781);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r24775, x, MPFR_RNDN);
        ;
        mpfr_add(r24777, r24775, r24776, MPFR_RNDN);
        mpfr_mul(r24778, r24777, r24775, MPFR_RNDN);
        mpfr_div(r24779, r24774, r24778, MPFR_RNDN);
        mpfr_sub(r24780, r24775, r24776, MPFR_RNDN);
        mpfr_div(r24781, r24779, r24780, MPFR_RNDN);
        return mpfr_get_d(r24781, MPFR_RNDN);
}

static mpfr_t r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r24782, "2", 10, MPFR_RNDN);
        mpfr_init(r24783);
        mpfr_init_set_str(r24784, "1", 10, MPFR_RNDN);
        mpfr_init(r24785);
        mpfr_init(r24786);
        mpfr_init(r24787);
        mpfr_init(r24788);
        mpfr_init(r24789);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r24783, x, MPFR_RNDN);
        ;
        mpfr_add(r24785, r24783, r24784, MPFR_RNDN);
        mpfr_mul(r24786, r24785, r24783, MPFR_RNDN);
        mpfr_div(r24787, r24782, r24786, MPFR_RNDN);
        mpfr_sub(r24788, r24783, r24784, MPFR_RNDN);
        mpfr_div(r24789, r24787, r24788, MPFR_RNDN);
        return mpfr_get_d(r24789, MPFR_RNDN);
}