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