#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 r7786 = 2;
        float r7787 = x;
        float r7788 = r7786 * r7787;
        float r7789 = exp(r7788);
        float r7790 = 1;
        float r7791 = r7789 - r7790;
        float r7792 = exp(r7787);
        float r7793 = r7792 - r7790;
        float r7794 = r7791 / r7793;
        float r7795 = sqrt(r7794);
        return r7795;
}

double f_id(double x) {
        double r7796 = 2;
        double r7797 = x;
        double r7798 = r7796 * r7797;
        double r7799 = exp(r7798);
        double r7800 = 1;
        double r7801 = r7799 - r7800;
        double r7802 = exp(r7797);
        double r7803 = r7802 - r7800;
        double r7804 = r7801 / r7803;
        double r7805 = sqrt(r7804);
        return r7805;
}


double f_of(float x) {
        float r7806 = x;
        float r7807 = exp(r7806);
        float r7808 = 1;
        float r7809 = r7807 + r7808;
        float r7810 = sqrt(r7809);
        return r7810;
}

double f_od(double x) {
        double r7811 = x;
        double r7812 = exp(r7811);
        double r7813 = 1;
        double r7814 = r7812 + r7813;
        double r7815 = sqrt(r7814);
        return r7815;
}

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 r7816, r7817, r7818, r7819, r7820, r7821, r7822, r7823, r7824, r7825;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r7816, "2", 10, MPFR_RNDN);
        mpfr_init(r7817);
        mpfr_init(r7818);
        mpfr_init(r7819);
        mpfr_init_set_str(r7820, "1", 10, MPFR_RNDN);
        mpfr_init(r7821);
        mpfr_init(r7822);
        mpfr_init(r7823);
        mpfr_init(r7824);
        mpfr_init(r7825);
}

double f_im(double x) {
        ;
        mpfr_set_d(r7817, x, MPFR_RNDN);
        mpfr_mul(r7818, r7816, r7817, MPFR_RNDN);
        mpfr_exp(r7819, r7818, MPFR_RNDN);
        ;
        mpfr_sub(r7821, r7819, r7820, MPFR_RNDN);
        mpfr_exp(r7822, r7817, MPFR_RNDN);
        mpfr_sub(r7823, r7822, r7820, MPFR_RNDN);
        mpfr_div(r7824, r7821, r7823, MPFR_RNDN);
        mpfr_sqrt(r7825, r7824, MPFR_RNDN);
        return mpfr_get_d(r7825, MPFR_RNDN);
}

static mpfr_t r7826, r7827, r7828, r7829, r7830;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7826);
        mpfr_init(r7827);
        mpfr_init_set_str(r7828, "1", 10, MPFR_RNDN);
        mpfr_init(r7829);
        mpfr_init(r7830);
}

double f_fm(double x) {
        mpfr_set_d(r7826, x, MPFR_RNDN);
        mpfr_exp(r7827, r7826, MPFR_RNDN);
        ;
        mpfr_add(r7829, r7827, r7828, MPFR_RNDN);
        mpfr_sqrt(r7830, r7829, MPFR_RNDN);
        return mpfr_get_d(r7830, MPFR_RNDN);
}

static mpfr_t r7831, r7832, r7833, r7834, r7835;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r7831);
        mpfr_init(r7832);
        mpfr_init_set_str(r7833, "1", 10, MPFR_RNDN);
        mpfr_init(r7834);
        mpfr_init(r7835);
}

double f_dm(double x) {
        mpfr_set_d(r7831, x, MPFR_RNDN);
        mpfr_exp(r7832, r7831, MPFR_RNDN);
        ;
        mpfr_add(r7834, r7832, r7833, MPFR_RNDN);
        mpfr_sqrt(r7835, r7834, MPFR_RNDN);
        return mpfr_get_d(r7835, MPFR_RNDN);
}