#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 r26530 = x;
        float r26531 = log(r26530);
        float r26532 = log(r26531);
        float r26533 = r26531 - r26532;
        return r26533;
}

double f_id(double x) {
        double r26534 = x;
        double r26535 = log(r26534);
        double r26536 = log(r26535);
        double r26537 = r26535 - r26536;
        return r26537;
}


double f_of(float x) {
        float r26538 = x;
        float r26539 = log(r26538);
        float r26540 = r26538 / r26539;
        float r26541 = log(r26540);
        return r26541;
}

double f_od(double x) {
        double r26542 = x;
        double r26543 = log(r26542);
        double r26544 = r26542 / r26543;
        double r26545 = log(r26544);
        return r26545;
}

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 r26546, r26547, r26548, r26549;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26546);
        mpfr_init(r26547);
        mpfr_init(r26548);
        mpfr_init(r26549);
}

double f_im(double x) {
        mpfr_set_d(r26546, x, MPFR_RNDN);
        mpfr_log(r26547, r26546, MPFR_RNDN);
        mpfr_log(r26548, r26547, MPFR_RNDN);
        mpfr_sub(r26549, r26547, r26548, MPFR_RNDN);
        return mpfr_get_d(r26549, MPFR_RNDN);
}

static mpfr_t r26550, r26551, r26552, r26553;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26550);
        mpfr_init(r26551);
        mpfr_init(r26552);
        mpfr_init(r26553);
}

double f_fm(double x) {
        mpfr_set_d(r26550, x, MPFR_RNDN);
        mpfr_log(r26551, r26550, MPFR_RNDN);
        mpfr_div(r26552, r26550, r26551, MPFR_RNDN);
        mpfr_log(r26553, r26552, MPFR_RNDN);
        return mpfr_get_d(r26553, MPFR_RNDN);
}

static mpfr_t r26554, r26555, r26556, r26557;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26554);
        mpfr_init(r26555);
        mpfr_init(r26556);
        mpfr_init(r26557);
}

double f_dm(double x) {
        mpfr_set_d(r26554, x, MPFR_RNDN);
        mpfr_log(r26555, r26554, MPFR_RNDN);
        mpfr_div(r26556, r26554, r26555, MPFR_RNDN);
        mpfr_log(r26557, r26556, MPFR_RNDN);
        return mpfr_get_d(r26557, MPFR_RNDN);
}