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

char *name = "2frac (problem 3.3.1)";

double f_if(float x) {
        float r9253 = 1.0f;
        float r9254 = x;
        float r9255 = r9254 + r9253;
        float r9256 = r9253 / r9255;
        float r9257 = r9253 / r9254;
        float r9258 = r9256 - r9257;
        return r9258;
}

double f_id(double x) {
        double r9259 = 1.0;
        double r9260 = x;
        double r9261 = r9260 + r9259;
        double r9262 = r9259 / r9261;
        double r9263 = r9259 / r9260;
        double r9264 = r9262 - r9263;
        return r9264;
}


double f_of(float x) {
        float r9265 = 1.0f;
        float r9266 = -r9265;
        float r9267 = x;
        float r9268 = fma(r9267, r9267, r9267);
        float r9269 = r9266 / r9268;
        return r9269;
}

double f_od(double x) {
        double r9270 = 1.0;
        double r9271 = -r9270;
        double r9272 = x;
        double r9273 = fma(r9272, r9272, r9272);
        double r9274 = r9271 / r9273;
        return r9274;
}

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 r9275, r9276, r9277, r9278, r9279, r9280;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r9275, "1", 10, MPFR_RNDN);
        mpfr_init(r9276);
        mpfr_init(r9277);
        mpfr_init(r9278);
        mpfr_init(r9279);
        mpfr_init(r9280);
}

double f_im(double x) {
        ;
        mpfr_set_d(r9276, x, MPFR_RNDN);
        mpfr_add(r9277, r9276, r9275, MPFR_RNDN);
        mpfr_div(r9278, r9275, r9277, MPFR_RNDN);
        mpfr_div(r9279, r9275, r9276, MPFR_RNDN);
        mpfr_sub(r9280, r9278, r9279, MPFR_RNDN);
        return mpfr_get_d(r9280, MPFR_RNDN);
}

static mpfr_t r9281, r9282, r9283, r9284, r9285;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r9281, "1", 10, MPFR_RNDN);
        mpfr_init(r9282);
        mpfr_init(r9283);
        mpfr_init(r9284);
        mpfr_init(r9285);
}

double f_fm(double x) {
        ;
        mpfr_neg(r9282, r9281, MPFR_RNDN);
        mpfr_set_d(r9283, x, MPFR_RNDN);
        mpfr_fma(r9284, r9283, r9283, r9283, MPFR_RNDN);
        mpfr_div(r9285, r9282, r9284, MPFR_RNDN);
        return mpfr_get_d(r9285, MPFR_RNDN);
}

static mpfr_t r9286, r9287, r9288, r9289, r9290;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init_set_str(r9286, "1", 10, MPFR_RNDN);
        mpfr_init(r9287);
        mpfr_init(r9288);
        mpfr_init(r9289);
        mpfr_init(r9290);
}

double f_dm(double x) {
        ;
        mpfr_neg(r9287, r9286, MPFR_RNDN);
        mpfr_set_d(r9288, x, MPFR_RNDN);
        mpfr_fma(r9289, r9288, r9288, r9288, MPFR_RNDN);
        mpfr_div(r9290, r9287, r9289, MPFR_RNDN);
        return mpfr_get_d(r9290, MPFR_RNDN);
}