#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 r27425 = 1;
        float r27426 = x;
        float r27427 = sqrt(r27426);
        float r27428 = r27425 / r27427;
        float r27429 = r27426 + r27425;
        float r27430 = sqrt(r27429);
        float r27431 = r27425 / r27430;
        float r27432 = r27428 - r27431;
        return r27432;
}

double f_id(double x) {
        double r27433 = 1;
        double r27434 = x;
        double r27435 = sqrt(r27434);
        double r27436 = r27433 / r27435;
        double r27437 = r27434 + r27433;
        double r27438 = sqrt(r27437);
        double r27439 = r27433 / r27438;
        double r27440 = r27436 - r27439;
        return r27440;
}


double f_of(float x) {
        float r27441 = 1;
        float r27442 = x;
        float r27443 = sqrt(r27442);
        float r27444 = r27441 / r27443;
        float r27445 = r27442 + r27441;
        float r27446 = sqrt(r27445);
        float r27447 = r27441 / r27446;
        float r27448 = r27444 - r27447;
        return r27448;
}

double f_od(double x) {
        double r27449 = 1;
        double r27450 = x;
        double r27451 = sqrt(r27450);
        double r27452 = r27449 / r27451;
        double r27453 = r27450 + r27449;
        double r27454 = sqrt(r27453);
        double r27455 = r27449 / r27454;
        double r27456 = r27452 - r27455;
        return r27456;
}

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 r27457, r27458, r27459, r27460, r27461, r27462, r27463, r27464;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27457, "1", 10, MPFR_RNDN);
        mpfr_init(r27458);
        mpfr_init(r27459);
        mpfr_init(r27460);
        mpfr_init(r27461);
        mpfr_init(r27462);
        mpfr_init(r27463);
        mpfr_init(r27464);
}

double f_im(double x) {
        ;
        mpfr_set_d(r27458, x, MPFR_RNDN);
        mpfr_sqrt(r27459, r27458, MPFR_RNDN);
        mpfr_div(r27460, r27457, r27459, MPFR_RNDN);
        mpfr_add(r27461, r27458, r27457, MPFR_RNDN);
        mpfr_sqrt(r27462, r27461, MPFR_RNDN);
        mpfr_div(r27463, r27457, r27462, MPFR_RNDN);
        mpfr_sub(r27464, r27460, r27463, MPFR_RNDN);
        return mpfr_get_d(r27464, MPFR_RNDN);
}

static mpfr_t r27465, r27466, r27467, r27468, r27469, r27470, r27471, r27472;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27465, "1", 10, MPFR_RNDN);
        mpfr_init(r27466);
        mpfr_init(r27467);
        mpfr_init(r27468);
        mpfr_init(r27469);
        mpfr_init(r27470);
        mpfr_init(r27471);
        mpfr_init(r27472);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27466, x, MPFR_RNDN);
        mpfr_sqrt(r27467, r27466, MPFR_RNDN);
        mpfr_div(r27468, r27465, r27467, MPFR_RNDN);
        mpfr_add(r27469, r27466, r27465, MPFR_RNDN);
        mpfr_sqrt(r27470, r27469, MPFR_RNDN);
        mpfr_div(r27471, r27465, r27470, MPFR_RNDN);
        mpfr_sub(r27472, r27468, r27471, MPFR_RNDN);
        return mpfr_get_d(r27472, MPFR_RNDN);
}

static mpfr_t r27473, r27474, r27475, r27476, r27477, r27478, r27479, r27480;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1104);
        mpfr_init_set_str(r27473, "1", 10, MPFR_RNDN);
        mpfr_init(r27474);
        mpfr_init(r27475);
        mpfr_init(r27476);
        mpfr_init(r27477);
        mpfr_init(r27478);
        mpfr_init(r27479);
        mpfr_init(r27480);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27474, x, MPFR_RNDN);
        mpfr_sqrt(r27475, r27474, MPFR_RNDN);
        mpfr_div(r27476, r27473, r27475, MPFR_RNDN);
        mpfr_add(r27477, r27474, r27473, MPFR_RNDN);
        mpfr_sqrt(r27478, r27477, MPFR_RNDN);
        mpfr_div(r27479, r27473, r27478, MPFR_RNDN);
        mpfr_sub(r27480, r27476, r27479, MPFR_RNDN);
        return mpfr_get_d(r27480, MPFR_RNDN);
}