#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r18463 = wj; float r18464 = exp(r18463); float r18465 = r18463 * r18464; float r18466 = x; float r18467 = r18465 - r18466; float r18468 = r18464 + r18465; float r18469 = r18467 / r18468; float r18470 = r18463 - r18469; return r18470; } double f_id(double wj, double x) { double r18471 = wj; double r18472 = exp(r18471); double r18473 = r18471 * r18472; double r18474 = x; double r18475 = r18473 - r18474; double r18476 = r18472 + r18473; double r18477 = r18475 / r18476; double r18478 = r18471 - r18477; return r18478; } double f_of(float wj, float x) { float r18479 = wj; float r18480 = 2.8588366508483887f; bool r18481 = r18479 <= r18480; float r18482 = exp(r18479); float r18483 = r18479 * r18482; float r18484 = x; float r18485 = r18483 - r18484; float r18486 = r18482 + r18483; float r18487 = r18485 / r18486; float r18488 = r18479 - r18487; float r18489 = 1.0f; float r18490 = r18479 + r18489; float r18491 = r18479 / r18490; float r18492 = fma(r18479, r18482, r18482); float r18493 = r18484 / r18492; float r18494 = r18491 - r18493; float r18495 = r18479 - r18494; float r18496 = r18481 ? r18488 : r18495; return r18496; } double f_od(double wj, double x) { double r18497 = wj; double r18498 = 2.8588366508483887; bool r18499 = r18497 <= r18498; double r18500 = exp(r18497); double r18501 = r18497 * r18500; double r18502 = x; double r18503 = r18501 - r18502; double r18504 = r18500 + r18501; double r18505 = r18503 / r18504; double r18506 = r18497 - r18505; double r18507 = 1.0; double r18508 = r18497 + r18507; double r18509 = r18497 / r18508; double r18510 = fma(r18497, r18500, r18500); double r18511 = r18502 / r18510; double r18512 = r18509 - r18511; double r18513 = r18497 - r18512; double r18514 = r18499 ? r18506 : r18513; return r18514; } 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 r18515, r18516, r18517, r18518, r18519, r18520, r18521, r18522; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18515); mpfr_init(r18516); mpfr_init(r18517); mpfr_init(r18518); mpfr_init(r18519); mpfr_init(r18520); mpfr_init(r18521); mpfr_init(r18522); } double f_im(double wj, double x) { mpfr_set_d(r18515, wj, MPFR_RNDN); mpfr_exp(r18516, r18515, MPFR_RNDN); mpfr_mul(r18517, r18515, r18516, MPFR_RNDN); mpfr_set_d(r18518, x, MPFR_RNDN); mpfr_sub(r18519, r18517, r18518, MPFR_RNDN); mpfr_add(r18520, r18516, r18517, MPFR_RNDN); mpfr_div(r18521, r18519, r18520, MPFR_RNDN); mpfr_sub(r18522, r18515, r18521, MPFR_RNDN); return mpfr_get_d(r18522, MPFR_RNDN); } static mpfr_t r18523, r18524, r18525, r18526, r18527, r18528, r18529, r18530, r18531, r18532, r18533, r18534, r18535, r18536, r18537, r18538, r18539, r18540; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18523); mpfr_init_set_str(r18524, "2.8588367f0", 10, MPFR_RNDN); mpfr_init(r18525); mpfr_init(r18526); mpfr_init(r18527); mpfr_init(r18528); mpfr_init(r18529); mpfr_init(r18530); mpfr_init(r18531); mpfr_init(r18532); mpfr_init_set_str(r18533, "1", 10, MPFR_RNDN); mpfr_init(r18534); mpfr_init(r18535); mpfr_init(r18536); mpfr_init(r18537); mpfr_init(r18538); mpfr_init(r18539); mpfr_init(r18540); } double f_fm(double wj, double x) { mpfr_set_d(r18523, wj, MPFR_RNDN); ; mpfr_set_si(r18525, mpfr_cmp(r18523, r18524) <= 0, MPFR_RNDN); mpfr_exp(r18526, r18523, MPFR_RNDN); mpfr_mul(r18527, r18523, r18526, MPFR_RNDN); mpfr_set_d(r18528, x, MPFR_RNDN); mpfr_sub(r18529, r18527, r18528, MPFR_RNDN); mpfr_add(r18530, r18526, r18527, MPFR_RNDN); mpfr_div(r18531, r18529, r18530, MPFR_RNDN); mpfr_sub(r18532, r18523, r18531, MPFR_RNDN); ; mpfr_add(r18534, r18523, r18533, MPFR_RNDN); mpfr_div(r18535, r18523, r18534, MPFR_RNDN); mpfr_fma(r18536, r18523, r18526, r18526, MPFR_RNDN); mpfr_div(r18537, r18528, r18536, MPFR_RNDN); mpfr_sub(r18538, r18535, r18537, MPFR_RNDN); mpfr_sub(r18539, r18523, r18538, MPFR_RNDN); if (mpfr_get_si(r18525, MPFR_RNDN)) { mpfr_set(r18540, r18532, MPFR_RNDN); } else { mpfr_set(r18540, r18539, MPFR_RNDN); }; return mpfr_get_d(r18540, MPFR_RNDN); } static mpfr_t r18541, r18542, r18543, r18544, r18545, r18546, r18547, r18548, r18549, r18550, r18551, r18552, r18553, r18554, r18555, r18556, r18557, r18558; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18541); mpfr_init_set_str(r18542, "2.8588367f0", 10, MPFR_RNDN); mpfr_init(r18543); mpfr_init(r18544); mpfr_init(r18545); mpfr_init(r18546); mpfr_init(r18547); mpfr_init(r18548); mpfr_init(r18549); mpfr_init(r18550); mpfr_init_set_str(r18551, "1", 10, MPFR_RNDN); mpfr_init(r18552); mpfr_init(r18553); mpfr_init(r18554); mpfr_init(r18555); mpfr_init(r18556); mpfr_init(r18557); mpfr_init(r18558); } double f_dm(double wj, double x) { mpfr_set_d(r18541, wj, MPFR_RNDN); ; mpfr_set_si(r18543, mpfr_cmp(r18541, r18542) <= 0, MPFR_RNDN); mpfr_exp(r18544, r18541, MPFR_RNDN); mpfr_mul(r18545, r18541, r18544, MPFR_RNDN); mpfr_set_d(r18546, x, MPFR_RNDN); mpfr_sub(r18547, r18545, r18546, MPFR_RNDN); mpfr_add(r18548, r18544, r18545, MPFR_RNDN); mpfr_div(r18549, r18547, r18548, MPFR_RNDN); mpfr_sub(r18550, r18541, r18549, MPFR_RNDN); ; mpfr_add(r18552, r18541, r18551, MPFR_RNDN); mpfr_div(r18553, r18541, r18552, MPFR_RNDN); mpfr_fma(r18554, r18541, r18544, r18544, MPFR_RNDN); mpfr_div(r18555, r18546, r18554, MPFR_RNDN); mpfr_sub(r18556, r18553, r18555, MPFR_RNDN); mpfr_sub(r18557, r18541, r18556, MPFR_RNDN); if (mpfr_get_si(r18543, MPFR_RNDN)) { mpfr_set(r18558, r18550, MPFR_RNDN); } else { mpfr_set(r18558, r18557, MPFR_RNDN); }; return mpfr_get_d(r18558, MPFR_RNDN); }