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

char *name = "Jmat.Real.lambertw, estimator";

double f_if(float x) {
        float r17775 = x;
        float r17776 = log(r17775);
        float r17777 = log(r17776);
        float r17778 = r17776 - r17777;
        return r17778;
}

double f_id(double x) {
        double r17779 = x;
        double r17780 = log(r17779);
        double r17781 = log(r17780);
        double r17782 = r17780 - r17781;
        return r17782;
}


double f_of(float x) {
        float r17783 = x;
        float r17784 = log(r17783);
        float r17785 = r17783 / r17784;
        float r17786 = log(r17785);
        return r17786;
}

double f_od(double x) {
        double r17787 = x;
        double r17788 = log(r17787);
        double r17789 = r17787 / r17788;
        double r17790 = log(r17789);
        return r17790;
}

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 r17791, r17792, r17793, r17794;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r17791);
        mpfr_init(r17792);
        mpfr_init(r17793);
        mpfr_init(r17794);
}

double f_im(double x) {
        mpfr_set_d(r17791, x, MPFR_RNDN);
        mpfr_log(r17792, r17791, MPFR_RNDN);
        mpfr_log(r17793, r17792, MPFR_RNDN);
        mpfr_sub(r17794, r17792, r17793, MPFR_RNDN);
        return mpfr_get_d(r17794, MPFR_RNDN);
}

static mpfr_t r17795, r17796, r17797, r17798;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17795);
        mpfr_init(r17796);
        mpfr_init(r17797);
        mpfr_init(r17798);
}

double f_fm(double x) {
        mpfr_set_d(r17795, x, MPFR_RNDN);
        mpfr_log(r17796, r17795, MPFR_RNDN);
        mpfr_div(r17797, r17795, r17796, MPFR_RNDN);
        mpfr_log(r17798, r17797, MPFR_RNDN);
        return mpfr_get_d(r17798, MPFR_RNDN);
}

static mpfr_t r17799, r17800, r17801, r17802;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r17799);
        mpfr_init(r17800);
        mpfr_init(r17801);
        mpfr_init(r17802);
}

double f_dm(double x) {
        mpfr_set_d(r17799, x, MPFR_RNDN);
        mpfr_log(r17800, r17799, MPFR_RNDN);
        mpfr_div(r17801, r17799, r17800, MPFR_RNDN);
        mpfr_log(r17802, r17801, MPFR_RNDN);
        return mpfr_get_d(r17802, MPFR_RNDN);
}