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

char *name = "2isqrt (example 3.6)";

double f_if(float x) {
        float r27381 = 1;
        float r27382 = x;
        float r27383 = sqrt(r27382);
        float r27384 = r27381 / r27383;
        float r27385 = r27382 + r27381;
        float r27386 = sqrt(r27385);
        float r27387 = r27381 / r27386;
        float r27388 = r27384 - r27387;
        return r27388;
}

double f_id(double x) {
        double r27389 = 1;
        double r27390 = x;
        double r27391 = sqrt(r27390);
        double r27392 = r27389 / r27391;
        double r27393 = r27390 + r27389;
        double r27394 = sqrt(r27393);
        double r27395 = r27389 / r27394;
        double r27396 = r27392 - r27395;
        return r27396;
}


double f_of(float x) {
        float r27397 = 1;
        float r27398 = x;
        float r27399 = sqrt(r27398);
        float r27400 = r27397 / r27399;
        float r27401 = r27398 + r27397;
        float r27402 = sqrt(r27401);
        float r27403 = r27397 / r27402;
        float r27404 = r27400 - r27403;
        return r27404;
}

double f_od(double x) {
        double r27405 = 1;
        double r27406 = x;
        double r27407 = sqrt(r27406);
        double r27408 = r27405 / r27407;
        double r27409 = r27406 + r27405;
        double r27410 = sqrt(r27409);
        double r27411 = r27405 / r27410;
        double r27412 = r27408 - r27411;
        return r27412;
}

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 r27413, r27414, r27415, r27416, r27417, r27418, r27419, r27420;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27413, "1", 10, MPFR_RNDN);
        mpfr_init(r27414);
        mpfr_init(r27415);
        mpfr_init(r27416);
        mpfr_init(r27417);
        mpfr_init(r27418);
        mpfr_init(r27419);
        mpfr_init(r27420);
}

double f_im(double x) {
        ;
        mpfr_set_d(r27414, x, MPFR_RNDN);
        mpfr_sqrt(r27415, r27414, MPFR_RNDN);
        mpfr_div(r27416, r27413, r27415, MPFR_RNDN);
        mpfr_add(r27417, r27414, r27413, MPFR_RNDN);
        mpfr_sqrt(r27418, r27417, MPFR_RNDN);
        mpfr_div(r27419, r27413, r27418, MPFR_RNDN);
        mpfr_sub(r27420, r27416, r27419, MPFR_RNDN);
        return mpfr_get_d(r27420, MPFR_RNDN);
}

static mpfr_t r27421, r27422, r27423, r27424, r27425, r27426, r27427, r27428;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27421, "1", 10, MPFR_RNDN);
        mpfr_init(r27422);
        mpfr_init(r27423);
        mpfr_init(r27424);
        mpfr_init(r27425);
        mpfr_init(r27426);
        mpfr_init(r27427);
        mpfr_init(r27428);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27422, x, MPFR_RNDN);
        mpfr_sqrt(r27423, r27422, MPFR_RNDN);
        mpfr_div(r27424, r27421, r27423, MPFR_RNDN);
        mpfr_add(r27425, r27422, r27421, MPFR_RNDN);
        mpfr_sqrt(r27426, r27425, MPFR_RNDN);
        mpfr_div(r27427, r27421, r27426, MPFR_RNDN);
        mpfr_sub(r27428, r27424, r27427, MPFR_RNDN);
        return mpfr_get_d(r27428, MPFR_RNDN);
}

static mpfr_t r27429, r27430, r27431, r27432, r27433, r27434, r27435, r27436;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27429, "1", 10, MPFR_RNDN);
        mpfr_init(r27430);
        mpfr_init(r27431);
        mpfr_init(r27432);
        mpfr_init(r27433);
        mpfr_init(r27434);
        mpfr_init(r27435);
        mpfr_init(r27436);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27430, x, MPFR_RNDN);
        mpfr_sqrt(r27431, r27430, MPFR_RNDN);
        mpfr_div(r27432, r27429, r27431, MPFR_RNDN);
        mpfr_add(r27433, r27430, r27429, MPFR_RNDN);
        mpfr_sqrt(r27434, r27433, MPFR_RNDN);
        mpfr_div(r27435, r27429, r27434, MPFR_RNDN);
        mpfr_sub(r27436, r27432, r27435, MPFR_RNDN);
        return mpfr_get_d(r27436, MPFR_RNDN);
}