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