#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 r10668 = 1.0f;
        float r10669 = x;
        float r10670 = r10669 + r10668;
        float r10671 = r10668 / r10670;
        float r10672 = 2.0f;
        float r10673 = r10672 / r10669;
        float r10674 = r10671 - r10673;
        float r10675 = r10669 - r10668;
        float r10676 = r10668 / r10675;
        float r10677 = r10674 + r10676;
        return r10677;
}

double f_id(double x) {
        double r10678 = 1.0;
        double r10679 = x;
        double r10680 = r10679 + r10678;
        double r10681 = r10678 / r10680;
        double r10682 = 2.0;
        double r10683 = r10682 / r10679;
        double r10684 = r10681 - r10683;
        double r10685 = r10679 - r10678;
        double r10686 = r10678 / r10685;
        double r10687 = r10684 + r10686;
        return r10687;
}


double f_of(float x) {
        float r10688 = 2.0f;
        float r10689 = x;
        float r10690 = 1.0f;
        float r10691 = r10689 + r10690;
        float r10692 = r10691 * r10689;
        float r10693 = r10688 / r10692;
        float r10694 = r10689 - r10690;
        float r10695 = r10693 / r10694;
        return r10695;
}

double f_od(double x) {
        double r10696 = 2.0;
        double r10697 = x;
        double r10698 = 1.0;
        double r10699 = r10697 + r10698;
        double r10700 = r10699 * r10697;
        double r10701 = r10696 / r10700;
        double r10702 = r10697 - r10698;
        double r10703 = r10701 / r10702;
        return r10703;
}

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 r10704, r10705, r10706, r10707, r10708, r10709, r10710, r10711, r10712, r10713;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r10704, "1", 10, MPFR_RNDN);
        mpfr_init(r10705);
        mpfr_init(r10706);
        mpfr_init(r10707);
        mpfr_init_set_str(r10708, "2", 10, MPFR_RNDN);
        mpfr_init(r10709);
        mpfr_init(r10710);
        mpfr_init(r10711);
        mpfr_init(r10712);
        mpfr_init(r10713);
}

double f_im(double x) {
        ;
        mpfr_set_d(r10705, x, MPFR_RNDN);
        mpfr_add(r10706, r10705, r10704, MPFR_RNDN);
        mpfr_div(r10707, r10704, r10706, MPFR_RNDN);
        ;
        mpfr_div(r10709, r10708, r10705, MPFR_RNDN);
        mpfr_sub(r10710, r10707, r10709, MPFR_RNDN);
        mpfr_sub(r10711, r10705, r10704, MPFR_RNDN);
        mpfr_div(r10712, r10704, r10711, MPFR_RNDN);
        mpfr_add(r10713, r10710, r10712, MPFR_RNDN);
        return mpfr_get_d(r10713, MPFR_RNDN);
}

static mpfr_t r10714, r10715, r10716, r10717, r10718, r10719, r10720, r10721;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r10714, "2", 10, MPFR_RNDN);
        mpfr_init(r10715);
        mpfr_init_set_str(r10716, "1", 10, MPFR_RNDN);
        mpfr_init(r10717);
        mpfr_init(r10718);
        mpfr_init(r10719);
        mpfr_init(r10720);
        mpfr_init(r10721);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10715, x, MPFR_RNDN);
        ;
        mpfr_add(r10717, r10715, r10716, MPFR_RNDN);
        mpfr_mul(r10718, r10717, r10715, MPFR_RNDN);
        mpfr_div(r10719, r10714, r10718, MPFR_RNDN);
        mpfr_sub(r10720, r10715, r10716, MPFR_RNDN);
        mpfr_div(r10721, r10719, r10720, MPFR_RNDN);
        return mpfr_get_d(r10721, MPFR_RNDN);
}

static mpfr_t r10722, r10723, r10724, r10725, r10726, r10727, r10728, r10729;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r10722, "2", 10, MPFR_RNDN);
        mpfr_init(r10723);
        mpfr_init_set_str(r10724, "1", 10, MPFR_RNDN);
        mpfr_init(r10725);
        mpfr_init(r10726);
        mpfr_init(r10727);
        mpfr_init(r10728);
        mpfr_init(r10729);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10723, x, MPFR_RNDN);
        ;
        mpfr_add(r10725, r10723, r10724, MPFR_RNDN);
        mpfr_mul(r10726, r10725, r10723, MPFR_RNDN);
        mpfr_div(r10727, r10722, r10726, MPFR_RNDN);
        mpfr_sub(r10728, r10723, r10724, MPFR_RNDN);
        mpfr_div(r10729, r10727, r10728, MPFR_RNDN);
        return mpfr_get_d(r10729, MPFR_RNDN);
}