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