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