#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 r18435 = wj;
        float r18436 = exp(r18435);
        float r18437 = r18435 * r18436;
        float r18438 = x;
        float r18439 = r18437 - r18438;
        float r18440 = r18436 + r18437;
        float r18441 = r18439 / r18440;
        float r18442 = r18435 - r18441;
        return r18442;
}

double f_id(double wj, double x) {
        double r18443 = wj;
        double r18444 = exp(r18443);
        double r18445 = r18443 * r18444;
        double r18446 = x;
        double r18447 = r18445 - r18446;
        double r18448 = r18444 + r18445;
        double r18449 = r18447 / r18448;
        double r18450 = r18443 - r18449;
        return r18450;
}


double f_of(float wj, float x) {
        float r18451 = wj;
        float r18452 = 1.0f;
        float r18453 = r18451 + r18452;
        float r18454 = r18451 / r18453;
        float r18455 = x;
        float r18456 = exp(r18451);
        float r18457 = fma(r18451, r18456, r18456);
        float r18458 = r18455 / r18457;
        float r18459 = r18454 - r18458;
        float r18460 = r18451 - r18459;
        return r18460;
}

double f_od(double wj, double x) {
        double r18461 = wj;
        double r18462 = 1.0;
        double r18463 = r18461 + r18462;
        double r18464 = r18461 / r18463;
        double r18465 = x;
        double r18466 = exp(r18461);
        double r18467 = fma(r18461, r18466, r18466);
        double r18468 = r18465 / r18467;
        double r18469 = r18464 - r18468;
        double r18470 = r18461 - r18469;
        return r18470;
}

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 r18471, r18472, r18473, r18474, r18475, r18476, r18477, r18478;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r18471);
        mpfr_init(r18472);
        mpfr_init(r18473);
        mpfr_init(r18474);
        mpfr_init(r18475);
        mpfr_init(r18476);
        mpfr_init(r18477);
        mpfr_init(r18478);
}

double f_im(double wj, double x) {
        mpfr_set_d(r18471, wj, MPFR_RNDN);
        mpfr_exp(r18472, r18471, MPFR_RNDN);
        mpfr_mul(r18473, r18471, r18472, MPFR_RNDN);
        mpfr_set_d(r18474, x, MPFR_RNDN);
        mpfr_sub(r18475, r18473, r18474, MPFR_RNDN);
        mpfr_add(r18476, r18472, r18473, MPFR_RNDN);
        mpfr_div(r18477, r18475, r18476, MPFR_RNDN);
        mpfr_sub(r18478, r18471, r18477, MPFR_RNDN);
        return mpfr_get_d(r18478, MPFR_RNDN);
}

static mpfr_t r18479, r18480, r18481, r18482, r18483, r18484, r18485, r18486, r18487, r18488;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18479);
        mpfr_init_set_str(r18480, "1", 10, MPFR_RNDN);
        mpfr_init(r18481);
        mpfr_init(r18482);
        mpfr_init(r18483);
        mpfr_init(r18484);
        mpfr_init(r18485);
        mpfr_init(r18486);
        mpfr_init(r18487);
        mpfr_init(r18488);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r18479, wj, MPFR_RNDN);
        ;
        mpfr_add(r18481, r18479, r18480, MPFR_RNDN);
        mpfr_div(r18482, r18479, r18481, MPFR_RNDN);
        mpfr_set_d(r18483, x, MPFR_RNDN);
        mpfr_exp(r18484, r18479, MPFR_RNDN);
        mpfr_fma(r18485, r18479, r18484, r18484, MPFR_RNDN);
        mpfr_div(r18486, r18483, r18485, MPFR_RNDN);
        mpfr_sub(r18487, r18482, r18486, MPFR_RNDN);
        mpfr_sub(r18488, r18479, r18487, MPFR_RNDN);
        return mpfr_get_d(r18488, MPFR_RNDN);
}

static mpfr_t r18489, r18490, r18491, r18492, r18493, r18494, r18495, r18496, r18497, r18498;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r18489);
        mpfr_init_set_str(r18490, "1", 10, MPFR_RNDN);
        mpfr_init(r18491);
        mpfr_init(r18492);
        mpfr_init(r18493);
        mpfr_init(r18494);
        mpfr_init(r18495);
        mpfr_init(r18496);
        mpfr_init(r18497);
        mpfr_init(r18498);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r18489, wj, MPFR_RNDN);
        ;
        mpfr_add(r18491, r18489, r18490, MPFR_RNDN);
        mpfr_div(r18492, r18489, r18491, MPFR_RNDN);
        mpfr_set_d(r18493, x, MPFR_RNDN);
        mpfr_exp(r18494, r18489, MPFR_RNDN);
        mpfr_fma(r18495, r18489, r18494, r18494, MPFR_RNDN);
        mpfr_div(r18496, r18493, r18495, MPFR_RNDN);
        mpfr_sub(r18497, r18492, r18496, MPFR_RNDN);
        mpfr_sub(r18498, r18489, r18497, MPFR_RNDN);
        return mpfr_get_d(r18498, MPFR_RNDN);
}