#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r25056 = wj; float r25057 = exp(r25056); float r25058 = r25056 * r25057; float r25059 = x; float r25060 = r25058 - r25059; float r25061 = r25057 + r25058; float r25062 = r25060 / r25061; float r25063 = r25056 - r25062; return r25063; } double f_id(double wj, double x) { double r25064 = wj; double r25065 = exp(r25064); double r25066 = r25064 * r25065; double r25067 = x; double r25068 = r25066 - r25067; double r25069 = r25065 + r25066; double r25070 = r25068 / r25069; double r25071 = r25064 - r25070; return r25071; } double f_of(float wj, float x) { float r25072 = wj; float r25073 = 1; float r25074 = r25072 + r25073; float r25075 = r25073 / r25074; float r25076 = exp(r25075); float r25077 = log(r25076); float r25078 = x; float r25079 = exp(r25072); float r25080 = r25078 / r25079; float r25081 = r25072 - r25080; float r25082 = r25077 * r25081; float r25083 = r25072 - r25082; float r25084 = 7.565976676171613e-13; bool r25085 = r25083 <= r25084; float r25086 = r25072 * r25072; float r25087 = r25086 - r25072; float r25088 = r25081 * r25087; float r25089 = r25080 - r25088; float r25090 = r25085 ? r25089 : r25083; return r25090; } double f_od(double wj, double x) { double r25091 = wj; double r25092 = 1; double r25093 = r25091 + r25092; double r25094 = r25092 / r25093; double r25095 = exp(r25094); double r25096 = log(r25095); double r25097 = x; double r25098 = exp(r25091); double r25099 = r25097 / r25098; double r25100 = r25091 - r25099; double r25101 = r25096 * r25100; double r25102 = r25091 - r25101; double r25103 = 7.565976676171613e-13; bool r25104 = r25102 <= r25103; double r25105 = r25091 * r25091; double r25106 = r25105 - r25091; double r25107 = r25100 * r25106; double r25108 = r25099 - r25107; double r25109 = r25104 ? r25108 : r25102; return r25109; } 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 r25110, r25111, r25112, r25113, r25114, r25115, r25116, r25117; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r25110); mpfr_init(r25111); mpfr_init(r25112); mpfr_init(r25113); mpfr_init(r25114); mpfr_init(r25115); mpfr_init(r25116); mpfr_init(r25117); } double f_im(double wj, double x) { mpfr_set_d(r25110, wj, MPFR_RNDN); mpfr_exp(r25111, r25110, MPFR_RNDN); mpfr_mul(r25112, r25110, r25111, MPFR_RNDN); mpfr_set_d(r25113, x, MPFR_RNDN); mpfr_sub(r25114, r25112, r25113, MPFR_RNDN); mpfr_add(r25115, r25111, r25112, MPFR_RNDN); mpfr_div(r25116, r25114, r25115, MPFR_RNDN); mpfr_sub(r25117, r25110, r25116, MPFR_RNDN); return mpfr_get_d(r25117, MPFR_RNDN); } static mpfr_t r25118, r25119, r25120, r25121, r25122, r25123, r25124, r25125, r25126, r25127, r25128, r25129, r25130, r25131, r25132, r25133, r25134, r25135, r25136; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r25118); mpfr_init_set_str(r25119, "1", 10, MPFR_RNDN); mpfr_init(r25120); mpfr_init(r25121); mpfr_init(r25122); mpfr_init(r25123); mpfr_init(r25124); mpfr_init(r25125); mpfr_init(r25126); mpfr_init(r25127); mpfr_init(r25128); mpfr_init(r25129); mpfr_init_set_str(r25130, "7.565976676171613e-13", 10, MPFR_RNDN); mpfr_init(r25131); mpfr_init(r25132); mpfr_init(r25133); mpfr_init(r25134); mpfr_init(r25135); mpfr_init(r25136); } double f_fm(double wj, double x) { mpfr_set_d(r25118, wj, MPFR_RNDN); ; mpfr_add(r25120, r25118, r25119, MPFR_RNDN); mpfr_div(r25121, r25119, r25120, MPFR_RNDN); mpfr_exp(r25122, r25121, MPFR_RNDN); mpfr_log(r25123, r25122, MPFR_RNDN); mpfr_set_d(r25124, x, MPFR_RNDN); mpfr_exp(r25125, r25118, MPFR_RNDN); mpfr_div(r25126, r25124, r25125, MPFR_RNDN); mpfr_sub(r25127, r25118, r25126, MPFR_RNDN); mpfr_mul(r25128, r25123, r25127, MPFR_RNDN); mpfr_sub(r25129, r25118, r25128, MPFR_RNDN); ; mpfr_set_si(r25131, mpfr_cmp(r25129, r25130) <= 0, MPFR_RNDN); mpfr_mul(r25132, r25118, r25118, MPFR_RNDN); mpfr_sub(r25133, r25132, r25118, MPFR_RNDN); mpfr_mul(r25134, r25127, r25133, MPFR_RNDN); mpfr_sub(r25135, r25126, r25134, MPFR_RNDN); if (mpfr_get_si(r25131, MPFR_RNDN)) { mpfr_set(r25136, r25135, MPFR_RNDN); } else { mpfr_set(r25136, r25129, MPFR_RNDN); }; return mpfr_get_d(r25136, MPFR_RNDN); } static mpfr_t r25137, r25138, r25139, r25140, r25141, r25142, r25143, r25144, r25145, r25146, r25147, r25148, r25149, r25150, r25151, r25152, r25153, r25154, r25155; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r25137); mpfr_init_set_str(r25138, "1", 10, MPFR_RNDN); mpfr_init(r25139); mpfr_init(r25140); mpfr_init(r25141); mpfr_init(r25142); mpfr_init(r25143); mpfr_init(r25144); mpfr_init(r25145); mpfr_init(r25146); mpfr_init(r25147); mpfr_init(r25148); mpfr_init_set_str(r25149, "7.565976676171613e-13", 10, MPFR_RNDN); mpfr_init(r25150); mpfr_init(r25151); mpfr_init(r25152); mpfr_init(r25153); mpfr_init(r25154); mpfr_init(r25155); } double f_dm(double wj, double x) { mpfr_set_d(r25137, wj, MPFR_RNDN); ; mpfr_add(r25139, r25137, r25138, MPFR_RNDN); mpfr_div(r25140, r25138, r25139, MPFR_RNDN); mpfr_exp(r25141, r25140, MPFR_RNDN); mpfr_log(r25142, r25141, MPFR_RNDN); mpfr_set_d(r25143, x, MPFR_RNDN); mpfr_exp(r25144, r25137, MPFR_RNDN); mpfr_div(r25145, r25143, r25144, MPFR_RNDN); mpfr_sub(r25146, r25137, r25145, MPFR_RNDN); mpfr_mul(r25147, r25142, r25146, MPFR_RNDN); mpfr_sub(r25148, r25137, r25147, MPFR_RNDN); ; mpfr_set_si(r25150, mpfr_cmp(r25148, r25149) <= 0, MPFR_RNDN); mpfr_mul(r25151, r25137, r25137, MPFR_RNDN); mpfr_sub(r25152, r25151, r25137, MPFR_RNDN); mpfr_mul(r25153, r25146, r25152, MPFR_RNDN); mpfr_sub(r25154, r25145, r25153, MPFR_RNDN); if (mpfr_get_si(r25150, MPFR_RNDN)) { mpfr_set(r25155, r25154, MPFR_RNDN); } else { mpfr_set(r25155, r25148, MPFR_RNDN); }; return mpfr_get_d(r25155, MPFR_RNDN); }