#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r18955 = wj; float r18956 = exp(r18955); float r18957 = r18955 * r18956; float r18958 = x; float r18959 = r18957 - r18958; float r18960 = r18956 + r18957; float r18961 = r18959 / r18960; float r18962 = r18955 - r18961; return r18962; } double f_id(double wj, double x) { double r18963 = wj; double r18964 = exp(r18963); double r18965 = r18963 * r18964; double r18966 = x; double r18967 = r18965 - r18966; double r18968 = r18964 + r18965; double r18969 = r18967 / r18968; double r18970 = r18963 - r18969; return r18970; } double f_of(float wj, float x) { float r18971 = wj; float r18972 = exp(r18971); float r18973 = r18971 * r18972; float r18974 = x; float r18975 = r18973 - r18974; float r18976 = r18972 + r18973; float r18977 = r18975 / r18976; float r18978 = 0.9145612716674805f; bool r18979 = r18977 <= r18978; float r18980 = r18971 - r18977; float r18981 = 1.0f; float r18982 = r18971 + r18981; float r18983 = r18971 / r18982; float r18984 = r18981 + r18971; float r18985 = r18974 / r18984; float r18986 = r18985 / r18972; float r18987 = r18983 - r18986; float r18988 = r18971 - r18987; float r18989 = r18979 ? r18980 : r18988; return r18989; } double f_od(double wj, double x) { double r18990 = wj; double r18991 = exp(r18990); double r18992 = r18990 * r18991; double r18993 = x; double r18994 = r18992 - r18993; double r18995 = r18991 + r18992; double r18996 = r18994 / r18995; double r18997 = 0.9145612716674805; bool r18998 = r18996 <= r18997; double r18999 = r18990 - r18996; double r19000 = 1.0; double r19001 = r18990 + r19000; double r19002 = r18990 / r19001; double r19003 = r19000 + r18990; double r19004 = r18993 / r19003; double r19005 = r19004 / r18991; double r19006 = r19002 - r19005; double r19007 = r18990 - r19006; double r19008 = r18998 ? r18999 : r19007; return r19008; } 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 r19009, r19010, r19011, r19012, r19013, r19014, r19015, r19016; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19009); mpfr_init(r19010); mpfr_init(r19011); mpfr_init(r19012); mpfr_init(r19013); mpfr_init(r19014); mpfr_init(r19015); mpfr_init(r19016); } double f_im(double wj, double x) { mpfr_set_d(r19009, wj, MPFR_RNDN); mpfr_exp(r19010, r19009, MPFR_RNDN); mpfr_mul(r19011, r19009, r19010, MPFR_RNDN); mpfr_set_d(r19012, x, MPFR_RNDN); mpfr_sub(r19013, r19011, r19012, MPFR_RNDN); mpfr_add(r19014, r19010, r19011, MPFR_RNDN); mpfr_div(r19015, r19013, r19014, MPFR_RNDN); mpfr_sub(r19016, r19009, r19015, MPFR_RNDN); return mpfr_get_d(r19016, MPFR_RNDN); } static mpfr_t r19017, r19018, r19019, r19020, r19021, r19022, r19023, r19024, r19025, r19026, r19027, r19028, r19029, r19030, r19031, r19032, r19033, r19034, r19035; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19017); mpfr_init(r19018); mpfr_init(r19019); mpfr_init(r19020); mpfr_init(r19021); mpfr_init(r19022); mpfr_init(r19023); mpfr_init_set_str(r19024, "0.9145613f0", 10, MPFR_RNDN); mpfr_init(r19025); mpfr_init(r19026); mpfr_init_set_str(r19027, "1", 10, MPFR_RNDN); mpfr_init(r19028); mpfr_init(r19029); mpfr_init(r19030); mpfr_init(r19031); mpfr_init(r19032); mpfr_init(r19033); mpfr_init(r19034); mpfr_init(r19035); } double f_fm(double wj, double x) { mpfr_set_d(r19017, wj, MPFR_RNDN); mpfr_exp(r19018, r19017, MPFR_RNDN); mpfr_mul(r19019, r19017, r19018, MPFR_RNDN); mpfr_set_d(r19020, x, MPFR_RNDN); mpfr_sub(r19021, r19019, r19020, MPFR_RNDN); mpfr_add(r19022, r19018, r19019, MPFR_RNDN); mpfr_div(r19023, r19021, r19022, MPFR_RNDN); ; mpfr_set_si(r19025, mpfr_cmp(r19023, r19024) <= 0, MPFR_RNDN); mpfr_sub(r19026, r19017, r19023, MPFR_RNDN); ; mpfr_add(r19028, r19017, r19027, MPFR_RNDN); mpfr_div(r19029, r19017, r19028, MPFR_RNDN); mpfr_add(r19030, r19027, r19017, MPFR_RNDN); mpfr_div(r19031, r19020, r19030, MPFR_RNDN); mpfr_div(r19032, r19031, r19018, MPFR_RNDN); mpfr_sub(r19033, r19029, r19032, MPFR_RNDN); mpfr_sub(r19034, r19017, r19033, MPFR_RNDN); if (mpfr_get_si(r19025, MPFR_RNDN)) { mpfr_set(r19035, r19026, MPFR_RNDN); } else { mpfr_set(r19035, r19034, MPFR_RNDN); }; return mpfr_get_d(r19035, MPFR_RNDN); } static mpfr_t r19036, r19037, r19038, r19039, r19040, r19041, r19042, r19043, r19044, r19045, r19046, r19047, r19048, r19049, r19050, r19051, r19052, r19053, r19054; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19036); mpfr_init(r19037); mpfr_init(r19038); mpfr_init(r19039); mpfr_init(r19040); mpfr_init(r19041); mpfr_init(r19042); mpfr_init_set_str(r19043, "0.9145613f0", 10, MPFR_RNDN); mpfr_init(r19044); mpfr_init(r19045); mpfr_init_set_str(r19046, "1", 10, MPFR_RNDN); mpfr_init(r19047); mpfr_init(r19048); mpfr_init(r19049); mpfr_init(r19050); mpfr_init(r19051); mpfr_init(r19052); mpfr_init(r19053); mpfr_init(r19054); } double f_dm(double wj, double x) { mpfr_set_d(r19036, wj, MPFR_RNDN); mpfr_exp(r19037, r19036, MPFR_RNDN); mpfr_mul(r19038, r19036, r19037, MPFR_RNDN); mpfr_set_d(r19039, x, MPFR_RNDN); mpfr_sub(r19040, r19038, r19039, MPFR_RNDN); mpfr_add(r19041, r19037, r19038, MPFR_RNDN); mpfr_div(r19042, r19040, r19041, MPFR_RNDN); ; mpfr_set_si(r19044, mpfr_cmp(r19042, r19043) <= 0, MPFR_RNDN); mpfr_sub(r19045, r19036, r19042, MPFR_RNDN); ; mpfr_add(r19047, r19036, r19046, MPFR_RNDN); mpfr_div(r19048, r19036, r19047, MPFR_RNDN); mpfr_add(r19049, r19046, r19036, MPFR_RNDN); mpfr_div(r19050, r19039, r19049, MPFR_RNDN); mpfr_div(r19051, r19050, r19037, MPFR_RNDN); mpfr_sub(r19052, r19048, r19051, MPFR_RNDN); mpfr_sub(r19053, r19036, r19052, MPFR_RNDN); if (mpfr_get_si(r19044, MPFR_RNDN)) { mpfr_set(r19054, r19045, MPFR_RNDN); } else { mpfr_set(r19054, r19053, MPFR_RNDN); }; return mpfr_get_d(r19054, MPFR_RNDN); }