#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 r20331 = wj;
        float r20332 = exp(r20331);
        float r20333 = r20331 * r20332;
        float r20334 = x;
        float r20335 = r20333 - r20334;
        float r20336 = r20332 + r20333;
        float r20337 = r20335 / r20336;
        float r20338 = r20331 - r20337;
        return r20338;
}

double f_id(double wj, double x) {
        double r20339 = wj;
        double r20340 = exp(r20339);
        double r20341 = r20339 * r20340;
        double r20342 = x;
        double r20343 = r20341 - r20342;
        double r20344 = r20340 + r20341;
        double r20345 = r20343 / r20344;
        double r20346 = r20339 - r20345;
        return r20346;
}


double f_of(float wj, float x) {
        float r20347 = wj;
        float r20348 = 4.0f;
        float r20349 = pow(r20347, r20348);
        float r20350 = r20347 * r20347;
        float r20351 = r20347 * (r20347 * r20347);
        float r20352 = r20350 - r20351;
        float r20353 = r20349 + r20352;
        float r20354 = x;
        float r20355 = exp(r20347);
        float r20356 = 1.0f;
        float r20357 = r20356 + r20347;
        float r20358 = r20355 * r20357;
        float r20359 = r20354 / r20358;
        float r20360 = r20353 + r20359;
        return r20360;
}

double f_od(double wj, double x) {
        double r20361 = wj;
        double r20362 = 4.0;
        double r20363 = pow(r20361, r20362);
        double r20364 = r20361 * r20361;
        double r20365 = r20361 * (r20361 * r20361);
        double r20366 = r20364 - r20365;
        double r20367 = r20363 + r20366;
        double r20368 = x;
        double r20369 = exp(r20361);
        double r20370 = 1.0;
        double r20371 = r20370 + r20361;
        double r20372 = r20369 * r20371;
        double r20373 = r20368 / r20372;
        double r20374 = r20367 + r20373;
        return r20374;
}

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 r20375, r20376, r20377, r20378, r20379, r20380, r20381, r20382;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20375);
        mpfr_init(r20376);
        mpfr_init(r20377);
        mpfr_init(r20378);
        mpfr_init(r20379);
        mpfr_init(r20380);
        mpfr_init(r20381);
        mpfr_init(r20382);
}

double f_im(double wj, double x) {
        mpfr_set_d(r20375, wj, MPFR_RNDN);
        mpfr_exp(r20376, r20375, MPFR_RNDN);
        mpfr_mul(r20377, r20375, r20376, MPFR_RNDN);
        mpfr_set_d(r20378, x, MPFR_RNDN);
        mpfr_sub(r20379, r20377, r20378, MPFR_RNDN);
        mpfr_add(r20380, r20376, r20377, MPFR_RNDN);
        mpfr_div(r20381, r20379, r20380, MPFR_RNDN);
        mpfr_sub(r20382, r20375, r20381, MPFR_RNDN);
        return mpfr_get_d(r20382, MPFR_RNDN);
}

static mpfr_t r20383, r20384, r20385, r20386, r20387, r20388, r20389, r20390, r20391, r20392, r20393, r20394, r20395, r20396;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20383);
        mpfr_init_set_str(r20384, "4", 10, MPFR_RNDN);
        mpfr_init(r20385);
        mpfr_init(r20386);
        mpfr_init(r20387);
        mpfr_init(r20388);
        mpfr_init(r20389);
        mpfr_init(r20390);
        mpfr_init(r20391);
        mpfr_init_set_str(r20392, "1", 10, MPFR_RNDN);
        mpfr_init(r20393);
        mpfr_init(r20394);
        mpfr_init(r20395);
        mpfr_init(r20396);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r20383, wj, MPFR_RNDN);
        ;
        mpfr_pow(r20385, r20383, r20384, MPFR_RNDN);
        mpfr_sqr(r20386, r20383, MPFR_RNDN);
        mpfr_mul(r20387, r20383, r20383, MPFR_RNDN); mpfr_mul(r20387, r20387, r20383, MPFR_RNDN);
        mpfr_sub(r20388, r20386, r20387, MPFR_RNDN);
        mpfr_add(r20389, r20385, r20388, MPFR_RNDN);
        mpfr_set_d(r20390, x, MPFR_RNDN);
        mpfr_exp(r20391, r20383, MPFR_RNDN);
        ;
        mpfr_add(r20393, r20392, r20383, MPFR_RNDN);
        mpfr_mul(r20394, r20391, r20393, MPFR_RNDN);
        mpfr_div(r20395, r20390, r20394, MPFR_RNDN);
        mpfr_add(r20396, r20389, r20395, MPFR_RNDN);
        return mpfr_get_d(r20396, MPFR_RNDN);
}

static mpfr_t r20397, r20398, r20399, r20400, r20401, r20402, r20403, r20404, r20405, r20406, r20407, r20408, r20409, r20410;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20397);
        mpfr_init_set_str(r20398, "4", 10, MPFR_RNDN);
        mpfr_init(r20399);
        mpfr_init(r20400);
        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);
}

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