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

char *name = "qlog (example 3.10)";

double f_if(float x) {
        float r19728 = 1.0f;
        float r19729 = x;
        float r19730 = r19728 - r19729;
        float r19731 = log(r19730);
        float r19732 = r19728 + r19729;
        float r19733 = log(r19732);
        float r19734 = r19731 / r19733;
        return r19734;
}

double f_id(double x) {
        double r19735 = 1.0;
        double r19736 = x;
        double r19737 = r19735 - r19736;
        double r19738 = log(r19737);
        double r19739 = r19735 + r19736;
        double r19740 = log(r19739);
        double r19741 = r19738 / r19740;
        return r19741;
}


double f_of(float x) {
        float r19742 = 0.5f;
        float r19743 = x;
        float r19744 = r19743 * r19743;
        float r19745 = r19742 * r19744;
        float r19746 = 1.0f;
        float r19747 = r19746 + r19743;
        float r19748 = r19745 + r19747;
        float r19749 = -r19748;
        return r19749;
}

double f_od(double x) {
        double r19750 = 0.5;
        double r19751 = x;
        double r19752 = r19751 * r19751;
        double r19753 = r19750 * r19752;
        double r19754 = 1.0;
        double r19755 = r19754 + r19751;
        double r19756 = r19753 + r19755;
        double r19757 = -r19756;
        return r19757;
}

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 r19758, r19759, r19760, r19761, r19762, r19763, r19764;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r19758, "1", 10, MPFR_RNDN);
        mpfr_init(r19759);
        mpfr_init(r19760);
        mpfr_init(r19761);
        mpfr_init(r19762);
        mpfr_init(r19763);
        mpfr_init(r19764);
}

double f_im(double x) {
        ;
        mpfr_set_d(r19759, x, MPFR_RNDN);
        mpfr_sub(r19760, r19758, r19759, MPFR_RNDN);
        mpfr_log(r19761, r19760, MPFR_RNDN);
        mpfr_add(r19762, r19758, r19759, MPFR_RNDN);
        mpfr_log(r19763, r19762, MPFR_RNDN);
        mpfr_div(r19764, r19761, r19763, MPFR_RNDN);
        return mpfr_get_d(r19764, MPFR_RNDN);
}

static mpfr_t r19765, r19766, r19767, r19768, r19769, r19770, r19771, r19772;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r19765, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19766);
        mpfr_init(r19767);
        mpfr_init(r19768);
        mpfr_init_set_str(r19769, "1", 10, MPFR_RNDN);
        mpfr_init(r19770);
        mpfr_init(r19771);
        mpfr_init(r19772);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r19766, x, MPFR_RNDN);
        mpfr_sqr(r19767, r19766, MPFR_RNDN);
        mpfr_mul(r19768, r19765, r19767, MPFR_RNDN);
        ;
        mpfr_add(r19770, r19769, r19766, MPFR_RNDN);
        mpfr_add(r19771, r19768, r19770, MPFR_RNDN);
        mpfr_neg(r19772, r19771, MPFR_RNDN);
        return mpfr_get_d(r19772, MPFR_RNDN);
}

static mpfr_t r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r19773, "1/2", 10, MPFR_RNDN);
        mpfr_init(r19774);
        mpfr_init(r19775);
        mpfr_init(r19776);
        mpfr_init_set_str(r19777, "1", 10, MPFR_RNDN);
        mpfr_init(r19778);
        mpfr_init(r19779);
        mpfr_init(r19780);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r19774, x, MPFR_RNDN);
        mpfr_sqr(r19775, r19774, MPFR_RNDN);
        mpfr_mul(r19776, r19773, r19775, MPFR_RNDN);
        ;
        mpfr_add(r19778, r19777, r19774, MPFR_RNDN);
        mpfr_add(r19779, r19776, r19778, MPFR_RNDN);
        mpfr_neg(r19780, r19779, MPFR_RNDN);
        return mpfr_get_d(r19780, MPFR_RNDN);
}