#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19464 = wj; float r19465 = exp(r19464); float r19466 = r19464 * r19465; float r19467 = x; float r19468 = r19466 - r19467; float r19469 = r19465 + r19466; float r19470 = r19468 / r19469; float r19471 = r19464 - r19470; return r19471; } double f_id(double wj, double x) { double r19472 = wj; double r19473 = exp(r19472); double r19474 = r19472 * r19473; double r19475 = x; double r19476 = r19474 - r19475; double r19477 = r19473 + r19474; double r19478 = r19476 / r19477; double r19479 = r19472 - r19478; return r19479; } double f_of(float wj, float x) { float r19480 = wj; float r19481 = 24.52827262878418f; bool r19482 = r19480 <= r19481; float r19483 = exp(r19480); float r19484 = r19480 * r19483; float r19485 = x; float r19486 = r19484 - r19485; float r19487 = r19483 + r19484; float r19488 = r19486 / r19487; float r19489 = r19480 - r19488; float r19490 = 1.0f; float r19491 = r19490 + r19480; float r19492 = r19480 / r19491; float r19493 = r19480 - r19492; float r19494 = r19485 / r19487; float r19495 = r19493 + r19494; float r19496 = r19482 ? r19489 : r19495; return r19496; } double f_od(double wj, double x) { double r19497 = wj; double r19498 = 24.52827262878418; bool r19499 = r19497 <= r19498; double r19500 = exp(r19497); double r19501 = r19497 * r19500; double r19502 = x; double r19503 = r19501 - r19502; double r19504 = r19500 + r19501; double r19505 = r19503 / r19504; double r19506 = r19497 - r19505; double r19507 = 1.0; double r19508 = r19507 + r19497; double r19509 = r19497 / r19508; double r19510 = r19497 - r19509; double r19511 = r19502 / r19504; double r19512 = r19510 + r19511; double r19513 = r19499 ? r19506 : r19512; return r19513; } 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 r19514, r19515, r19516, r19517, r19518, r19519, r19520, r19521; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19514); mpfr_init(r19515); mpfr_init(r19516); mpfr_init(r19517); mpfr_init(r19518); mpfr_init(r19519); mpfr_init(r19520); mpfr_init(r19521); } double f_im(double wj, double x) { mpfr_set_d(r19514, wj, MPFR_RNDN); mpfr_exp(r19515, r19514, MPFR_RNDN); mpfr_mul(r19516, r19514, r19515, MPFR_RNDN); mpfr_set_d(r19517, x, MPFR_RNDN); mpfr_sub(r19518, r19516, r19517, MPFR_RNDN); mpfr_add(r19519, r19515, r19516, MPFR_RNDN); mpfr_div(r19520, r19518, r19519, MPFR_RNDN); mpfr_sub(r19521, r19514, r19520, MPFR_RNDN); return mpfr_get_d(r19521, MPFR_RNDN); } static mpfr_t r19522, r19523, r19524, r19525, r19526, r19527, r19528, r19529, r19530, r19531, r19532, r19533, r19534, r19535, r19536, r19537, r19538; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19522); mpfr_init_set_str(r19523, "24.528273f0", 10, MPFR_RNDN); mpfr_init(r19524); mpfr_init(r19525); mpfr_init(r19526); mpfr_init(r19527); mpfr_init(r19528); mpfr_init(r19529); mpfr_init(r19530); mpfr_init(r19531); mpfr_init_set_str(r19532, "1", 10, MPFR_RNDN); mpfr_init(r19533); mpfr_init(r19534); mpfr_init(r19535); mpfr_init(r19536); mpfr_init(r19537); mpfr_init(r19538); } double f_fm(double wj, double x) { mpfr_set_d(r19522, wj, MPFR_RNDN); ; mpfr_set_si(r19524, mpfr_cmp(r19522, r19523) <= 0, MPFR_RNDN); mpfr_exp(r19525, r19522, MPFR_RNDN); mpfr_mul(r19526, r19522, r19525, MPFR_RNDN); mpfr_set_d(r19527, x, MPFR_RNDN); mpfr_sub(r19528, r19526, r19527, MPFR_RNDN); mpfr_add(r19529, r19525, r19526, MPFR_RNDN); mpfr_div(r19530, r19528, r19529, MPFR_RNDN); mpfr_sub(r19531, r19522, r19530, MPFR_RNDN); ; mpfr_add(r19533, r19532, r19522, MPFR_RNDN); mpfr_div(r19534, r19522, r19533, MPFR_RNDN); mpfr_sub(r19535, r19522, r19534, MPFR_RNDN); mpfr_div(r19536, r19527, r19529, MPFR_RNDN); mpfr_add(r19537, r19535, r19536, MPFR_RNDN); if (mpfr_get_si(r19524, MPFR_RNDN)) { mpfr_set(r19538, r19531, MPFR_RNDN); } else { mpfr_set(r19538, r19537, MPFR_RNDN); }; return mpfr_get_d(r19538, MPFR_RNDN); } static mpfr_t r19539, r19540, r19541, r19542, r19543, r19544, r19545, r19546, r19547, r19548, r19549, r19550, r19551, r19552, r19553, r19554, r19555; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19539); mpfr_init_set_str(r19540, "24.528273f0", 10, MPFR_RNDN); mpfr_init(r19541); mpfr_init(r19542); mpfr_init(r19543); mpfr_init(r19544); mpfr_init(r19545); mpfr_init(r19546); mpfr_init(r19547); mpfr_init(r19548); mpfr_init_set_str(r19549, "1", 10, MPFR_RNDN); mpfr_init(r19550); mpfr_init(r19551); mpfr_init(r19552); mpfr_init(r19553); mpfr_init(r19554); mpfr_init(r19555); } double f_dm(double wj, double x) { mpfr_set_d(r19539, wj, MPFR_RNDN); ; mpfr_set_si(r19541, mpfr_cmp(r19539, r19540) <= 0, MPFR_RNDN); mpfr_exp(r19542, r19539, MPFR_RNDN); mpfr_mul(r19543, r19539, r19542, MPFR_RNDN); mpfr_set_d(r19544, x, MPFR_RNDN); mpfr_sub(r19545, r19543, r19544, MPFR_RNDN); mpfr_add(r19546, r19542, r19543, MPFR_RNDN); mpfr_div(r19547, r19545, r19546, MPFR_RNDN); mpfr_sub(r19548, r19539, r19547, MPFR_RNDN); ; mpfr_add(r19550, r19549, r19539, MPFR_RNDN); mpfr_div(r19551, r19539, r19550, MPFR_RNDN); mpfr_sub(r19552, r19539, r19551, MPFR_RNDN); mpfr_div(r19553, r19544, r19546, MPFR_RNDN); mpfr_add(r19554, r19552, r19553, MPFR_RNDN); if (mpfr_get_si(r19541, MPFR_RNDN)) { mpfr_set(r19555, r19548, MPFR_RNDN); } else { mpfr_set(r19555, r19554, MPFR_RNDN); }; return mpfr_get_d(r19555, MPFR_RNDN); }