#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 = 1;
        float r7807 = x;
        float r7808 = exp(r7807);
        float r7809 = 3;
        float r7810 = pow(r7808, r7809);
        float r7811 = r7806 + r7810;
        float r7812 = sqrt(r7811);
        float r7813 = r7807 + r7807;
        float r7814 = exp(r7813);
        float r7815 = r7806 - r7808;
        float r7816 = r7814 + r7815;
        float r7817 = sqrt(r7816);
        float r7818 = r7812 / r7817;
        return r7818;
}

double f_od(double x) {
        double r7819 = 1;
        double r7820 = x;
        double r7821 = exp(r7820);
        double r7822 = 3;
        double r7823 = pow(r7821, r7822);
        double r7824 = r7819 + r7823;
        double r7825 = sqrt(r7824);
        double r7826 = r7820 + r7820;
        double r7827 = exp(r7826);
        double r7828 = r7819 - r7821;
        double r7829 = r7827 + r7828;
        double r7830 = sqrt(r7829);
        double r7831 = r7825 / r7830;
        return r7831;
}

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 r7832, r7833, r7834, r7835, r7836, r7837, r7838, r7839, r7840, r7841;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r7832, "2", 10, MPFR_RNDN);
        mpfr_init(r7833);
        mpfr_init(r7834);
        mpfr_init(r7835);
        mpfr_init_set_str(r7836, "1", 10, MPFR_RNDN);
        mpfr_init(r7837);
        mpfr_init(r7838);
        mpfr_init(r7839);
        mpfr_init(r7840);
        mpfr_init(r7841);
}

double f_im(double x) {
        ;
        mpfr_set_d(r7833, x, MPFR_RNDN);
        mpfr_mul(r7834, r7832, r7833, MPFR_RNDN);
        mpfr_exp(r7835, r7834, MPFR_RNDN);
        ;
        mpfr_sub(r7837, r7835, r7836, MPFR_RNDN);
        mpfr_exp(r7838, r7833, MPFR_RNDN);
        mpfr_sub(r7839, r7838, r7836, MPFR_RNDN);
        mpfr_div(r7840, r7837, r7839, MPFR_RNDN);
        mpfr_sqrt(r7841, r7840, MPFR_RNDN);
        return mpfr_get_d(r7841, MPFR_RNDN);
}

static mpfr_t r7842, r7843, r7844, r7845, r7846, r7847, r7848, r7849, r7850, r7851, r7852, r7853, r7854;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r7842, "1", 10, MPFR_RNDN);
        mpfr_init(r7843);
        mpfr_init(r7844);
        mpfr_init_set_str(r7845, "3", 10, MPFR_RNDN);
        mpfr_init(r7846);
        mpfr_init(r7847);
        mpfr_init(r7848);
        mpfr_init(r7849);
        mpfr_init(r7850);
        mpfr_init(r7851);
        mpfr_init(r7852);
        mpfr_init(r7853);
        mpfr_init(r7854);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r7843, x, MPFR_RNDN);
        mpfr_exp(r7844, r7843, MPFR_RNDN);
        ;
        mpfr_pow(r7846, r7844, r7845, MPFR_RNDN);
        mpfr_add(r7847, r7842, r7846, MPFR_RNDN);
        mpfr_sqrt(r7848, r7847, MPFR_RNDN);
        mpfr_add(r7849, r7843, r7843, MPFR_RNDN);
        mpfr_exp(r7850, r7849, MPFR_RNDN);
        mpfr_sub(r7851, r7842, r7844, MPFR_RNDN);
        mpfr_add(r7852, r7850, r7851, MPFR_RNDN);
        mpfr_sqrt(r7853, r7852, MPFR_RNDN);
        mpfr_div(r7854, r7848, r7853, MPFR_RNDN);
        return mpfr_get_d(r7854, MPFR_RNDN);
}

static mpfr_t r7855, r7856, r7857, r7858, r7859, r7860, r7861, r7862, r7863, r7864, r7865, r7866, r7867;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r7855, "1", 10, MPFR_RNDN);
        mpfr_init(r7856);
        mpfr_init(r7857);
        mpfr_init_set_str(r7858, "3", 10, MPFR_RNDN);
        mpfr_init(r7859);
        mpfr_init(r7860);
        mpfr_init(r7861);
        mpfr_init(r7862);
        mpfr_init(r7863);
        mpfr_init(r7864);
        mpfr_init(r7865);
        mpfr_init(r7866);
        mpfr_init(r7867);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r7856, x, MPFR_RNDN);
        mpfr_exp(r7857, r7856, MPFR_RNDN);
        ;
        mpfr_pow(r7859, r7857, r7858, MPFR_RNDN);
        mpfr_add(r7860, r7855, r7859, MPFR_RNDN);
        mpfr_sqrt(r7861, r7860, MPFR_RNDN);
        mpfr_add(r7862, r7856, r7856, MPFR_RNDN);
        mpfr_exp(r7863, r7862, MPFR_RNDN);
        mpfr_sub(r7864, r7855, r7857, MPFR_RNDN);
        mpfr_add(r7865, r7863, r7864, MPFR_RNDN);
        mpfr_sqrt(r7866, r7865, MPFR_RNDN);
        mpfr_div(r7867, r7861, r7866, MPFR_RNDN);
        return mpfr_get_d(r7867, MPFR_RNDN);
}