#include #include #include #include #include char *name = "VandenBroeck and Keller, Equation (20)"; double f_if(float f) { float r19907 = 1; float r19908 = atan2(1.0, 0.0); float r19909 = 4; float r19910 = r19908 / r19909; float r19911 = r19907 / r19910; float r19912 = f; float r19913 = r19910 * r19912; float r19914 = exp(r19913); float r19915 = -r19913; float r19916 = exp(r19915); float r19917 = r19914 + r19916; float r19918 = r19914 - r19916; float r19919 = r19917 / r19918; float r19920 = log(r19919); float r19921 = r19911 * r19920; float r19922 = -r19921; return r19922; } double f_id(double f) { double r19923 = 1; double r19924 = atan2(1.0, 0.0); double r19925 = 4; double r19926 = r19924 / r19925; double r19927 = r19923 / r19926; double r19928 = f; double r19929 = r19926 * r19928; double r19930 = exp(r19929); double r19931 = -r19929; double r19932 = exp(r19931); double r19933 = r19930 + r19932; double r19934 = r19930 - r19932; double r19935 = r19933 / r19934; double r19936 = log(r19935); double r19937 = r19927 * r19936; double r19938 = -r19937; return r19938; } double f_of(float f) { float r19939 = 4; float r19940 = r19939 + r19939; float r19941 = atan2(1.0, 0.0); float r19942 = r19940 / r19941; float r19943 = r19939 / r19941; float r19944 = cbrt(r19943); float r19945 = log(r19944); float r19946 = r19942 * r19945; float r19947 = f; float r19948 = log(r19947); float r19949 = r19945 - r19948; float r19950 = r19943 * r19949; float r19951 = r19946 + r19950; float r19952 = r19947 * r19947; float r19953 = 1/12; float r19954 = r19953 * r19941; float r19955 = r19952 * r19954; float r19956 = 7/5760; float r19957 = r19956 * r19941; float r19958 = pow(r19947, r19939); float r19959 = r19957 * r19958; float r19960 = r19941 * r19941; float r19961 = r19959 * r19960; float r19962 = r19955 - r19961; float r19963 = r19951 + r19962; float r19964 = -r19963; return r19964; } double f_od(double f) { double r19965 = 4; double r19966 = r19965 + r19965; double r19967 = atan2(1.0, 0.0); double r19968 = r19966 / r19967; double r19969 = r19965 / r19967; double r19970 = cbrt(r19969); double r19971 = log(r19970); double r19972 = r19968 * r19971; double r19973 = f; double r19974 = log(r19973); double r19975 = r19971 - r19974; double r19976 = r19969 * r19975; double r19977 = r19972 + r19976; double r19978 = r19973 * r19973; double r19979 = 1/12; double r19980 = r19979 * r19967; double r19981 = r19978 * r19980; double r19982 = 7/5760; double r19983 = r19982 * r19967; double r19984 = pow(r19973, r19965); double r19985 = r19983 * r19984; double r19986 = r19967 * r19967; double r19987 = r19985 * r19986; double r19988 = r19981 - r19987; double r19989 = r19977 + r19988; double r19990 = -r19989; return r19990; } 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 r19991, r19992, r19993, r19994, r19995, r19996, r19997, r19998, r19999, r20000, r20001, r20002, r20003, r20004, r20005, r20006; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init_set_str(r19991, "1", 10, MPFR_RNDN); mpfr_init(r19992); mpfr_init_set_str(r19993, "4", 10, MPFR_RNDN); mpfr_init(r19994); mpfr_init(r19995); mpfr_init(r19996); mpfr_init(r19997); mpfr_init(r19998); mpfr_init(r19999); mpfr_init(r20000); mpfr_init(r20001); mpfr_init(r20002); mpfr_init(r20003); mpfr_init(r20004); mpfr_init(r20005); mpfr_init(r20006); } double f_im(double f) { ; mpfr_const_pi(r19992, MPFR_RNDN); ; mpfr_div(r19994, r19992, r19993, MPFR_RNDN); mpfr_div(r19995, r19991, r19994, MPFR_RNDN); mpfr_set_d(r19996, f, MPFR_RNDN); mpfr_mul(r19997, r19994, r19996, MPFR_RNDN); mpfr_exp(r19998, r19997, MPFR_RNDN); mpfr_neg(r19999, r19997, MPFR_RNDN); mpfr_exp(r20000, r19999, MPFR_RNDN); mpfr_add(r20001, r19998, r20000, MPFR_RNDN); mpfr_sub(r20002, r19998, r20000, MPFR_RNDN); mpfr_div(r20003, r20001, r20002, MPFR_RNDN); mpfr_log(r20004, r20003, MPFR_RNDN); mpfr_mul(r20005, r19995, r20004, MPFR_RNDN); mpfr_neg(r20006, r20005, MPFR_RNDN); return mpfr_get_d(r20006, MPFR_RNDN); } static mpfr_t r20007, r20008, r20009, r20010, r20011, r20012, r20013, r20014, r20015, r20016, r20017, r20018, r20019, r20020, r20021, r20022, r20023, r20024, r20025, r20026, r20027, r20028, r20029, r20030, r20031, r20032; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20007, "4", 10, MPFR_RNDN); mpfr_init(r20008); mpfr_init(r20009); mpfr_init(r20010); mpfr_init(r20011); mpfr_init(r20012); mpfr_init(r20013); mpfr_init(r20014); mpfr_init(r20015); mpfr_init(r20016); mpfr_init(r20017); mpfr_init(r20018); mpfr_init(r20019); mpfr_init(r20020); mpfr_init_set_str(r20021, "1/12", 10, MPFR_RNDN); mpfr_init(r20022); mpfr_init(r20023); mpfr_init_set_str(r20024, "7/5760", 10, MPFR_RNDN); mpfr_init(r20025); mpfr_init(r20026); mpfr_init(r20027); mpfr_init(r20028); mpfr_init(r20029); mpfr_init(r20030); mpfr_init(r20031); mpfr_init(r20032); } double f_fm(double f) { ; mpfr_add(r20008, r20007, r20007, MPFR_RNDN); mpfr_const_pi(r20009, MPFR_RNDN); mpfr_div(r20010, r20008, r20009, MPFR_RNDN); mpfr_div(r20011, r20007, r20009, MPFR_RNDN); mpfr_cbrt(r20012, r20011, MPFR_RNDN); mpfr_log(r20013, r20012, MPFR_RNDN); mpfr_mul(r20014, r20010, r20013, MPFR_RNDN); mpfr_set_d(r20015, f, MPFR_RNDN); mpfr_log(r20016, r20015, MPFR_RNDN); mpfr_sub(r20017, r20013, r20016, MPFR_RNDN); mpfr_mul(r20018, r20011, r20017, MPFR_RNDN); mpfr_add(r20019, r20014, r20018, MPFR_RNDN); mpfr_mul(r20020, r20015, r20015, MPFR_RNDN); ; mpfr_mul(r20022, r20021, r20009, MPFR_RNDN); mpfr_mul(r20023, r20020, r20022, MPFR_RNDN); ; mpfr_mul(r20025, r20024, r20009, MPFR_RNDN); mpfr_pow(r20026, r20015, r20007, MPFR_RNDN); mpfr_mul(r20027, r20025, r20026, MPFR_RNDN); mpfr_mul(r20028, r20009, r20009, MPFR_RNDN); mpfr_mul(r20029, r20027, r20028, MPFR_RNDN); mpfr_sub(r20030, r20023, r20029, MPFR_RNDN); mpfr_add(r20031, r20019, r20030, MPFR_RNDN); mpfr_neg(r20032, r20031, MPFR_RNDN); return mpfr_get_d(r20032, MPFR_RNDN); } static mpfr_t r20033, r20034, r20035, r20036, r20037, r20038, r20039, r20040, r20041, r20042, r20043, r20044, r20045, r20046, r20047, r20048, r20049, r20050, r20051, r20052, r20053, r20054, r20055, r20056, r20057, r20058; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init_set_str(r20033, "4", 10, MPFR_RNDN); mpfr_init(r20034); mpfr_init(r20035); mpfr_init(r20036); mpfr_init(r20037); mpfr_init(r20038); mpfr_init(r20039); mpfr_init(r20040); mpfr_init(r20041); mpfr_init(r20042); mpfr_init(r20043); mpfr_init(r20044); mpfr_init(r20045); mpfr_init(r20046); mpfr_init_set_str(r20047, "1/12", 10, MPFR_RNDN); mpfr_init(r20048); mpfr_init(r20049); mpfr_init_set_str(r20050, "7/5760", 10, MPFR_RNDN); mpfr_init(r20051); mpfr_init(r20052); mpfr_init(r20053); mpfr_init(r20054); mpfr_init(r20055); mpfr_init(r20056); mpfr_init(r20057); mpfr_init(r20058); } double f_dm(double f) { ; mpfr_add(r20034, r20033, r20033, MPFR_RNDN); mpfr_const_pi(r20035, MPFR_RNDN); mpfr_div(r20036, r20034, r20035, MPFR_RNDN); mpfr_div(r20037, r20033, r20035, MPFR_RNDN); mpfr_cbrt(r20038, r20037, MPFR_RNDN); mpfr_log(r20039, r20038, MPFR_RNDN); mpfr_mul(r20040, r20036, r20039, MPFR_RNDN); mpfr_set_d(r20041, f, MPFR_RNDN); mpfr_log(r20042, r20041, MPFR_RNDN); mpfr_sub(r20043, r20039, r20042, MPFR_RNDN); mpfr_mul(r20044, r20037, r20043, MPFR_RNDN); mpfr_add(r20045, r20040, r20044, MPFR_RNDN); mpfr_mul(r20046, r20041, r20041, MPFR_RNDN); ; mpfr_mul(r20048, r20047, r20035, MPFR_RNDN); mpfr_mul(r20049, r20046, r20048, MPFR_RNDN); ; mpfr_mul(r20051, r20050, r20035, MPFR_RNDN); mpfr_pow(r20052, r20041, r20033, MPFR_RNDN); mpfr_mul(r20053, r20051, r20052, MPFR_RNDN); mpfr_mul(r20054, r20035, r20035, MPFR_RNDN); mpfr_mul(r20055, r20053, r20054, MPFR_RNDN); mpfr_sub(r20056, r20049, r20055, MPFR_RNDN); mpfr_add(r20057, r20045, r20056, MPFR_RNDN); mpfr_neg(r20058, r20057, MPFR_RNDN); return mpfr_get_d(r20058, MPFR_RNDN); }