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

char *name = "NMSE example 3.1";

double f_if(float x) {
        float r17743 = x;
        float r17744 = 1.0f;
        float r17745 = r17743 + r17744;
        float r17746 = sqrt(r17745);
        float r17747 = sqrt(r17743);
        float r17748 = r17746 - r17747;
        return r17748;
}

double f_id(double x) {
        double r17749 = x;
        double r17750 = 1.0;
        double r17751 = r17749 + r17750;
        double r17752 = sqrt(r17751);
        double r17753 = sqrt(r17749);
        double r17754 = r17752 - r17753;
        return r17754;
}


double f_of(float x) {
        float r17755 = 1.0f;
        float r17756 = x;
        float r17757 = r17756 + r17755;
        float r17758 = sqrt(r17757);
        float r17759 = sqrt(r17756);
        float r17760 = r17758 + r17759;
        float r17761 = r17755 / r17760;
        return r17761;
}

double f_od(double x) {
        double r17762 = 1.0;
        double r17763 = x;
        double r17764 = r17763 + r17762;
        double r17765 = sqrt(r17764);
        double r17766 = sqrt(r17763);
        double r17767 = r17765 + r17766;
        double r17768 = r17762 / r17767;
        return r17768;
}

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 r17769, r17770, r17771, r17772, r17773, r17774;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17769);
        mpfr_init_set_str(r17770, "1", 10, MPFR_RNDN);
        mpfr_init(r17771);
        mpfr_init(r17772);
        mpfr_init(r17773);
        mpfr_init(r17774);
}

double f_im(double x) {
        mpfr_set_d(r17769, x, MPFR_RNDN);
        ;
        mpfr_add(r17771, r17769, r17770, MPFR_RNDN);
        mpfr_sqrt(r17772, r17771, MPFR_RNDN);
        mpfr_sqrt(r17773, r17769, MPFR_RNDN);
        mpfr_sub(r17774, r17772, r17773, MPFR_RNDN);
        return mpfr_get_d(r17774, MPFR_RNDN);
}

static mpfr_t r17775, r17776, r17777, r17778, r17779, r17780, r17781;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17775, "1", 10, MPFR_RNDN);
        mpfr_init(r17776);
        mpfr_init(r17777);
        mpfr_init(r17778);
        mpfr_init(r17779);
        mpfr_init(r17780);
        mpfr_init(r17781);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r17776, x, MPFR_RNDN);
        mpfr_add(r17777, r17776, r17775, MPFR_RNDN);
        mpfr_sqrt(r17778, r17777, MPFR_RNDN);
        mpfr_sqrt(r17779, r17776, MPFR_RNDN);
        mpfr_add(r17780, r17778, r17779, MPFR_RNDN);
        mpfr_div(r17781, r17775, r17780, MPFR_RNDN);
        return mpfr_get_d(r17781, MPFR_RNDN);
}

static mpfr_t r17782, r17783, r17784, r17785, r17786, r17787, r17788;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17782, "1", 10, MPFR_RNDN);
        mpfr_init(r17783);
        mpfr_init(r17784);
        mpfr_init(r17785);
        mpfr_init(r17786);
        mpfr_init(r17787);
        mpfr_init(r17788);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r17783, x, MPFR_RNDN);
        mpfr_add(r17784, r17783, r17782, MPFR_RNDN);
        mpfr_sqrt(r17785, r17784, MPFR_RNDN);
        mpfr_sqrt(r17786, r17783, MPFR_RNDN);
        mpfr_add(r17787, r17785, r17786, MPFR_RNDN);
        mpfr_div(r17788, r17782, r17787, MPFR_RNDN);
        return mpfr_get_d(r17788, MPFR_RNDN);
}