#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r24863 = wj; float r24864 = exp(r24863); float r24865 = r24863 * r24864; float r24866 = x; float r24867 = r24865 - r24866; float r24868 = r24864 + r24865; float r24869 = r24867 / r24868; float r24870 = r24863 - r24869; return r24870; } double f_id(double wj, double x) { double r24871 = wj; double r24872 = exp(r24871); double r24873 = r24871 * r24872; double r24874 = x; double r24875 = r24873 - r24874; double r24876 = r24872 + r24873; double r24877 = r24875 / r24876; double r24878 = r24871 - r24877; return r24878; } double f_of(float wj, float x) { float r24879 = wj; float r24880 = 1; float r24881 = r24880 + r24879; float r24882 = r24879 / r24881; float r24883 = r24879 - r24882; float r24884 = x; float r24885 = exp(r24879); float r24886 = r24879 * r24885; float r24887 = r24885 + r24886; float r24888 = r24884 / r24887; float r24889 = r24883 + r24888; float r24890 = 3.639973196389132e-18; bool r24891 = r24889 <= r24890; float r24892 = 2; float r24893 = pow(r24879, r24892); float r24894 = r24893 + r24884; float r24895 = r24879 * r24884; float r24896 = r24892 * r24895; float r24897 = r24894 - r24896; float r24898 = r24891 ? r24897 : r24889; return r24898; } double f_od(double wj, double x) { double r24899 = wj; double r24900 = 1; double r24901 = r24900 + r24899; double r24902 = r24899 / r24901; double r24903 = r24899 - r24902; double r24904 = x; double r24905 = exp(r24899); double r24906 = r24899 * r24905; double r24907 = r24905 + r24906; double r24908 = r24904 / r24907; double r24909 = r24903 + r24908; double r24910 = 3.639973196389132e-18; bool r24911 = r24909 <= r24910; double r24912 = 2; double r24913 = pow(r24899, r24912); double r24914 = r24913 + r24904; double r24915 = r24899 * r24904; double r24916 = r24912 * r24915; double r24917 = r24914 - r24916; double r24918 = r24911 ? r24917 : r24909; return r24918; } 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 r24919, r24920, r24921, r24922, r24923, r24924, r24925, r24926; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r24919); mpfr_init(r24920); mpfr_init(r24921); mpfr_init(r24922); mpfr_init(r24923); mpfr_init(r24924); mpfr_init(r24925); mpfr_init(r24926); } double f_im(double wj, double x) { mpfr_set_d(r24919, wj, MPFR_RNDN); mpfr_exp(r24920, r24919, MPFR_RNDN); mpfr_mul(r24921, r24919, r24920, MPFR_RNDN); mpfr_set_d(r24922, x, MPFR_RNDN); mpfr_sub(r24923, r24921, r24922, MPFR_RNDN); mpfr_add(r24924, r24920, r24921, MPFR_RNDN); mpfr_div(r24925, r24923, r24924, MPFR_RNDN); mpfr_sub(r24926, r24919, r24925, MPFR_RNDN); return mpfr_get_d(r24926, MPFR_RNDN); } static mpfr_t r24927, r24928, r24929, r24930, r24931, r24932, r24933, r24934, r24935, r24936, r24937, r24938, r24939, r24940, r24941, r24942, r24943, r24944, r24945, r24946; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r24927); mpfr_init_set_str(r24928, "1", 10, MPFR_RNDN); mpfr_init(r24929); mpfr_init(r24930); mpfr_init(r24931); mpfr_init(r24932); mpfr_init(r24933); mpfr_init(r24934); mpfr_init(r24935); mpfr_init(r24936); mpfr_init(r24937); mpfr_init_set_str(r24938, "3.639973196389132e-18", 10, MPFR_RNDN); mpfr_init(r24939); mpfr_init_set_str(r24940, "2", 10, MPFR_RNDN); mpfr_init(r24941); mpfr_init(r24942); mpfr_init(r24943); mpfr_init(r24944); mpfr_init(r24945); mpfr_init(r24946); } double f_fm(double wj, double x) { mpfr_set_d(r24927, wj, MPFR_RNDN); ; mpfr_add(r24929, r24928, r24927, MPFR_RNDN); mpfr_div(r24930, r24927, r24929, MPFR_RNDN); mpfr_sub(r24931, r24927, r24930, MPFR_RNDN); mpfr_set_d(r24932, x, MPFR_RNDN); mpfr_exp(r24933, r24927, MPFR_RNDN); mpfr_mul(r24934, r24927, r24933, MPFR_RNDN); mpfr_add(r24935, r24933, r24934, MPFR_RNDN); mpfr_div(r24936, r24932, r24935, MPFR_RNDN); mpfr_add(r24937, r24931, r24936, MPFR_RNDN); ; mpfr_set_si(r24939, mpfr_cmp(r24937, r24938) <= 0, MPFR_RNDN); ; mpfr_pow(r24941, r24927, r24940, MPFR_RNDN); mpfr_add(r24942, r24941, r24932, MPFR_RNDN); mpfr_mul(r24943, r24927, r24932, MPFR_RNDN); mpfr_mul(r24944, r24940, r24943, MPFR_RNDN); mpfr_sub(r24945, r24942, r24944, MPFR_RNDN); if (mpfr_get_si(r24939, MPFR_RNDN)) { mpfr_set(r24946, r24945, MPFR_RNDN); } else { mpfr_set(r24946, r24937, MPFR_RNDN); }; return mpfr_get_d(r24946, MPFR_RNDN); } static mpfr_t r24947, r24948, r24949, r24950, r24951, r24952, r24953, r24954, r24955, r24956, r24957, r24958, r24959, r24960, r24961, r24962, r24963, r24964, r24965, r24966; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r24947); mpfr_init_set_str(r24948, "1", 10, MPFR_RNDN); mpfr_init(r24949); mpfr_init(r24950); mpfr_init(r24951); mpfr_init(r24952); mpfr_init(r24953); mpfr_init(r24954); mpfr_init(r24955); mpfr_init(r24956); mpfr_init(r24957); mpfr_init_set_str(r24958, "3.639973196389132e-18", 10, MPFR_RNDN); mpfr_init(r24959); mpfr_init_set_str(r24960, "2", 10, MPFR_RNDN); mpfr_init(r24961); mpfr_init(r24962); mpfr_init(r24963); mpfr_init(r24964); mpfr_init(r24965); mpfr_init(r24966); } double f_dm(double wj, double x) { mpfr_set_d(r24947, wj, MPFR_RNDN); ; mpfr_add(r24949, r24948, r24947, MPFR_RNDN); mpfr_div(r24950, r24947, r24949, MPFR_RNDN); mpfr_sub(r24951, r24947, r24950, MPFR_RNDN); mpfr_set_d(r24952, x, MPFR_RNDN); mpfr_exp(r24953, r24947, MPFR_RNDN); mpfr_mul(r24954, r24947, r24953, MPFR_RNDN); mpfr_add(r24955, r24953, r24954, MPFR_RNDN); mpfr_div(r24956, r24952, r24955, MPFR_RNDN); mpfr_add(r24957, r24951, r24956, MPFR_RNDN); ; mpfr_set_si(r24959, mpfr_cmp(r24957, r24958) <= 0, MPFR_RNDN); ; mpfr_pow(r24961, r24947, r24960, MPFR_RNDN); mpfr_add(r24962, r24961, r24952, MPFR_RNDN); mpfr_mul(r24963, r24947, r24952, MPFR_RNDN); mpfr_mul(r24964, r24960, r24963, MPFR_RNDN); mpfr_sub(r24965, r24962, r24964, MPFR_RNDN); if (mpfr_get_si(r24959, MPFR_RNDN)) { mpfr_set(r24966, r24965, MPFR_RNDN); } else { mpfr_set(r24966, r24957, MPFR_RNDN); }; return mpfr_get_d(r24966, MPFR_RNDN); }