#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 r20344 = wj;
        float r20345 = exp(r20344);
        float r20346 = r20344 * r20345;
        float r20347 = x;
        float r20348 = r20346 - r20347;
        float r20349 = r20345 + r20346;
        float r20350 = r20348 / r20349;
        float r20351 = r20344 - r20350;
        return r20351;
}

double f_id(double wj, double x) {
        double r20352 = wj;
        double r20353 = exp(r20352);
        double r20354 = r20352 * r20353;
        double r20355 = x;
        double r20356 = r20354 - r20355;
        double r20357 = r20353 + r20354;
        double r20358 = r20356 / r20357;
        double r20359 = r20352 - r20358;
        return r20359;
}


double f_of(float wj, float x) {
        float r20360 = wj;
        float r20361 = r20360 * r20360;
        float r20362 = r20360 * (r20360 * r20360);
        float r20363 = r20361 - r20362;
        float r20364 = x;
        float r20365 = 1.0f;
        float r20366 = r20360 + r20365;
        float r20367 = exp(r20360);
        float r20368 = r20366 * r20367;
        float r20369 = r20364 / r20368;
        float r20370 = r20363 + r20369;
        return r20370;
}

double f_od(double wj, double x) {
        double r20371 = wj;
        double r20372 = r20371 * r20371;
        double r20373 = r20371 * (r20371 * r20371);
        double r20374 = r20372 - r20373;
        double r20375 = x;
        double r20376 = 1.0;
        double r20377 = r20371 + r20376;
        double r20378 = exp(r20371);
        double r20379 = r20377 * r20378;
        double r20380 = r20375 / r20379;
        double r20381 = r20374 + r20380;
        return r20381;
}

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 r20382, r20383, r20384, r20385, r20386, r20387, r20388, r20389;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20382);
        mpfr_init(r20383);
        mpfr_init(r20384);
        mpfr_init(r20385);
        mpfr_init(r20386);
        mpfr_init(r20387);
        mpfr_init(r20388);
        mpfr_init(r20389);
}

double f_im(double wj, double x) {
        mpfr_set_d(r20382, wj, MPFR_RNDN);
        mpfr_exp(r20383, r20382, MPFR_RNDN);
        mpfr_mul(r20384, r20382, r20383, MPFR_RNDN);
        mpfr_set_d(r20385, x, MPFR_RNDN);
        mpfr_sub(r20386, r20384, r20385, MPFR_RNDN);
        mpfr_add(r20387, r20383, r20384, MPFR_RNDN);
        mpfr_div(r20388, r20386, r20387, MPFR_RNDN);
        mpfr_sub(r20389, r20382, r20388, MPFR_RNDN);
        return mpfr_get_d(r20389, MPFR_RNDN);
}

static mpfr_t r20390, r20391, r20392, r20393, r20394, r20395, r20396, r20397, r20398, r20399, r20400;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20390);
        mpfr_init(r20391);
        mpfr_init(r20392);
        mpfr_init(r20393);
        mpfr_init(r20394);
        mpfr_init_set_str(r20395, "1", 10, MPFR_RNDN);
        mpfr_init(r20396);
        mpfr_init(r20397);
        mpfr_init(r20398);
        mpfr_init(r20399);
        mpfr_init(r20400);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r20390, wj, MPFR_RNDN);
        mpfr_sqr(r20391, r20390, MPFR_RNDN);
        mpfr_mul(r20392, r20390, r20390, MPFR_RNDN); mpfr_mul(r20392, r20392, r20390, MPFR_RNDN);
        mpfr_sub(r20393, r20391, r20392, MPFR_RNDN);
        mpfr_set_d(r20394, x, MPFR_RNDN);
        ;
        mpfr_add(r20396, r20390, r20395, MPFR_RNDN);
        mpfr_exp(r20397, r20390, MPFR_RNDN);
        mpfr_mul(r20398, r20396, r20397, MPFR_RNDN);
        mpfr_div(r20399, r20394, r20398, MPFR_RNDN);
        mpfr_add(r20400, r20393, r20399, MPFR_RNDN);
        return mpfr_get_d(r20400, MPFR_RNDN);
}

static mpfr_t r20401, r20402, r20403, r20404, r20405, r20406, r20407, r20408, r20409, r20410, r20411;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20401);
        mpfr_init(r20402);
        mpfr_init(r20403);
        mpfr_init(r20404);
        mpfr_init(r20405);
        mpfr_init_set_str(r20406, "1", 10, MPFR_RNDN);
        mpfr_init(r20407);
        mpfr_init(r20408);
        mpfr_init(r20409);
        mpfr_init(r20410);
        mpfr_init(r20411);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r20401, wj, MPFR_RNDN);
        mpfr_sqr(r20402, r20401, MPFR_RNDN);
        mpfr_mul(r20403, r20401, r20401, MPFR_RNDN); mpfr_mul(r20403, r20403, r20401, MPFR_RNDN);
        mpfr_sub(r20404, r20402, r20403, MPFR_RNDN);
        mpfr_set_d(r20405, x, MPFR_RNDN);
        ;
        mpfr_add(r20407, r20401, r20406, MPFR_RNDN);
        mpfr_exp(r20408, r20401, MPFR_RNDN);
        mpfr_mul(r20409, r20407, r20408, MPFR_RNDN);
        mpfr_div(r20410, r20405, r20409, MPFR_RNDN);
        mpfr_add(r20411, r20404, r20410, MPFR_RNDN);
        return mpfr_get_d(r20411, MPFR_RNDN);
}