#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 r26430 = x;
        float r26431 = log(r26430);
        float r26432 = log(r26431);
        float r26433 = r26431 - r26432;
        return r26433;
}

double f_id(double x) {
        double r26434 = x;
        double r26435 = log(r26434);
        double r26436 = log(r26435);
        double r26437 = r26435 - r26436;
        return r26437;
}


double f_of(float x) {
        float r26438 = x;
        float r26439 = log(r26438);
        float r26440 = r26438 / r26439;
        float r26441 = log(r26440);
        return r26441;
}

double f_od(double x) {
        double r26442 = x;
        double r26443 = log(r26442);
        double r26444 = r26442 / r26443;
        double r26445 = log(r26444);
        return r26445;
}

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 r26446, r26447, r26448, r26449;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26446);
        mpfr_init(r26447);
        mpfr_init(r26448);
        mpfr_init(r26449);
}

double f_im(double x) {
        mpfr_set_d(r26446, x, MPFR_RNDN);
        mpfr_log(r26447, r26446, MPFR_RNDN);
        mpfr_log(r26448, r26447, MPFR_RNDN);
        mpfr_sub(r26449, r26447, r26448, MPFR_RNDN);
        return mpfr_get_d(r26449, MPFR_RNDN);
}

static mpfr_t r26450, r26451, r26452, r26453;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26450);
        mpfr_init(r26451);
        mpfr_init(r26452);
        mpfr_init(r26453);
}

double f_fm(double x) {
        mpfr_set_d(r26450, x, MPFR_RNDN);
        mpfr_log(r26451, r26450, MPFR_RNDN);
        mpfr_div(r26452, r26450, r26451, MPFR_RNDN);
        mpfr_log(r26453, r26452, MPFR_RNDN);
        return mpfr_get_d(r26453, MPFR_RNDN);
}

static mpfr_t r26454, r26455, r26456, r26457;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26454);
        mpfr_init(r26455);
        mpfr_init(r26456);
        mpfr_init(r26457);
}

double f_dm(double x) {
        mpfr_set_d(r26454, x, MPFR_RNDN);
        mpfr_log(r26455, r26454, MPFR_RNDN);
        mpfr_div(r26456, r26454, r26455, MPFR_RNDN);
        mpfr_log(r26457, r26456, MPFR_RNDN);
        return mpfr_get_d(r26457, MPFR_RNDN);
}