#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 r27590 = x;
        float r27591 = 1;
        float r27592 = r27590 + r27591;
        float r27593 = sqrt(r27592);
        float r27594 = sqrt(r27590);
        float r27595 = r27593 - r27594;
        return r27595;
}

double f_id(double x) {
        double r27596 = x;
        double r27597 = 1;
        double r27598 = r27596 + r27597;
        double r27599 = sqrt(r27598);
        double r27600 = sqrt(r27596);
        double r27601 = r27599 - r27600;
        return r27601;
}


double f_of(float x) {
        float r27602 = 1;
        float r27603 = x;
        float r27604 = r27603 + r27602;
        float r27605 = sqrt(r27604);
        float r27606 = sqrt(r27603);
        float r27607 = r27605 + r27606;
        float r27608 = r27607 * r27607;
        float r27609 = r27602 / r27608;
        float r27610 = 1/2;
        float r27611 = pow(r27609, r27610);
        return r27611;
}

double f_od(double x) {
        double r27612 = 1;
        double r27613 = x;
        double r27614 = r27613 + r27612;
        double r27615 = sqrt(r27614);
        double r27616 = sqrt(r27613);
        double r27617 = r27615 + r27616;
        double r27618 = r27617 * r27617;
        double r27619 = r27612 / r27618;
        double r27620 = 1/2;
        double r27621 = pow(r27619, r27620);
        return r27621;
}

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 r27622, r27623, r27624, r27625, r27626, r27627;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27622);
        mpfr_init_set_str(r27623, "1", 10, MPFR_RNDN);
        mpfr_init(r27624);
        mpfr_init(r27625);
        mpfr_init(r27626);
        mpfr_init(r27627);
}

double f_im(double x) {
        mpfr_set_d(r27622, x, MPFR_RNDN);
        ;
        mpfr_add(r27624, r27622, r27623, MPFR_RNDN);
        mpfr_sqrt(r27625, r27624, MPFR_RNDN);
        mpfr_sqrt(r27626, r27622, MPFR_RNDN);
        mpfr_sub(r27627, r27625, r27626, MPFR_RNDN);
        return mpfr_get_d(r27627, MPFR_RNDN);
}

static mpfr_t r27628, r27629, r27630, r27631, r27632, r27633, r27634, r27635, r27636, r27637;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27628, "1", 10, MPFR_RNDN);
        mpfr_init(r27629);
        mpfr_init(r27630);
        mpfr_init(r27631);
        mpfr_init(r27632);
        mpfr_init(r27633);
        mpfr_init(r27634);
        mpfr_init(r27635);
        mpfr_init_set_str(r27636, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27637);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27629, x, MPFR_RNDN);
        mpfr_add(r27630, r27629, r27628, MPFR_RNDN);
        mpfr_sqrt(r27631, r27630, MPFR_RNDN);
        mpfr_sqrt(r27632, r27629, MPFR_RNDN);
        mpfr_add(r27633, r27631, r27632, MPFR_RNDN);
        mpfr_mul(r27634, r27633, r27633, MPFR_RNDN);
        mpfr_div(r27635, r27628, r27634, MPFR_RNDN);
        ;
        mpfr_pow(r27637, r27635, r27636, MPFR_RNDN);
        return mpfr_get_d(r27637, MPFR_RNDN);
}

static mpfr_t r27638, r27639, r27640, r27641, r27642, r27643, r27644, r27645, r27646, r27647;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27638, "1", 10, MPFR_RNDN);
        mpfr_init(r27639);
        mpfr_init(r27640);
        mpfr_init(r27641);
        mpfr_init(r27642);
        mpfr_init(r27643);
        mpfr_init(r27644);
        mpfr_init(r27645);
        mpfr_init_set_str(r27646, "1/2", 10, MPFR_RNDN);
        mpfr_init(r27647);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27639, x, MPFR_RNDN);
        mpfr_add(r27640, r27639, r27638, MPFR_RNDN);
        mpfr_sqrt(r27641, r27640, MPFR_RNDN);
        mpfr_sqrt(r27642, r27639, MPFR_RNDN);
        mpfr_add(r27643, r27641, r27642, MPFR_RNDN);
        mpfr_mul(r27644, r27643, r27643, MPFR_RNDN);
        mpfr_div(r27645, r27638, r27644, MPFR_RNDN);
        ;
        mpfr_pow(r27647, r27645, r27646, MPFR_RNDN);
        return mpfr_get_d(r27647, MPFR_RNDN);
}