#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 r27047 = 1;
        float r27048 = x;
        float r27049 = r27048 + r27047;
        float r27050 = r27047 / r27049;
        float r27051 = 2;
        float r27052 = r27051 / r27048;
        float r27053 = r27050 - r27052;
        float r27054 = r27048 - r27047;
        float r27055 = r27047 / r27054;
        float r27056 = r27053 + r27055;
        return r27056;
}

double f_id(double x) {
        double r27057 = 1;
        double r27058 = x;
        double r27059 = r27058 + r27057;
        double r27060 = r27057 / r27059;
        double r27061 = 2;
        double r27062 = r27061 / r27058;
        double r27063 = r27060 - r27062;
        double r27064 = r27058 - r27057;
        double r27065 = r27057 / r27064;
        double r27066 = r27063 + r27065;
        return r27066;
}


double f_of(float x) {
        float r27067 = 2;
        float r27068 = x;
        float r27069 = 1;
        float r27070 = r27068 + r27069;
        float r27071 = r27070 * r27068;
        float r27072 = r27067 / r27071;
        float r27073 = r27068 - r27069;
        float r27074 = r27072 / r27073;
        return r27074;
}

double f_od(double x) {
        double r27075 = 2;
        double r27076 = x;
        double r27077 = 1;
        double r27078 = r27076 + r27077;
        double r27079 = r27078 * r27076;
        double r27080 = r27075 / r27079;
        double r27081 = r27076 - r27077;
        double r27082 = r27080 / r27081;
        return r27082;
}

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 r27083, r27084, r27085, r27086, r27087, r27088, r27089, r27090, r27091, r27092;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27083, "1", 10, MPFR_RNDN);
        mpfr_init(r27084);
        mpfr_init(r27085);
        mpfr_init(r27086);
        mpfr_init_set_str(r27087, "2", 10, MPFR_RNDN);
        mpfr_init(r27088);
        mpfr_init(r27089);
        mpfr_init(r27090);
        mpfr_init(r27091);
        mpfr_init(r27092);
}

double f_im(double x) {
        ;
        mpfr_set_d(r27084, x, MPFR_RNDN);
        mpfr_add(r27085, r27084, r27083, MPFR_RNDN);
        mpfr_div(r27086, r27083, r27085, MPFR_RNDN);
        ;
        mpfr_div(r27088, r27087, r27084, MPFR_RNDN);
        mpfr_sub(r27089, r27086, r27088, MPFR_RNDN);
        mpfr_sub(r27090, r27084, r27083, MPFR_RNDN);
        mpfr_div(r27091, r27083, r27090, MPFR_RNDN);
        mpfr_add(r27092, r27089, r27091, MPFR_RNDN);
        return mpfr_get_d(r27092, MPFR_RNDN);
}

static mpfr_t r27093, r27094, r27095, r27096, r27097, r27098, r27099, r27100;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27093, "2", 10, MPFR_RNDN);
        mpfr_init(r27094);
        mpfr_init_set_str(r27095, "1", 10, MPFR_RNDN);
        mpfr_init(r27096);
        mpfr_init(r27097);
        mpfr_init(r27098);
        mpfr_init(r27099);
        mpfr_init(r27100);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27094, x, MPFR_RNDN);
        ;
        mpfr_add(r27096, r27094, r27095, MPFR_RNDN);
        mpfr_mul(r27097, r27096, r27094, MPFR_RNDN);
        mpfr_div(r27098, r27093, r27097, MPFR_RNDN);
        mpfr_sub(r27099, r27094, r27095, MPFR_RNDN);
        mpfr_div(r27100, r27098, r27099, MPFR_RNDN);
        return mpfr_get_d(r27100, MPFR_RNDN);
}

static mpfr_t r27101, r27102, r27103, r27104, r27105, r27106, r27107, r27108;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27101, "2", 10, MPFR_RNDN);
        mpfr_init(r27102);
        mpfr_init_set_str(r27103, "1", 10, MPFR_RNDN);
        mpfr_init(r27104);
        mpfr_init(r27105);
        mpfr_init(r27106);
        mpfr_init(r27107);
        mpfr_init(r27108);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27102, x, MPFR_RNDN);
        ;
        mpfr_add(r27104, r27102, r27103, MPFR_RNDN);
        mpfr_mul(r27105, r27104, r27102, MPFR_RNDN);
        mpfr_div(r27106, r27101, r27105, MPFR_RNDN);
        mpfr_sub(r27107, r27102, r27103, MPFR_RNDN);
        mpfr_div(r27108, r27106, r27107, MPFR_RNDN);
        return mpfr_get_d(r27108, MPFR_RNDN);
}