#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 r25489 = wj;
        float r25490 = exp(r25489);
        float r25491 = r25489 * r25490;
        float r25492 = x;
        float r25493 = r25491 - r25492;
        float r25494 = r25490 + r25491;
        float r25495 = r25493 / r25494;
        float r25496 = r25489 - r25495;
        return r25496;
}

double f_id(double wj, double x) {
        double r25497 = wj;
        double r25498 = exp(r25497);
        double r25499 = r25497 * r25498;
        double r25500 = x;
        double r25501 = r25499 - r25500;
        double r25502 = r25498 + r25499;
        double r25503 = r25501 / r25502;
        double r25504 = r25497 - r25503;
        return r25504;
}


double f_of(float wj, float x) {
        float r25505 = x;
        float r25506 = wj;
        float r25507 = exp(r25506);
        float r25508 = fma(r25506, r25507, r25507);
        float r25509 = r25505 / r25508;
        float r25510 = r25506 * r25506;
        float r25511 = 1;
        float r25512 = r25511 - r25506;
        float r25513 = 4;
        float r25514 = pow(r25506, r25513);
        float r25515 = fma(r25510, r25512, r25514);
        float r25516 = r25509 + r25515;
        return r25516;
}

double f_od(double wj, double x) {
        double r25517 = x;
        double r25518 = wj;
        double r25519 = exp(r25518);
        double r25520 = fma(r25518, r25519, r25519);
        double r25521 = r25517 / r25520;
        double r25522 = r25518 * r25518;
        double r25523 = 1;
        double r25524 = r25523 - r25518;
        double r25525 = 4;
        double r25526 = pow(r25518, r25525);
        double r25527 = fma(r25522, r25524, r25526);
        double r25528 = r25521 + r25527;
        return r25528;
}

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 r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(912);
        mpfr_init(r25529);
        mpfr_init(r25530);
        mpfr_init(r25531);
        mpfr_init(r25532);
        mpfr_init(r25533);
        mpfr_init(r25534);
        mpfr_init(r25535);
        mpfr_init(r25536);
}

double f_im(double wj, double x) {
        mpfr_set_d(r25529, wj, MPFR_RNDN);
        mpfr_exp(r25530, r25529, MPFR_RNDN);
        mpfr_mul(r25531, r25529, r25530, MPFR_RNDN);
        mpfr_set_d(r25532, x, MPFR_RNDN);
        mpfr_sub(r25533, r25531, r25532, MPFR_RNDN);
        mpfr_add(r25534, r25530, r25531, MPFR_RNDN);
        mpfr_div(r25535, r25533, r25534, MPFR_RNDN);
        mpfr_sub(r25536, r25529, r25535, MPFR_RNDN);
        return mpfr_get_d(r25536, MPFR_RNDN);
}

static mpfr_t r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(912);
        mpfr_init(r25537);
        mpfr_init(r25538);
        mpfr_init(r25539);
        mpfr_init(r25540);
        mpfr_init(r25541);
        mpfr_init(r25542);
        mpfr_init_set_str(r25543, "1", 10, MPFR_RNDN);
        mpfr_init(r25544);
        mpfr_init_set_str(r25545, "4", 10, MPFR_RNDN);
        mpfr_init(r25546);
        mpfr_init(r25547);
        mpfr_init(r25548);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r25537, x, MPFR_RNDN);
        mpfr_set_d(r25538, wj, MPFR_RNDN);
        mpfr_exp(r25539, r25538, MPFR_RNDN);
        mpfr_fma(r25540, r25538, r25539, r25539, MPFR_RNDN);
        mpfr_div(r25541, r25537, r25540, MPFR_RNDN);
        mpfr_mul(r25542, r25538, r25538, MPFR_RNDN);
        ;
        mpfr_sub(r25544, r25543, r25538, MPFR_RNDN);
        ;
        mpfr_pow(r25546, r25538, r25545, MPFR_RNDN);
        mpfr_fma(r25547, r25542, r25544, r25546, MPFR_RNDN);
        mpfr_add(r25548, r25541, r25547, MPFR_RNDN);
        return mpfr_get_d(r25548, MPFR_RNDN);
}

static mpfr_t r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(912);
        mpfr_init(r25549);
        mpfr_init(r25550);
        mpfr_init(r25551);
        mpfr_init(r25552);
        mpfr_init(r25553);
        mpfr_init(r25554);
        mpfr_init_set_str(r25555, "1", 10, MPFR_RNDN);
        mpfr_init(r25556);
        mpfr_init_set_str(r25557, "4", 10, MPFR_RNDN);
        mpfr_init(r25558);
        mpfr_init(r25559);
        mpfr_init(r25560);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r25549, x, MPFR_RNDN);
        mpfr_set_d(r25550, wj, MPFR_RNDN);
        mpfr_exp(r25551, r25550, MPFR_RNDN);
        mpfr_fma(r25552, r25550, r25551, r25551, MPFR_RNDN);
        mpfr_div(r25553, r25549, r25552, MPFR_RNDN);
        mpfr_mul(r25554, r25550, r25550, MPFR_RNDN);
        ;
        mpfr_sub(r25556, r25555, r25550, MPFR_RNDN);
        ;
        mpfr_pow(r25558, r25550, r25557, MPFR_RNDN);
        mpfr_fma(r25559, r25554, r25556, r25558, MPFR_RNDN);
        mpfr_add(r25560, r25553, r25559, MPFR_RNDN);
        return mpfr_get_d(r25560, MPFR_RNDN);
}