#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 r11091 = 1;
        float r11092 = x;
        float r11093 = r11092 + r11091;
        float r11094 = r11091 / r11093;
        float r11095 = 2;
        float r11096 = r11095 / r11092;
        float r11097 = r11094 - r11096;
        float r11098 = r11092 - r11091;
        float r11099 = r11091 / r11098;
        float r11100 = r11097 + r11099;
        return r11100;
}

double f_id(double x) {
        double r11101 = 1;
        double r11102 = x;
        double r11103 = r11102 + r11101;
        double r11104 = r11101 / r11103;
        double r11105 = 2;
        double r11106 = r11105 / r11102;
        double r11107 = r11104 - r11106;
        double r11108 = r11102 - r11101;
        double r11109 = r11101 / r11108;
        double r11110 = r11107 + r11109;
        return r11110;
}


double f_of(float x) {
        float r11111 = 2;
        float r11112 = x;
        float r11113 = 1;
        float r11114 = r11112 - r11113;
        float r11115 = fma(r11112, r11112, r11112);
        float r11116 = r11114 * r11115;
        float r11117 = r11111 / r11116;
        return r11117;
}

double f_od(double x) {
        double r11118 = 2;
        double r11119 = x;
        double r11120 = 1;
        double r11121 = r11119 - r11120;
        double r11122 = fma(r11119, r11119, r11119);
        double r11123 = r11121 * r11122;
        double r11124 = r11118 / r11123;
        return r11124;
}

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 r11125, r11126, r11127, r11128, r11129, r11130, r11131, r11132, r11133, r11134;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r11125, "1", 10, MPFR_RNDN);
        mpfr_init(r11126);
        mpfr_init(r11127);
        mpfr_init(r11128);
        mpfr_init_set_str(r11129, "2", 10, MPFR_RNDN);
        mpfr_init(r11130);
        mpfr_init(r11131);
        mpfr_init(r11132);
        mpfr_init(r11133);
        mpfr_init(r11134);
}

double f_im(double x) {
        ;
        mpfr_set_d(r11126, x, MPFR_RNDN);
        mpfr_add(r11127, r11126, r11125, MPFR_RNDN);
        mpfr_div(r11128, r11125, r11127, MPFR_RNDN);
        ;
        mpfr_div(r11130, r11129, r11126, MPFR_RNDN);
        mpfr_sub(r11131, r11128, r11130, MPFR_RNDN);
        mpfr_sub(r11132, r11126, r11125, MPFR_RNDN);
        mpfr_div(r11133, r11125, r11132, MPFR_RNDN);
        mpfr_add(r11134, r11131, r11133, MPFR_RNDN);
        return mpfr_get_d(r11134, MPFR_RNDN);
}

static mpfr_t r11135, r11136, r11137, r11138, r11139, r11140, r11141;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r11135, "2", 10, MPFR_RNDN);
        mpfr_init(r11136);
        mpfr_init_set_str(r11137, "1", 10, MPFR_RNDN);
        mpfr_init(r11138);
        mpfr_init(r11139);
        mpfr_init(r11140);
        mpfr_init(r11141);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r11136, x, MPFR_RNDN);
        ;
        mpfr_sub(r11138, r11136, r11137, MPFR_RNDN);
        mpfr_fma(r11139, r11136, r11136, r11136, MPFR_RNDN);
        mpfr_mul(r11140, r11138, r11139, MPFR_RNDN);
        mpfr_div(r11141, r11135, r11140, MPFR_RNDN);
        return mpfr_get_d(r11141, MPFR_RNDN);
}

static mpfr_t r11142, r11143, r11144, r11145, r11146, r11147, r11148;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r11142, "2", 10, MPFR_RNDN);
        mpfr_init(r11143);
        mpfr_init_set_str(r11144, "1", 10, MPFR_RNDN);
        mpfr_init(r11145);
        mpfr_init(r11146);
        mpfr_init(r11147);
        mpfr_init(r11148);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r11143, x, MPFR_RNDN);
        ;
        mpfr_sub(r11145, r11143, r11144, MPFR_RNDN);
        mpfr_fma(r11146, r11143, r11143, r11143, MPFR_RNDN);
        mpfr_mul(r11147, r11145, r11146, MPFR_RNDN);
        mpfr_div(r11148, r11142, r11147, MPFR_RNDN);
        return mpfr_get_d(r11148, MPFR_RNDN);
}