#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19127 = wj; float r19128 = exp(r19127); float r19129 = r19127 * r19128; float r19130 = x; float r19131 = r19129 - r19130; float r19132 = r19128 + r19129; float r19133 = r19131 / r19132; float r19134 = r19127 - r19133; return r19134; } double f_id(double wj, double x) { double r19135 = wj; double r19136 = exp(r19135); double r19137 = r19135 * r19136; double r19138 = x; double r19139 = r19137 - r19138; double r19140 = r19136 + r19137; double r19141 = r19139 / r19140; double r19142 = r19135 - r19141; return r19142; } double f_of(float wj, float x) { float r19143 = wj; float r19144 = 0.4936105012893677f; bool r19145 = r19143 <= r19144; float r19146 = exp(r19143); float r19147 = r19143 * r19146; float r19148 = x; float r19149 = r19147 - r19148; float r19150 = r19146 + r19147; float r19151 = r19149 / r19150; float r19152 = r19143 - r19151; float r19153 = 1.0f; float r19154 = r19153 + r19143; float r19155 = r19143 / r19154; float r19156 = r19143 - r19155; float r19157 = r19148 / r19150; float r19158 = r19156 + r19157; float r19159 = r19145 ? r19152 : r19158; return r19159; } double f_od(double wj, double x) { double r19160 = wj; double r19161 = 0.4936105012893677; bool r19162 = r19160 <= r19161; double r19163 = exp(r19160); double r19164 = r19160 * r19163; double r19165 = x; double r19166 = r19164 - r19165; double r19167 = r19163 + r19164; double r19168 = r19166 / r19167; double r19169 = r19160 - r19168; double r19170 = 1.0; double r19171 = r19170 + r19160; double r19172 = r19160 / r19171; double r19173 = r19160 - r19172; double r19174 = r19165 / r19167; double r19175 = r19173 + r19174; double r19176 = r19162 ? r19169 : r19175; return r19176; } 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 r19177, r19178, r19179, r19180, r19181, r19182, r19183, r19184; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19177); mpfr_init(r19178); mpfr_init(r19179); mpfr_init(r19180); mpfr_init(r19181); mpfr_init(r19182); mpfr_init(r19183); mpfr_init(r19184); } double f_im(double wj, double x) { mpfr_set_d(r19177, wj, MPFR_RNDN); mpfr_exp(r19178, r19177, MPFR_RNDN); mpfr_mul(r19179, r19177, r19178, MPFR_RNDN); mpfr_set_d(r19180, x, MPFR_RNDN); mpfr_sub(r19181, r19179, r19180, MPFR_RNDN); mpfr_add(r19182, r19178, r19179, MPFR_RNDN); mpfr_div(r19183, r19181, r19182, MPFR_RNDN); mpfr_sub(r19184, r19177, r19183, MPFR_RNDN); return mpfr_get_d(r19184, MPFR_RNDN); } static mpfr_t r19185, r19186, r19187, r19188, r19189, r19190, r19191, r19192, r19193, r19194, r19195, r19196, r19197, r19198, r19199, r19200, r19201; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19185); mpfr_init_set_str(r19186, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19187); mpfr_init(r19188); mpfr_init(r19189); mpfr_init(r19190); mpfr_init(r19191); mpfr_init(r19192); mpfr_init(r19193); mpfr_init(r19194); mpfr_init_set_str(r19195, "1", 10, MPFR_RNDN); mpfr_init(r19196); mpfr_init(r19197); mpfr_init(r19198); mpfr_init(r19199); mpfr_init(r19200); mpfr_init(r19201); } double f_fm(double wj, double x) { mpfr_set_d(r19185, wj, MPFR_RNDN); ; mpfr_set_si(r19187, mpfr_cmp(r19185, r19186) <= 0, MPFR_RNDN); mpfr_exp(r19188, r19185, MPFR_RNDN); mpfr_mul(r19189, r19185, r19188, MPFR_RNDN); mpfr_set_d(r19190, x, MPFR_RNDN); mpfr_sub(r19191, r19189, r19190, MPFR_RNDN); mpfr_add(r19192, r19188, r19189, MPFR_RNDN); mpfr_div(r19193, r19191, r19192, MPFR_RNDN); mpfr_sub(r19194, r19185, r19193, MPFR_RNDN); ; mpfr_add(r19196, r19195, r19185, MPFR_RNDN); mpfr_div(r19197, r19185, r19196, MPFR_RNDN); mpfr_sub(r19198, r19185, r19197, MPFR_RNDN); mpfr_div(r19199, r19190, r19192, MPFR_RNDN); mpfr_add(r19200, r19198, r19199, MPFR_RNDN); if (mpfr_get_si(r19187, MPFR_RNDN)) { mpfr_set(r19201, r19194, MPFR_RNDN); } else { mpfr_set(r19201, r19200, MPFR_RNDN); }; return mpfr_get_d(r19201, MPFR_RNDN); } static mpfr_t r19202, r19203, r19204, r19205, r19206, r19207, r19208, r19209, r19210, r19211, r19212, r19213, r19214, r19215, r19216, r19217, r19218; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19202); mpfr_init_set_str(r19203, "0.4936105f0", 10, MPFR_RNDN); mpfr_init(r19204); mpfr_init(r19205); mpfr_init(r19206); mpfr_init(r19207); mpfr_init(r19208); mpfr_init(r19209); mpfr_init(r19210); mpfr_init(r19211); mpfr_init_set_str(r19212, "1", 10, MPFR_RNDN); mpfr_init(r19213); mpfr_init(r19214); mpfr_init(r19215); mpfr_init(r19216); mpfr_init(r19217); mpfr_init(r19218); } double f_dm(double wj, double x) { mpfr_set_d(r19202, wj, MPFR_RNDN); ; mpfr_set_si(r19204, mpfr_cmp(r19202, r19203) <= 0, MPFR_RNDN); mpfr_exp(r19205, r19202, MPFR_RNDN); mpfr_mul(r19206, r19202, r19205, MPFR_RNDN); mpfr_set_d(r19207, x, MPFR_RNDN); mpfr_sub(r19208, r19206, r19207, MPFR_RNDN); mpfr_add(r19209, r19205, r19206, MPFR_RNDN); mpfr_div(r19210, r19208, r19209, MPFR_RNDN); mpfr_sub(r19211, r19202, r19210, MPFR_RNDN); ; mpfr_add(r19213, r19212, r19202, MPFR_RNDN); mpfr_div(r19214, r19202, r19213, MPFR_RNDN); mpfr_sub(r19215, r19202, r19214, MPFR_RNDN); mpfr_div(r19216, r19207, r19209, MPFR_RNDN); mpfr_add(r19217, r19215, r19216, MPFR_RNDN); if (mpfr_get_si(r19204, MPFR_RNDN)) { mpfr_set(r19218, r19211, MPFR_RNDN); } else { mpfr_set(r19218, r19217, MPFR_RNDN); }; return mpfr_get_d(r19218, MPFR_RNDN); }