#include #include #include #include #include char *name = "Jmat.Real.lambertw, newton loop step"; double f_if(float wj, float x) { float r24775 = wj; float r24776 = exp(r24775); float r24777 = r24775 * r24776; float r24778 = x; float r24779 = r24777 - r24778; float r24780 = r24776 + r24777; float r24781 = r24779 / r24780; float r24782 = r24775 - r24781; return r24782; } double f_id(double wj, double x) { double r24783 = wj; double r24784 = exp(r24783); double r24785 = r24783 * r24784; double r24786 = x; double r24787 = r24785 - r24786; double r24788 = r24784 + r24785; double r24789 = r24787 / r24788; double r24790 = r24783 - r24789; return r24790; } double f_of(float wj, float x) { float r24791 = wj; float r24792 = 1; float r24793 = r24792 + r24791; float r24794 = r24791 / r24793; float r24795 = r24791 - r24794; float r24796 = x; float r24797 = exp(r24791); float r24798 = r24791 * r24797; float r24799 = r24797 + r24798; float r24800 = r24796 / r24799; float r24801 = r24795 + r24800; float r24802 = 3.639973196389132e-18; bool r24803 = r24801 <= r24802; float r24804 = 2; float r24805 = pow(r24791, r24804); float r24806 = r24805 + r24796; float r24807 = r24791 * r24796; float r24808 = r24804 * r24807; float r24809 = r24806 - r24808; float r24810 = r24803 ? r24809 : r24801; return r24810; } double f_od(double wj, double x) { double r24811 = wj; double r24812 = 1; double r24813 = r24812 + r24811; double r24814 = r24811 / r24813; double r24815 = r24811 - r24814; double r24816 = x; double r24817 = exp(r24811); double r24818 = r24811 * r24817; double r24819 = r24817 + r24818; double r24820 = r24816 / r24819; double r24821 = r24815 + r24820; double r24822 = 3.639973196389132e-18; bool r24823 = r24821 <= r24822; double r24824 = 2; double r24825 = pow(r24811, r24824); double r24826 = r24825 + r24816; double r24827 = r24811 * r24816; double r24828 = r24824 * r24827; double r24829 = r24826 - r24828; double r24830 = r24823 ? r24829 : r24821; return r24830; } 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 r24831, r24832, r24833, r24834, r24835, r24836, r24837, r24838; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r24831); mpfr_init(r24832); mpfr_init(r24833); mpfr_init(r24834); mpfr_init(r24835); mpfr_init(r24836); mpfr_init(r24837); mpfr_init(r24838); } double f_im(double wj, double x) { mpfr_set_d(r24831, wj, MPFR_RNDN); mpfr_exp(r24832, r24831, MPFR_RNDN); mpfr_mul(r24833, r24831, r24832, MPFR_RNDN); mpfr_set_d(r24834, x, MPFR_RNDN); mpfr_sub(r24835, r24833, r24834, MPFR_RNDN); mpfr_add(r24836, r24832, r24833, MPFR_RNDN); mpfr_div(r24837, r24835, r24836, MPFR_RNDN); mpfr_sub(r24838, r24831, r24837, MPFR_RNDN); return mpfr_get_d(r24838, MPFR_RNDN); } static mpfr_t r24839, r24840, r24841, r24842, r24843, r24844, r24845, r24846, r24847, r24848, r24849, r24850, r24851, r24852, r24853, r24854, r24855, r24856, r24857, r24858; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r24839); mpfr_init_set_str(r24840, "1", 10, MPFR_RNDN); mpfr_init(r24841); mpfr_init(r24842); mpfr_init(r24843); mpfr_init(r24844); mpfr_init(r24845); mpfr_init(r24846); mpfr_init(r24847); mpfr_init(r24848); mpfr_init(r24849); mpfr_init_set_str(r24850, "3.639973196389132e-18", 10, MPFR_RNDN); mpfr_init(r24851); mpfr_init_set_str(r24852, "2", 10, MPFR_RNDN); mpfr_init(r24853); mpfr_init(r24854); mpfr_init(r24855); mpfr_init(r24856); mpfr_init(r24857); mpfr_init(r24858); } double f_fm(double wj, double x) { mpfr_set_d(r24839, wj, MPFR_RNDN); ; mpfr_add(r24841, r24840, r24839, MPFR_RNDN); mpfr_div(r24842, r24839, r24841, MPFR_RNDN); mpfr_sub(r24843, r24839, r24842, MPFR_RNDN); mpfr_set_d(r24844, x, MPFR_RNDN); mpfr_exp(r24845, r24839, MPFR_RNDN); mpfr_mul(r24846, r24839, r24845, MPFR_RNDN); mpfr_add(r24847, r24845, r24846, MPFR_RNDN); mpfr_div(r24848, r24844, r24847, MPFR_RNDN); mpfr_add(r24849, r24843, r24848, MPFR_RNDN); ; mpfr_set_si(r24851, mpfr_cmp(r24849, r24850) <= 0, MPFR_RNDN); ; mpfr_pow(r24853, r24839, r24852, MPFR_RNDN); mpfr_add(r24854, r24853, r24844, MPFR_RNDN); mpfr_mul(r24855, r24839, r24844, MPFR_RNDN); mpfr_mul(r24856, r24852, r24855, MPFR_RNDN); mpfr_sub(r24857, r24854, r24856, MPFR_RNDN); if (mpfr_get_si(r24851, MPFR_RNDN)) { mpfr_set(r24858, r24857, MPFR_RNDN); } else { mpfr_set(r24858, r24849, MPFR_RNDN); }; return mpfr_get_d(r24858, MPFR_RNDN); } static mpfr_t r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868, r24869, r24870, r24871, r24872, r24873, r24874, r24875, r24876, r24877, r24878; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r24859); mpfr_init_set_str(r24860, "1", 10, MPFR_RNDN); mpfr_init(r24861); mpfr_init(r24862); mpfr_init(r24863); mpfr_init(r24864); mpfr_init(r24865); mpfr_init(r24866); mpfr_init(r24867); mpfr_init(r24868); mpfr_init(r24869); mpfr_init_set_str(r24870, "3.639973196389132e-18", 10, MPFR_RNDN); mpfr_init(r24871); mpfr_init_set_str(r24872, "2", 10, MPFR_RNDN); mpfr_init(r24873); mpfr_init(r24874); mpfr_init(r24875); mpfr_init(r24876); mpfr_init(r24877); mpfr_init(r24878); } double f_dm(double wj, double x) { mpfr_set_d(r24859, wj, MPFR_RNDN); ; mpfr_add(r24861, r24860, r24859, MPFR_RNDN); mpfr_div(r24862, r24859, r24861, MPFR_RNDN); mpfr_sub(r24863, r24859, r24862, MPFR_RNDN); mpfr_set_d(r24864, x, MPFR_RNDN); mpfr_exp(r24865, r24859, MPFR_RNDN); mpfr_mul(r24866, r24859, r24865, MPFR_RNDN); mpfr_add(r24867, r24865, r24866, MPFR_RNDN); mpfr_div(r24868, r24864, r24867, MPFR_RNDN); mpfr_add(r24869, r24863, r24868, MPFR_RNDN); ; mpfr_set_si(r24871, mpfr_cmp(r24869, r24870) <= 0, MPFR_RNDN); ; mpfr_pow(r24873, r24859, r24872, MPFR_RNDN); mpfr_add(r24874, r24873, r24864, MPFR_RNDN); mpfr_mul(r24875, r24859, r24864, MPFR_RNDN); mpfr_mul(r24876, r24872, r24875, MPFR_RNDN); mpfr_sub(r24877, r24874, r24876, MPFR_RNDN); if (mpfr_get_si(r24871, MPFR_RNDN)) { mpfr_set(r24878, r24877, MPFR_RNDN); } else { mpfr_set(r24878, r24869, MPFR_RNDN); }; return mpfr_get_d(r24878, MPFR_RNDN); }