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

char *name = "2sqrt (example 3.1)";

double f_if(float x) {
        float r10678 = x;
        float r10679 = 1.0f;
        float r10680 = r10678 + r10679;
        float r10681 = sqrt(r10680);
        float r10682 = sqrt(r10678);
        float r10683 = r10681 - r10682;
        return r10683;
}

double f_id(double x) {
        double r10684 = x;
        double r10685 = 1.0;
        double r10686 = r10684 + r10685;
        double r10687 = sqrt(r10686);
        double r10688 = sqrt(r10684);
        double r10689 = r10687 - r10688;
        return r10689;
}


double f_of(float x) {
        float r10690 = 1.0f;
        float r10691 = x;
        float r10692 = r10691 + r10690;
        float r10693 = sqrt(r10692);
        float r10694 = sqrt(r10691);
        float r10695 = r10693 + r10694;
        float r10696 = r10690 / r10695;
        return r10696;
}

double f_od(double x) {
        double r10697 = 1.0;
        double r10698 = x;
        double r10699 = r10698 + r10697;
        double r10700 = sqrt(r10699);
        double r10701 = sqrt(r10698);
        double r10702 = r10700 + r10701;
        double r10703 = r10697 / 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;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r10704);
        mpfr_init_set_str(r10705, "1", 10, MPFR_RNDN);
        mpfr_init(r10706);
        mpfr_init(r10707);
        mpfr_init(r10708);
        mpfr_init(r10709);
}

double f_im(double x) {
        mpfr_set_d(r10704, x, MPFR_RNDN);
        ;
        mpfr_add(r10706, r10704, r10705, MPFR_RNDN);
        mpfr_sqrt(r10707, r10706, MPFR_RNDN);
        mpfr_sqrt(r10708, r10704, MPFR_RNDN);
        mpfr_sub(r10709, r10707, r10708, MPFR_RNDN);
        return mpfr_get_d(r10709, MPFR_RNDN);
}

static mpfr_t r10710, r10711, r10712, r10713, r10714, r10715, r10716;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10710, "1", 10, MPFR_RNDN);
        mpfr_init(r10711);
        mpfr_init(r10712);
        mpfr_init(r10713);
        mpfr_init(r10714);
        mpfr_init(r10715);
        mpfr_init(r10716);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10711, x, MPFR_RNDN);
        mpfr_add(r10712, r10711, r10710, MPFR_RNDN);
        mpfr_sqrt(r10713, r10712, MPFR_RNDN);
        mpfr_sqrt(r10714, r10711, MPFR_RNDN);
        mpfr_add(r10715, r10713, r10714, MPFR_RNDN);
        mpfr_div(r10716, r10710, r10715, MPFR_RNDN);
        return mpfr_get_d(r10716, MPFR_RNDN);
}

static mpfr_t r10717, r10718, r10719, r10720, r10721, r10722, r10723;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r10717, "1", 10, MPFR_RNDN);
        mpfr_init(r10718);
        mpfr_init(r10719);
        mpfr_init(r10720);
        mpfr_init(r10721);
        mpfr_init(r10722);
        mpfr_init(r10723);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10718, x, MPFR_RNDN);
        mpfr_add(r10719, r10718, r10717, MPFR_RNDN);
        mpfr_sqrt(r10720, r10719, MPFR_RNDN);
        mpfr_sqrt(r10721, r10718, MPFR_RNDN);
        mpfr_add(r10722, r10720, r10721, MPFR_RNDN);
        mpfr_div(r10723, r10717, r10722, MPFR_RNDN);
        return mpfr_get_d(r10723, MPFR_RNDN);
}