#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 r19971 = wj;
        float r19972 = exp(r19971);
        float r19973 = r19971 * r19972;
        float r19974 = x;
        float r19975 = r19973 - r19974;
        float r19976 = r19972 + r19973;
        float r19977 = r19975 / r19976;
        float r19978 = r19971 - r19977;
        return r19978;
}

double f_id(double wj, double x) {
        double r19979 = wj;
        double r19980 = exp(r19979);
        double r19981 = r19979 * r19980;
        double r19982 = x;
        double r19983 = r19981 - r19982;
        double r19984 = r19980 + r19981;
        double r19985 = r19983 / r19984;
        double r19986 = r19979 - r19985;
        return r19986;
}


double f_of(float wj, float x) {
        float r19987 = wj;
        float r19988 = exp(r19987);
        float r19989 = r19987 * r19988;
        float r19990 = x;
        float r19991 = r19989 - r19990;
        float r19992 = r19988 + r19989;
        float r19993 = r19991 / r19992;
        float r19994 = r19987 - r19993;
        float r19995 = 1.247241741003307e-15f;
        bool r19996 = r19994 <= r19995;
        float r19997 = r19987 * r19987;
        float r19998 = r19997 + r19990;
        float r19999 = 2.0f;
        float r20000 = r19987 * r19990;
        float r20001 = r19999 * r20000;
        float r20002 = r19998 - r20001;
        float r20003 = 1.0f;
        float r20004 = r19987 + r20003;
        float r20005 = r19987 / r20004;
        float r20006 = r20003 + r19987;
        float r20007 = r19990 / r20006;
        float r20008 = r20007 / r19988;
        float r20009 = r20005 - r20008;
        float r20010 = r19987 - r20009;
        float r20011 = r19996 ? r20002 : r20010;
        return r20011;
}

double f_od(double wj, double x) {
        double r20012 = wj;
        double r20013 = exp(r20012);
        double r20014 = r20012 * r20013;
        double r20015 = x;
        double r20016 = r20014 - r20015;
        double r20017 = r20013 + r20014;
        double r20018 = r20016 / r20017;
        double r20019 = r20012 - r20018;
        double r20020 = 1.247241741003307e-15;
        bool r20021 = r20019 <= r20020;
        double r20022 = r20012 * r20012;
        double r20023 = r20022 + r20015;
        double r20024 = 2.0;
        double r20025 = r20012 * r20015;
        double r20026 = r20024 * r20025;
        double r20027 = r20023 - r20026;
        double r20028 = 1.0;
        double r20029 = r20012 + r20028;
        double r20030 = r20012 / r20029;
        double r20031 = r20028 + r20012;
        double r20032 = r20015 / r20031;
        double r20033 = r20032 / r20013;
        double r20034 = r20030 - r20033;
        double r20035 = r20012 - r20034;
        double r20036 = r20021 ? r20027 : r20035;
        return r20036;
}

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 r20037, r20038, r20039, r20040, r20041, r20042, r20043, r20044;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20037);
        mpfr_init(r20038);
        mpfr_init(r20039);
        mpfr_init(r20040);
        mpfr_init(r20041);
        mpfr_init(r20042);
        mpfr_init(r20043);
        mpfr_init(r20044);
}

double f_im(double wj, double x) {
        mpfr_set_d(r20037, wj, MPFR_RNDN);
        mpfr_exp(r20038, r20037, MPFR_RNDN);
        mpfr_mul(r20039, r20037, r20038, MPFR_RNDN);
        mpfr_set_d(r20040, x, MPFR_RNDN);
        mpfr_sub(r20041, r20039, r20040, MPFR_RNDN);
        mpfr_add(r20042, r20038, r20039, MPFR_RNDN);
        mpfr_div(r20043, r20041, r20042, MPFR_RNDN);
        mpfr_sub(r20044, r20037, r20043, MPFR_RNDN);
        return mpfr_get_d(r20044, MPFR_RNDN);
}

static mpfr_t r20045, r20046, r20047, r20048, r20049, r20050, r20051, r20052, r20053, r20054, r20055, r20056, r20057, r20058, r20059, r20060, r20061, r20062, r20063, r20064, r20065, r20066, r20067, r20068, r20069;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20045);
        mpfr_init(r20046);
        mpfr_init(r20047);
        mpfr_init(r20048);
        mpfr_init(r20049);
        mpfr_init(r20050);
        mpfr_init(r20051);
        mpfr_init(r20052);
        mpfr_init_set_str(r20053, "1.247241741003307e-15", 10, MPFR_RNDN);
        mpfr_init(r20054);
        mpfr_init(r20055);
        mpfr_init(r20056);
        mpfr_init_set_str(r20057, "2", 10, MPFR_RNDN);
        mpfr_init(r20058);
        mpfr_init(r20059);
        mpfr_init(r20060);
        mpfr_init_set_str(r20061, "1", 10, MPFR_RNDN);
        mpfr_init(r20062);
        mpfr_init(r20063);
        mpfr_init(r20064);
        mpfr_init(r20065);
        mpfr_init(r20066);
        mpfr_init(r20067);
        mpfr_init(r20068);
        mpfr_init(r20069);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r20045, wj, MPFR_RNDN);
        mpfr_exp(r20046, r20045, MPFR_RNDN);
        mpfr_mul(r20047, r20045, r20046, MPFR_RNDN);
        mpfr_set_d(r20048, x, MPFR_RNDN);
        mpfr_sub(r20049, r20047, r20048, MPFR_RNDN);
        mpfr_add(r20050, r20046, r20047, MPFR_RNDN);
        mpfr_div(r20051, r20049, r20050, MPFR_RNDN);
        mpfr_sub(r20052, r20045, r20051, MPFR_RNDN);
        ;
        mpfr_set_si(r20054, mpfr_cmp(r20052, r20053) <= 0, MPFR_RNDN);
        mpfr_sqr(r20055, r20045, MPFR_RNDN);
        mpfr_add(r20056, r20055, r20048, MPFR_RNDN);
        ;
        mpfr_mul(r20058, r20045, r20048, MPFR_RNDN);
        mpfr_mul(r20059, r20057, r20058, MPFR_RNDN);
        mpfr_sub(r20060, r20056, r20059, MPFR_RNDN);
        ;
        mpfr_add(r20062, r20045, r20061, MPFR_RNDN);
        mpfr_div(r20063, r20045, r20062, MPFR_RNDN);
        mpfr_add(r20064, r20061, r20045, MPFR_RNDN);
        mpfr_div(r20065, r20048, r20064, MPFR_RNDN);
        mpfr_div(r20066, r20065, r20046, MPFR_RNDN);
        mpfr_sub(r20067, r20063, r20066, MPFR_RNDN);
        mpfr_sub(r20068, r20045, r20067, MPFR_RNDN);
        if (mpfr_get_si(r20054, MPFR_RNDN)) { mpfr_set(r20069, r20060, MPFR_RNDN); } else { mpfr_set(r20069, r20068, MPFR_RNDN); };
        return mpfr_get_d(r20069, MPFR_RNDN);
}

static mpfr_t r20070, r20071, r20072, r20073, r20074, r20075, r20076, r20077, r20078, r20079, r20080, r20081, r20082, r20083, r20084, r20085, r20086, r20087, r20088, r20089, r20090, r20091, r20092, r20093, r20094;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20070);
        mpfr_init(r20071);
        mpfr_init(r20072);
        mpfr_init(r20073);
        mpfr_init(r20074);
        mpfr_init(r20075);
        mpfr_init(r20076);
        mpfr_init(r20077);
        mpfr_init_set_str(r20078, "1.247241741003307e-15", 10, MPFR_RNDN);
        mpfr_init(r20079);
        mpfr_init(r20080);
        mpfr_init(r20081);
        mpfr_init_set_str(r20082, "2", 10, MPFR_RNDN);
        mpfr_init(r20083);
        mpfr_init(r20084);
        mpfr_init(r20085);
        mpfr_init_set_str(r20086, "1", 10, MPFR_RNDN);
        mpfr_init(r20087);
        mpfr_init(r20088);
        mpfr_init(r20089);
        mpfr_init(r20090);
        mpfr_init(r20091);
        mpfr_init(r20092);
        mpfr_init(r20093);
        mpfr_init(r20094);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r20070, wj, MPFR_RNDN);
        mpfr_exp(r20071, r20070, MPFR_RNDN);
        mpfr_mul(r20072, r20070, r20071, MPFR_RNDN);
        mpfr_set_d(r20073, x, MPFR_RNDN);
        mpfr_sub(r20074, r20072, r20073, MPFR_RNDN);
        mpfr_add(r20075, r20071, r20072, MPFR_RNDN);
        mpfr_div(r20076, r20074, r20075, MPFR_RNDN);
        mpfr_sub(r20077, r20070, r20076, MPFR_RNDN);
        ;
        mpfr_set_si(r20079, mpfr_cmp(r20077, r20078) <= 0, MPFR_RNDN);
        mpfr_sqr(r20080, r20070, MPFR_RNDN);
        mpfr_add(r20081, r20080, r20073, MPFR_RNDN);
        ;
        mpfr_mul(r20083, r20070, r20073, MPFR_RNDN);
        mpfr_mul(r20084, r20082, r20083, MPFR_RNDN);
        mpfr_sub(r20085, r20081, r20084, MPFR_RNDN);
        ;
        mpfr_add(r20087, r20070, r20086, MPFR_RNDN);
        mpfr_div(r20088, r20070, r20087, MPFR_RNDN);
        mpfr_add(r20089, r20086, r20070, MPFR_RNDN);
        mpfr_div(r20090, r20073, r20089, MPFR_RNDN);
        mpfr_div(r20091, r20090, r20071, MPFR_RNDN);
        mpfr_sub(r20092, r20088, r20091, MPFR_RNDN);
        mpfr_sub(r20093, r20070, r20092, MPFR_RNDN);
        if (mpfr_get_si(r20079, MPFR_RNDN)) { mpfr_set(r20094, r20085, MPFR_RNDN); } else { mpfr_set(r20094, r20093, MPFR_RNDN); };
        return mpfr_get_d(r20094, MPFR_RNDN);
}