#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19891 = wj; float r19892 = exp(r19891); float r19893 = r19891 * r19892; float r19894 = x; float r19895 = r19893 - r19894; float r19896 = r19892 + r19893; float r19897 = r19895 / r19896; float r19898 = r19891 - r19897; return r19898; } double f_id(double wj, double x) { double r19899 = wj; double r19900 = exp(r19899); double r19901 = r19899 * r19900; double r19902 = x; double r19903 = r19901 - r19902; double r19904 = r19900 + r19901; double r19905 = r19903 / r19904; double r19906 = r19899 - r19905; return r19906; } double f_of(float wj, float x) { float r19907 = wj; float r19908 = exp(r19907); float r19909 = r19907 * r19908; float r19910 = x; float r19911 = r19909 - r19910; float r19912 = r19908 + r19909; float r19913 = r19911 / r19912; float r19914 = 2.4270278086576645e-11f; bool r19915 = r19913 <= r19914; float r19916 = r19907 * r19907; float r19917 = r19910 + r19916; float r19918 = 1.0f; float r19919 = r19907 + r19918; float r19920 = r19907 / r19919; float r19921 = r19918 + r19907; float r19922 = r19910 / r19921; float r19923 = r19922 / r19908; float r19924 = r19920 - r19923; float r19925 = r19907 - r19924; float r19926 = r19915 ? r19917 : r19925; return r19926; } double f_od(double wj, double x) { double r19927 = wj; double r19928 = exp(r19927); double r19929 = r19927 * r19928; double r19930 = x; double r19931 = r19929 - r19930; double r19932 = r19928 + r19929; double r19933 = r19931 / r19932; double r19934 = 2.4270278086576645e-11; bool r19935 = r19933 <= r19934; double r19936 = r19927 * r19927; double r19937 = r19930 + r19936; double r19938 = 1.0; double r19939 = r19927 + r19938; double r19940 = r19927 / r19939; double r19941 = r19938 + r19927; double r19942 = r19930 / r19941; double r19943 = r19942 / r19928; double r19944 = r19940 - r19943; double r19945 = r19927 - r19944; double r19946 = r19935 ? r19937 : r19945; return r19946; } 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 r19947, r19948, r19949, r19950, r19951, r19952, r19953, r19954; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19947); mpfr_init(r19948); mpfr_init(r19949); mpfr_init(r19950); mpfr_init(r19951); mpfr_init(r19952); mpfr_init(r19953); mpfr_init(r19954); } double f_im(double wj, double x) { mpfr_set_d(r19947, wj, MPFR_RNDN); mpfr_exp(r19948, r19947, MPFR_RNDN); mpfr_mul(r19949, r19947, r19948, MPFR_RNDN); mpfr_set_d(r19950, x, MPFR_RNDN); mpfr_sub(r19951, r19949, r19950, MPFR_RNDN); mpfr_add(r19952, r19948, r19949, MPFR_RNDN); mpfr_div(r19953, r19951, r19952, MPFR_RNDN); mpfr_sub(r19954, r19947, r19953, MPFR_RNDN); return mpfr_get_d(r19954, MPFR_RNDN); } static mpfr_t r19955, r19956, r19957, r19958, r19959, r19960, r19961, r19962, r19963, r19964, r19965, r19966, r19967, r19968, r19969, r19970, r19971, r19972, r19973, r19974; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19955); mpfr_init(r19956); mpfr_init(r19957); mpfr_init(r19958); mpfr_init(r19959); mpfr_init(r19960); mpfr_init(r19961); mpfr_init_set_str(r19962, "2.4270278086576645e-11", 10, MPFR_RNDN); mpfr_init(r19963); mpfr_init(r19964); mpfr_init(r19965); mpfr_init_set_str(r19966, "1", 10, MPFR_RNDN); mpfr_init(r19967); mpfr_init(r19968); mpfr_init(r19969); mpfr_init(r19970); mpfr_init(r19971); mpfr_init(r19972); mpfr_init(r19973); mpfr_init(r19974); } double f_fm(double wj, double x) { mpfr_set_d(r19955, wj, MPFR_RNDN); mpfr_exp(r19956, r19955, MPFR_RNDN); mpfr_mul(r19957, r19955, r19956, MPFR_RNDN); mpfr_set_d(r19958, x, MPFR_RNDN); mpfr_sub(r19959, r19957, r19958, MPFR_RNDN); mpfr_add(r19960, r19956, r19957, MPFR_RNDN); mpfr_div(r19961, r19959, r19960, MPFR_RNDN); ; mpfr_set_si(r19963, mpfr_cmp(r19961, r19962) <= 0, MPFR_RNDN); mpfr_sqr(r19964, r19955, MPFR_RNDN); mpfr_add(r19965, r19958, r19964, MPFR_RNDN); ; mpfr_add(r19967, r19955, r19966, MPFR_RNDN); mpfr_div(r19968, r19955, r19967, MPFR_RNDN); mpfr_add(r19969, r19966, r19955, MPFR_RNDN); mpfr_div(r19970, r19958, r19969, MPFR_RNDN); mpfr_div(r19971, r19970, r19956, MPFR_RNDN); mpfr_sub(r19972, r19968, r19971, MPFR_RNDN); mpfr_sub(r19973, r19955, r19972, MPFR_RNDN); if (mpfr_get_si(r19963, MPFR_RNDN)) { mpfr_set(r19974, r19965, MPFR_RNDN); } else { mpfr_set(r19974, r19973, MPFR_RNDN); }; return mpfr_get_d(r19974, MPFR_RNDN); } static mpfr_t r19975, r19976, r19977, r19978, r19979, r19980, r19981, r19982, r19983, r19984, r19985, r19986, r19987, r19988, r19989, r19990, r19991, r19992, r19993, r19994; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19975); mpfr_init(r19976); mpfr_init(r19977); mpfr_init(r19978); mpfr_init(r19979); mpfr_init(r19980); mpfr_init(r19981); mpfr_init_set_str(r19982, "2.4270278086576645e-11", 10, MPFR_RNDN); mpfr_init(r19983); mpfr_init(r19984); mpfr_init(r19985); mpfr_init_set_str(r19986, "1", 10, MPFR_RNDN); mpfr_init(r19987); mpfr_init(r19988); mpfr_init(r19989); mpfr_init(r19990); mpfr_init(r19991); mpfr_init(r19992); mpfr_init(r19993); mpfr_init(r19994); } double f_dm(double wj, double x) { mpfr_set_d(r19975, wj, MPFR_RNDN); mpfr_exp(r19976, r19975, MPFR_RNDN); mpfr_mul(r19977, r19975, r19976, MPFR_RNDN); mpfr_set_d(r19978, x, MPFR_RNDN); mpfr_sub(r19979, r19977, r19978, MPFR_RNDN); mpfr_add(r19980, r19976, r19977, MPFR_RNDN); mpfr_div(r19981, r19979, r19980, MPFR_RNDN); ; mpfr_set_si(r19983, mpfr_cmp(r19981, r19982) <= 0, MPFR_RNDN); mpfr_sqr(r19984, r19975, MPFR_RNDN); mpfr_add(r19985, r19978, r19984, MPFR_RNDN); ; mpfr_add(r19987, r19975, r19986, MPFR_RNDN); mpfr_div(r19988, r19975, r19987, MPFR_RNDN); mpfr_add(r19989, r19986, r19975, MPFR_RNDN); mpfr_div(r19990, r19978, r19989, MPFR_RNDN); mpfr_div(r19991, r19990, r19976, MPFR_RNDN); mpfr_sub(r19992, r19988, r19991, MPFR_RNDN); mpfr_sub(r19993, r19975, r19992, MPFR_RNDN); if (mpfr_get_si(r19983, MPFR_RNDN)) { mpfr_set(r19994, r19985, MPFR_RNDN); } else { mpfr_set(r19994, r19993, MPFR_RNDN); }; return mpfr_get_d(r19994, MPFR_RNDN); }