#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r25709 = wj; float r25710 = exp(r25709); float r25711 = r25709 * r25710; float r25712 = x; float r25713 = r25711 - r25712; float r25714 = r25710 + r25711; float r25715 = r25713 / r25714; float r25716 = r25709 - r25715; return r25716; } double f_id(double wj, double x) { double r25717 = wj; double r25718 = exp(r25717); double r25719 = r25717 * r25718; double r25720 = x; double r25721 = r25719 - r25720; double r25722 = r25718 + r25719; double r25723 = r25721 / r25722; double r25724 = r25717 - r25723; return r25724; } double f_of(float wj, float x) { float r25725 = x; float r25726 = 2; float r25727 = r25725 * r25726; float r25728 = wj; float r25729 = -r25728; float r25730 = fma(r25728, r25728, r25725); float r25731 = fma(r25727, r25729, r25730); float r25732 = 1.0428080602965551e-26; bool r25733 = r25731 <= r25732; float r25734 = 1; float r25735 = r25728 + r25734; float r25736 = r25734 / r25735; float r25737 = exp(r25728); float r25738 = r25725 / r25737; float r25739 = r25728 - r25738; float r25740 = r25736 * r25739; float r25741 = r25728 - r25740; float r25742 = r25733 ? r25731 : r25741; return r25742; } double f_od(double wj, double x) { double r25743 = x; double r25744 = 2; double r25745 = r25743 * r25744; double r25746 = wj; double r25747 = -r25746; double r25748 = fma(r25746, r25746, r25743); double r25749 = fma(r25745, r25747, r25748); double r25750 = 1.0428080602965551e-26; bool r25751 = r25749 <= r25750; double r25752 = 1; double r25753 = r25746 + r25752; double r25754 = r25752 / r25753; double r25755 = exp(r25746); double r25756 = r25743 / r25755; double r25757 = r25746 - r25756; double r25758 = r25754 * r25757; double r25759 = r25746 - r25758; double r25760 = r25751 ? r25749 : r25759; return r25760; } 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 r25761, r25762, r25763, r25764, r25765, r25766, r25767, r25768; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r25761); mpfr_init(r25762); mpfr_init(r25763); mpfr_init(r25764); mpfr_init(r25765); mpfr_init(r25766); mpfr_init(r25767); mpfr_init(r25768); } double f_im(double wj, double x) { mpfr_set_d(r25761, wj, MPFR_RNDN); mpfr_exp(r25762, r25761, MPFR_RNDN); mpfr_mul(r25763, r25761, r25762, MPFR_RNDN); mpfr_set_d(r25764, x, MPFR_RNDN); mpfr_sub(r25765, r25763, r25764, MPFR_RNDN); mpfr_add(r25766, r25762, r25763, MPFR_RNDN); mpfr_div(r25767, r25765, r25766, MPFR_RNDN); mpfr_sub(r25768, r25761, r25767, MPFR_RNDN); return mpfr_get_d(r25768, MPFR_RNDN); } static mpfr_t r25769, r25770, r25771, r25772, r25773, r25774, r25775, r25776, r25777, r25778, r25779, r25780, r25781, r25782, r25783, r25784, r25785, r25786; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r25769); mpfr_init_set_str(r25770, "2", 10, MPFR_RNDN); mpfr_init(r25771); mpfr_init(r25772); mpfr_init(r25773); mpfr_init(r25774); mpfr_init(r25775); mpfr_init_set_str(r25776, "1.0428080602965551e-26", 10, MPFR_RNDN); mpfr_init(r25777); mpfr_init_set_str(r25778, "1", 10, MPFR_RNDN); mpfr_init(r25779); mpfr_init(r25780); mpfr_init(r25781); mpfr_init(r25782); mpfr_init(r25783); mpfr_init(r25784); mpfr_init(r25785); mpfr_init(r25786); } double f_fm(double wj, double x) { mpfr_set_d(r25769, x, MPFR_RNDN); ; mpfr_mul(r25771, r25769, r25770, MPFR_RNDN); mpfr_set_d(r25772, wj, MPFR_RNDN); mpfr_neg(r25773, r25772, MPFR_RNDN); mpfr_fma(r25774, r25772, r25772, r25769, MPFR_RNDN); mpfr_fma(r25775, r25771, r25773, r25774, MPFR_RNDN); ; mpfr_set_si(r25777, mpfr_cmp(r25775, r25776) <= 0, MPFR_RNDN); ; mpfr_add(r25779, r25772, r25778, MPFR_RNDN); mpfr_div(r25780, r25778, r25779, MPFR_RNDN); mpfr_exp(r25781, r25772, MPFR_RNDN); mpfr_div(r25782, r25769, r25781, MPFR_RNDN); mpfr_sub(r25783, r25772, r25782, MPFR_RNDN); mpfr_mul(r25784, r25780, r25783, MPFR_RNDN); mpfr_sub(r25785, r25772, r25784, MPFR_RNDN); if (mpfr_get_si(r25777, MPFR_RNDN)) { mpfr_set(r25786, r25775, MPFR_RNDN); } else { mpfr_set(r25786, r25785, MPFR_RNDN); }; return mpfr_get_d(r25786, MPFR_RNDN); } static mpfr_t r25787, r25788, r25789, r25790, r25791, r25792, r25793, r25794, r25795, r25796, r25797, r25798, r25799, r25800, r25801, r25802, r25803, r25804; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r25787); mpfr_init_set_str(r25788, "2", 10, MPFR_RNDN); mpfr_init(r25789); mpfr_init(r25790); mpfr_init(r25791); mpfr_init(r25792); mpfr_init(r25793); mpfr_init_set_str(r25794, "1.0428080602965551e-26", 10, MPFR_RNDN); mpfr_init(r25795); mpfr_init_set_str(r25796, "1", 10, MPFR_RNDN); mpfr_init(r25797); mpfr_init(r25798); mpfr_init(r25799); mpfr_init(r25800); mpfr_init(r25801); mpfr_init(r25802); mpfr_init(r25803); mpfr_init(r25804); } double f_dm(double wj, double x) { mpfr_set_d(r25787, x, MPFR_RNDN); ; mpfr_mul(r25789, r25787, r25788, MPFR_RNDN); mpfr_set_d(r25790, wj, MPFR_RNDN); mpfr_neg(r25791, r25790, MPFR_RNDN); mpfr_fma(r25792, r25790, r25790, r25787, MPFR_RNDN); mpfr_fma(r25793, r25789, r25791, r25792, MPFR_RNDN); ; mpfr_set_si(r25795, mpfr_cmp(r25793, r25794) <= 0, MPFR_RNDN); ; mpfr_add(r25797, r25790, r25796, MPFR_RNDN); mpfr_div(r25798, r25796, r25797, MPFR_RNDN); mpfr_exp(r25799, r25790, MPFR_RNDN); mpfr_div(r25800, r25787, r25799, MPFR_RNDN); mpfr_sub(r25801, r25790, r25800, MPFR_RNDN); mpfr_mul(r25802, r25798, r25801, MPFR_RNDN); mpfr_sub(r25803, r25790, r25802, MPFR_RNDN); if (mpfr_get_si(r25795, MPFR_RNDN)) { mpfr_set(r25804, r25793, MPFR_RNDN); } else { mpfr_set(r25804, r25803, MPFR_RNDN); }; return mpfr_get_d(r25804, MPFR_RNDN); }