#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 r26418 = wj;
        float r26419 = exp(r26418);
        float r26420 = r26418 * r26419;
        float r26421 = x;
        float r26422 = r26420 - r26421;
        float r26423 = r26419 + r26420;
        float r26424 = r26422 / r26423;
        float r26425 = r26418 - r26424;
        return r26425;
}

double f_id(double wj, double x) {
        double r26426 = wj;
        double r26427 = exp(r26426);
        double r26428 = r26426 * r26427;
        double r26429 = x;
        double r26430 = r26428 - r26429;
        double r26431 = r26427 + r26428;
        double r26432 = r26430 / r26431;
        double r26433 = r26426 - r26432;
        return r26433;
}


double f_of(float wj, float x) {
        float r26434 = wj;
        float r26435 = 2;
        float r26436 = pow(r26434, r26435);
        float r26437 = x;
        float r26438 = r26436 + r26437;
        float r26439 = r26434 * r26437;
        float r26440 = r26435 * r26439;
        float r26441 = r26438 - r26440;
        float r26442 = 1.1273399280310297e-25;
        bool r26443 = r26441 <= r26442;
        float r26444 = 1;
        float r26445 = r26444 + r26434;
        float r26446 = r26434 / r26445;
        float r26447 = r26434 - r26446;
        float r26448 = cbrt(r26447);
        float r26449 = r26448 * r26448;
        float r26450 = cbrt(r26448);
        float r26451 = r26450 * r26450;
        float r26452 = r26451 * r26450;
        float r26453 = r26449 * r26452;
        float r26454 = exp(r26434);
        float r26455 = r26434 * r26454;
        float r26456 = r26454 + r26455;
        float r26457 = r26437 / r26456;
        float r26458 = r26453 + r26457;
        float r26459 = r26443 ? r26441 : r26458;
        return r26459;
}

double f_od(double wj, double x) {
        double r26460 = wj;
        double r26461 = 2;
        double r26462 = pow(r26460, r26461);
        double r26463 = x;
        double r26464 = r26462 + r26463;
        double r26465 = r26460 * r26463;
        double r26466 = r26461 * r26465;
        double r26467 = r26464 - r26466;
        double r26468 = 1.1273399280310297e-25;
        bool r26469 = r26467 <= r26468;
        double r26470 = 1;
        double r26471 = r26470 + r26460;
        double r26472 = r26460 / r26471;
        double r26473 = r26460 - r26472;
        double r26474 = cbrt(r26473);
        double r26475 = r26474 * r26474;
        double r26476 = cbrt(r26474);
        double r26477 = r26476 * r26476;
        double r26478 = r26477 * r26476;
        double r26479 = r26475 * r26478;
        double r26480 = exp(r26460);
        double r26481 = r26460 * r26480;
        double r26482 = r26480 + r26481;
        double r26483 = r26463 / r26482;
        double r26484 = r26479 + r26483;
        double r26485 = r26469 ? r26467 : r26484;
        return r26485;
}

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 r26486, r26487, r26488, r26489, r26490, r26491, r26492, r26493;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r26486);
        mpfr_init(r26487);
        mpfr_init(r26488);
        mpfr_init(r26489);
        mpfr_init(r26490);
        mpfr_init(r26491);
        mpfr_init(r26492);
        mpfr_init(r26493);
}

double f_im(double wj, double x) {
        mpfr_set_d(r26486, wj, MPFR_RNDN);
        mpfr_exp(r26487, r26486, MPFR_RNDN);
        mpfr_mul(r26488, r26486, r26487, MPFR_RNDN);
        mpfr_set_d(r26489, x, MPFR_RNDN);
        mpfr_sub(r26490, r26488, r26489, MPFR_RNDN);
        mpfr_add(r26491, r26487, r26488, MPFR_RNDN);
        mpfr_div(r26492, r26490, r26491, MPFR_RNDN);
        mpfr_sub(r26493, r26486, r26492, MPFR_RNDN);
        return mpfr_get_d(r26493, MPFR_RNDN);
}

static mpfr_t r26494, r26495, r26496, r26497, r26498, r26499, r26500, r26501, r26502, r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511, r26512, r26513, r26514, r26515, r26516, r26517, r26518, r26519;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26494);
        mpfr_init_set_str(r26495, "2", 10, MPFR_RNDN);
        mpfr_init(r26496);
        mpfr_init(r26497);
        mpfr_init(r26498);
        mpfr_init(r26499);
        mpfr_init(r26500);
        mpfr_init(r26501);
        mpfr_init_set_str(r26502, "1.1273399280310297e-25", 10, MPFR_RNDN);
        mpfr_init(r26503);
        mpfr_init_set_str(r26504, "1", 10, MPFR_RNDN);
        mpfr_init(r26505);
        mpfr_init(r26506);
        mpfr_init(r26507);
        mpfr_init(r26508);
        mpfr_init(r26509);
        mpfr_init(r26510);
        mpfr_init(r26511);
        mpfr_init(r26512);
        mpfr_init(r26513);
        mpfr_init(r26514);
        mpfr_init(r26515);
        mpfr_init(r26516);
        mpfr_init(r26517);
        mpfr_init(r26518);
        mpfr_init(r26519);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r26494, wj, MPFR_RNDN);
        ;
        mpfr_pow(r26496, r26494, r26495, MPFR_RNDN);
        mpfr_set_d(r26497, x, MPFR_RNDN);
        mpfr_add(r26498, r26496, r26497, MPFR_RNDN);
        mpfr_mul(r26499, r26494, r26497, MPFR_RNDN);
        mpfr_mul(r26500, r26495, r26499, MPFR_RNDN);
        mpfr_sub(r26501, r26498, r26500, MPFR_RNDN);
        ;
        mpfr_set_si(r26503, mpfr_cmp(r26501, r26502) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r26505, r26504, r26494, MPFR_RNDN);
        mpfr_div(r26506, r26494, r26505, MPFR_RNDN);
        mpfr_sub(r26507, r26494, r26506, MPFR_RNDN);
        mpfr_cbrt(r26508, r26507, MPFR_RNDN);
        mpfr_mul(r26509, r26508, r26508, MPFR_RNDN);
        mpfr_cbrt(r26510, r26508, MPFR_RNDN);
        mpfr_mul(r26511, r26510, r26510, MPFR_RNDN);
        mpfr_mul(r26512, r26511, r26510, MPFR_RNDN);
        mpfr_mul(r26513, r26509, r26512, MPFR_RNDN);
        mpfr_exp(r26514, r26494, MPFR_RNDN);
        mpfr_mul(r26515, r26494, r26514, MPFR_RNDN);
        mpfr_add(r26516, r26514, r26515, MPFR_RNDN);
        mpfr_div(r26517, r26497, r26516, MPFR_RNDN);
        mpfr_add(r26518, r26513, r26517, MPFR_RNDN);
        if (mpfr_get_si(r26503, MPFR_RNDN)) { mpfr_set(r26519, r26501, MPFR_RNDN); } else { mpfr_set(r26519, r26518, MPFR_RNDN); };
        return mpfr_get_d(r26519, MPFR_RNDN);
}

static mpfr_t r26520, r26521, r26522, r26523, r26524, r26525, r26526, r26527, r26528, r26529, r26530, r26531, r26532, r26533, r26534, r26535, r26536, r26537, r26538, r26539, r26540, r26541, r26542, r26543, r26544, r26545;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26520);
        mpfr_init_set_str(r26521, "2", 10, MPFR_RNDN);
        mpfr_init(r26522);
        mpfr_init(r26523);
        mpfr_init(r26524);
        mpfr_init(r26525);
        mpfr_init(r26526);
        mpfr_init(r26527);
        mpfr_init_set_str(r26528, "1.1273399280310297e-25", 10, MPFR_RNDN);
        mpfr_init(r26529);
        mpfr_init_set_str(r26530, "1", 10, MPFR_RNDN);
        mpfr_init(r26531);
        mpfr_init(r26532);
        mpfr_init(r26533);
        mpfr_init(r26534);
        mpfr_init(r26535);
        mpfr_init(r26536);
        mpfr_init(r26537);
        mpfr_init(r26538);
        mpfr_init(r26539);
        mpfr_init(r26540);
        mpfr_init(r26541);
        mpfr_init(r26542);
        mpfr_init(r26543);
        mpfr_init(r26544);
        mpfr_init(r26545);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r26520, wj, MPFR_RNDN);
        ;
        mpfr_pow(r26522, r26520, r26521, MPFR_RNDN);
        mpfr_set_d(r26523, x, MPFR_RNDN);
        mpfr_add(r26524, r26522, r26523, MPFR_RNDN);
        mpfr_mul(r26525, r26520, r26523, MPFR_RNDN);
        mpfr_mul(r26526, r26521, r26525, MPFR_RNDN);
        mpfr_sub(r26527, r26524, r26526, MPFR_RNDN);
        ;
        mpfr_set_si(r26529, mpfr_cmp(r26527, r26528) <= 0, MPFR_RNDN);
        ;
        mpfr_add(r26531, r26530, r26520, MPFR_RNDN);
        mpfr_div(r26532, r26520, r26531, MPFR_RNDN);
        mpfr_sub(r26533, r26520, r26532, MPFR_RNDN);
        mpfr_cbrt(r26534, r26533, MPFR_RNDN);
        mpfr_mul(r26535, r26534, r26534, MPFR_RNDN);
        mpfr_cbrt(r26536, r26534, MPFR_RNDN);
        mpfr_mul(r26537, r26536, r26536, MPFR_RNDN);
        mpfr_mul(r26538, r26537, r26536, MPFR_RNDN);
        mpfr_mul(r26539, r26535, r26538, MPFR_RNDN);
        mpfr_exp(r26540, r26520, MPFR_RNDN);
        mpfr_mul(r26541, r26520, r26540, MPFR_RNDN);
        mpfr_add(r26542, r26540, r26541, MPFR_RNDN);
        mpfr_div(r26543, r26523, r26542, MPFR_RNDN);
        mpfr_add(r26544, r26539, r26543, MPFR_RNDN);
        if (mpfr_get_si(r26529, MPFR_RNDN)) { mpfr_set(r26545, r26527, MPFR_RNDN); } else { mpfr_set(r26545, r26544, MPFR_RNDN); };
        return mpfr_get_d(r26545, MPFR_RNDN);
}