#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19143 = wj; float r19144 = exp(r19143); float r19145 = r19143 * r19144; float r19146 = x; float r19147 = r19145 - r19146; float r19148 = r19144 + r19145; float r19149 = r19147 / r19148; float r19150 = r19143 - r19149; return r19150; } double f_id(double wj, double x) { double r19151 = wj; double r19152 = exp(r19151); double r19153 = r19151 * r19152; double r19154 = x; double r19155 = r19153 - r19154; double r19156 = r19152 + r19153; double r19157 = r19155 / r19156; double r19158 = r19151 - r19157; return r19158; } double f_of(float wj, float x) { float r19159 = wj; float r19160 = 0.9588147401809692f; bool r19161 = r19159 <= r19160; float r19162 = exp(r19159); float r19163 = r19159 * r19162; float r19164 = x; float r19165 = r19163 - r19164; float r19166 = r19162 + r19163; float r19167 = r19165 / r19166; float r19168 = r19159 - r19167; float r19169 = 1.0f; float r19170 = r19159 + r19169; float r19171 = r19159 / r19170; float r19172 = r19169 + r19159; float r19173 = r19164 / r19172; float r19174 = r19173 / r19162; float r19175 = r19171 - r19174; float r19176 = r19159 - r19175; float r19177 = r19161 ? r19168 : r19176; return r19177; } double f_od(double wj, double x) { double r19178 = wj; double r19179 = 0.9588147401809692; bool r19180 = r19178 <= r19179; double r19181 = exp(r19178); double r19182 = r19178 * r19181; double r19183 = x; double r19184 = r19182 - r19183; double r19185 = r19181 + r19182; double r19186 = r19184 / r19185; double r19187 = r19178 - r19186; double r19188 = 1.0; double r19189 = r19178 + r19188; double r19190 = r19178 / r19189; double r19191 = r19188 + r19178; double r19192 = r19183 / r19191; double r19193 = r19192 / r19181; double r19194 = r19190 - r19193; double r19195 = r19178 - r19194; double r19196 = r19180 ? r19187 : r19195; return r19196; } 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 r19197, r19198, r19199, r19200, r19201, r19202, r19203, r19204; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19197); mpfr_init(r19198); mpfr_init(r19199); mpfr_init(r19200); mpfr_init(r19201); mpfr_init(r19202); mpfr_init(r19203); mpfr_init(r19204); } double f_im(double wj, double x) { mpfr_set_d(r19197, wj, MPFR_RNDN); mpfr_exp(r19198, r19197, MPFR_RNDN); mpfr_mul(r19199, r19197, r19198, MPFR_RNDN); mpfr_set_d(r19200, x, MPFR_RNDN); mpfr_sub(r19201, r19199, r19200, MPFR_RNDN); mpfr_add(r19202, r19198, r19199, MPFR_RNDN); mpfr_div(r19203, r19201, r19202, MPFR_RNDN); mpfr_sub(r19204, r19197, r19203, MPFR_RNDN); return mpfr_get_d(r19204, MPFR_RNDN); } static mpfr_t r19205, r19206, r19207, r19208, r19209, r19210, r19211, r19212, r19213, r19214, r19215, r19216, r19217, r19218, r19219, r19220, r19221, r19222, r19223; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19205); mpfr_init_set_str(r19206, "0.95881474f0", 10, MPFR_RNDN); mpfr_init(r19207); mpfr_init(r19208); mpfr_init(r19209); mpfr_init(r19210); mpfr_init(r19211); mpfr_init(r19212); mpfr_init(r19213); mpfr_init(r19214); mpfr_init_set_str(r19215, "1", 10, MPFR_RNDN); mpfr_init(r19216); mpfr_init(r19217); mpfr_init(r19218); mpfr_init(r19219); mpfr_init(r19220); mpfr_init(r19221); mpfr_init(r19222); mpfr_init(r19223); } double f_fm(double wj, double x) { mpfr_set_d(r19205, wj, MPFR_RNDN); ; mpfr_set_si(r19207, mpfr_cmp(r19205, r19206) <= 0, MPFR_RNDN); mpfr_exp(r19208, r19205, MPFR_RNDN); mpfr_mul(r19209, r19205, r19208, MPFR_RNDN); mpfr_set_d(r19210, x, MPFR_RNDN); mpfr_sub(r19211, r19209, r19210, MPFR_RNDN); mpfr_add(r19212, r19208, r19209, MPFR_RNDN); mpfr_div(r19213, r19211, r19212, MPFR_RNDN); mpfr_sub(r19214, r19205, r19213, MPFR_RNDN); ; mpfr_add(r19216, r19205, r19215, MPFR_RNDN); mpfr_div(r19217, r19205, r19216, MPFR_RNDN); mpfr_add(r19218, r19215, r19205, MPFR_RNDN); mpfr_div(r19219, r19210, r19218, MPFR_RNDN); mpfr_div(r19220, r19219, r19208, MPFR_RNDN); mpfr_sub(r19221, r19217, r19220, MPFR_RNDN); mpfr_sub(r19222, r19205, r19221, MPFR_RNDN); if (mpfr_get_si(r19207, MPFR_RNDN)) { mpfr_set(r19223, r19214, MPFR_RNDN); } else { mpfr_set(r19223, r19222, MPFR_RNDN); }; return mpfr_get_d(r19223, MPFR_RNDN); } static mpfr_t r19224, r19225, r19226, r19227, r19228, r19229, r19230, r19231, r19232, r19233, r19234, r19235, r19236, r19237, r19238, r19239, r19240, r19241, r19242; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19224); mpfr_init_set_str(r19225, "0.95881474f0", 10, MPFR_RNDN); mpfr_init(r19226); mpfr_init(r19227); mpfr_init(r19228); mpfr_init(r19229); mpfr_init(r19230); mpfr_init(r19231); mpfr_init(r19232); mpfr_init(r19233); mpfr_init_set_str(r19234, "1", 10, MPFR_RNDN); mpfr_init(r19235); mpfr_init(r19236); mpfr_init(r19237); mpfr_init(r19238); mpfr_init(r19239); mpfr_init(r19240); mpfr_init(r19241); mpfr_init(r19242); } double f_dm(double wj, double x) { mpfr_set_d(r19224, wj, MPFR_RNDN); ; mpfr_set_si(r19226, mpfr_cmp(r19224, r19225) <= 0, MPFR_RNDN); mpfr_exp(r19227, r19224, MPFR_RNDN); mpfr_mul(r19228, r19224, r19227, MPFR_RNDN); mpfr_set_d(r19229, x, MPFR_RNDN); mpfr_sub(r19230, r19228, r19229, MPFR_RNDN); mpfr_add(r19231, r19227, r19228, MPFR_RNDN); mpfr_div(r19232, r19230, r19231, MPFR_RNDN); mpfr_sub(r19233, r19224, r19232, MPFR_RNDN); ; mpfr_add(r19235, r19224, r19234, MPFR_RNDN); mpfr_div(r19236, r19224, r19235, MPFR_RNDN); mpfr_add(r19237, r19234, r19224, MPFR_RNDN); mpfr_div(r19238, r19229, r19237, MPFR_RNDN); mpfr_div(r19239, r19238, r19227, MPFR_RNDN); mpfr_sub(r19240, r19236, r19239, MPFR_RNDN); mpfr_sub(r19241, r19224, r19240, MPFR_RNDN); if (mpfr_get_si(r19226, MPFR_RNDN)) { mpfr_set(r19242, r19233, MPFR_RNDN); } else { mpfr_set(r19242, r19241, MPFR_RNDN); }; return mpfr_get_d(r19242, MPFR_RNDN); }