#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19227 = wj; float r19228 = exp(r19227); float r19229 = r19227 * r19228; float r19230 = x; float r19231 = r19229 - r19230; float r19232 = r19228 + r19229; float r19233 = r19231 / r19232; float r19234 = r19227 - r19233; return r19234; } double f_id(double wj, double x) { double r19235 = wj; double r19236 = exp(r19235); double r19237 = r19235 * r19236; double r19238 = x; double r19239 = r19237 - r19238; double r19240 = r19236 + r19237; double r19241 = r19239 / r19240; double r19242 = r19235 - r19241; return r19242; } double f_of(float wj, float x) { float r19243 = wj; float r19244 = 0.4936105012893677f; bool r19245 = r19243 <= r19244; float r19246 = exp(r19243); float r19247 = r19243 * r19246; float r19248 = x; float r19249 = r19247 - r19248; float r19250 = r19246 + r19247; float r19251 = r19249 / r19250; float r19252 = r19243 - r19251; float r19253 = 1.0f; float r19254 = r19253 + r19243; float r19255 = r19243 / r19254; float r19256 = r19243 - r19255; float r19257 = r19248 / r19250; float r19258 = r19256 + r19257; float r19259 = r19245 ? r19252 : r19258; return r19259; } double f_od(double wj, double x) { double r19260 = wj; double r19261 = 0.4936105012893677; bool r19262 = r19260 <= r19261; double r19263 = exp(r19260); double r19264 = r19260 * r19263; double r19265 = x; double r19266 = r19264 - r19265; double r19267 = r19263 + r19264; double r19268 = r19266 / r19267; double r19269 = r19260 - r19268; double r19270 = 1.0; double r19271 = r19270 + r19260; double r19272 = r19260 / r19271; double r19273 = r19260 - r19272; double r19274 = r19265 / r19267; double r19275 = r19273 + r19274; double r19276 = r19262 ? r19269 : r19275; return r19276; } 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 r19277, r19278, r19279, r19280, r19281, r19282, r19283, r19284; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19277); mpfr_init(r19278); mpfr_init(r19279); mpfr_init(r19280); mpfr_init(r19281); mpfr_init(r19282); mpfr_init(r19283); mpfr_init(r19284); } double f_im(double wj, double x) { mpfr_set_d(r19277, wj, MPFR_RNDN); mpfr_exp(r19278, r19277, MPFR_RNDN); mpfr_mul(r19279, r19277, r19278, MPFR_RNDN); mpfr_set_d(r19280, x, MPFR_RNDN); mpfr_sub(r19281, r19279, r19280, MPFR_RNDN); mpfr_add(r19282, r19278, r19279, MPFR_RNDN); mpfr_div(r19283, r19281, r19282, MPFR_RNDN); mpfr_sub(r19284, r19277, r19283, MPFR_RNDN); return mpfr_get_d(r19284, MPFR_RNDN); } static mpfr_t r19285, r19286, r19287, r19288, r19289, r19290, r19291, r19292, r19293, r19294, r19295, r19296, r19297, r19298, r19299, r19300, r19301; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19285); mpfr_init_set_str(r19286, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19287); mpfr_init(r19288); mpfr_init(r19289); mpfr_init(r19290); mpfr_init(r19291); mpfr_init(r19292); mpfr_init(r19293); mpfr_init(r19294); mpfr_init_set_str(r19295, "1", 10, MPFR_RNDN); mpfr_init(r19296); mpfr_init(r19297); mpfr_init(r19298); mpfr_init(r19299); mpfr_init(r19300); mpfr_init(r19301); } double f_fm(double wj, double x) { mpfr_set_d(r19285, wj, MPFR_RNDN); ; mpfr_set_si(r19287, mpfr_cmp(r19285, r19286) <= 0, MPFR_RNDN); mpfr_exp(r19288, r19285, MPFR_RNDN); mpfr_mul(r19289, r19285, r19288, MPFR_RNDN); mpfr_set_d(r19290, x, MPFR_RNDN); mpfr_sub(r19291, r19289, r19290, MPFR_RNDN); mpfr_add(r19292, r19288, r19289, MPFR_RNDN); mpfr_div(r19293, r19291, r19292, MPFR_RNDN); mpfr_sub(r19294, r19285, r19293, MPFR_RNDN); ; mpfr_add(r19296, r19295, r19285, MPFR_RNDN); mpfr_div(r19297, r19285, r19296, MPFR_RNDN); mpfr_sub(r19298, r19285, r19297, MPFR_RNDN); mpfr_div(r19299, r19290, r19292, MPFR_RNDN); mpfr_add(r19300, r19298, r19299, MPFR_RNDN); if (mpfr_get_si(r19287, MPFR_RNDN)) { mpfr_set(r19301, r19294, MPFR_RNDN); } else { mpfr_set(r19301, r19300, MPFR_RNDN); }; return mpfr_get_d(r19301, MPFR_RNDN); } static mpfr_t r19302, r19303, r19304, r19305, r19306, r19307, r19308, r19309, r19310, r19311, r19312, r19313, r19314, r19315, r19316, r19317, r19318; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19302); mpfr_init_set_str(r19303, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19304); mpfr_init(r19305); mpfr_init(r19306); mpfr_init(r19307); mpfr_init(r19308); mpfr_init(r19309); mpfr_init(r19310); mpfr_init(r19311); mpfr_init_set_str(r19312, "1", 10, MPFR_RNDN); mpfr_init(r19313); mpfr_init(r19314); mpfr_init(r19315); mpfr_init(r19316); mpfr_init(r19317); mpfr_init(r19318); } double f_dm(double wj, double x) { mpfr_set_d(r19302, wj, MPFR_RNDN); ; mpfr_set_si(r19304, mpfr_cmp(r19302, r19303) <= 0, MPFR_RNDN); mpfr_exp(r19305, r19302, MPFR_RNDN); mpfr_mul(r19306, r19302, r19305, MPFR_RNDN); mpfr_set_d(r19307, x, MPFR_RNDN); mpfr_sub(r19308, r19306, r19307, MPFR_RNDN); mpfr_add(r19309, r19305, r19306, MPFR_RNDN); mpfr_div(r19310, r19308, r19309, MPFR_RNDN); mpfr_sub(r19311, r19302, r19310, MPFR_RNDN); ; mpfr_add(r19313, r19312, r19302, MPFR_RNDN); mpfr_div(r19314, r19302, r19313, MPFR_RNDN); mpfr_sub(r19315, r19302, r19314, MPFR_RNDN); mpfr_div(r19316, r19307, r19309, MPFR_RNDN); mpfr_add(r19317, r19315, r19316, MPFR_RNDN); if (mpfr_get_si(r19304, MPFR_RNDN)) { mpfr_set(r19318, r19311, MPFR_RNDN); } else { mpfr_set(r19318, r19317, MPFR_RNDN); }; return mpfr_get_d(r19318, MPFR_RNDN); }