#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19319 = wj; float r19320 = exp(r19319); float r19321 = r19319 * r19320; float r19322 = x; float r19323 = r19321 - r19322; float r19324 = r19320 + r19321; float r19325 = r19323 / r19324; float r19326 = r19319 - r19325; return r19326; } double f_id(double wj, double x) { double r19327 = wj; double r19328 = exp(r19327); double r19329 = r19327 * r19328; double r19330 = x; double r19331 = r19329 - r19330; double r19332 = r19328 + r19329; double r19333 = r19331 / r19332; double r19334 = r19327 - r19333; return r19334; } double f_of(float wj, float x) { float r19335 = wj; float r19336 = exp(r19335); float r19337 = r19335 * r19336; float r19338 = x; float r19339 = r19337 - r19338; float r19340 = r19336 + r19337; float r19341 = r19339 / r19340; float r19342 = r19335 - r19341; float r19343 = 4.477717099564361e-09f; bool r19344 = r19342 <= r19343; float r19345 = r19335 * r19335; float r19346 = r19335 - r19345; float r19347 = 4.0f; float r19348 = pow(r19335, r19347); float r19349 = fma(r19335, r19346, r19348); float r19350 = fma(r19335, r19336, r19336); float r19351 = r19338 / r19350; float r19352 = r19349 + r19351; float r19353 = 1.0f; float r19354 = r19353 + r19335; float r19355 = r19335 / r19354; float r19356 = r19335 - r19355; float r19357 = r19338 / r19340; float r19358 = r19356 + r19357; float r19359 = r19344 ? r19352 : r19358; return r19359; } double f_od(double wj, double x) { double r19360 = wj; double r19361 = exp(r19360); double r19362 = r19360 * r19361; double r19363 = x; double r19364 = r19362 - r19363; double r19365 = r19361 + r19362; double r19366 = r19364 / r19365; double r19367 = r19360 - r19366; double r19368 = 4.477717099564361e-09; bool r19369 = r19367 <= r19368; double r19370 = r19360 * r19360; double r19371 = r19360 - r19370; double r19372 = 4.0; double r19373 = pow(r19360, r19372); double r19374 = fma(r19360, r19371, r19373); double r19375 = fma(r19360, r19361, r19361); double r19376 = r19363 / r19375; double r19377 = r19374 + r19376; double r19378 = 1.0; double r19379 = r19378 + r19360; double r19380 = r19360 / r19379; double r19381 = r19360 - r19380; double r19382 = r19363 / r19365; double r19383 = r19381 + r19382; double r19384 = r19369 ? r19377 : r19383; return r19384; } 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 r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19385); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); } double f_im(double wj, double x) { mpfr_set_d(r19385, wj, MPFR_RNDN); mpfr_exp(r19386, r19385, MPFR_RNDN); mpfr_mul(r19387, r19385, r19386, MPFR_RNDN); mpfr_set_d(r19388, x, MPFR_RNDN); mpfr_sub(r19389, r19387, r19388, MPFR_RNDN); mpfr_add(r19390, r19386, r19387, MPFR_RNDN); mpfr_div(r19391, r19389, r19390, MPFR_RNDN); mpfr_sub(r19392, r19385, r19391, MPFR_RNDN); return mpfr_get_d(r19392, MPFR_RNDN); } static mpfr_t r19393, r19394, r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414, r19415, r19416, r19417; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19393); mpfr_init(r19394); mpfr_init(r19395); mpfr_init(r19396); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init_set_str(r19401, "4.477717099564361e-09", 10, MPFR_RNDN); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init_set_str(r19405, "4", 10, MPFR_RNDN); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init_set_str(r19411, "1", 10, MPFR_RNDN); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); mpfr_init(r19415); mpfr_init(r19416); mpfr_init(r19417); } double f_fm(double wj, double x) { mpfr_set_d(r19393, wj, MPFR_RNDN); mpfr_exp(r19394, r19393, MPFR_RNDN); mpfr_mul(r19395, r19393, r19394, MPFR_RNDN); mpfr_set_d(r19396, x, MPFR_RNDN); mpfr_sub(r19397, r19395, r19396, MPFR_RNDN); mpfr_add(r19398, r19394, r19395, MPFR_RNDN); mpfr_div(r19399, r19397, r19398, MPFR_RNDN); mpfr_sub(r19400, r19393, r19399, MPFR_RNDN); ; mpfr_set_si(r19402, mpfr_cmp(r19400, r19401) <= 0, MPFR_RNDN); mpfr_sqr(r19403, r19393, MPFR_RNDN); mpfr_sub(r19404, r19393, r19403, MPFR_RNDN); ; mpfr_pow(r19406, r19393, r19405, MPFR_RNDN); mpfr_fma(r19407, r19393, r19404, r19406, MPFR_RNDN); mpfr_fma(r19408, r19393, r19394, r19394, MPFR_RNDN); mpfr_div(r19409, r19396, r19408, MPFR_RNDN); mpfr_add(r19410, r19407, r19409, MPFR_RNDN); ; mpfr_add(r19412, r19411, r19393, MPFR_RNDN); mpfr_div(r19413, r19393, r19412, MPFR_RNDN); mpfr_sub(r19414, r19393, r19413, MPFR_RNDN); mpfr_div(r19415, r19396, r19398, MPFR_RNDN); mpfr_add(r19416, r19414, r19415, MPFR_RNDN); if (mpfr_get_si(r19402, MPFR_RNDN)) { mpfr_set(r19417, r19410, MPFR_RNDN); } else { mpfr_set(r19417, r19416, MPFR_RNDN); }; return mpfr_get_d(r19417, MPFR_RNDN); } static mpfr_t r19418, r19419, r19420, r19421, r19422, r19423, r19424, r19425, r19426, r19427, r19428, r19429, r19430, r19431, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440, r19441, r19442; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19418); mpfr_init(r19419); mpfr_init(r19420); mpfr_init(r19421); mpfr_init(r19422); mpfr_init(r19423); mpfr_init(r19424); mpfr_init(r19425); mpfr_init_set_str(r19426, "4.477717099564361e-09", 10, MPFR_RNDN); mpfr_init(r19427); mpfr_init(r19428); mpfr_init(r19429); mpfr_init_set_str(r19430, "4", 10, MPFR_RNDN); mpfr_init(r19431); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); mpfr_init(r19435); mpfr_init_set_str(r19436, "1", 10, MPFR_RNDN); mpfr_init(r19437); mpfr_init(r19438); mpfr_init(r19439); mpfr_init(r19440); mpfr_init(r19441); mpfr_init(r19442); } double f_dm(double wj, double x) { mpfr_set_d(r19418, wj, MPFR_RNDN); mpfr_exp(r19419, r19418, MPFR_RNDN); mpfr_mul(r19420, r19418, r19419, MPFR_RNDN); mpfr_set_d(r19421, x, MPFR_RNDN); mpfr_sub(r19422, r19420, r19421, MPFR_RNDN); mpfr_add(r19423, r19419, r19420, MPFR_RNDN); mpfr_div(r19424, r19422, r19423, MPFR_RNDN); mpfr_sub(r19425, r19418, r19424, MPFR_RNDN); ; mpfr_set_si(r19427, mpfr_cmp(r19425, r19426) <= 0, MPFR_RNDN); mpfr_sqr(r19428, r19418, MPFR_RNDN); mpfr_sub(r19429, r19418, r19428, MPFR_RNDN); ; mpfr_pow(r19431, r19418, r19430, MPFR_RNDN); mpfr_fma(r19432, r19418, r19429, r19431, MPFR_RNDN); mpfr_fma(r19433, r19418, r19419, r19419, MPFR_RNDN); mpfr_div(r19434, r19421, r19433, MPFR_RNDN); mpfr_add(r19435, r19432, r19434, MPFR_RNDN); ; mpfr_add(r19437, r19436, r19418, MPFR_RNDN); mpfr_div(r19438, r19418, r19437, MPFR_RNDN); mpfr_sub(r19439, r19418, r19438, MPFR_RNDN); mpfr_div(r19440, r19421, r19423, MPFR_RNDN); mpfr_add(r19441, r19439, r19440, MPFR_RNDN); if (mpfr_get_si(r19427, MPFR_RNDN)) { mpfr_set(r19442, r19435, MPFR_RNDN); } else { mpfr_set(r19442, r19441, MPFR_RNDN); }; return mpfr_get_d(r19442, MPFR_RNDN); }