#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 r27402 = x;
        float r27403 = 1;
        float r27404 = r27402 + r27403;
        float r27405 = sqrt(r27404);
        float r27406 = sqrt(r27402);
        float r27407 = r27405 - r27406;
        return r27407;
}

double f_id(double x) {
        double r27408 = x;
        double r27409 = 1;
        double r27410 = r27408 + r27409;
        double r27411 = sqrt(r27410);
        double r27412 = sqrt(r27408);
        double r27413 = r27411 - r27412;
        return r27413;
}


double f_of(float x) {
        float r27414 = 1;
        float r27415 = x;
        float r27416 = r27415 + r27414;
        float r27417 = sqrt(r27416);
        float r27418 = sqrt(r27415);
        float r27419 = r27417 + r27418;
        float r27420 = r27414 / r27419;
        return r27420;
}

double f_od(double x) {
        double r27421 = 1;
        double r27422 = x;
        double r27423 = r27422 + r27421;
        double r27424 = sqrt(r27423);
        double r27425 = sqrt(r27422);
        double r27426 = r27424 + r27425;
        double r27427 = r27421 / r27426;
        return r27427;
}

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 r27428, r27429, r27430, r27431, r27432, r27433;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r27428);
        mpfr_init_set_str(r27429, "1", 10, MPFR_RNDN);
        mpfr_init(r27430);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
}

double f_im(double x) {
        mpfr_set_d(r27428, x, MPFR_RNDN);
        ;
        mpfr_add(r27430, r27428, r27429, MPFR_RNDN);
        mpfr_sqrt(r27431, r27430, MPFR_RNDN);
        mpfr_sqrt(r27432, r27428, MPFR_RNDN);
        mpfr_sub(r27433, r27431, r27432, MPFR_RNDN);
        return mpfr_get_d(r27433, MPFR_RNDN);
}

static mpfr_t r27434, r27435, r27436, r27437, r27438, r27439, r27440;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27434, "1", 10, MPFR_RNDN);
        mpfr_init(r27435);
        mpfr_init(r27436);
        mpfr_init(r27437);
        mpfr_init(r27438);
        mpfr_init(r27439);
        mpfr_init(r27440);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27435, x, MPFR_RNDN);
        mpfr_add(r27436, r27435, r27434, MPFR_RNDN);
        mpfr_sqrt(r27437, r27436, MPFR_RNDN);
        mpfr_sqrt(r27438, r27435, MPFR_RNDN);
        mpfr_add(r27439, r27437, r27438, MPFR_RNDN);
        mpfr_div(r27440, r27434, r27439, MPFR_RNDN);
        return mpfr_get_d(r27440, MPFR_RNDN);
}

static mpfr_t r27441, r27442, r27443, r27444, r27445, r27446, r27447;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r27441, "1", 10, MPFR_RNDN);
        mpfr_init(r27442);
        mpfr_init(r27443);
        mpfr_init(r27444);
        mpfr_init(r27445);
        mpfr_init(r27446);
        mpfr_init(r27447);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27442, x, MPFR_RNDN);
        mpfr_add(r27443, r27442, r27441, MPFR_RNDN);
        mpfr_sqrt(r27444, r27443, MPFR_RNDN);
        mpfr_sqrt(r27445, r27442, MPFR_RNDN);
        mpfr_add(r27446, r27444, r27445, MPFR_RNDN);
        mpfr_div(r27447, r27441, r27446, MPFR_RNDN);
        return mpfr_get_d(r27447, MPFR_RNDN);
}