#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19039 = wj; float r19040 = exp(r19039); float r19041 = r19039 * r19040; float r19042 = x; float r19043 = r19041 - r19042; float r19044 = r19040 + r19041; float r19045 = r19043 / r19044; float r19046 = r19039 - r19045; return r19046; } double f_id(double wj, double x) { double r19047 = wj; double r19048 = exp(r19047); double r19049 = r19047 * r19048; double r19050 = x; double r19051 = r19049 - r19050; double r19052 = r19048 + r19049; double r19053 = r19051 / r19052; double r19054 = r19047 - r19053; return r19054; } double f_of(float wj, float x) { float r19055 = wj; float r19056 = 11.522058486938477f; bool r19057 = r19055 <= r19056; float r19058 = exp(r19055); float r19059 = r19055 * r19058; float r19060 = x; float r19061 = r19059 - r19060; float r19062 = r19058 + r19059; float r19063 = r19061 / r19062; float r19064 = r19055 - r19063; float r19065 = 1.0f; float r19066 = r19055 + r19065; float r19067 = r19055 / r19066; float r19068 = fma(r19055, r19058, r19058); float r19069 = r19060 / r19068; float r19070 = r19067 - r19069; float r19071 = r19055 - r19070; float r19072 = r19057 ? r19064 : r19071; return r19072; } double f_od(double wj, double x) { double r19073 = wj; double r19074 = 11.522058486938477; bool r19075 = r19073 <= r19074; double r19076 = exp(r19073); double r19077 = r19073 * r19076; double r19078 = x; double r19079 = r19077 - r19078; double r19080 = r19076 + r19077; double r19081 = r19079 / r19080; double r19082 = r19073 - r19081; double r19083 = 1.0; double r19084 = r19073 + r19083; double r19085 = r19073 / r19084; double r19086 = fma(r19073, r19076, r19076); double r19087 = r19078 / r19086; double r19088 = r19085 - r19087; double r19089 = r19073 - r19088; double r19090 = r19075 ? r19082 : r19089; return r19090; } 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 r19091, r19092, r19093, r19094, r19095, r19096, r19097, r19098; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19091); mpfr_init(r19092); mpfr_init(r19093); mpfr_init(r19094); mpfr_init(r19095); mpfr_init(r19096); mpfr_init(r19097); mpfr_init(r19098); } double f_im(double wj, double x) { mpfr_set_d(r19091, wj, MPFR_RNDN); mpfr_exp(r19092, r19091, MPFR_RNDN); mpfr_mul(r19093, r19091, r19092, MPFR_RNDN); mpfr_set_d(r19094, x, MPFR_RNDN); mpfr_sub(r19095, r19093, r19094, MPFR_RNDN); mpfr_add(r19096, r19092, r19093, MPFR_RNDN); mpfr_div(r19097, r19095, r19096, MPFR_RNDN); mpfr_sub(r19098, r19091, r19097, MPFR_RNDN); return mpfr_get_d(r19098, MPFR_RNDN); } static mpfr_t r19099, r19100, r19101, r19102, r19103, r19104, r19105, r19106, r19107, r19108, r19109, r19110, r19111, r19112, r19113, r19114, r19115, r19116; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19099); mpfr_init_set_str(r19100, "11.5220585f0", 10, MPFR_RNDN); mpfr_init(r19101); mpfr_init(r19102); mpfr_init(r19103); mpfr_init(r19104); mpfr_init(r19105); mpfr_init(r19106); mpfr_init(r19107); mpfr_init(r19108); mpfr_init_set_str(r19109, "1", 10, MPFR_RNDN); mpfr_init(r19110); mpfr_init(r19111); mpfr_init(r19112); mpfr_init(r19113); mpfr_init(r19114); mpfr_init(r19115); mpfr_init(r19116); } double f_fm(double wj, double x) { mpfr_set_d(r19099, wj, MPFR_RNDN); ; mpfr_set_si(r19101, mpfr_cmp(r19099, r19100) <= 0, MPFR_RNDN); mpfr_exp(r19102, r19099, MPFR_RNDN); mpfr_mul(r19103, r19099, r19102, MPFR_RNDN); mpfr_set_d(r19104, x, MPFR_RNDN); mpfr_sub(r19105, r19103, r19104, MPFR_RNDN); mpfr_add(r19106, r19102, r19103, MPFR_RNDN); mpfr_div(r19107, r19105, r19106, MPFR_RNDN); mpfr_sub(r19108, r19099, r19107, MPFR_RNDN); ; mpfr_add(r19110, r19099, r19109, MPFR_RNDN); mpfr_div(r19111, r19099, r19110, MPFR_RNDN); mpfr_fma(r19112, r19099, r19102, r19102, MPFR_RNDN); mpfr_div(r19113, r19104, r19112, MPFR_RNDN); mpfr_sub(r19114, r19111, r19113, MPFR_RNDN); mpfr_sub(r19115, r19099, r19114, MPFR_RNDN); if (mpfr_get_si(r19101, MPFR_RNDN)) { mpfr_set(r19116, r19108, MPFR_RNDN); } else { mpfr_set(r19116, r19115, MPFR_RNDN); }; return mpfr_get_d(r19116, MPFR_RNDN); } static mpfr_t r19117, r19118, r19119, r19120, r19121, r19122, r19123, r19124, r19125, r19126, r19127, r19128, r19129, r19130, r19131, r19132, r19133, r19134; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19117); mpfr_init_set_str(r19118, "11.5220585f0", 10, MPFR_RNDN); mpfr_init(r19119); mpfr_init(r19120); mpfr_init(r19121); mpfr_init(r19122); mpfr_init(r19123); mpfr_init(r19124); mpfr_init(r19125); mpfr_init(r19126); mpfr_init_set_str(r19127, "1", 10, MPFR_RNDN); mpfr_init(r19128); mpfr_init(r19129); mpfr_init(r19130); mpfr_init(r19131); mpfr_init(r19132); mpfr_init(r19133); mpfr_init(r19134); } double f_dm(double wj, double x) { mpfr_set_d(r19117, wj, MPFR_RNDN); ; mpfr_set_si(r19119, mpfr_cmp(r19117, r19118) <= 0, MPFR_RNDN); mpfr_exp(r19120, r19117, MPFR_RNDN); mpfr_mul(r19121, r19117, r19120, MPFR_RNDN); mpfr_set_d(r19122, x, MPFR_RNDN); mpfr_sub(r19123, r19121, r19122, MPFR_RNDN); mpfr_add(r19124, r19120, r19121, MPFR_RNDN); mpfr_div(r19125, r19123, r19124, MPFR_RNDN); mpfr_sub(r19126, r19117, r19125, MPFR_RNDN); ; mpfr_add(r19128, r19117, r19127, MPFR_RNDN); mpfr_div(r19129, r19117, r19128, MPFR_RNDN); mpfr_fma(r19130, r19117, r19120, r19120, MPFR_RNDN); mpfr_div(r19131, r19122, r19130, MPFR_RNDN); mpfr_sub(r19132, r19129, r19131, MPFR_RNDN); mpfr_sub(r19133, r19117, r19132, MPFR_RNDN); if (mpfr_get_si(r19119, MPFR_RNDN)) { mpfr_set(r19134, r19126, MPFR_RNDN); } else { mpfr_set(r19134, r19133, MPFR_RNDN); }; return mpfr_get_d(r19134, MPFR_RNDN); }