#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 r26424 = wj;
        float r26425 = exp(r26424);
        float r26426 = r26424 * r26425;
        float r26427 = x;
        float r26428 = r26426 - r26427;
        float r26429 = r26425 + r26426;
        float r26430 = r26428 / r26429;
        float r26431 = r26424 - r26430;
        return r26431;
}

double f_id(double wj, double x) {
        double r26432 = wj;
        double r26433 = exp(r26432);
        double r26434 = r26432 * r26433;
        double r26435 = x;
        double r26436 = r26434 - r26435;
        double r26437 = r26433 + r26434;
        double r26438 = r26436 / r26437;
        double r26439 = r26432 - r26438;
        return r26439;
}


double f_of(float wj, float x) {
        float r26440 = wj;
        float r26441 = r26440 * r26440;
        float r26442 = 2.2027417560231758e-30;
        bool r26443 = r26441 <= r26442;
        float r26444 = x;
        float r26445 = r26444 + r26441;
        float r26446 = 1;
        float r26447 = r26440 + r26446;
        float r26448 = r26440 / r26447;
        float r26449 = r26446 + r26440;
        float r26450 = r26444 / r26449;
        float r26451 = exp(r26440);
        float r26452 = r26450 / r26451;
        float r26453 = r26448 - r26452;
        float r26454 = r26440 - r26453;
        float r26455 = r26443 ? r26445 : r26454;
        return r26455;
}

double f_od(double wj, double x) {
        double r26456 = wj;
        double r26457 = r26456 * r26456;
        double r26458 = 2.2027417560231758e-30;
        bool r26459 = r26457 <= r26458;
        double r26460 = x;
        double r26461 = r26460 + r26457;
        double r26462 = 1;
        double r26463 = r26456 + r26462;
        double r26464 = r26456 / r26463;
        double r26465 = r26462 + r26456;
        double r26466 = r26460 / r26465;
        double r26467 = exp(r26456);
        double r26468 = r26466 / r26467;
        double r26469 = r26464 - r26468;
        double r26470 = r26456 - r26469;
        double r26471 = r26459 ? r26461 : r26470;
        return r26471;
}

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 r26472, r26473, r26474, r26475, r26476, r26477, r26478, r26479;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r26472);
        mpfr_init(r26473);
        mpfr_init(r26474);
        mpfr_init(r26475);
        mpfr_init(r26476);
        mpfr_init(r26477);
        mpfr_init(r26478);
        mpfr_init(r26479);
}

double f_im(double wj, double x) {
        mpfr_set_d(r26472, wj, MPFR_RNDN);
        mpfr_exp(r26473, r26472, MPFR_RNDN);
        mpfr_mul(r26474, r26472, r26473, MPFR_RNDN);
        mpfr_set_d(r26475, x, MPFR_RNDN);
        mpfr_sub(r26476, r26474, r26475, MPFR_RNDN);
        mpfr_add(r26477, r26473, r26474, MPFR_RNDN);
        mpfr_div(r26478, r26476, r26477, MPFR_RNDN);
        mpfr_sub(r26479, r26472, r26478, MPFR_RNDN);
        return mpfr_get_d(r26479, MPFR_RNDN);
}

static mpfr_t r26480, r26481, r26482, r26483, r26484, r26485, r26486, r26487, r26488, r26489, r26490, r26491, r26492, r26493, r26494, r26495;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26480);
        mpfr_init(r26481);
        mpfr_init_set_str(r26482, "2.2027417560231758e-30", 10, MPFR_RNDN);
        mpfr_init(r26483);
        mpfr_init(r26484);
        mpfr_init(r26485);
        mpfr_init_set_str(r26486, "1", 10, MPFR_RNDN);
        mpfr_init(r26487);
        mpfr_init(r26488);
        mpfr_init(r26489);
        mpfr_init(r26490);
        mpfr_init(r26491);
        mpfr_init(r26492);
        mpfr_init(r26493);
        mpfr_init(r26494);
        mpfr_init(r26495);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r26480, wj, MPFR_RNDN);
        mpfr_mul(r26481, r26480, r26480, MPFR_RNDN);
        ;
        mpfr_set_si(r26483, mpfr_cmp(r26481, r26482) <= 0, MPFR_RNDN);
        mpfr_set_d(r26484, x, MPFR_RNDN);
        mpfr_add(r26485, r26484, r26481, MPFR_RNDN);
        ;
        mpfr_add(r26487, r26480, r26486, MPFR_RNDN);
        mpfr_div(r26488, r26480, r26487, MPFR_RNDN);
        mpfr_add(r26489, r26486, r26480, MPFR_RNDN);
        mpfr_div(r26490, r26484, r26489, MPFR_RNDN);
        mpfr_exp(r26491, r26480, MPFR_RNDN);
        mpfr_div(r26492, r26490, r26491, MPFR_RNDN);
        mpfr_sub(r26493, r26488, r26492, MPFR_RNDN);
        mpfr_sub(r26494, r26480, r26493, MPFR_RNDN);
        if (mpfr_get_si(r26483, MPFR_RNDN)) { mpfr_set(r26495, r26485, MPFR_RNDN); } else { mpfr_set(r26495, r26494, MPFR_RNDN); };
        return mpfr_get_d(r26495, MPFR_RNDN);
}

static mpfr_t r26496, r26497, r26498, r26499, r26500, r26501, r26502, r26503, r26504, r26505, r26506, r26507, r26508, r26509, r26510, r26511;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26496);
        mpfr_init(r26497);
        mpfr_init_set_str(r26498, "2.2027417560231758e-30", 10, MPFR_RNDN);
        mpfr_init(r26499);
        mpfr_init(r26500);
        mpfr_init(r26501);
        mpfr_init_set_str(r26502, "1", 10, MPFR_RNDN);
        mpfr_init(r26503);
        mpfr_init(r26504);
        mpfr_init(r26505);
        mpfr_init(r26506);
        mpfr_init(r26507);
        mpfr_init(r26508);
        mpfr_init(r26509);
        mpfr_init(r26510);
        mpfr_init(r26511);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r26496, wj, MPFR_RNDN);
        mpfr_mul(r26497, r26496, r26496, MPFR_RNDN);
        ;
        mpfr_set_si(r26499, mpfr_cmp(r26497, r26498) <= 0, MPFR_RNDN);
        mpfr_set_d(r26500, x, MPFR_RNDN);
        mpfr_add(r26501, r26500, r26497, MPFR_RNDN);
        ;
        mpfr_add(r26503, r26496, r26502, MPFR_RNDN);
        mpfr_div(r26504, r26496, r26503, MPFR_RNDN);
        mpfr_add(r26505, r26502, r26496, MPFR_RNDN);
        mpfr_div(r26506, r26500, r26505, MPFR_RNDN);
        mpfr_exp(r26507, r26496, MPFR_RNDN);
        mpfr_div(r26508, r26506, r26507, MPFR_RNDN);
        mpfr_sub(r26509, r26504, r26508, MPFR_RNDN);
        mpfr_sub(r26510, r26496, r26509, MPFR_RNDN);
        if (mpfr_get_si(r26499, MPFR_RNDN)) { mpfr_set(r26511, r26501, MPFR_RNDN); } else { mpfr_set(r26511, r26510, MPFR_RNDN); };
        return mpfr_get_d(r26511, MPFR_RNDN);
}