#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r20812 = 1; float r20813 = atan2(1.0, 0.0); float r20814 = 4; float r20815 = r20813 / r20814; float r20816 = r20812 / r20815; float r20817 = f; float r20818 = r20815 * r20817; float r20819 = exp(r20818); float r20820 = -r20818; float r20821 = exp(r20820); float r20822 = r20819 + r20821; float r20823 = r20819 - r20821; float r20824 = r20822 / r20823; float r20825 = log(r20824); float r20826 = r20816 * r20825; float r20827 = -r20826; return r20827; } double f_id(double f) { double r20828 = 1; double r20829 = atan2(1.0, 0.0); double r20830 = 4; double r20831 = r20829 / r20830; double r20832 = r20828 / r20831; double r20833 = f; double r20834 = r20831 * r20833; double r20835 = exp(r20834); double r20836 = -r20834; double r20837 = exp(r20836); double r20838 = r20835 + r20837; double r20839 = r20835 - r20837; double r20840 = r20838 / r20839; double r20841 = log(r20840); double r20842 = r20832 * r20841; double r20843 = -r20842; return r20843; } double f_of(float f) { float r20844 = 4; float r20845 = atan2(1.0, 0.0); float r20846 = r20844 / r20845; float r20847 = f; float r20848 = log(r20847); float r20849 = log(r20844); float r20850 = r20848 - r20849; float r20851 = r20846 * r20850; float r20852 = log(r20845); float r20853 = r20846 * r20852; float r20854 = r20851 + r20853; float r20855 = pow(r20847, r20844); float r20856 = 7/5760; float r20857 = r20845 * r20856; float r20858 = r20845 * r20857; float r20859 = 1/12; float r20860 = -r20859; float r20861 = r20847 * r20847; float r20862 = r20860 * r20861; float r20863 = fma(r20855, r20858, r20862); float r20864 = r20845 * r20863; float r20865 = r20854 + r20864; return r20865; } double f_od(double f) { double r20866 = 4; double r20867 = atan2(1.0, 0.0); double r20868 = r20866 / r20867; double r20869 = f; double r20870 = log(r20869); double r20871 = log(r20866); double r20872 = r20870 - r20871; double r20873 = r20868 * r20872; double r20874 = log(r20867); double r20875 = r20868 * r20874; double r20876 = r20873 + r20875; double r20877 = pow(r20869, r20866); double r20878 = 7/5760; double r20879 = r20867 * r20878; double r20880 = r20867 * r20879; double r20881 = 1/12; double r20882 = -r20881; double r20883 = r20869 * r20869; double r20884 = r20882 * r20883; double r20885 = fma(r20877, r20880, r20884); double r20886 = r20867 * r20885; double r20887 = r20876 + r20886; return r20887; } 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 r20888, r20889, r20890, r20891, r20892, r20893, r20894, r20895, r20896, r20897, r20898, r20899, r20900, r20901, r20902, r20903; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20888, "1", 10, MPFR_RNDN); mpfr_init(r20889); mpfr_init_set_str(r20890, "4", 10, MPFR_RNDN); mpfr_init(r20891); mpfr_init(r20892); mpfr_init(r20893); mpfr_init(r20894); mpfr_init(r20895); mpfr_init(r20896); mpfr_init(r20897); mpfr_init(r20898); mpfr_init(r20899); mpfr_init(r20900); mpfr_init(r20901); mpfr_init(r20902); mpfr_init(r20903); } double f_im(double f) { ; mpfr_const_pi(r20889, MPFR_RNDN); ; mpfr_div(r20891, r20889, r20890, MPFR_RNDN); mpfr_div(r20892, r20888, r20891, MPFR_RNDN); mpfr_set_d(r20893, f, MPFR_RNDN); mpfr_mul(r20894, r20891, r20893, MPFR_RNDN); mpfr_exp(r20895, r20894, MPFR_RNDN); mpfr_neg(r20896, r20894, MPFR_RNDN); mpfr_exp(r20897, r20896, MPFR_RNDN); mpfr_add(r20898, r20895, r20897, MPFR_RNDN); mpfr_sub(r20899, r20895, r20897, MPFR_RNDN); mpfr_div(r20900, r20898, r20899, MPFR_RNDN); mpfr_log(r20901, r20900, MPFR_RNDN); mpfr_mul(r20902, r20892, r20901, MPFR_RNDN); mpfr_neg(r20903, r20902, MPFR_RNDN); return mpfr_get_d(r20903, MPFR_RNDN); } static mpfr_t r20904, r20905, r20906, r20907, r20908, r20909, r20910, r20911, r20912, r20913, r20914, r20915, r20916, r20917, r20918, r20919, r20920, r20921, r20922, r20923, r20924, r20925; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20904, "4", 10, MPFR_RNDN); mpfr_init(r20905); mpfr_init(r20906); mpfr_init(r20907); mpfr_init(r20908); mpfr_init(r20909); mpfr_init(r20910); mpfr_init(r20911); mpfr_init(r20912); mpfr_init(r20913); mpfr_init(r20914); mpfr_init(r20915); mpfr_init_set_str(r20916, "7/5760", 10, MPFR_RNDN); mpfr_init(r20917); mpfr_init(r20918); mpfr_init_set_str(r20919, "1/12", 10, MPFR_RNDN); mpfr_init(r20920); mpfr_init(r20921); mpfr_init(r20922); mpfr_init(r20923); mpfr_init(r20924); mpfr_init(r20925); } double f_fm(double f) { ; mpfr_const_pi(r20905, MPFR_RNDN); mpfr_div(r20906, r20904, r20905, MPFR_RNDN); mpfr_set_d(r20907, f, MPFR_RNDN); mpfr_log(r20908, r20907, MPFR_RNDN); mpfr_log(r20909, r20904, MPFR_RNDN); mpfr_sub(r20910, r20908, r20909, MPFR_RNDN); mpfr_mul(r20911, r20906, r20910, MPFR_RNDN); mpfr_log(r20912, r20905, MPFR_RNDN); mpfr_mul(r20913, r20906, r20912, MPFR_RNDN); mpfr_add(r20914, r20911, r20913, MPFR_RNDN); mpfr_pow(r20915, r20907, r20904, MPFR_RNDN); ; mpfr_mul(r20917, r20905, r20916, MPFR_RNDN); mpfr_mul(r20918, r20905, r20917, MPFR_RNDN); ; mpfr_neg(r20920, r20919, MPFR_RNDN); mpfr_mul(r20921, r20907, r20907, MPFR_RNDN); mpfr_mul(r20922, r20920, r20921, MPFR_RNDN); mpfr_fma(r20923, r20915, r20918, r20922, MPFR_RNDN); mpfr_mul(r20924, r20905, r20923, MPFR_RNDN); mpfr_add(r20925, r20914, r20924, MPFR_RNDN); return mpfr_get_d(r20925, MPFR_RNDN); } static mpfr_t r20926, r20927, r20928, r20929, r20930, r20931, r20932, r20933, r20934, r20935, r20936, r20937, r20938, r20939, r20940, r20941, r20942, r20943, r20944, r20945, r20946, r20947; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20926, "4", 10, MPFR_RNDN); mpfr_init(r20927); mpfr_init(r20928); mpfr_init(r20929); mpfr_init(r20930); mpfr_init(r20931); mpfr_init(r20932); mpfr_init(r20933); mpfr_init(r20934); mpfr_init(r20935); mpfr_init(r20936); mpfr_init(r20937); mpfr_init_set_str(r20938, "7/5760", 10, MPFR_RNDN); mpfr_init(r20939); mpfr_init(r20940); mpfr_init_set_str(r20941, "1/12", 10, MPFR_RNDN); mpfr_init(r20942); mpfr_init(r20943); mpfr_init(r20944); mpfr_init(r20945); mpfr_init(r20946); mpfr_init(r20947); } double f_dm(double f) { ; mpfr_const_pi(r20927, MPFR_RNDN); mpfr_div(r20928, r20926, r20927, MPFR_RNDN); mpfr_set_d(r20929, f, MPFR_RNDN); mpfr_log(r20930, r20929, MPFR_RNDN); mpfr_log(r20931, r20926, MPFR_RNDN); mpfr_sub(r20932, r20930, r20931, MPFR_RNDN); mpfr_mul(r20933, r20928, r20932, MPFR_RNDN); mpfr_log(r20934, r20927, MPFR_RNDN); mpfr_mul(r20935, r20928, r20934, MPFR_RNDN); mpfr_add(r20936, r20933, r20935, MPFR_RNDN); mpfr_pow(r20937, r20929, r20926, MPFR_RNDN); ; mpfr_mul(r20939, r20927, r20938, MPFR_RNDN); mpfr_mul(r20940, r20927, r20939, MPFR_RNDN); ; mpfr_neg(r20942, r20941, MPFR_RNDN); mpfr_mul(r20943, r20929, r20929, MPFR_RNDN); mpfr_mul(r20944, r20942, r20943, MPFR_RNDN); mpfr_fma(r20945, r20937, r20940, r20944, MPFR_RNDN); mpfr_mul(r20946, r20927, r20945, MPFR_RNDN); mpfr_add(r20947, r20936, r20946, MPFR_RNDN); return mpfr_get_d(r20947, MPFR_RNDN); }