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

char *name = "NMSE example 3.10";

double f_if(float x) {
        float r17750 = 1.0f;
        float r17751 = x;
        float r17752 = r17750 - r17751;
        float r17753 = log(r17752);
        float r17754 = r17750 + r17751;
        float r17755 = log(r17754);
        float r17756 = r17753 / r17755;
        return r17756;
}

double f_id(double x) {
        double r17757 = 1.0;
        double r17758 = x;
        double r17759 = r17757 - r17758;
        double r17760 = log(r17759);
        double r17761 = r17757 + r17758;
        double r17762 = log(r17761);
        double r17763 = r17760 / r17762;
        return r17763;
}


double f_of(float x) {
        float r17764 = 0.5f;
        float r17765 = x;
        float r17766 = r17765 * r17765;
        float r17767 = r17764 * r17766;
        float r17768 = 1.0f;
        float r17769 = r17768 + r17765;
        float r17770 = r17767 + r17769;
        float r17771 = -r17770;
        return r17771;
}

double f_od(double x) {
        double r17772 = 0.5;
        double r17773 = x;
        double r17774 = r17773 * r17773;
        double r17775 = r17772 * r17774;
        double r17776 = 1.0;
        double r17777 = r17776 + r17773;
        double r17778 = r17775 + r17777;
        double r17779 = -r17778;
        return r17779;
}

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 r17780, r17781, r17782, r17783, r17784, r17785, r17786;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17780, "1", 10, MPFR_RNDN);
        mpfr_init(r17781);
        mpfr_init(r17782);
        mpfr_init(r17783);
        mpfr_init(r17784);
        mpfr_init(r17785);
        mpfr_init(r17786);
}

double f_im(double x) {
        ;
        mpfr_set_d(r17781, x, MPFR_RNDN);
        mpfr_sub(r17782, r17780, r17781, MPFR_RNDN);
        mpfr_log(r17783, r17782, MPFR_RNDN);
        mpfr_add(r17784, r17780, r17781, MPFR_RNDN);
        mpfr_log(r17785, r17784, MPFR_RNDN);
        mpfr_div(r17786, r17783, r17785, MPFR_RNDN);
        return mpfr_get_d(r17786, MPFR_RNDN);
}

static mpfr_t r17787, r17788, r17789, r17790, r17791, r17792, r17793, r17794;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17787, "1/2", 10, MPFR_RNDN);
        mpfr_init(r17788);
        mpfr_init(r17789);
        mpfr_init(r17790);
        mpfr_init_set_str(r17791, "1", 10, MPFR_RNDN);
        mpfr_init(r17792);
        mpfr_init(r17793);
        mpfr_init(r17794);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r17788, x, MPFR_RNDN);
        mpfr_sqr(r17789, r17788, MPFR_RNDN);
        mpfr_mul(r17790, r17787, r17789, MPFR_RNDN);
        ;
        mpfr_add(r17792, r17791, r17788, MPFR_RNDN);
        mpfr_add(r17793, r17790, r17792, MPFR_RNDN);
        mpfr_neg(r17794, r17793, MPFR_RNDN);
        return mpfr_get_d(r17794, MPFR_RNDN);
}

static mpfr_t r17795, r17796, r17797, r17798, r17799, r17800, r17801, r17802;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init_set_str(r17795, "1/2", 10, MPFR_RNDN);
        mpfr_init(r17796);
        mpfr_init(r17797);
        mpfr_init(r17798);
        mpfr_init_set_str(r17799, "1", 10, MPFR_RNDN);
        mpfr_init(r17800);
        mpfr_init(r17801);
        mpfr_init(r17802);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r17796, x, MPFR_RNDN);
        mpfr_sqr(r17797, r17796, MPFR_RNDN);
        mpfr_mul(r17798, r17795, r17797, MPFR_RNDN);
        ;
        mpfr_add(r17800, r17799, r17796, MPFR_RNDN);
        mpfr_add(r17801, r17798, r17800, MPFR_RNDN);
        mpfr_neg(r17802, r17801, MPFR_RNDN);
        return mpfr_get_d(r17802, MPFR_RNDN);
}