#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r20407 = wj; float r20408 = exp(r20407); float r20409 = r20407 * r20408; float r20410 = x; float r20411 = r20409 - r20410; float r20412 = r20408 + r20409; float r20413 = r20411 / r20412; float r20414 = r20407 - r20413; return r20414; } double f_id(double wj, double x) { double r20415 = wj; double r20416 = exp(r20415); double r20417 = r20415 * r20416; double r20418 = x; double r20419 = r20417 - r20418; double r20420 = r20416 + r20417; double r20421 = r20419 / r20420; double r20422 = r20415 - r20421; return r20422; } double f_of(float wj, float x) { float r20423 = wj; float r20424 = 0.6019543409347534f; bool r20425 = r20423 <= r20424; float r20426 = exp(r20423); float r20427 = r20423 * r20426; float r20428 = x; float r20429 = r20427 - r20428; float r20430 = r20426 + r20427; float r20431 = r20429 / r20430; float r20432 = r20423 - r20431; float r20433 = 1.0f; float r20434 = r20423 + r20433; float r20435 = r20423 / r20434; float r20436 = r20433 + r20423; float r20437 = r20428 / r20436; float r20438 = r20437 / r20426; float r20439 = r20435 - r20438; float r20440 = r20423 - r20439; float r20441 = r20425 ? r20432 : r20440; return r20441; } double f_od(double wj, double x) { double r20442 = wj; double r20443 = 0.6019543409347534; bool r20444 = r20442 <= r20443; double r20445 = exp(r20442); double r20446 = r20442 * r20445; double r20447 = x; double r20448 = r20446 - r20447; double r20449 = r20445 + r20446; double r20450 = r20448 / r20449; double r20451 = r20442 - r20450; double r20452 = 1.0; double r20453 = r20442 + r20452; double r20454 = r20442 / r20453; double r20455 = r20452 + r20442; double r20456 = r20447 / r20455; double r20457 = r20456 / r20445; double r20458 = r20454 - r20457; double r20459 = r20442 - r20458; double r20460 = r20444 ? r20451 : r20459; return r20460; } 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 r20461, r20462, r20463, r20464, r20465, r20466, r20467, r20468; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r20461); mpfr_init(r20462); mpfr_init(r20463); mpfr_init(r20464); mpfr_init(r20465); mpfr_init(r20466); mpfr_init(r20467); mpfr_init(r20468); } double f_im(double wj, double x) { mpfr_set_d(r20461, wj, MPFR_RNDN); mpfr_exp(r20462, r20461, MPFR_RNDN); mpfr_mul(r20463, r20461, r20462, MPFR_RNDN); mpfr_set_d(r20464, x, MPFR_RNDN); mpfr_sub(r20465, r20463, r20464, MPFR_RNDN); mpfr_add(r20466, r20462, r20463, MPFR_RNDN); mpfr_div(r20467, r20465, r20466, MPFR_RNDN); mpfr_sub(r20468, r20461, r20467, MPFR_RNDN); return mpfr_get_d(r20468, MPFR_RNDN); } static mpfr_t r20469, r20470, r20471, r20472, r20473, r20474, r20475, r20476, r20477, r20478, r20479, r20480, r20481, r20482, r20483, r20484, r20485, r20486, r20487; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r20469); mpfr_init_set_str(r20470, "0.60195434f0", 10, MPFR_RNDN); mpfr_init(r20471); mpfr_init(r20472); mpfr_init(r20473); mpfr_init(r20474); mpfr_init(r20475); mpfr_init(r20476); mpfr_init(r20477); mpfr_init(r20478); mpfr_init_set_str(r20479, "1", 10, MPFR_RNDN); mpfr_init(r20480); mpfr_init(r20481); mpfr_init(r20482); mpfr_init(r20483); mpfr_init(r20484); mpfr_init(r20485); mpfr_init(r20486); mpfr_init(r20487); } double f_fm(double wj, double x) { mpfr_set_d(r20469, wj, MPFR_RNDN); ; mpfr_set_si(r20471, mpfr_cmp(r20469, r20470) <= 0, MPFR_RNDN); mpfr_exp(r20472, r20469, MPFR_RNDN); mpfr_mul(r20473, r20469, r20472, MPFR_RNDN); mpfr_set_d(r20474, x, MPFR_RNDN); mpfr_sub(r20475, r20473, r20474, MPFR_RNDN); mpfr_add(r20476, r20472, r20473, MPFR_RNDN); mpfr_div(r20477, r20475, r20476, MPFR_RNDN); mpfr_sub(r20478, r20469, r20477, MPFR_RNDN); ; mpfr_add(r20480, r20469, r20479, MPFR_RNDN); mpfr_div(r20481, r20469, r20480, MPFR_RNDN); mpfr_add(r20482, r20479, r20469, MPFR_RNDN); mpfr_div(r20483, r20474, r20482, MPFR_RNDN); mpfr_div(r20484, r20483, r20472, MPFR_RNDN); mpfr_sub(r20485, r20481, r20484, MPFR_RNDN); mpfr_sub(r20486, r20469, r20485, MPFR_RNDN); if (mpfr_get_si(r20471, MPFR_RNDN)) { mpfr_set(r20487, r20478, MPFR_RNDN); } else { mpfr_set(r20487, r20486, MPFR_RNDN); }; return mpfr_get_d(r20487, MPFR_RNDN); } static mpfr_t r20488, r20489, r20490, r20491, r20492, r20493, r20494, r20495, r20496, r20497, r20498, r20499, r20500, r20501, r20502, r20503, r20504, r20505, r20506; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r20488); mpfr_init_set_str(r20489, "0.60195434f0", 10, MPFR_RNDN); mpfr_init(r20490); mpfr_init(r20491); mpfr_init(r20492); mpfr_init(r20493); mpfr_init(r20494); mpfr_init(r20495); mpfr_init(r20496); mpfr_init(r20497); mpfr_init_set_str(r20498, "1", 10, MPFR_RNDN); mpfr_init(r20499); mpfr_init(r20500); mpfr_init(r20501); mpfr_init(r20502); mpfr_init(r20503); mpfr_init(r20504); mpfr_init(r20505); mpfr_init(r20506); } double f_dm(double wj, double x) { mpfr_set_d(r20488, wj, MPFR_RNDN); ; mpfr_set_si(r20490, mpfr_cmp(r20488, r20489) <= 0, MPFR_RNDN); mpfr_exp(r20491, r20488, MPFR_RNDN); mpfr_mul(r20492, r20488, r20491, MPFR_RNDN); mpfr_set_d(r20493, x, MPFR_RNDN); mpfr_sub(r20494, r20492, r20493, MPFR_RNDN); mpfr_add(r20495, r20491, r20492, MPFR_RNDN); mpfr_div(r20496, r20494, r20495, MPFR_RNDN); mpfr_sub(r20497, r20488, r20496, MPFR_RNDN); ; mpfr_add(r20499, r20488, r20498, MPFR_RNDN); mpfr_div(r20500, r20488, r20499, MPFR_RNDN); mpfr_add(r20501, r20498, r20488, MPFR_RNDN); mpfr_div(r20502, r20493, r20501, MPFR_RNDN); mpfr_div(r20503, r20502, r20491, MPFR_RNDN); mpfr_sub(r20504, r20500, r20503, MPFR_RNDN); mpfr_sub(r20505, r20488, r20504, MPFR_RNDN); if (mpfr_get_si(r20490, MPFR_RNDN)) { mpfr_set(r20506, r20497, MPFR_RNDN); } else { mpfr_set(r20506, r20505, MPFR_RNDN); }; return mpfr_get_d(r20506, MPFR_RNDN); }