#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19695 = wj; float r19696 = exp(r19695); float r19697 = r19695 * r19696; float r19698 = x; float r19699 = r19697 - r19698; float r19700 = r19696 + r19697; float r19701 = r19699 / r19700; float r19702 = r19695 - r19701; return r19702; } double f_id(double wj, double x) { double r19703 = wj; double r19704 = exp(r19703); double r19705 = r19703 * r19704; double r19706 = x; double r19707 = r19705 - r19706; double r19708 = r19704 + r19705; double r19709 = r19707 / r19708; double r19710 = r19703 - r19709; return r19710; } double f_of(float wj, float x) { float r19711 = wj; float r19712 = 0.4936105012893677f; bool r19713 = r19711 <= r19712; float r19714 = exp(r19711); float r19715 = r19711 * r19714; float r19716 = x; float r19717 = r19715 - r19716; float r19718 = r19714 + r19715; float r19719 = r19717 / r19718; float r19720 = r19711 - r19719; float r19721 = 1.0f; float r19722 = r19721 + r19711; float r19723 = r19711 / r19722; float r19724 = r19711 - r19723; float r19725 = r19716 / r19718; float r19726 = r19724 + r19725; float r19727 = r19713 ? r19720 : r19726; return r19727; } double f_od(double wj, double x) { double r19728 = wj; double r19729 = 0.4936105012893677; bool r19730 = r19728 <= r19729; double r19731 = exp(r19728); double r19732 = r19728 * r19731; double r19733 = x; double r19734 = r19732 - r19733; double r19735 = r19731 + r19732; double r19736 = r19734 / r19735; double r19737 = r19728 - r19736; double r19738 = 1.0; double r19739 = r19738 + r19728; double r19740 = r19728 / r19739; double r19741 = r19728 - r19740; double r19742 = r19733 / r19735; double r19743 = r19741 + r19742; double r19744 = r19730 ? r19737 : r19743; return r19744; } 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 r19745, r19746, r19747, r19748, r19749, r19750, r19751, r19752; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19745); mpfr_init(r19746); mpfr_init(r19747); mpfr_init(r19748); mpfr_init(r19749); mpfr_init(r19750); mpfr_init(r19751); mpfr_init(r19752); } double f_im(double wj, double x) { mpfr_set_d(r19745, wj, MPFR_RNDN); mpfr_exp(r19746, r19745, MPFR_RNDN); mpfr_mul(r19747, r19745, r19746, MPFR_RNDN); mpfr_set_d(r19748, x, MPFR_RNDN); mpfr_sub(r19749, r19747, r19748, MPFR_RNDN); mpfr_add(r19750, r19746, r19747, MPFR_RNDN); mpfr_div(r19751, r19749, r19750, MPFR_RNDN); mpfr_sub(r19752, r19745, r19751, MPFR_RNDN); return mpfr_get_d(r19752, MPFR_RNDN); } static mpfr_t r19753, r19754, r19755, r19756, r19757, r19758, r19759, r19760, r19761, r19762, r19763, r19764, r19765, r19766, r19767, r19768, r19769; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19753); mpfr_init_set_str(r19754, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19755); mpfr_init(r19756); mpfr_init(r19757); mpfr_init(r19758); mpfr_init(r19759); mpfr_init(r19760); mpfr_init(r19761); mpfr_init(r19762); mpfr_init_set_str(r19763, "1", 10, MPFR_RNDN); mpfr_init(r19764); mpfr_init(r19765); mpfr_init(r19766); mpfr_init(r19767); mpfr_init(r19768); mpfr_init(r19769); } double f_fm(double wj, double x) { mpfr_set_d(r19753, wj, MPFR_RNDN); ; mpfr_set_si(r19755, mpfr_cmp(r19753, r19754) <= 0, MPFR_RNDN); mpfr_exp(r19756, r19753, MPFR_RNDN); mpfr_mul(r19757, r19753, r19756, MPFR_RNDN); mpfr_set_d(r19758, x, MPFR_RNDN); mpfr_sub(r19759, r19757, r19758, MPFR_RNDN); mpfr_add(r19760, r19756, r19757, MPFR_RNDN); mpfr_div(r19761, r19759, r19760, MPFR_RNDN); mpfr_sub(r19762, r19753, r19761, MPFR_RNDN); ; mpfr_add(r19764, r19763, r19753, MPFR_RNDN); mpfr_div(r19765, r19753, r19764, MPFR_RNDN); mpfr_sub(r19766, r19753, r19765, MPFR_RNDN); mpfr_div(r19767, r19758, r19760, MPFR_RNDN); mpfr_add(r19768, r19766, r19767, MPFR_RNDN); if (mpfr_get_si(r19755, MPFR_RNDN)) { mpfr_set(r19769, r19762, MPFR_RNDN); } else { mpfr_set(r19769, r19768, MPFR_RNDN); }; return mpfr_get_d(r19769, MPFR_RNDN); } static mpfr_t r19770, r19771, r19772, r19773, r19774, r19775, r19776, r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784, r19785, r19786; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19770); mpfr_init_set_str(r19771, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19772); mpfr_init(r19773); mpfr_init(r19774); mpfr_init(r19775); mpfr_init(r19776); mpfr_init(r19777); mpfr_init(r19778); mpfr_init(r19779); mpfr_init_set_str(r19780, "1", 10, MPFR_RNDN); mpfr_init(r19781); mpfr_init(r19782); mpfr_init(r19783); mpfr_init(r19784); mpfr_init(r19785); mpfr_init(r19786); } double f_dm(double wj, double x) { mpfr_set_d(r19770, wj, MPFR_RNDN); ; mpfr_set_si(r19772, mpfr_cmp(r19770, r19771) <= 0, MPFR_RNDN); mpfr_exp(r19773, r19770, MPFR_RNDN); mpfr_mul(r19774, r19770, r19773, MPFR_RNDN); mpfr_set_d(r19775, x, MPFR_RNDN); mpfr_sub(r19776, r19774, r19775, MPFR_RNDN); mpfr_add(r19777, r19773, r19774, MPFR_RNDN); mpfr_div(r19778, r19776, r19777, MPFR_RNDN); mpfr_sub(r19779, r19770, r19778, MPFR_RNDN); ; mpfr_add(r19781, r19780, r19770, MPFR_RNDN); mpfr_div(r19782, r19770, r19781, MPFR_RNDN); mpfr_sub(r19783, r19770, r19782, MPFR_RNDN); mpfr_div(r19784, r19775, r19777, MPFR_RNDN); mpfr_add(r19785, r19783, r19784, MPFR_RNDN); if (mpfr_get_si(r19772, MPFR_RNDN)) { mpfr_set(r19786, r19779, MPFR_RNDN); } else { mpfr_set(r19786, r19785, MPFR_RNDN); }; return mpfr_get_d(r19786, MPFR_RNDN); }