#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19915 = wj; float r19916 = exp(r19915); float r19917 = r19915 * r19916; float r19918 = x; float r19919 = r19917 - r19918; float r19920 = r19916 + r19917; float r19921 = r19919 / r19920; float r19922 = r19915 - r19921; return r19922; } double f_id(double wj, double x) { double r19923 = wj; double r19924 = exp(r19923); double r19925 = r19923 * r19924; double r19926 = x; double r19927 = r19925 - r19926; double r19928 = r19924 + r19925; double r19929 = r19927 / r19928; double r19930 = r19923 - r19929; return r19930; } double f_of(float wj, float x) { float r19931 = wj; float r19932 = 0.6019543409347534f; bool r19933 = r19931 <= r19932; float r19934 = exp(r19931); float r19935 = r19931 * r19934; float r19936 = x; float r19937 = r19935 - r19936; float r19938 = r19934 + r19935; float r19939 = r19937 / r19938; float r19940 = r19931 - r19939; float r19941 = 1.0f; float r19942 = r19931 + r19941; float r19943 = r19931 / r19942; float r19944 = r19941 + r19931; float r19945 = r19936 / r19944; float r19946 = r19945 / r19934; float r19947 = r19943 - r19946; float r19948 = r19931 - r19947; float r19949 = r19933 ? r19940 : r19948; return r19949; } double f_od(double wj, double x) { double r19950 = wj; double r19951 = 0.6019543409347534; bool r19952 = r19950 <= r19951; double r19953 = exp(r19950); double r19954 = r19950 * r19953; double r19955 = x; double r19956 = r19954 - r19955; double r19957 = r19953 + r19954; double r19958 = r19956 / r19957; double r19959 = r19950 - r19958; double r19960 = 1.0; double r19961 = r19950 + r19960; double r19962 = r19950 / r19961; double r19963 = r19960 + r19950; double r19964 = r19955 / r19963; double r19965 = r19964 / r19953; double r19966 = r19962 - r19965; double r19967 = r19950 - r19966; double r19968 = r19952 ? r19959 : r19967; return r19968; } 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 r19969, r19970, r19971, r19972, r19973, r19974, r19975, r19976; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19969); mpfr_init(r19970); mpfr_init(r19971); mpfr_init(r19972); mpfr_init(r19973); mpfr_init(r19974); mpfr_init(r19975); mpfr_init(r19976); } double f_im(double wj, double x) { mpfr_set_d(r19969, wj, MPFR_RNDN); mpfr_exp(r19970, r19969, MPFR_RNDN); mpfr_mul(r19971, r19969, r19970, MPFR_RNDN); mpfr_set_d(r19972, x, MPFR_RNDN); mpfr_sub(r19973, r19971, r19972, MPFR_RNDN); mpfr_add(r19974, r19970, r19971, MPFR_RNDN); mpfr_div(r19975, r19973, r19974, MPFR_RNDN); mpfr_sub(r19976, r19969, r19975, MPFR_RNDN); return mpfr_get_d(r19976, MPFR_RNDN); } static mpfr_t r19977, r19978, r19979, r19980, r19981, r19982, r19983, r19984, r19985, r19986, r19987, r19988, r19989, r19990, r19991, r19992, r19993, r19994, r19995; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19977); mpfr_init_set_str(r19978, "0.60195434f0", 10, MPFR_RNDN); mpfr_init(r19979); mpfr_init(r19980); mpfr_init(r19981); mpfr_init(r19982); mpfr_init(r19983); mpfr_init(r19984); mpfr_init(r19985); mpfr_init(r19986); mpfr_init_set_str(r19987, "1", 10, MPFR_RNDN); mpfr_init(r19988); mpfr_init(r19989); mpfr_init(r19990); mpfr_init(r19991); mpfr_init(r19992); mpfr_init(r19993); mpfr_init(r19994); mpfr_init(r19995); } double f_fm(double wj, double x) { mpfr_set_d(r19977, wj, MPFR_RNDN); ; mpfr_set_si(r19979, mpfr_cmp(r19977, r19978) <= 0, MPFR_RNDN); mpfr_exp(r19980, r19977, MPFR_RNDN); mpfr_mul(r19981, r19977, r19980, MPFR_RNDN); mpfr_set_d(r19982, x, MPFR_RNDN); mpfr_sub(r19983, r19981, r19982, MPFR_RNDN); mpfr_add(r19984, r19980, r19981, MPFR_RNDN); mpfr_div(r19985, r19983, r19984, MPFR_RNDN); mpfr_sub(r19986, r19977, r19985, MPFR_RNDN); ; mpfr_add(r19988, r19977, r19987, MPFR_RNDN); mpfr_div(r19989, r19977, r19988, MPFR_RNDN); mpfr_add(r19990, r19987, r19977, MPFR_RNDN); mpfr_div(r19991, r19982, r19990, MPFR_RNDN); mpfr_div(r19992, r19991, r19980, MPFR_RNDN); mpfr_sub(r19993, r19989, r19992, MPFR_RNDN); mpfr_sub(r19994, r19977, r19993, MPFR_RNDN); if (mpfr_get_si(r19979, MPFR_RNDN)) { mpfr_set(r19995, r19986, MPFR_RNDN); } else { mpfr_set(r19995, r19994, MPFR_RNDN); }; return mpfr_get_d(r19995, MPFR_RNDN); } static mpfr_t r19996, r19997, r19998, r19999, r20000, r20001, r20002, r20003, r20004, r20005, r20006, r20007, r20008, r20009, r20010, r20011, r20012, r20013, r20014; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19996); mpfr_init_set_str(r19997, "0.60195434f0", 10, MPFR_RNDN); mpfr_init(r19998); mpfr_init(r19999); mpfr_init(r20000); mpfr_init(r20001); mpfr_init(r20002); mpfr_init(r20003); mpfr_init(r20004); mpfr_init(r20005); mpfr_init_set_str(r20006, "1", 10, MPFR_RNDN); mpfr_init(r20007); mpfr_init(r20008); mpfr_init(r20009); mpfr_init(r20010); mpfr_init(r20011); mpfr_init(r20012); mpfr_init(r20013); mpfr_init(r20014); } double f_dm(double wj, double x) { mpfr_set_d(r19996, wj, MPFR_RNDN); ; mpfr_set_si(r19998, mpfr_cmp(r19996, r19997) <= 0, MPFR_RNDN); mpfr_exp(r19999, r19996, MPFR_RNDN); mpfr_mul(r20000, r19996, r19999, MPFR_RNDN); mpfr_set_d(r20001, x, MPFR_RNDN); mpfr_sub(r20002, r20000, r20001, MPFR_RNDN); mpfr_add(r20003, r19999, r20000, MPFR_RNDN); mpfr_div(r20004, r20002, r20003, MPFR_RNDN); mpfr_sub(r20005, r19996, r20004, MPFR_RNDN); ; mpfr_add(r20007, r19996, r20006, MPFR_RNDN); mpfr_div(r20008, r19996, r20007, MPFR_RNDN); mpfr_add(r20009, r20006, r19996, MPFR_RNDN); mpfr_div(r20010, r20001, r20009, MPFR_RNDN); mpfr_div(r20011, r20010, r19999, MPFR_RNDN); mpfr_sub(r20012, r20008, r20011, MPFR_RNDN); mpfr_sub(r20013, r19996, r20012, MPFR_RNDN); if (mpfr_get_si(r19998, MPFR_RNDN)) { mpfr_set(r20014, r20005, MPFR_RNDN); } else { mpfr_set(r20014, r20013, MPFR_RNDN); }; return mpfr_get_d(r20014, MPFR_RNDN); }