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

char *name = "logs (example 3.8)";

double f_if(float n) {
        float r20749 = n;
        float r20750 = 1;
        float r20751 = r20749 + r20750;
        float r20752 = log(r20751);
        float r20753 = r20751 * r20752;
        float r20754 = log(r20749);
        float r20755 = r20749 * r20754;
        float r20756 = r20753 - r20755;
        float r20757 = r20756 - r20750;
        return r20757;
}

double f_id(double n) {
        double r20758 = n;
        double r20759 = 1;
        double r20760 = r20758 + r20759;
        double r20761 = log(r20760);
        double r20762 = r20760 * r20761;
        double r20763 = log(r20758);
        double r20764 = r20758 * r20763;
        double r20765 = r20762 - r20764;
        double r20766 = r20765 - r20759;
        return r20766;
}


double f_of(float n) {
        float r20767 = 1;
        float r20768 = n;
        float r20769 = r20767 + r20768;
        float r20770 = log(r20769);
        float r20771 = r20770 * r20769;
        float r20772 = -r20768;
        float r20773 = log(r20768);
        float r20774 = r20772 * r20773;
        float r20775 = r20774 + r20767;
        float r20776 = r20771 - r20775;
        return r20776;
}

double f_od(double n) {
        double r20777 = 1;
        double r20778 = n;
        double r20779 = r20777 + r20778;
        double r20780 = log(r20779);
        double r20781 = r20780 * r20779;
        double r20782 = -r20778;
        double r20783 = log(r20778);
        double r20784 = r20782 * r20783;
        double r20785 = r20784 + r20777;
        double r20786 = r20781 - r20785;
        return r20786;
}

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 r20787, r20788, r20789, r20790, r20791, r20792, r20793, r20794, r20795;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20787);
        mpfr_init_set_str(r20788, "1", 10, MPFR_RNDN);
        mpfr_init(r20789);
        mpfr_init(r20790);
        mpfr_init(r20791);
        mpfr_init(r20792);
        mpfr_init(r20793);
        mpfr_init(r20794);
        mpfr_init(r20795);
}

double f_im(double n) {
        mpfr_set_d(r20787, n, MPFR_RNDN);
        ;
        mpfr_add(r20789, r20787, r20788, MPFR_RNDN);
        mpfr_log(r20790, r20789, MPFR_RNDN);
        mpfr_mul(r20791, r20789, r20790, MPFR_RNDN);
        mpfr_log(r20792, r20787, MPFR_RNDN);
        mpfr_mul(r20793, r20787, r20792, MPFR_RNDN);
        mpfr_sub(r20794, r20791, r20793, MPFR_RNDN);
        mpfr_sub(r20795, r20794, r20788, MPFR_RNDN);
        return mpfr_get_d(r20795, MPFR_RNDN);
}

static mpfr_t r20796, r20797, r20798, r20799, r20800, r20801, r20802, r20803, r20804, r20805;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20796, "1", 10, MPFR_RNDN);
        mpfr_init(r20797);
        mpfr_init(r20798);
        mpfr_init(r20799);
        mpfr_init(r20800);
        mpfr_init(r20801);
        mpfr_init(r20802);
        mpfr_init(r20803);
        mpfr_init(r20804);
        mpfr_init(r20805);
}

double f_fm(double n) {
        ;
        mpfr_set_d(r20797, n, MPFR_RNDN);
        mpfr_add(r20798, r20796, r20797, MPFR_RNDN);
        mpfr_log(r20799, r20798, MPFR_RNDN);
        mpfr_mul(r20800, r20799, r20798, MPFR_RNDN);
        mpfr_neg(r20801, r20797, MPFR_RNDN);
        mpfr_log(r20802, r20797, MPFR_RNDN);
        mpfr_mul(r20803, r20801, r20802, MPFR_RNDN);
        mpfr_add(r20804, r20803, r20796, MPFR_RNDN);
        mpfr_sub(r20805, r20800, r20804, MPFR_RNDN);
        return mpfr_get_d(r20805, MPFR_RNDN);
}

static mpfr_t r20806, r20807, r20808, r20809, r20810, r20811, r20812, r20813, r20814, r20815;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r20806, "1", 10, MPFR_RNDN);
        mpfr_init(r20807);
        mpfr_init(r20808);
        mpfr_init(r20809);
        mpfr_init(r20810);
        mpfr_init(r20811);
        mpfr_init(r20812);
        mpfr_init(r20813);
        mpfr_init(r20814);
        mpfr_init(r20815);
}

double f_dm(double n) {
        ;
        mpfr_set_d(r20807, n, MPFR_RNDN);
        mpfr_add(r20808, r20806, r20807, MPFR_RNDN);
        mpfr_log(r20809, r20808, MPFR_RNDN);
        mpfr_mul(r20810, r20809, r20808, MPFR_RNDN);
        mpfr_neg(r20811, r20807, MPFR_RNDN);
        mpfr_log(r20812, r20807, MPFR_RNDN);
        mpfr_mul(r20813, r20811, r20812, MPFR_RNDN);
        mpfr_add(r20814, r20813, r20806, MPFR_RNDN);
        mpfr_sub(r20815, r20810, r20814, MPFR_RNDN);
        return mpfr_get_d(r20815, MPFR_RNDN);
}