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