#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r18711 = wj; float r18712 = exp(r18711); float r18713 = r18711 * r18712; float r18714 = x; float r18715 = r18713 - r18714; float r18716 = r18712 + r18713; float r18717 = r18715 / r18716; float r18718 = r18711 - r18717; return r18718; } double f_id(double wj, double x) { double r18719 = wj; double r18720 = exp(r18719); double r18721 = r18719 * r18720; double r18722 = x; double r18723 = r18721 - r18722; double r18724 = r18720 + r18721; double r18725 = r18723 / r18724; double r18726 = r18719 - r18725; return r18726; } double f_of(float wj, float x) { float r18727 = wj; float r18728 = 11.522058486938477f; bool r18729 = r18727 <= r18728; float r18730 = exp(r18727); float r18731 = r18727 * r18730; float r18732 = x; float r18733 = r18731 - r18732; float r18734 = r18730 + r18731; float r18735 = r18733 / r18734; float r18736 = r18727 - r18735; float r18737 = 1.0f; float r18738 = r18727 + r18737; float r18739 = r18727 / r18738; float r18740 = fma(r18727, r18730, r18730); float r18741 = r18732 / r18740; float r18742 = r18739 - r18741; float r18743 = r18727 - r18742; float r18744 = r18729 ? r18736 : r18743; return r18744; } double f_od(double wj, double x) { double r18745 = wj; double r18746 = 11.522058486938477; bool r18747 = r18745 <= r18746; double r18748 = exp(r18745); double r18749 = r18745 * r18748; double r18750 = x; double r18751 = r18749 - r18750; double r18752 = r18748 + r18749; double r18753 = r18751 / r18752; double r18754 = r18745 - r18753; double r18755 = 1.0; double r18756 = r18745 + r18755; double r18757 = r18745 / r18756; double r18758 = fma(r18745, r18748, r18748); double r18759 = r18750 / r18758; double r18760 = r18757 - r18759; double r18761 = r18745 - r18760; double r18762 = r18747 ? r18754 : r18761; return r18762; } 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 r18763, r18764, r18765, r18766, r18767, r18768, r18769, r18770; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18763); mpfr_init(r18764); mpfr_init(r18765); mpfr_init(r18766); mpfr_init(r18767); mpfr_init(r18768); mpfr_init(r18769); mpfr_init(r18770); } double f_im(double wj, double x) { mpfr_set_d(r18763, wj, MPFR_RNDN); mpfr_exp(r18764, r18763, MPFR_RNDN); mpfr_mul(r18765, r18763, r18764, MPFR_RNDN); mpfr_set_d(r18766, x, MPFR_RNDN); mpfr_sub(r18767, r18765, r18766, MPFR_RNDN); mpfr_add(r18768, r18764, r18765, MPFR_RNDN); mpfr_div(r18769, r18767, r18768, MPFR_RNDN); mpfr_sub(r18770, r18763, r18769, MPFR_RNDN); return mpfr_get_d(r18770, MPFR_RNDN); } static mpfr_t r18771, r18772, r18773, r18774, r18775, r18776, r18777, r18778, r18779, r18780, r18781, r18782, r18783, r18784, r18785, r18786, r18787, r18788; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18771); mpfr_init_set_str(r18772, "11.5220585f0", 10, MPFR_RNDN); mpfr_init(r18773); mpfr_init(r18774); mpfr_init(r18775); mpfr_init(r18776); mpfr_init(r18777); mpfr_init(r18778); mpfr_init(r18779); mpfr_init(r18780); mpfr_init_set_str(r18781, "1", 10, MPFR_RNDN); mpfr_init(r18782); mpfr_init(r18783); mpfr_init(r18784); mpfr_init(r18785); mpfr_init(r18786); mpfr_init(r18787); mpfr_init(r18788); } double f_fm(double wj, double x) { mpfr_set_d(r18771, wj, MPFR_RNDN); ; mpfr_set_si(r18773, mpfr_cmp(r18771, r18772) <= 0, MPFR_RNDN); mpfr_exp(r18774, r18771, MPFR_RNDN); mpfr_mul(r18775, r18771, r18774, MPFR_RNDN); mpfr_set_d(r18776, x, MPFR_RNDN); mpfr_sub(r18777, r18775, r18776, MPFR_RNDN); mpfr_add(r18778, r18774, r18775, MPFR_RNDN); mpfr_div(r18779, r18777, r18778, MPFR_RNDN); mpfr_sub(r18780, r18771, r18779, MPFR_RNDN); ; mpfr_add(r18782, r18771, r18781, MPFR_RNDN); mpfr_div(r18783, r18771, r18782, MPFR_RNDN); mpfr_fma(r18784, r18771, r18774, r18774, MPFR_RNDN); mpfr_div(r18785, r18776, r18784, MPFR_RNDN); mpfr_sub(r18786, r18783, r18785, MPFR_RNDN); mpfr_sub(r18787, r18771, r18786, MPFR_RNDN); if (mpfr_get_si(r18773, MPFR_RNDN)) { mpfr_set(r18788, r18780, MPFR_RNDN); } else { mpfr_set(r18788, r18787, MPFR_RNDN); }; return mpfr_get_d(r18788, MPFR_RNDN); } static mpfr_t r18789, r18790, r18791, r18792, r18793, r18794, r18795, r18796, r18797, r18798, r18799, r18800, r18801, r18802, r18803, r18804, r18805, r18806; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18789); mpfr_init_set_str(r18790, "11.5220585f0", 10, MPFR_RNDN); mpfr_init(r18791); mpfr_init(r18792); mpfr_init(r18793); mpfr_init(r18794); mpfr_init(r18795); mpfr_init(r18796); mpfr_init(r18797); mpfr_init(r18798); mpfr_init_set_str(r18799, "1", 10, MPFR_RNDN); mpfr_init(r18800); mpfr_init(r18801); mpfr_init(r18802); mpfr_init(r18803); mpfr_init(r18804); mpfr_init(r18805); mpfr_init(r18806); } double f_dm(double wj, double x) { mpfr_set_d(r18789, wj, MPFR_RNDN); ; mpfr_set_si(r18791, mpfr_cmp(r18789, r18790) <= 0, MPFR_RNDN); mpfr_exp(r18792, r18789, MPFR_RNDN); mpfr_mul(r18793, r18789, r18792, MPFR_RNDN); mpfr_set_d(r18794, x, MPFR_RNDN); mpfr_sub(r18795, r18793, r18794, MPFR_RNDN); mpfr_add(r18796, r18792, r18793, MPFR_RNDN); mpfr_div(r18797, r18795, r18796, MPFR_RNDN); mpfr_sub(r18798, r18789, r18797, MPFR_RNDN); ; mpfr_add(r18800, r18789, r18799, MPFR_RNDN); mpfr_div(r18801, r18789, r18800, MPFR_RNDN); mpfr_fma(r18802, r18789, r18792, r18792, MPFR_RNDN); mpfr_div(r18803, r18794, r18802, MPFR_RNDN); mpfr_sub(r18804, r18801, r18803, MPFR_RNDN); mpfr_sub(r18805, r18789, r18804, MPFR_RNDN); if (mpfr_get_si(r18791, MPFR_RNDN)) { mpfr_set(r18806, r18798, MPFR_RNDN); } else { mpfr_set(r18806, r18805, MPFR_RNDN); }; return mpfr_get_d(r18806, MPFR_RNDN); }