#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r20171 = wj; float r20172 = exp(r20171); float r20173 = r20171 * r20172; float r20174 = x; float r20175 = r20173 - r20174; float r20176 = r20172 + r20173; float r20177 = r20175 / r20176; float r20178 = r20171 - r20177; return r20178; } double f_id(double wj, double x) { double r20179 = wj; double r20180 = exp(r20179); double r20181 = r20179 * r20180; double r20182 = x; double r20183 = r20181 - r20182; double r20184 = r20180 + r20181; double r20185 = r20183 / r20184; double r20186 = r20179 - r20185; return r20186; } double f_of(float wj, float x) { float r20187 = wj; float r20188 = exp(r20187); float r20189 = r20187 * r20188; float r20190 = x; float r20191 = r20189 - r20190; float r20192 = r20188 + r20189; float r20193 = r20191 / r20192; float r20194 = r20187 - r20193; float r20195 = 0.02884620428085327f; bool r20196 = r20194 <= r20195; float r20197 = 1.0f; float r20198 = r20197 + r20187; float r20199 = r20187 / r20198; float r20200 = r20187 - r20199; float r20201 = r20190 / r20192; float r20202 = r20200 + r20201; float r20203 = r20196 ? r20194 : r20202; return r20203; } double f_od(double wj, double x) { double r20204 = wj; double r20205 = exp(r20204); double r20206 = r20204 * r20205; double r20207 = x; double r20208 = r20206 - r20207; double r20209 = r20205 + r20206; double r20210 = r20208 / r20209; double r20211 = r20204 - r20210; double r20212 = 0.02884620428085327; bool r20213 = r20211 <= r20212; double r20214 = 1.0; double r20215 = r20214 + r20204; double r20216 = r20204 / r20215; double r20217 = r20204 - r20216; double r20218 = r20207 / r20209; double r20219 = r20217 + r20218; double r20220 = r20213 ? r20211 : r20219; return r20220; } 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 r20221, r20222, r20223, r20224, r20225, r20226, r20227, r20228; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r20221); mpfr_init(r20222); mpfr_init(r20223); mpfr_init(r20224); mpfr_init(r20225); mpfr_init(r20226); mpfr_init(r20227); mpfr_init(r20228); } double f_im(double wj, double x) { mpfr_set_d(r20221, wj, MPFR_RNDN); mpfr_exp(r20222, r20221, MPFR_RNDN); mpfr_mul(r20223, r20221, r20222, MPFR_RNDN); mpfr_set_d(r20224, x, MPFR_RNDN); mpfr_sub(r20225, r20223, r20224, MPFR_RNDN); mpfr_add(r20226, r20222, r20223, MPFR_RNDN); mpfr_div(r20227, r20225, r20226, MPFR_RNDN); mpfr_sub(r20228, r20221, r20227, MPFR_RNDN); return mpfr_get_d(r20228, MPFR_RNDN); } static mpfr_t r20229, r20230, r20231, r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242, r20243, r20244, r20245; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r20229); mpfr_init(r20230); mpfr_init(r20231); mpfr_init(r20232); mpfr_init(r20233); mpfr_init(r20234); mpfr_init(r20235); mpfr_init(r20236); mpfr_init_set_str(r20237, "0.028846204f0", 10, MPFR_RNDN); mpfr_init(r20238); mpfr_init_set_str(r20239, "1", 10, MPFR_RNDN); mpfr_init(r20240); mpfr_init(r20241); mpfr_init(r20242); mpfr_init(r20243); mpfr_init(r20244); mpfr_init(r20245); } double f_fm(double wj, double x) { mpfr_set_d(r20229, wj, MPFR_RNDN); mpfr_exp(r20230, r20229, MPFR_RNDN); mpfr_mul(r20231, r20229, r20230, MPFR_RNDN); mpfr_set_d(r20232, x, MPFR_RNDN); mpfr_sub(r20233, r20231, r20232, MPFR_RNDN); mpfr_add(r20234, r20230, r20231, MPFR_RNDN); mpfr_div(r20235, r20233, r20234, MPFR_RNDN); mpfr_sub(r20236, r20229, r20235, MPFR_RNDN); ; mpfr_set_si(r20238, mpfr_cmp(r20236, r20237) <= 0, MPFR_RNDN); ; mpfr_add(r20240, r20239, r20229, MPFR_RNDN); mpfr_div(r20241, r20229, r20240, MPFR_RNDN); mpfr_sub(r20242, r20229, r20241, MPFR_RNDN); mpfr_div(r20243, r20232, r20234, MPFR_RNDN); mpfr_add(r20244, r20242, r20243, MPFR_RNDN); if (mpfr_get_si(r20238, MPFR_RNDN)) { mpfr_set(r20245, r20236, MPFR_RNDN); } else { mpfr_set(r20245, r20244, MPFR_RNDN); }; return mpfr_get_d(r20245, MPFR_RNDN); } static mpfr_t r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255, r20256, r20257, r20258, r20259, r20260, r20261, r20262; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r20246); mpfr_init(r20247); mpfr_init(r20248); mpfr_init(r20249); mpfr_init(r20250); mpfr_init(r20251); mpfr_init(r20252); mpfr_init(r20253); mpfr_init_set_str(r20254, "0.028846204f0", 10, MPFR_RNDN); mpfr_init(r20255); mpfr_init_set_str(r20256, "1", 10, MPFR_RNDN); mpfr_init(r20257); mpfr_init(r20258); mpfr_init(r20259); mpfr_init(r20260); mpfr_init(r20261); mpfr_init(r20262); } double f_dm(double wj, double x) { mpfr_set_d(r20246, wj, MPFR_RNDN); mpfr_exp(r20247, r20246, MPFR_RNDN); mpfr_mul(r20248, r20246, r20247, MPFR_RNDN); mpfr_set_d(r20249, x, MPFR_RNDN); mpfr_sub(r20250, r20248, r20249, MPFR_RNDN); mpfr_add(r20251, r20247, r20248, MPFR_RNDN); mpfr_div(r20252, r20250, r20251, MPFR_RNDN); mpfr_sub(r20253, r20246, r20252, MPFR_RNDN); ; mpfr_set_si(r20255, mpfr_cmp(r20253, r20254) <= 0, MPFR_RNDN); ; mpfr_add(r20257, r20256, r20246, MPFR_RNDN); mpfr_div(r20258, r20246, r20257, MPFR_RNDN); mpfr_sub(r20259, r20246, r20258, MPFR_RNDN); mpfr_div(r20260, r20249, r20251, MPFR_RNDN); mpfr_add(r20261, r20259, r20260, MPFR_RNDN); if (mpfr_get_si(r20255, MPFR_RNDN)) { mpfr_set(r20262, r20253, MPFR_RNDN); } else { mpfr_set(r20262, r20261, MPFR_RNDN); }; return mpfr_get_d(r20262, MPFR_RNDN); }