#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r25490 = wj; float r25491 = exp(r25490); float r25492 = r25490 * r25491; float r25493 = x; float r25494 = r25492 - r25493; float r25495 = r25491 + r25492; float r25496 = r25494 / r25495; float r25497 = r25490 - r25496; return r25497; } double f_id(double wj, double x) { double r25498 = wj; double r25499 = exp(r25498); double r25500 = r25498 * r25499; double r25501 = x; double r25502 = r25500 - r25501; double r25503 = r25499 + r25500; double r25504 = r25502 / r25503; double r25505 = r25498 - r25504; return r25505; } double f_of(float wj, float x) { float r25506 = wj; float r25507 = exp(r25506); float r25508 = r25506 * r25507; float r25509 = x; float r25510 = r25508 - r25509; float r25511 = r25507 + r25508; float r25512 = r25510 / r25511; float r25513 = r25506 - r25512; float r25514 = 2.2390702063761374e-15; bool r25515 = r25513 <= r25514; float r25516 = fma(r25506, r25506, r25509); float r25517 = 1; float r25518 = r25506 + r25517; float r25519 = r25506 / r25518; float r25520 = fma(r25507, r25506, r25507); float r25521 = r25509 / r25520; float r25522 = r25519 - r25521; float r25523 = r25506 - r25522; float r25524 = r25515 ? r25516 : r25523; return r25524; } double f_od(double wj, double x) { double r25525 = wj; double r25526 = exp(r25525); double r25527 = r25525 * r25526; double r25528 = x; double r25529 = r25527 - r25528; double r25530 = r25526 + r25527; double r25531 = r25529 / r25530; double r25532 = r25525 - r25531; double r25533 = 2.2390702063761374e-15; bool r25534 = r25532 <= r25533; double r25535 = fma(r25525, r25525, r25528); double r25536 = 1; double r25537 = r25525 + r25536; double r25538 = r25525 / r25537; double r25539 = fma(r25526, r25525, r25526); double r25540 = r25528 / r25539; double r25541 = r25538 - r25540; double r25542 = r25525 - r25541; double r25543 = r25534 ? r25535 : r25542; return r25543; } 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 r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r25544); mpfr_init(r25545); mpfr_init(r25546); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); } double f_im(double wj, double x) { mpfr_set_d(r25544, wj, MPFR_RNDN); mpfr_exp(r25545, r25544, MPFR_RNDN); mpfr_mul(r25546, r25544, r25545, MPFR_RNDN); mpfr_set_d(r25547, x, MPFR_RNDN); mpfr_sub(r25548, r25546, r25547, MPFR_RNDN); mpfr_add(r25549, r25545, r25546, MPFR_RNDN); mpfr_div(r25550, r25548, r25549, MPFR_RNDN); mpfr_sub(r25551, r25544, r25550, MPFR_RNDN); return mpfr_get_d(r25551, MPFR_RNDN); } static mpfr_t r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567, r25568, r25569, r25570; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init_set_str(r25560, "2.2390702063761374e-15", 10, MPFR_RNDN); mpfr_init(r25561); mpfr_init(r25562); mpfr_init_set_str(r25563, "1", 10, MPFR_RNDN); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); } double f_fm(double wj, double x) { mpfr_set_d(r25552, wj, MPFR_RNDN); mpfr_exp(r25553, r25552, MPFR_RNDN); mpfr_mul(r25554, r25552, r25553, MPFR_RNDN); mpfr_set_d(r25555, x, MPFR_RNDN); mpfr_sub(r25556, r25554, r25555, MPFR_RNDN); mpfr_add(r25557, r25553, r25554, MPFR_RNDN); mpfr_div(r25558, r25556, r25557, MPFR_RNDN); mpfr_sub(r25559, r25552, r25558, MPFR_RNDN); ; mpfr_set_si(r25561, mpfr_cmp(r25559, r25560) <= 0, MPFR_RNDN); mpfr_fma(r25562, r25552, r25552, r25555, MPFR_RNDN); ; mpfr_add(r25564, r25552, r25563, MPFR_RNDN); mpfr_div(r25565, r25552, r25564, MPFR_RNDN); mpfr_fma(r25566, r25553, r25552, r25553, MPFR_RNDN); mpfr_div(r25567, r25555, r25566, MPFR_RNDN); mpfr_sub(r25568, r25565, r25567, MPFR_RNDN); mpfr_sub(r25569, r25552, r25568, MPFR_RNDN); if (mpfr_get_si(r25561, MPFR_RNDN)) { mpfr_set(r25570, r25562, MPFR_RNDN); } else { mpfr_set(r25570, r25569, MPFR_RNDN); }; return mpfr_get_d(r25570, MPFR_RNDN); } static mpfr_t r25571, r25572, r25573, r25574, r25575, r25576, r25577, r25578, r25579, r25580, r25581, r25582, r25583, r25584, r25585, r25586, r25587, r25588, r25589; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); mpfr_init(r25574); mpfr_init(r25575); mpfr_init(r25576); mpfr_init(r25577); mpfr_init(r25578); mpfr_init_set_str(r25579, "2.2390702063761374e-15", 10, MPFR_RNDN); mpfr_init(r25580); mpfr_init(r25581); mpfr_init_set_str(r25582, "1", 10, MPFR_RNDN); mpfr_init(r25583); mpfr_init(r25584); mpfr_init(r25585); mpfr_init(r25586); mpfr_init(r25587); mpfr_init(r25588); mpfr_init(r25589); } double f_dm(double wj, double x) { mpfr_set_d(r25571, wj, MPFR_RNDN); mpfr_exp(r25572, r25571, MPFR_RNDN); mpfr_mul(r25573, r25571, r25572, MPFR_RNDN); mpfr_set_d(r25574, x, MPFR_RNDN); mpfr_sub(r25575, r25573, r25574, MPFR_RNDN); mpfr_add(r25576, r25572, r25573, MPFR_RNDN); mpfr_div(r25577, r25575, r25576, MPFR_RNDN); mpfr_sub(r25578, r25571, r25577, MPFR_RNDN); ; mpfr_set_si(r25580, mpfr_cmp(r25578, r25579) <= 0, MPFR_RNDN); mpfr_fma(r25581, r25571, r25571, r25574, MPFR_RNDN); ; mpfr_add(r25583, r25571, r25582, MPFR_RNDN); mpfr_div(r25584, r25571, r25583, MPFR_RNDN); mpfr_fma(r25585, r25572, r25571, r25572, MPFR_RNDN); mpfr_div(r25586, r25574, r25585, MPFR_RNDN); mpfr_sub(r25587, r25584, r25586, MPFR_RNDN); mpfr_sub(r25588, r25571, r25587, MPFR_RNDN); if (mpfr_get_si(r25580, MPFR_RNDN)) { mpfr_set(r25589, r25581, MPFR_RNDN); } else { mpfr_set(r25589, r25588, MPFR_RNDN); }; return mpfr_get_d(r25589, MPFR_RNDN); }