#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19311 = wj; float r19312 = exp(r19311); float r19313 = r19311 * r19312; float r19314 = x; float r19315 = r19313 - r19314; float r19316 = r19312 + r19313; float r19317 = r19315 / r19316; float r19318 = r19311 - r19317; return r19318; } double f_id(double wj, double x) { double r19319 = wj; double r19320 = exp(r19319); double r19321 = r19319 * r19320; double r19322 = x; double r19323 = r19321 - r19322; double r19324 = r19320 + r19321; double r19325 = r19323 / r19324; double r19326 = r19319 - r19325; return r19326; } double f_of(float wj, float x) { float r19327 = wj; float r19328 = 7.279548168182373f; bool r19329 = r19327 <= r19328; float r19330 = exp(r19327); float r19331 = r19327 * r19330; float r19332 = x; float r19333 = r19331 - r19332; float r19334 = r19330 + r19331; float r19335 = r19333 / r19334; float r19336 = r19327 - r19335; float r19337 = 1.0f; float r19338 = r19337 + r19327; float r19339 = r19327 / r19338; float r19340 = r19327 - r19339; float r19341 = fma(r19327, r19330, r19330); float r19342 = r19332 / r19341; float r19343 = log1p(r19342); float r19344 = expm1(r19343); float r19345 = r19340 + r19344; float r19346 = r19329 ? r19336 : r19345; return r19346; } double f_od(double wj, double x) { double r19347 = wj; double r19348 = 7.279548168182373; bool r19349 = r19347 <= r19348; double r19350 = exp(r19347); double r19351 = r19347 * r19350; double r19352 = x; double r19353 = r19351 - r19352; double r19354 = r19350 + r19351; double r19355 = r19353 / r19354; double r19356 = r19347 - r19355; double r19357 = 1.0; double r19358 = r19357 + r19347; double r19359 = r19347 / r19358; double r19360 = r19347 - r19359; double r19361 = fma(r19347, r19350, r19350); double r19362 = r19352 / r19361; double r19363 = log1p(r19362); double r19364 = expm1(r19363); double r19365 = r19360 + r19364; double r19366 = r19349 ? r19356 : r19365; return r19366; } 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 r19367, r19368, r19369, r19370, r19371, r19372, r19373, r19374; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19367); mpfr_init(r19368); mpfr_init(r19369); mpfr_init(r19370); mpfr_init(r19371); mpfr_init(r19372); mpfr_init(r19373); mpfr_init(r19374); } double f_im(double wj, double x) { mpfr_set_d(r19367, wj, MPFR_RNDN); mpfr_exp(r19368, r19367, MPFR_RNDN); mpfr_mul(r19369, r19367, r19368, MPFR_RNDN); mpfr_set_d(r19370, x, MPFR_RNDN); mpfr_sub(r19371, r19369, r19370, MPFR_RNDN); mpfr_add(r19372, r19368, r19369, MPFR_RNDN); mpfr_div(r19373, r19371, r19372, MPFR_RNDN); mpfr_sub(r19374, r19367, r19373, MPFR_RNDN); return mpfr_get_d(r19374, MPFR_RNDN); } static mpfr_t r19375, r19376, r19377, r19378, r19379, r19380, r19381, r19382, r19383, r19384, r19385, r19386, r19387, r19388, r19389, r19390, r19391, r19392, r19393, r19394; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19375); mpfr_init_set_str(r19376, "7.279548f0", 10, MPFR_RNDN); mpfr_init(r19377); mpfr_init(r19378); mpfr_init(r19379); mpfr_init(r19380); mpfr_init(r19381); mpfr_init(r19382); mpfr_init(r19383); mpfr_init(r19384); mpfr_init_set_str(r19385, "1", 10, MPFR_RNDN); mpfr_init(r19386); mpfr_init(r19387); mpfr_init(r19388); mpfr_init(r19389); mpfr_init(r19390); mpfr_init(r19391); mpfr_init(r19392); mpfr_init(r19393); mpfr_init(r19394); } double f_fm(double wj, double x) { mpfr_set_d(r19375, wj, MPFR_RNDN); ; mpfr_set_si(r19377, mpfr_cmp(r19375, r19376) <= 0, MPFR_RNDN); mpfr_exp(r19378, r19375, MPFR_RNDN); mpfr_mul(r19379, r19375, r19378, MPFR_RNDN); mpfr_set_d(r19380, x, MPFR_RNDN); mpfr_sub(r19381, r19379, r19380, MPFR_RNDN); mpfr_add(r19382, r19378, r19379, MPFR_RNDN); mpfr_div(r19383, r19381, r19382, MPFR_RNDN); mpfr_sub(r19384, r19375, r19383, MPFR_RNDN); ; mpfr_add(r19386, r19385, r19375, MPFR_RNDN); mpfr_div(r19387, r19375, r19386, MPFR_RNDN); mpfr_sub(r19388, r19375, r19387, MPFR_RNDN); mpfr_fma(r19389, r19375, r19378, r19378, MPFR_RNDN); mpfr_div(r19390, r19380, r19389, MPFR_RNDN); mpfr_log1p(r19391, r19390, MPFR_RNDN); mpfr_expm1(r19392, r19391, MPFR_RNDN); mpfr_add(r19393, r19388, r19392, MPFR_RNDN); if (mpfr_get_si(r19377, MPFR_RNDN)) { mpfr_set(r19394, r19384, MPFR_RNDN); } else { mpfr_set(r19394, r19393, MPFR_RNDN); }; return mpfr_get_d(r19394, MPFR_RNDN); } static mpfr_t r19395, r19396, r19397, r19398, r19399, r19400, r19401, r19402, r19403, r19404, r19405, r19406, r19407, r19408, r19409, r19410, r19411, r19412, r19413, r19414; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19395); mpfr_init_set_str(r19396, "7.279548f0", 10, MPFR_RNDN); mpfr_init(r19397); mpfr_init(r19398); mpfr_init(r19399); mpfr_init(r19400); mpfr_init(r19401); mpfr_init(r19402); mpfr_init(r19403); mpfr_init(r19404); mpfr_init_set_str(r19405, "1", 10, MPFR_RNDN); mpfr_init(r19406); mpfr_init(r19407); mpfr_init(r19408); mpfr_init(r19409); mpfr_init(r19410); mpfr_init(r19411); mpfr_init(r19412); mpfr_init(r19413); mpfr_init(r19414); } double f_dm(double wj, double x) { mpfr_set_d(r19395, wj, MPFR_RNDN); ; mpfr_set_si(r19397, mpfr_cmp(r19395, r19396) <= 0, MPFR_RNDN); mpfr_exp(r19398, r19395, MPFR_RNDN); mpfr_mul(r19399, r19395, r19398, MPFR_RNDN); mpfr_set_d(r19400, x, MPFR_RNDN); mpfr_sub(r19401, r19399, r19400, MPFR_RNDN); mpfr_add(r19402, r19398, r19399, MPFR_RNDN); mpfr_div(r19403, r19401, r19402, MPFR_RNDN); mpfr_sub(r19404, r19395, r19403, MPFR_RNDN); ; mpfr_add(r19406, r19405, r19395, MPFR_RNDN); mpfr_div(r19407, r19395, r19406, MPFR_RNDN); mpfr_sub(r19408, r19395, r19407, MPFR_RNDN); mpfr_fma(r19409, r19395, r19398, r19398, MPFR_RNDN); mpfr_div(r19410, r19400, r19409, MPFR_RNDN); mpfr_log1p(r19411, r19410, MPFR_RNDN); mpfr_expm1(r19412, r19411, MPFR_RNDN); mpfr_add(r19413, r19408, r19412, MPFR_RNDN); if (mpfr_get_si(r19397, MPFR_RNDN)) { mpfr_set(r19414, r19404, MPFR_RNDN); } else { mpfr_set(r19414, r19413, MPFR_RNDN); }; return mpfr_get_d(r19414, MPFR_RNDN); }