#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 r26388 = wj;
        float r26389 = exp(r26388);
        float r26390 = r26388 * r26389;
        float r26391 = x;
        float r26392 = r26390 - r26391;
        float r26393 = r26389 + r26390;
        float r26394 = r26392 / r26393;
        float r26395 = r26388 - r26394;
        return r26395;
}

double f_id(double wj, double x) {
        double r26396 = wj;
        double r26397 = exp(r26396);
        double r26398 = r26396 * r26397;
        double r26399 = x;
        double r26400 = r26398 - r26399;
        double r26401 = r26397 + r26398;
        double r26402 = r26400 / r26401;
        double r26403 = r26396 - r26402;
        return r26403;
}


double f_of(float wj, float x) {
        float r26404 = wj;
        float r26405 = x;
        float r26406 = exp(r26404);
        float r26407 = r26405 / r26406;
        float r26408 = r26404 - r26407;
        float r26409 = r26404 * r26404;
        float r26410 = r26409 * r26406;
        float r26411 = r26406 - r26410;
        float r26412 = r26408 / r26411;
        float r26413 = 1;
        float r26414 = -r26413;
        float r26415 = 1/3;
        float r26416 = 1/2;
        float r26417 = fma(r26415, r26404, r26416);
        float r26418 = r26417 * r26409;
        float r26419 = r26414 + r26418;
        float r26420 = fma(r26412, r26419, r26404);
        float r26421 = 2.608735788713234e-17;
        bool r26422 = r26420 <= r26421;
        float r26423 = fma(r26404, r26404, r26405);
        float r26424 = r26404 + r26413;
        float r26425 = r26404 / r26424;
        float r26426 = fma(r26406, r26404, r26406);
        float r26427 = r26405 / r26426;
        float r26428 = r26425 - r26427;
        float r26429 = r26404 - r26428;
        float r26430 = r26422 ? r26423 : r26429;
        return r26430;
}

double f_od(double wj, double x) {
        double r26431 = wj;
        double r26432 = x;
        double r26433 = exp(r26431);
        double r26434 = r26432 / r26433;
        double r26435 = r26431 - r26434;
        double r26436 = r26431 * r26431;
        double r26437 = r26436 * r26433;
        double r26438 = r26433 - r26437;
        double r26439 = r26435 / r26438;
        double r26440 = 1;
        double r26441 = -r26440;
        double r26442 = 1/3;
        double r26443 = 1/2;
        double r26444 = fma(r26442, r26431, r26443);
        double r26445 = r26444 * r26436;
        double r26446 = r26441 + r26445;
        double r26447 = fma(r26439, r26446, r26431);
        double r26448 = 2.608735788713234e-17;
        bool r26449 = r26447 <= r26448;
        double r26450 = fma(r26431, r26431, r26432);
        double r26451 = r26431 + r26440;
        double r26452 = r26431 / r26451;
        double r26453 = fma(r26433, r26431, r26433);
        double r26454 = r26432 / r26453;
        double r26455 = r26452 - r26454;
        double r26456 = r26431 - r26455;
        double r26457 = r26449 ? r26450 : r26456;
        return r26457;
}

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 r26458, r26459, r26460, r26461, r26462, r26463, r26464, r26465;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(848);
        mpfr_init(r26458);
        mpfr_init(r26459);
        mpfr_init(r26460);
        mpfr_init(r26461);
        mpfr_init(r26462);
        mpfr_init(r26463);
        mpfr_init(r26464);
        mpfr_init(r26465);
}

double f_im(double wj, double x) {
        mpfr_set_d(r26458, wj, MPFR_RNDN);
        mpfr_exp(r26459, r26458, MPFR_RNDN);
        mpfr_mul(r26460, r26458, r26459, MPFR_RNDN);
        mpfr_set_d(r26461, x, MPFR_RNDN);
        mpfr_sub(r26462, r26460, r26461, MPFR_RNDN);
        mpfr_add(r26463, r26459, r26460, MPFR_RNDN);
        mpfr_div(r26464, r26462, r26463, MPFR_RNDN);
        mpfr_sub(r26465, r26458, r26464, MPFR_RNDN);
        return mpfr_get_d(r26465, MPFR_RNDN);
}

static mpfr_t r26466, r26467, r26468, r26469, r26470, r26471, r26472, r26473, r26474, r26475, r26476, r26477, r26478, r26479, r26480, r26481, r26482, r26483, r26484, r26485, r26486, r26487, r26488, r26489, r26490, r26491, r26492;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26466);
        mpfr_init(r26467);
        mpfr_init(r26468);
        mpfr_init(r26469);
        mpfr_init(r26470);
        mpfr_init(r26471);
        mpfr_init(r26472);
        mpfr_init(r26473);
        mpfr_init(r26474);
        mpfr_init_set_str(r26475, "1", 10, MPFR_RNDN);
        mpfr_init(r26476);
        mpfr_init_set_str(r26477, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26478, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26479);
        mpfr_init(r26480);
        mpfr_init(r26481);
        mpfr_init(r26482);
        mpfr_init_set_str(r26483, "2.608735788713234e-17", 10, MPFR_RNDN);
        mpfr_init(r26484);
        mpfr_init(r26485);
        mpfr_init(r26486);
        mpfr_init(r26487);
        mpfr_init(r26488);
        mpfr_init(r26489);
        mpfr_init(r26490);
        mpfr_init(r26491);
        mpfr_init(r26492);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r26466, wj, MPFR_RNDN);
        mpfr_set_d(r26467, x, MPFR_RNDN);
        mpfr_exp(r26468, r26466, MPFR_RNDN);
        mpfr_div(r26469, r26467, r26468, MPFR_RNDN);
        mpfr_sub(r26470, r26466, r26469, MPFR_RNDN);
        mpfr_mul(r26471, r26466, r26466, MPFR_RNDN);
        mpfr_mul(r26472, r26471, r26468, MPFR_RNDN);
        mpfr_sub(r26473, r26468, r26472, MPFR_RNDN);
        mpfr_div(r26474, r26470, r26473, MPFR_RNDN);
        ;
        mpfr_neg(r26476, r26475, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r26479, r26477, r26466, r26478, MPFR_RNDN);
        mpfr_mul(r26480, r26479, r26471, MPFR_RNDN);
        mpfr_add(r26481, r26476, r26480, MPFR_RNDN);
        mpfr_fma(r26482, r26474, r26481, r26466, MPFR_RNDN);
        ;
        mpfr_set_si(r26484, mpfr_cmp(r26482, r26483) <= 0, MPFR_RNDN);
        mpfr_fma(r26485, r26466, r26466, r26467, MPFR_RNDN);
        mpfr_add(r26486, r26466, r26475, MPFR_RNDN);
        mpfr_div(r26487, r26466, r26486, MPFR_RNDN);
        mpfr_fma(r26488, r26468, r26466, r26468, MPFR_RNDN);
        mpfr_div(r26489, r26467, r26488, MPFR_RNDN);
        mpfr_sub(r26490, r26487, r26489, MPFR_RNDN);
        mpfr_sub(r26491, r26466, r26490, MPFR_RNDN);
        if (mpfr_get_si(r26484, MPFR_RNDN)) { mpfr_set(r26492, r26485, MPFR_RNDN); } else { mpfr_set(r26492, r26491, MPFR_RNDN); };
        return mpfr_get_d(r26492, MPFR_RNDN);
}

static mpfr_t r26493, 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_dm() {
        mpfr_set_default_prec(848);
        mpfr_init(r26493);
        mpfr_init(r26494);
        mpfr_init(r26495);
        mpfr_init(r26496);
        mpfr_init(r26497);
        mpfr_init(r26498);
        mpfr_init(r26499);
        mpfr_init(r26500);
        mpfr_init(r26501);
        mpfr_init_set_str(r26502, "1", 10, MPFR_RNDN);
        mpfr_init(r26503);
        mpfr_init_set_str(r26504, "1/3", 10, MPFR_RNDN);
        mpfr_init_set_str(r26505, "1/2", 10, MPFR_RNDN);
        mpfr_init(r26506);
        mpfr_init(r26507);
        mpfr_init(r26508);
        mpfr_init(r26509);
        mpfr_init_set_str(r26510, "2.608735788713234e-17", 10, MPFR_RNDN);
        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_dm(double wj, double x) {
        mpfr_set_d(r26493, wj, MPFR_RNDN);
        mpfr_set_d(r26494, x, MPFR_RNDN);
        mpfr_exp(r26495, r26493, MPFR_RNDN);
        mpfr_div(r26496, r26494, r26495, MPFR_RNDN);
        mpfr_sub(r26497, r26493, r26496, MPFR_RNDN);
        mpfr_mul(r26498, r26493, r26493, MPFR_RNDN);
        mpfr_mul(r26499, r26498, r26495, MPFR_RNDN);
        mpfr_sub(r26500, r26495, r26499, MPFR_RNDN);
        mpfr_div(r26501, r26497, r26500, MPFR_RNDN);
        ;
        mpfr_neg(r26503, r26502, MPFR_RNDN);
        ;
        ;
        mpfr_fma(r26506, r26504, r26493, r26505, MPFR_RNDN);
        mpfr_mul(r26507, r26506, r26498, MPFR_RNDN);
        mpfr_add(r26508, r26503, r26507, MPFR_RNDN);
        mpfr_fma(r26509, r26501, r26508, r26493, MPFR_RNDN);
        ;
        mpfr_set_si(r26511, mpfr_cmp(r26509, r26510) <= 0, MPFR_RNDN);
        mpfr_fma(r26512, r26493, r26493, r26494, MPFR_RNDN);
        mpfr_add(r26513, r26493, r26502, MPFR_RNDN);
        mpfr_div(r26514, r26493, r26513, MPFR_RNDN);
        mpfr_fma(r26515, r26495, r26493, r26495, MPFR_RNDN);
        mpfr_div(r26516, r26494, r26515, MPFR_RNDN);
        mpfr_sub(r26517, r26514, r26516, MPFR_RNDN);
        mpfr_sub(r26518, r26493, r26517, MPFR_RNDN);
        if (mpfr_get_si(r26511, MPFR_RNDN)) { mpfr_set(r26519, r26512, MPFR_RNDN); } else { mpfr_set(r26519, r26518, MPFR_RNDN); };
        return mpfr_get_d(r26519, MPFR_RNDN);
}