#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19723 = wj; float r19724 = exp(r19723); float r19725 = r19723 * r19724; float r19726 = x; float r19727 = r19725 - r19726; float r19728 = r19724 + r19725; float r19729 = r19727 / r19728; float r19730 = r19723 - r19729; return r19730; } double f_id(double wj, double x) { double r19731 = wj; double r19732 = exp(r19731); double r19733 = r19731 * r19732; double r19734 = x; double r19735 = r19733 - r19734; double r19736 = r19732 + r19733; double r19737 = r19735 / r19736; double r19738 = r19731 - r19737; return r19738; } double f_of(float wj, float x) { float r19739 = wj; float r19740 = exp(r19739); float r19741 = r19739 * r19740; float r19742 = x; float r19743 = r19741 - r19742; float r19744 = r19740 + r19741; float r19745 = r19743 / r19744; float r19746 = r19739 - r19745; float r19747 = 0.0018067137571051717f; bool r19748 = r19746 <= r19747; float r19749 = cbrt(r19746); float r19750 = r19749 * (r19749 * r19749); float r19751 = 1.0f; float r19752 = r19751 + r19739; float r19753 = r19739 / r19752; float r19754 = r19739 - r19753; float r19755 = r19742 / r19744; float r19756 = r19754 + r19755; float r19757 = r19748 ? r19750 : r19756; return r19757; } double f_od(double wj, double x) { double r19758 = wj; double r19759 = exp(r19758); double r19760 = r19758 * r19759; double r19761 = x; double r19762 = r19760 - r19761; double r19763 = r19759 + r19760; double r19764 = r19762 / r19763; double r19765 = r19758 - r19764; double r19766 = 0.0018067137571051717; bool r19767 = r19765 <= r19766; double r19768 = cbrt(r19765); double r19769 = r19768 * (r19768 * r19768); double r19770 = 1.0; double r19771 = r19770 + r19758; double r19772 = r19758 / r19771; double r19773 = r19758 - r19772; double r19774 = r19761 / r19763; double r19775 = r19773 + r19774; double r19776 = r19767 ? r19769 : r19775; return r19776; } 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 r19777, r19778, r19779, r19780, r19781, r19782, r19783, r19784; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19777); mpfr_init(r19778); mpfr_init(r19779); mpfr_init(r19780); mpfr_init(r19781); mpfr_init(r19782); mpfr_init(r19783); mpfr_init(r19784); } double f_im(double wj, double x) { mpfr_set_d(r19777, wj, MPFR_RNDN); mpfr_exp(r19778, r19777, MPFR_RNDN); mpfr_mul(r19779, r19777, r19778, MPFR_RNDN); mpfr_set_d(r19780, x, MPFR_RNDN); mpfr_sub(r19781, r19779, r19780, MPFR_RNDN); mpfr_add(r19782, r19778, r19779, MPFR_RNDN); mpfr_div(r19783, r19781, r19782, MPFR_RNDN); mpfr_sub(r19784, r19777, r19783, MPFR_RNDN); return mpfr_get_d(r19784, MPFR_RNDN); } static mpfr_t r19785, r19786, r19787, r19788, r19789, r19790, r19791, r19792, r19793, r19794, r19795, r19796, r19797, r19798, r19799, r19800, r19801, r19802, r19803; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19785); mpfr_init(r19786); mpfr_init(r19787); mpfr_init(r19788); mpfr_init(r19789); mpfr_init(r19790); mpfr_init(r19791); mpfr_init(r19792); mpfr_init_set_str(r19793, "0.0018067138f0", 10, MPFR_RNDN); mpfr_init(r19794); mpfr_init(r19795); mpfr_init(r19796); mpfr_init_set_str(r19797, "1", 10, MPFR_RNDN); mpfr_init(r19798); mpfr_init(r19799); mpfr_init(r19800); mpfr_init(r19801); mpfr_init(r19802); mpfr_init(r19803); } double f_fm(double wj, double x) { mpfr_set_d(r19785, wj, MPFR_RNDN); mpfr_exp(r19786, r19785, MPFR_RNDN); mpfr_mul(r19787, r19785, r19786, MPFR_RNDN); mpfr_set_d(r19788, x, MPFR_RNDN); mpfr_sub(r19789, r19787, r19788, MPFR_RNDN); mpfr_add(r19790, r19786, r19787, MPFR_RNDN); mpfr_div(r19791, r19789, r19790, MPFR_RNDN); mpfr_sub(r19792, r19785, r19791, MPFR_RNDN); ; mpfr_set_si(r19794, mpfr_cmp(r19792, r19793) <= 0, MPFR_RNDN); mpfr_cbrt(r19795, r19792, MPFR_RNDN); mpfr_mul(r19796, r19795, r19795, MPFR_RNDN); mpfr_mul(r19796, r19796, r19795, MPFR_RNDN); ; mpfr_add(r19798, r19797, r19785, MPFR_RNDN); mpfr_div(r19799, r19785, r19798, MPFR_RNDN); mpfr_sub(r19800, r19785, r19799, MPFR_RNDN); mpfr_div(r19801, r19788, r19790, MPFR_RNDN); mpfr_add(r19802, r19800, r19801, MPFR_RNDN); if (mpfr_get_si(r19794, MPFR_RNDN)) { mpfr_set(r19803, r19796, MPFR_RNDN); } else { mpfr_set(r19803, r19802, MPFR_RNDN); }; return mpfr_get_d(r19803, MPFR_RNDN); } static mpfr_t r19804, r19805, r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19804); mpfr_init(r19805); mpfr_init(r19806); mpfr_init(r19807); mpfr_init(r19808); mpfr_init(r19809); mpfr_init(r19810); mpfr_init(r19811); mpfr_init_set_str(r19812, "0.0018067138f0", 10, MPFR_RNDN); mpfr_init(r19813); mpfr_init(r19814); mpfr_init(r19815); mpfr_init_set_str(r19816, "1", 10, MPFR_RNDN); mpfr_init(r19817); mpfr_init(r19818); mpfr_init(r19819); mpfr_init(r19820); mpfr_init(r19821); mpfr_init(r19822); } double f_dm(double wj, double x) { mpfr_set_d(r19804, wj, MPFR_RNDN); mpfr_exp(r19805, r19804, MPFR_RNDN); mpfr_mul(r19806, r19804, r19805, MPFR_RNDN); mpfr_set_d(r19807, x, MPFR_RNDN); mpfr_sub(r19808, r19806, r19807, MPFR_RNDN); mpfr_add(r19809, r19805, r19806, MPFR_RNDN); mpfr_div(r19810, r19808, r19809, MPFR_RNDN); mpfr_sub(r19811, r19804, r19810, MPFR_RNDN); ; mpfr_set_si(r19813, mpfr_cmp(r19811, r19812) <= 0, MPFR_RNDN); mpfr_cbrt(r19814, r19811, MPFR_RNDN); mpfr_mul(r19815, r19814, r19814, MPFR_RNDN); mpfr_mul(r19815, r19815, r19814, MPFR_RNDN); ; mpfr_add(r19817, r19816, r19804, MPFR_RNDN); mpfr_div(r19818, r19804, r19817, MPFR_RNDN); mpfr_sub(r19819, r19804, r19818, MPFR_RNDN); mpfr_div(r19820, r19807, r19809, MPFR_RNDN); mpfr_add(r19821, r19819, r19820, MPFR_RNDN); if (mpfr_get_si(r19813, MPFR_RNDN)) { mpfr_set(r19822, r19815, MPFR_RNDN); } else { mpfr_set(r19822, r19821, MPFR_RNDN); }; return mpfr_get_d(r19822, MPFR_RNDN); }