#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(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; return r20158; } double f_id(double wj, double x) { double r20159 = wj; double r20160 = exp(r20159); double r20161 = r20159 * r20160; double r20162 = x; double r20163 = r20161 - r20162; double r20164 = r20160 + r20161; double r20165 = r20163 / r20164; double r20166 = r20159 - r20165; return r20166; } double f_of(float wj, float x) { float r20167 = wj; float r20168 = exp(r20167); float r20169 = r20167 * r20168; float r20170 = x; float r20171 = r20169 - r20170; float r20172 = r20168 + r20169; float r20173 = r20171 / r20172; float r20174 = r20167 - r20173; float r20175 = 0.02884620428085327f; bool r20176 = r20174 <= r20175; float r20177 = 1.0f; float r20178 = r20177 + r20167; float r20179 = r20167 / r20178; float r20180 = r20167 - r20179; float r20181 = r20170 / r20172; float r20182 = r20180 + r20181; float r20183 = r20176 ? r20174 : r20182; return r20183; } double f_od(double wj, double x) { double r20184 = wj; double r20185 = exp(r20184); double r20186 = r20184 * r20185; double r20187 = x; double r20188 = r20186 - r20187; double r20189 = r20185 + r20186; double r20190 = r20188 / r20189; double r20191 = r20184 - r20190; double r20192 = 0.02884620428085327; bool r20193 = r20191 <= r20192; double r20194 = 1.0; double r20195 = r20194 + r20184; double r20196 = r20184 / r20195; double r20197 = r20184 - r20196; double r20198 = r20187 / r20189; double r20199 = r20197 + r20198; double r20200 = r20193 ? r20191 : 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; 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, "0.028846204f0", 10, MPFR_RNDN); mpfr_init(r20218); mpfr_init_set_str(r20219, "1", 10, MPFR_RNDN); mpfr_init(r20220); mpfr_init(r20221); mpfr_init(r20222); mpfr_init(r20223); mpfr_init(r20224); mpfr_init(r20225); } 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_add(r20220, r20219, r20209, MPFR_RNDN); mpfr_div(r20221, r20209, r20220, MPFR_RNDN); mpfr_sub(r20222, r20209, r20221, MPFR_RNDN); mpfr_div(r20223, r20212, r20214, MPFR_RNDN); mpfr_add(r20224, r20222, r20223, MPFR_RNDN); if (mpfr_get_si(r20218, MPFR_RNDN)) { mpfr_set(r20225, r20216, MPFR_RNDN); } else { mpfr_set(r20225, r20224, MPFR_RNDN); }; return mpfr_get_d(r20225, MPFR_RNDN); } static mpfr_t r20226, r20227, r20228, r20229, r20230, r20231, r20232, r20233, r20234, r20235, r20236, r20237, r20238, r20239, r20240, r20241, r20242; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r20226); mpfr_init(r20227); mpfr_init(r20228); mpfr_init(r20229); mpfr_init(r20230); mpfr_init(r20231); mpfr_init(r20232); mpfr_init(r20233); mpfr_init_set_str(r20234, "0.028846204f0", 10, MPFR_RNDN); mpfr_init(r20235); mpfr_init_set_str(r20236, "1", 10, MPFR_RNDN); mpfr_init(r20237); mpfr_init(r20238); mpfr_init(r20239); mpfr_init(r20240); mpfr_init(r20241); mpfr_init(r20242); } double f_dm(double wj, double x) { mpfr_set_d(r20226, wj, MPFR_RNDN); mpfr_exp(r20227, r20226, MPFR_RNDN); mpfr_mul(r20228, r20226, r20227, MPFR_RNDN); mpfr_set_d(r20229, x, MPFR_RNDN); mpfr_sub(r20230, r20228, r20229, MPFR_RNDN); mpfr_add(r20231, r20227, r20228, MPFR_RNDN); mpfr_div(r20232, r20230, r20231, MPFR_RNDN); mpfr_sub(r20233, r20226, r20232, MPFR_RNDN); ; mpfr_set_si(r20235, mpfr_cmp(r20233, r20234) <= 0, MPFR_RNDN); ; mpfr_add(r20237, r20236, r20226, MPFR_RNDN); mpfr_div(r20238, r20226, r20237, MPFR_RNDN); mpfr_sub(r20239, r20226, r20238, MPFR_RNDN); mpfr_div(r20240, r20229, r20231, MPFR_RNDN); mpfr_add(r20241, r20239, r20240, MPFR_RNDN); if (mpfr_get_si(r20235, MPFR_RNDN)) { mpfr_set(r20242, r20233, MPFR_RNDN); } else { mpfr_set(r20242, r20241, MPFR_RNDN); }; return mpfr_get_d(r20242, MPFR_RNDN); }