#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19191 = wj; float r19192 = exp(r19191); float r19193 = r19191 * r19192; float r19194 = x; float r19195 = r19193 - r19194; float r19196 = r19192 + r19193; float r19197 = r19195 / r19196; float r19198 = r19191 - r19197; return r19198; } double f_id(double wj, double x) { double r19199 = wj; double r19200 = exp(r19199); double r19201 = r19199 * r19200; double r19202 = x; double r19203 = r19201 - r19202; double r19204 = r19200 + r19201; double r19205 = r19203 / r19204; double r19206 = r19199 - r19205; return r19206; } double f_of(float wj, float x) { float r19207 = wj; float r19208 = exp(r19207); float r19209 = r19207 * r19208; float r19210 = x; float r19211 = r19209 - r19210; float r19212 = r19208 + r19209; float r19213 = r19211 / r19212; float r19214 = r19207 - r19213; float r19215 = 0.07101290673017502f; bool r19216 = r19214 <= r19215; float r19217 = 1.0f; float r19218 = r19217 + r19207; float r19219 = r19207 / r19218; float r19220 = r19207 - r19219; float r19221 = r19210 / r19212; float r19222 = r19220 + r19221; float r19223 = r19216 ? r19214 : r19222; return r19223; } double f_od(double wj, double x) { double r19224 = wj; double r19225 = exp(r19224); double r19226 = r19224 * r19225; double r19227 = x; double r19228 = r19226 - r19227; double r19229 = r19225 + r19226; double r19230 = r19228 / r19229; double r19231 = r19224 - r19230; double r19232 = 0.07101290673017502; bool r19233 = r19231 <= r19232; double r19234 = 1.0; double r19235 = r19234 + r19224; double r19236 = r19224 / r19235; double r19237 = r19224 - r19236; double r19238 = r19227 / r19229; double r19239 = r19237 + r19238; double r19240 = r19233 ? r19231 : r19239; return r19240; } 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 r19241, r19242, r19243, r19244, r19245, r19246, r19247, r19248; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19241); mpfr_init(r19242); mpfr_init(r19243); mpfr_init(r19244); mpfr_init(r19245); mpfr_init(r19246); mpfr_init(r19247); mpfr_init(r19248); } double f_im(double wj, double x) { mpfr_set_d(r19241, wj, MPFR_RNDN); mpfr_exp(r19242, r19241, MPFR_RNDN); mpfr_mul(r19243, r19241, r19242, MPFR_RNDN); mpfr_set_d(r19244, x, MPFR_RNDN); mpfr_sub(r19245, r19243, r19244, MPFR_RNDN); mpfr_add(r19246, r19242, r19243, MPFR_RNDN); mpfr_div(r19247, r19245, r19246, MPFR_RNDN); mpfr_sub(r19248, r19241, r19247, MPFR_RNDN); return mpfr_get_d(r19248, MPFR_RNDN); } static mpfr_t r19249, r19250, r19251, r19252, r19253, r19254, r19255, r19256, r19257, r19258, r19259, r19260, r19261, r19262, r19263, r19264, r19265; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19249); mpfr_init(r19250); mpfr_init(r19251); mpfr_init(r19252); mpfr_init(r19253); mpfr_init(r19254); mpfr_init(r19255); mpfr_init(r19256); mpfr_init_set_str(r19257, "0.07101291f0", 10, MPFR_RNDN); mpfr_init(r19258); mpfr_init_set_str(r19259, "1", 10, MPFR_RNDN); mpfr_init(r19260); mpfr_init(r19261); mpfr_init(r19262); mpfr_init(r19263); mpfr_init(r19264); mpfr_init(r19265); } double f_fm(double wj, double x) { mpfr_set_d(r19249, wj, MPFR_RNDN); mpfr_exp(r19250, r19249, MPFR_RNDN); mpfr_mul(r19251, r19249, r19250, MPFR_RNDN); mpfr_set_d(r19252, x, MPFR_RNDN); mpfr_sub(r19253, r19251, r19252, MPFR_RNDN); mpfr_add(r19254, r19250, r19251, MPFR_RNDN); mpfr_div(r19255, r19253, r19254, MPFR_RNDN); mpfr_sub(r19256, r19249, r19255, MPFR_RNDN); ; mpfr_set_si(r19258, mpfr_cmp(r19256, r19257) <= 0, MPFR_RNDN); ; mpfr_add(r19260, r19259, r19249, MPFR_RNDN); mpfr_div(r19261, r19249, r19260, MPFR_RNDN); mpfr_sub(r19262, r19249, r19261, MPFR_RNDN); mpfr_div(r19263, r19252, r19254, MPFR_RNDN); mpfr_add(r19264, r19262, r19263, MPFR_RNDN); if (mpfr_get_si(r19258, MPFR_RNDN)) { mpfr_set(r19265, r19256, MPFR_RNDN); } else { mpfr_set(r19265, r19264, MPFR_RNDN); }; return mpfr_get_d(r19265, MPFR_RNDN); } static mpfr_t r19266, r19267, r19268, r19269, r19270, r19271, r19272, r19273, r19274, r19275, r19276, r19277, r19278, r19279, r19280, r19281, r19282; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19266); mpfr_init(r19267); mpfr_init(r19268); mpfr_init(r19269); mpfr_init(r19270); mpfr_init(r19271); mpfr_init(r19272); mpfr_init(r19273); mpfr_init_set_str(r19274, "0.07101291f0", 10, MPFR_RNDN); mpfr_init(r19275); mpfr_init_set_str(r19276, "1", 10, MPFR_RNDN); mpfr_init(r19277); mpfr_init(r19278); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init(r19282); } double f_dm(double wj, double x) { mpfr_set_d(r19266, wj, MPFR_RNDN); mpfr_exp(r19267, r19266, MPFR_RNDN); mpfr_mul(r19268, r19266, r19267, MPFR_RNDN); mpfr_set_d(r19269, x, MPFR_RNDN); mpfr_sub(r19270, r19268, r19269, MPFR_RNDN); mpfr_add(r19271, r19267, r19268, MPFR_RNDN); mpfr_div(r19272, r19270, r19271, MPFR_RNDN); mpfr_sub(r19273, r19266, r19272, MPFR_RNDN); ; mpfr_set_si(r19275, mpfr_cmp(r19273, r19274) <= 0, MPFR_RNDN); ; mpfr_add(r19277, r19276, r19266, MPFR_RNDN); mpfr_div(r19278, r19266, r19277, MPFR_RNDN); mpfr_sub(r19279, r19266, r19278, MPFR_RNDN); mpfr_div(r19280, r19269, r19271, MPFR_RNDN); mpfr_add(r19281, r19279, r19280, MPFR_RNDN); if (mpfr_get_si(r19275, MPFR_RNDN)) { mpfr_set(r19282, r19273, MPFR_RNDN); } else { mpfr_set(r19282, r19281, MPFR_RNDN); }; return mpfr_get_d(r19282, MPFR_RNDN); }