#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19559 = wj; float r19560 = exp(r19559); float r19561 = r19559 * r19560; float r19562 = x; float r19563 = r19561 - r19562; float r19564 = r19560 + r19561; float r19565 = r19563 / r19564; float r19566 = r19559 - r19565; return r19566; } double f_id(double wj, double x) { double r19567 = wj; double r19568 = exp(r19567); double r19569 = r19567 * r19568; double r19570 = x; double r19571 = r19569 - r19570; double r19572 = r19568 + r19569; double r19573 = r19571 / r19572; double r19574 = r19567 - r19573; return r19574; } double f_of(float wj, float x) { float r19575 = wj; float r19576 = 0.03272881358861923f; bool r19577 = r19575 <= r19576; float r19578 = exp(r19575); float r19579 = r19575 * r19578; float r19580 = x; float r19581 = r19579 - r19580; float r19582 = r19578 + r19579; float r19583 = r19581 / r19582; float r19584 = r19575 - r19583; float r19585 = 1.0f; float r19586 = r19585 + r19575; float r19587 = r19575 / r19586; float r19588 = r19575 - r19587; float r19589 = pow(r19582, r19585); float r19590 = r19580 / r19589; float r19591 = r19588 + r19590; float r19592 = r19577 ? r19584 : r19591; return r19592; } double f_od(double wj, double x) { double r19593 = wj; double r19594 = 0.03272881358861923; bool r19595 = r19593 <= r19594; double r19596 = exp(r19593); double r19597 = r19593 * r19596; double r19598 = x; double r19599 = r19597 - r19598; double r19600 = r19596 + r19597; double r19601 = r19599 / r19600; double r19602 = r19593 - r19601; double r19603 = 1.0; double r19604 = r19603 + r19593; double r19605 = r19593 / r19604; double r19606 = r19593 - r19605; double r19607 = pow(r19600, r19603); double r19608 = r19598 / r19607; double r19609 = r19606 + r19608; double r19610 = r19595 ? r19602 : r19609; return r19610; } 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 r19611, r19612, r19613, r19614, r19615, r19616, r19617, r19618; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19611); mpfr_init(r19612); mpfr_init(r19613); mpfr_init(r19614); mpfr_init(r19615); mpfr_init(r19616); mpfr_init(r19617); mpfr_init(r19618); } double f_im(double wj, double x) { mpfr_set_d(r19611, wj, MPFR_RNDN); mpfr_exp(r19612, r19611, MPFR_RNDN); mpfr_mul(r19613, r19611, r19612, MPFR_RNDN); mpfr_set_d(r19614, x, MPFR_RNDN); mpfr_sub(r19615, r19613, r19614, MPFR_RNDN); mpfr_add(r19616, r19612, r19613, MPFR_RNDN); mpfr_div(r19617, r19615, r19616, MPFR_RNDN); mpfr_sub(r19618, r19611, r19617, MPFR_RNDN); return mpfr_get_d(r19618, MPFR_RNDN); } static mpfr_t r19619, r19620, r19621, r19622, r19623, r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19619); mpfr_init_set_str(r19620, "0.032728814f0", 10, MPFR_RNDN); mpfr_init(r19621); mpfr_init(r19622); mpfr_init(r19623); mpfr_init(r19624); mpfr_init(r19625); mpfr_init(r19626); mpfr_init(r19627); mpfr_init(r19628); mpfr_init_set_str(r19629, "1", 10, MPFR_RNDN); mpfr_init(r19630); mpfr_init(r19631); mpfr_init(r19632); mpfr_init(r19633); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); } double f_fm(double wj, double x) { mpfr_set_d(r19619, wj, MPFR_RNDN); ; mpfr_set_si(r19621, mpfr_cmp(r19619, r19620) <= 0, MPFR_RNDN); mpfr_exp(r19622, r19619, MPFR_RNDN); mpfr_mul(r19623, r19619, r19622, MPFR_RNDN); mpfr_set_d(r19624, x, MPFR_RNDN); mpfr_sub(r19625, r19623, r19624, MPFR_RNDN); mpfr_add(r19626, r19622, r19623, MPFR_RNDN); mpfr_div(r19627, r19625, r19626, MPFR_RNDN); mpfr_sub(r19628, r19619, r19627, MPFR_RNDN); ; mpfr_add(r19630, r19629, r19619, MPFR_RNDN); mpfr_div(r19631, r19619, r19630, MPFR_RNDN); mpfr_sub(r19632, r19619, r19631, MPFR_RNDN); mpfr_pow(r19633, r19626, r19629, MPFR_RNDN); mpfr_div(r19634, r19624, r19633, MPFR_RNDN); mpfr_add(r19635, r19632, r19634, MPFR_RNDN); if (mpfr_get_si(r19621, MPFR_RNDN)) { mpfr_set(r19636, r19628, MPFR_RNDN); } else { mpfr_set(r19636, r19635, MPFR_RNDN); }; return mpfr_get_d(r19636, MPFR_RNDN); } static mpfr_t r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650, r19651, r19652, r19653, r19654; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19637); mpfr_init_set_str(r19638, "0.032728814f0", 10, MPFR_RNDN); mpfr_init(r19639); mpfr_init(r19640); mpfr_init(r19641); mpfr_init(r19642); mpfr_init(r19643); mpfr_init(r19644); mpfr_init(r19645); mpfr_init(r19646); mpfr_init_set_str(r19647, "1", 10, MPFR_RNDN); mpfr_init(r19648); mpfr_init(r19649); mpfr_init(r19650); mpfr_init(r19651); mpfr_init(r19652); mpfr_init(r19653); mpfr_init(r19654); } double f_dm(double wj, double x) { mpfr_set_d(r19637, wj, MPFR_RNDN); ; mpfr_set_si(r19639, mpfr_cmp(r19637, r19638) <= 0, MPFR_RNDN); mpfr_exp(r19640, r19637, MPFR_RNDN); mpfr_mul(r19641, r19637, r19640, MPFR_RNDN); mpfr_set_d(r19642, x, MPFR_RNDN); mpfr_sub(r19643, r19641, r19642, MPFR_RNDN); mpfr_add(r19644, r19640, r19641, MPFR_RNDN); mpfr_div(r19645, r19643, r19644, MPFR_RNDN); mpfr_sub(r19646, r19637, r19645, MPFR_RNDN); ; mpfr_add(r19648, r19647, r19637, MPFR_RNDN); mpfr_div(r19649, r19637, r19648, MPFR_RNDN); mpfr_sub(r19650, r19637, r19649, MPFR_RNDN); mpfr_pow(r19651, r19644, r19647, MPFR_RNDN); mpfr_div(r19652, r19642, r19651, MPFR_RNDN); mpfr_add(r19653, r19650, r19652, MPFR_RNDN); if (mpfr_get_si(r19639, MPFR_RNDN)) { mpfr_set(r19654, r19646, MPFR_RNDN); } else { mpfr_set(r19654, r19653, MPFR_RNDN); }; return mpfr_get_d(r19654, MPFR_RNDN); }