#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r19803 = wj; float r19804 = exp(r19803); float r19805 = r19803 * r19804; float r19806 = x; float r19807 = r19805 - r19806; float r19808 = r19804 + r19805; float r19809 = r19807 / r19808; float r19810 = r19803 - r19809; return r19810; } double f_id(double wj, double x) { double r19811 = wj; double r19812 = exp(r19811); double r19813 = r19811 * r19812; double r19814 = x; double r19815 = r19813 - r19814; double r19816 = r19812 + r19813; double r19817 = r19815 / r19816; double r19818 = r19811 - r19817; return r19818; } double f_of(float wj, float x) { float r19819 = wj; float r19820 = -2.63216061572775e-13f; bool r19821 = r19819 <= r19820; float r19822 = exp(r19819); float r19823 = r19822 * r19819; float r19824 = x; float r19825 = r19823 - r19824; float r19826 = -r19825; float r19827 = 1.0f; float r19828 = r19827 + r19819; float r19829 = r19828 * r19822; float r19830 = r19826 / r19829; float r19831 = r19819 + r19830; float r19832 = 1.522317513311434e-06f; bool r19833 = r19819 <= r19832; float r19834 = r19819 * r19819; float r19835 = r19834 + r19824; float r19836 = 2.0f; float r19837 = r19819 * r19824; float r19838 = r19836 * r19837; float r19839 = r19835 - r19838; float r19840 = r19819 / r19828; float r19841 = r19819 - r19840; float r19842 = r19819 * r19822; float r19843 = r19822 + r19842; float r19844 = r19824 / r19843; float r19845 = r19841 + r19844; float r19846 = r19833 ? r19839 : r19845; float r19847 = r19821 ? r19831 : r19846; return r19847; } double f_od(double wj, double x) { double r19848 = wj; double r19849 = -2.63216061572775e-13; bool r19850 = r19848 <= r19849; double r19851 = exp(r19848); double r19852 = r19851 * r19848; double r19853 = x; double r19854 = r19852 - r19853; double r19855 = -r19854; double r19856 = 1.0; double r19857 = r19856 + r19848; double r19858 = r19857 * r19851; double r19859 = r19855 / r19858; double r19860 = r19848 + r19859; double r19861 = 1.522317513311434e-06; bool r19862 = r19848 <= r19861; double r19863 = r19848 * r19848; double r19864 = r19863 + r19853; double r19865 = 2.0; double r19866 = r19848 * r19853; double r19867 = r19865 * r19866; double r19868 = r19864 - r19867; double r19869 = r19848 / r19857; double r19870 = r19848 - r19869; double r19871 = r19848 * r19851; double r19872 = r19851 + r19871; double r19873 = r19853 / r19872; double r19874 = r19870 + r19873; double r19875 = r19862 ? r19868 : r19874; double r19876 = r19850 ? r19860 : r19875; return r19876; } 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 r19877, r19878, r19879, r19880, r19881, r19882, r19883, r19884; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r19877); mpfr_init(r19878); mpfr_init(r19879); mpfr_init(r19880); mpfr_init(r19881); mpfr_init(r19882); mpfr_init(r19883); mpfr_init(r19884); } double f_im(double wj, double x) { mpfr_set_d(r19877, wj, MPFR_RNDN); mpfr_exp(r19878, r19877, MPFR_RNDN); mpfr_mul(r19879, r19877, r19878, MPFR_RNDN); mpfr_set_d(r19880, x, MPFR_RNDN); mpfr_sub(r19881, r19879, r19880, MPFR_RNDN); mpfr_add(r19882, r19878, r19879, MPFR_RNDN); mpfr_div(r19883, r19881, r19882, MPFR_RNDN); mpfr_sub(r19884, r19877, r19883, MPFR_RNDN); return mpfr_get_d(r19884, MPFR_RNDN); } static mpfr_t r19885, r19886, r19887, r19888, r19889, r19890, r19891, r19892, r19893, r19894, r19895, r19896, r19897, r19898, r19899, r19900, r19901, r19902, r19903, r19904, r19905, r19906, r19907, r19908, r19909, r19910, r19911, r19912, r19913; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r19885); mpfr_init_set_str(r19886, "-2.63216061572775e-13", 10, MPFR_RNDN); mpfr_init(r19887); mpfr_init(r19888); mpfr_init(r19889); mpfr_init(r19890); mpfr_init(r19891); mpfr_init(r19892); mpfr_init_set_str(r19893, "1", 10, MPFR_RNDN); mpfr_init(r19894); mpfr_init(r19895); mpfr_init(r19896); mpfr_init(r19897); mpfr_init_set_str(r19898, "1.522317513311434e-06", 10, MPFR_RNDN); mpfr_init(r19899); mpfr_init(r19900); mpfr_init(r19901); mpfr_init_set_str(r19902, "2", 10, MPFR_RNDN); mpfr_init(r19903); mpfr_init(r19904); mpfr_init(r19905); mpfr_init(r19906); mpfr_init(r19907); mpfr_init(r19908); mpfr_init(r19909); mpfr_init(r19910); mpfr_init(r19911); mpfr_init(r19912); mpfr_init(r19913); } double f_fm(double wj, double x) { mpfr_set_d(r19885, wj, MPFR_RNDN); ; mpfr_set_si(r19887, mpfr_cmp(r19885, r19886) <= 0, MPFR_RNDN); mpfr_exp(r19888, r19885, MPFR_RNDN); mpfr_mul(r19889, r19888, r19885, MPFR_RNDN); mpfr_set_d(r19890, x, MPFR_RNDN); mpfr_sub(r19891, r19889, r19890, MPFR_RNDN); mpfr_neg(r19892, r19891, MPFR_RNDN); ; mpfr_add(r19894, r19893, r19885, MPFR_RNDN); mpfr_mul(r19895, r19894, r19888, MPFR_RNDN); mpfr_div(r19896, r19892, r19895, MPFR_RNDN); mpfr_add(r19897, r19885, r19896, MPFR_RNDN); ; mpfr_set_si(r19899, mpfr_cmp(r19885, r19898) <= 0, MPFR_RNDN); mpfr_sqr(r19900, r19885, MPFR_RNDN); mpfr_add(r19901, r19900, r19890, MPFR_RNDN); ; mpfr_mul(r19903, r19885, r19890, MPFR_RNDN); mpfr_mul(r19904, r19902, r19903, MPFR_RNDN); mpfr_sub(r19905, r19901, r19904, MPFR_RNDN); mpfr_div(r19906, r19885, r19894, MPFR_RNDN); mpfr_sub(r19907, r19885, r19906, MPFR_RNDN); mpfr_mul(r19908, r19885, r19888, MPFR_RNDN); mpfr_add(r19909, r19888, r19908, MPFR_RNDN); mpfr_div(r19910, r19890, r19909, MPFR_RNDN); mpfr_add(r19911, r19907, r19910, MPFR_RNDN); if (mpfr_get_si(r19899, MPFR_RNDN)) { mpfr_set(r19912, r19905, MPFR_RNDN); } else { mpfr_set(r19912, r19911, MPFR_RNDN); }; if (mpfr_get_si(r19887, MPFR_RNDN)) { mpfr_set(r19913, r19897, MPFR_RNDN); } else { mpfr_set(r19913, r19912, MPFR_RNDN); }; return mpfr_get_d(r19913, MPFR_RNDN); } static mpfr_t r19914, r19915, r19916, r19917, r19918, r19919, r19920, r19921, r19922, r19923, r19924, r19925, r19926, r19927, r19928, r19929, r19930, r19931, r19932, r19933, r19934, r19935, r19936, r19937, r19938, r19939, r19940, r19941, r19942; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r19914); mpfr_init_set_str(r19915, "-2.63216061572775e-13", 10, MPFR_RNDN); mpfr_init(r19916); mpfr_init(r19917); mpfr_init(r19918); mpfr_init(r19919); mpfr_init(r19920); mpfr_init(r19921); mpfr_init_set_str(r19922, "1", 10, MPFR_RNDN); mpfr_init(r19923); mpfr_init(r19924); mpfr_init(r19925); mpfr_init(r19926); mpfr_init_set_str(r19927, "1.522317513311434e-06", 10, MPFR_RNDN); mpfr_init(r19928); mpfr_init(r19929); mpfr_init(r19930); mpfr_init_set_str(r19931, "2", 10, MPFR_RNDN); mpfr_init(r19932); mpfr_init(r19933); mpfr_init(r19934); mpfr_init(r19935); mpfr_init(r19936); mpfr_init(r19937); mpfr_init(r19938); mpfr_init(r19939); mpfr_init(r19940); mpfr_init(r19941); mpfr_init(r19942); } double f_dm(double wj, double x) { mpfr_set_d(r19914, wj, MPFR_RNDN); ; mpfr_set_si(r19916, mpfr_cmp(r19914, r19915) <= 0, MPFR_RNDN); mpfr_exp(r19917, r19914, MPFR_RNDN); mpfr_mul(r19918, r19917, r19914, MPFR_RNDN); mpfr_set_d(r19919, x, MPFR_RNDN); mpfr_sub(r19920, r19918, r19919, MPFR_RNDN); mpfr_neg(r19921, r19920, MPFR_RNDN); ; mpfr_add(r19923, r19922, r19914, MPFR_RNDN); mpfr_mul(r19924, r19923, r19917, MPFR_RNDN); mpfr_div(r19925, r19921, r19924, MPFR_RNDN); mpfr_add(r19926, r19914, r19925, MPFR_RNDN); ; mpfr_set_si(r19928, mpfr_cmp(r19914, r19927) <= 0, MPFR_RNDN); mpfr_sqr(r19929, r19914, MPFR_RNDN); mpfr_add(r19930, r19929, r19919, MPFR_RNDN); ; mpfr_mul(r19932, r19914, r19919, MPFR_RNDN); mpfr_mul(r19933, r19931, r19932, MPFR_RNDN); mpfr_sub(r19934, r19930, r19933, MPFR_RNDN); mpfr_div(r19935, r19914, r19923, MPFR_RNDN); mpfr_sub(r19936, r19914, r19935, MPFR_RNDN); mpfr_mul(r19937, r19914, r19917, MPFR_RNDN); mpfr_add(r19938, r19917, r19937, MPFR_RNDN); mpfr_div(r19939, r19919, r19938, MPFR_RNDN); mpfr_add(r19940, r19936, r19939, MPFR_RNDN); if (mpfr_get_si(r19928, MPFR_RNDN)) { mpfr_set(r19941, r19934, MPFR_RNDN); } else { mpfr_set(r19941, r19940, MPFR_RNDN); }; if (mpfr_get_si(r19916, MPFR_RNDN)) { mpfr_set(r19942, r19926, MPFR_RNDN); } else { mpfr_set(r19942, r19941, MPFR_RNDN); }; return mpfr_get_d(r19942, MPFR_RNDN); }