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

char *name = "sqrtexp (problem 3.4.4)";

double f_if(float x) {
        float r18362 = 2.0f;
        float r18363 = x;
        float r18364 = r18362 * r18363;
        float r18365 = exp(r18364);
        float r18366 = 1.0f;
        float r18367 = r18365 - r18366;
        float r18368 = exp(r18363);
        float r18369 = r18368 - r18366;
        float r18370 = r18367 / r18369;
        float r18371 = sqrt(r18370);
        return r18371;
}

double f_id(double x) {
        double r18372 = 2.0;
        double r18373 = x;
        double r18374 = r18372 * r18373;
        double r18375 = exp(r18374);
        double r18376 = 1.0;
        double r18377 = r18375 - r18376;
        double r18378 = exp(r18373);
        double r18379 = r18378 - r18376;
        double r18380 = r18377 / r18379;
        double r18381 = sqrt(r18380);
        return r18381;
}


double f_of(float x) {
        float r18382 = 1.0f;
        float r18383 = x;
        float r18384 = exp(r18383);
        float r18385 = r18382 + r18384;
        float r18386 = r18385 / r18382;
        float r18387 = sqrt(r18386);
        return r18387;
}

double f_od(double x) {
        double r18388 = 1.0;
        double r18389 = x;
        double r18390 = exp(r18389);
        double r18391 = r18388 + r18390;
        double r18392 = r18391 / r18388;
        double r18393 = sqrt(r18392);
        return r18393;
}

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 r18394, r18395, r18396, r18397, r18398, r18399, r18400, r18401, r18402, r18403;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r18394, "2", 10, MPFR_RNDN);
        mpfr_init(r18395);
        mpfr_init(r18396);
        mpfr_init(r18397);
        mpfr_init_set_str(r18398, "1", 10, MPFR_RNDN);
        mpfr_init(r18399);
        mpfr_init(r18400);
        mpfr_init(r18401);
        mpfr_init(r18402);
        mpfr_init(r18403);
}

double f_im(double x) {
        ;
        mpfr_set_d(r18395, x, MPFR_RNDN);
        mpfr_mul(r18396, r18394, r18395, MPFR_RNDN);
        mpfr_exp(r18397, r18396, MPFR_RNDN);
        ;
        mpfr_sub(r18399, r18397, r18398, MPFR_RNDN);
        mpfr_exp(r18400, r18395, MPFR_RNDN);
        mpfr_sub(r18401, r18400, r18398, MPFR_RNDN);
        mpfr_div(r18402, r18399, r18401, MPFR_RNDN);
        mpfr_sqrt(r18403, r18402, MPFR_RNDN);
        return mpfr_get_d(r18403, MPFR_RNDN);
}

static mpfr_t r18404, r18405, r18406, r18407, r18408, r18409;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r18404, "1", 10, MPFR_RNDN);
        mpfr_init(r18405);
        mpfr_init(r18406);
        mpfr_init(r18407);
        mpfr_init(r18408);
        mpfr_init(r18409);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r18405, x, MPFR_RNDN);
        mpfr_exp(r18406, r18405, MPFR_RNDN);
        mpfr_add(r18407, r18404, r18406, MPFR_RNDN);
        mpfr_div(r18408, r18407, r18404, MPFR_RNDN);
        mpfr_sqrt(r18409, r18408, MPFR_RNDN);
        return mpfr_get_d(r18409, MPFR_RNDN);
}

static mpfr_t r18410, r18411, r18412, r18413, r18414, r18415;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r18410, "1", 10, MPFR_RNDN);
        mpfr_init(r18411);
        mpfr_init(r18412);
        mpfr_init(r18413);
        mpfr_init(r18414);
        mpfr_init(r18415);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r18411, x, MPFR_RNDN);
        mpfr_exp(r18412, r18411, MPFR_RNDN);
        mpfr_add(r18413, r18410, r18412, MPFR_RNDN);
        mpfr_div(r18414, r18413, r18410, MPFR_RNDN);
        mpfr_sqrt(r18415, r18414, MPFR_RNDN);
        return mpfr_get_d(r18415, MPFR_RNDN);
}