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

char *name = "Jmat.Real.lambertw, newton loop step";

double f_if(float wj, float x) {
        float r18515 = wj;
        float r18516 = exp(r18515);
        float r18517 = r18515 * r18516;
        float r18518 = x;
        float r18519 = r18517 - r18518;
        float r18520 = r18516 + r18517;
        float r18521 = r18519 / r18520;
        float r18522 = r18515 - r18521;
        return r18522;
}

double f_id(double wj, double x) {
        double r18523 = wj;
        double r18524 = exp(r18523);
        double r18525 = r18523 * r18524;
        double r18526 = x;
        double r18527 = r18525 - r18526;
        double r18528 = r18524 + r18525;
        double r18529 = r18527 / r18528;
        double r18530 = r18523 - r18529;
        return r18530;
}


double f_of(float wj, float x) {
        float r18531 = wj;
        float r18532 = 1.0f;
        float r18533 = r18531 + r18532;
        float r18534 = r18531 / r18533;
        float r18535 = x;
        float r18536 = exp(r18531);
        float r18537 = fma(r18531, r18536, r18536);
        float r18538 = r18535 / r18537;
        float r18539 = r18534 - r18538;
        float r18540 = r18531 - r18539;
        return r18540;
}

double f_od(double wj, double x) {
        double r18541 = wj;
        double r18542 = 1.0;
        double r18543 = r18541 + r18542;
        double r18544 = r18541 / r18543;
        double r18545 = x;
        double r18546 = exp(r18541);
        double r18547 = fma(r18541, r18546, r18546);
        double r18548 = r18545 / r18547;
        double r18549 = r18544 - r18548;
        double r18550 = r18541 - r18549;
        return r18550;
}

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 r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18551);
        mpfr_init(r18552);
        mpfr_init(r18553);
        mpfr_init(r18554);
        mpfr_init(r18555);
        mpfr_init(r18556);
        mpfr_init(r18557);
        mpfr_init(r18558);
}

double f_im(double wj, double x) {
        mpfr_set_d(r18551, wj, MPFR_RNDN);
        mpfr_exp(r18552, r18551, MPFR_RNDN);
        mpfr_mul(r18553, r18551, r18552, MPFR_RNDN);
        mpfr_set_d(r18554, x, MPFR_RNDN);
        mpfr_sub(r18555, r18553, r18554, MPFR_RNDN);
        mpfr_add(r18556, r18552, r18553, MPFR_RNDN);
        mpfr_div(r18557, r18555, r18556, MPFR_RNDN);
        mpfr_sub(r18558, r18551, r18557, MPFR_RNDN);
        return mpfr_get_d(r18558, MPFR_RNDN);
}

static mpfr_t r18559, r18560, r18561, r18562, r18563, r18564, r18565, r18566, r18567, r18568;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18559);
        mpfr_init_set_str(r18560, "1", 10, MPFR_RNDN);
        mpfr_init(r18561);
        mpfr_init(r18562);
        mpfr_init(r18563);
        mpfr_init(r18564);
        mpfr_init(r18565);
        mpfr_init(r18566);
        mpfr_init(r18567);
        mpfr_init(r18568);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r18559, wj, MPFR_RNDN);
        ;
        mpfr_add(r18561, r18559, r18560, MPFR_RNDN);
        mpfr_div(r18562, r18559, r18561, MPFR_RNDN);
        mpfr_set_d(r18563, x, MPFR_RNDN);
        mpfr_exp(r18564, r18559, MPFR_RNDN);
        mpfr_fma(r18565, r18559, r18564, r18564, MPFR_RNDN);
        mpfr_div(r18566, r18563, r18565, MPFR_RNDN);
        mpfr_sub(r18567, r18562, r18566, MPFR_RNDN);
        mpfr_sub(r18568, r18559, r18567, MPFR_RNDN);
        return mpfr_get_d(r18568, MPFR_RNDN);
}

static mpfr_t r18569, r18570, r18571, r18572, r18573, r18574, r18575, r18576, r18577, r18578;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18569);
        mpfr_init_set_str(r18570, "1", 10, MPFR_RNDN);
        mpfr_init(r18571);
        mpfr_init(r18572);
        mpfr_init(r18573);
        mpfr_init(r18574);
        mpfr_init(r18575);
        mpfr_init(r18576);
        mpfr_init(r18577);
        mpfr_init(r18578);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r18569, wj, MPFR_RNDN);
        ;
        mpfr_add(r18571, r18569, r18570, MPFR_RNDN);
        mpfr_div(r18572, r18569, r18571, MPFR_RNDN);
        mpfr_set_d(r18573, x, MPFR_RNDN);
        mpfr_exp(r18574, r18569, MPFR_RNDN);
        mpfr_fma(r18575, r18569, r18574, r18574, MPFR_RNDN);
        mpfr_div(r18576, r18573, r18575, MPFR_RNDN);
        mpfr_sub(r18577, r18572, r18576, MPFR_RNDN);
        mpfr_sub(r18578, r18569, r18577, MPFR_RNDN);
        return mpfr_get_d(r18578, MPFR_RNDN);
}