#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 r25484 = wj;
        float r25485 = exp(r25484);
        float r25486 = r25484 * r25485;
        float r25487 = x;
        float r25488 = r25486 - r25487;
        float r25489 = r25485 + r25486;
        float r25490 = r25488 / r25489;
        float r25491 = r25484 - r25490;
        return r25491;
}

double f_id(double wj, double x) {
        double r25492 = wj;
        double r25493 = exp(r25492);
        double r25494 = r25492 * r25493;
        double r25495 = x;
        double r25496 = r25494 - r25495;
        double r25497 = r25493 + r25494;
        double r25498 = r25496 / r25497;
        double r25499 = r25492 - r25498;
        return r25499;
}


double f_of(float wj, float x) {
        float r25500 = x;
        float r25501 = wj;
        float r25502 = exp(r25501);
        float r25503 = r25500 / r25502;
        float r25504 = 1;
        float r25505 = r25504 + r25501;
        float r25506 = r25503 / r25505;
        float r25507 = 4;
        float r25508 = pow(r25501, r25507);
        float r25509 = r25504 - r25501;
        float r25510 = r25501 * r25501;
        float r25511 = r25509 * r25510;
        float r25512 = r25508 + r25511;
        float r25513 = r25506 + r25512;
        return r25513;
}

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

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

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

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

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

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

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

static mpfr_t r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563;

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

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