#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 r20135 = wj;
        float r20136 = exp(r20135);
        float r20137 = r20135 * r20136;
        float r20138 = x;
        float r20139 = r20137 - r20138;
        float r20140 = r20136 + r20137;
        float r20141 = r20139 / r20140;
        float r20142 = r20135 - r20141;
        return r20142;
}

double f_id(double wj, double x) {
        double r20143 = wj;
        double r20144 = exp(r20143);
        double r20145 = r20143 * r20144;
        double r20146 = x;
        double r20147 = r20145 - r20146;
        double r20148 = r20144 + r20145;
        double r20149 = r20147 / r20148;
        double r20150 = r20143 - r20149;
        return r20150;
}


double f_of(float wj, float x) {
        float r20151 = wj;
        float r20152 = exp(r20151);
        float r20153 = r20151 * r20152;
        float r20154 = x;
        float r20155 = r20153 - r20154;
        float r20156 = r20152 + r20153;
        float r20157 = r20155 / r20156;
        float r20158 = r20151 - r20157;
        float r20159 = 3.8050401405631674e-17f;
        bool r20160 = r20158 <= r20159;
        float r20161 = r20151 * r20151;
        float r20162 = r20151 * (r20151 * r20151);
        float r20163 = r20161 - r20162;
        float r20164 = 1.0f;
        float r20165 = r20151 + r20164;
        float r20166 = r20165 * r20152;
        float r20167 = r20154 / r20166;
        float r20168 = r20163 + r20167;
        float r20169 = r20151 / r20165;
        float r20170 = r20164 + r20151;
        float r20171 = r20154 / r20170;
        float r20172 = r20171 / r20152;
        float r20173 = r20169 - r20172;
        float r20174 = r20151 - r20173;
        float r20175 = r20160 ? r20168 : r20174;
        return r20175;
}

double f_od(double wj, double x) {
        double r20176 = wj;
        double r20177 = exp(r20176);
        double r20178 = r20176 * r20177;
        double r20179 = x;
        double r20180 = r20178 - r20179;
        double r20181 = r20177 + r20178;
        double r20182 = r20180 / r20181;
        double r20183 = r20176 - r20182;
        double r20184 = 3.8050401405631674e-17;
        bool r20185 = r20183 <= r20184;
        double r20186 = r20176 * r20176;
        double r20187 = r20176 * (r20176 * r20176);
        double r20188 = r20186 - r20187;
        double r20189 = 1.0;
        double r20190 = r20176 + r20189;
        double r20191 = r20190 * r20177;
        double r20192 = r20179 / r20191;
        double r20193 = r20188 + r20192;
        double r20194 = r20176 / r20190;
        double r20195 = r20189 + r20176;
        double r20196 = r20179 / r20195;
        double r20197 = r20196 / r20177;
        double r20198 = r20194 - r20197;
        double r20199 = r20176 - r20198;
        double r20200 = r20185 ? r20193 : r20199;
        return r20200;
}

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 r20201, r20202, r20203, r20204, r20205, r20206, r20207, r20208;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(144);
        mpfr_init(r20201);
        mpfr_init(r20202);
        mpfr_init(r20203);
        mpfr_init(r20204);
        mpfr_init(r20205);
        mpfr_init(r20206);
        mpfr_init(r20207);
        mpfr_init(r20208);
}

double f_im(double wj, double x) {
        mpfr_set_d(r20201, wj, MPFR_RNDN);
        mpfr_exp(r20202, r20201, MPFR_RNDN);
        mpfr_mul(r20203, r20201, r20202, MPFR_RNDN);
        mpfr_set_d(r20204, x, MPFR_RNDN);
        mpfr_sub(r20205, r20203, r20204, MPFR_RNDN);
        mpfr_add(r20206, r20202, r20203, MPFR_RNDN);
        mpfr_div(r20207, r20205, r20206, MPFR_RNDN);
        mpfr_sub(r20208, r20201, r20207, MPFR_RNDN);
        return mpfr_get_d(r20208, MPFR_RNDN);
}

static mpfr_t r20209, r20210, r20211, r20212, r20213, r20214, r20215, r20216, r20217, r20218, r20219, r20220, r20221, r20222, r20223, r20224, r20225, r20226, r20227, r20228, r20229, r20230, r20231, r20232, r20233;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20209);
        mpfr_init(r20210);
        mpfr_init(r20211);
        mpfr_init(r20212);
        mpfr_init(r20213);
        mpfr_init(r20214);
        mpfr_init(r20215);
        mpfr_init(r20216);
        mpfr_init_set_str(r20217, "3.8050401405631674e-17", 10, MPFR_RNDN);
        mpfr_init(r20218);
        mpfr_init(r20219);
        mpfr_init(r20220);
        mpfr_init(r20221);
        mpfr_init_set_str(r20222, "1", 10, MPFR_RNDN);
        mpfr_init(r20223);
        mpfr_init(r20224);
        mpfr_init(r20225);
        mpfr_init(r20226);
        mpfr_init(r20227);
        mpfr_init(r20228);
        mpfr_init(r20229);
        mpfr_init(r20230);
        mpfr_init(r20231);
        mpfr_init(r20232);
        mpfr_init(r20233);
}

double f_fm(double wj, double x) {
        mpfr_set_d(r20209, wj, MPFR_RNDN);
        mpfr_exp(r20210, r20209, MPFR_RNDN);
        mpfr_mul(r20211, r20209, r20210, MPFR_RNDN);
        mpfr_set_d(r20212, x, MPFR_RNDN);
        mpfr_sub(r20213, r20211, r20212, MPFR_RNDN);
        mpfr_add(r20214, r20210, r20211, MPFR_RNDN);
        mpfr_div(r20215, r20213, r20214, MPFR_RNDN);
        mpfr_sub(r20216, r20209, r20215, MPFR_RNDN);
        ;
        mpfr_set_si(r20218, mpfr_cmp(r20216, r20217) <= 0, MPFR_RNDN);
        mpfr_sqr(r20219, r20209, MPFR_RNDN);
        mpfr_mul(r20220, r20209, r20209, MPFR_RNDN); mpfr_mul(r20220, r20220, r20209, MPFR_RNDN);
        mpfr_sub(r20221, r20219, r20220, MPFR_RNDN);
        ;
        mpfr_add(r20223, r20209, r20222, MPFR_RNDN);
        mpfr_mul(r20224, r20223, r20210, MPFR_RNDN);
        mpfr_div(r20225, r20212, r20224, MPFR_RNDN);
        mpfr_add(r20226, r20221, r20225, MPFR_RNDN);
        mpfr_div(r20227, r20209, r20223, MPFR_RNDN);
        mpfr_add(r20228, r20222, r20209, MPFR_RNDN);
        mpfr_div(r20229, r20212, r20228, MPFR_RNDN);
        mpfr_div(r20230, r20229, r20210, MPFR_RNDN);
        mpfr_sub(r20231, r20227, r20230, MPFR_RNDN);
        mpfr_sub(r20232, r20209, r20231, MPFR_RNDN);
        if (mpfr_get_si(r20218, MPFR_RNDN)) { mpfr_set(r20233, r20226, MPFR_RNDN); } else { mpfr_set(r20233, r20232, MPFR_RNDN); };
        return mpfr_get_d(r20233, MPFR_RNDN);
}

static mpfr_t r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242, r20243, r20244, r20245, r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(144);
        mpfr_init(r20234);
        mpfr_init(r20235);
        mpfr_init(r20236);
        mpfr_init(r20237);
        mpfr_init(r20238);
        mpfr_init(r20239);
        mpfr_init(r20240);
        mpfr_init(r20241);
        mpfr_init_set_str(r20242, "3.8050401405631674e-17", 10, MPFR_RNDN);
        mpfr_init(r20243);
        mpfr_init(r20244);
        mpfr_init(r20245);
        mpfr_init(r20246);
        mpfr_init_set_str(r20247, "1", 10, MPFR_RNDN);
        mpfr_init(r20248);
        mpfr_init(r20249);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init(r20252);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
        mpfr_init(r20256);
        mpfr_init(r20257);
        mpfr_init(r20258);
}

double f_dm(double wj, double x) {
        mpfr_set_d(r20234, wj, MPFR_RNDN);
        mpfr_exp(r20235, r20234, MPFR_RNDN);
        mpfr_mul(r20236, r20234, r20235, MPFR_RNDN);
        mpfr_set_d(r20237, x, MPFR_RNDN);
        mpfr_sub(r20238, r20236, r20237, MPFR_RNDN);
        mpfr_add(r20239, r20235, r20236, MPFR_RNDN);
        mpfr_div(r20240, r20238, r20239, MPFR_RNDN);
        mpfr_sub(r20241, r20234, r20240, MPFR_RNDN);
        ;
        mpfr_set_si(r20243, mpfr_cmp(r20241, r20242) <= 0, MPFR_RNDN);
        mpfr_sqr(r20244, r20234, MPFR_RNDN);
        mpfr_mul(r20245, r20234, r20234, MPFR_RNDN); mpfr_mul(r20245, r20245, r20234, MPFR_RNDN);
        mpfr_sub(r20246, r20244, r20245, MPFR_RNDN);
        ;
        mpfr_add(r20248, r20234, r20247, MPFR_RNDN);
        mpfr_mul(r20249, r20248, r20235, MPFR_RNDN);
        mpfr_div(r20250, r20237, r20249, MPFR_RNDN);
        mpfr_add(r20251, r20246, r20250, MPFR_RNDN);
        mpfr_div(r20252, r20234, r20248, MPFR_RNDN);
        mpfr_add(r20253, r20247, r20234, MPFR_RNDN);
        mpfr_div(r20254, r20237, r20253, MPFR_RNDN);
        mpfr_div(r20255, r20254, r20235, MPFR_RNDN);
        mpfr_sub(r20256, r20252, r20255, MPFR_RNDN);
        mpfr_sub(r20257, r20234, r20256, MPFR_RNDN);
        if (mpfr_get_si(r20243, MPFR_RNDN)) { mpfr_set(r20258, r20251, MPFR_RNDN); } else { mpfr_set(r20258, r20257, MPFR_RNDN); };
        return mpfr_get_d(r20258, MPFR_RNDN);
}