#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r22901 = N; float r22902 = 1; float r22903 = r22901 + r22902; float r22904 = log(r22903); float r22905 = log(r22901); float r22906 = r22904 - r22905; return r22906; } double f_id(double N) { double r22907 = N; double r22908 = 1; double r22909 = r22907 + r22908; double r22910 = log(r22909); double r22911 = log(r22907); double r22912 = r22910 - r22911; return r22912; } double f_of(float N) { float r22913 = N; float r22914 = 1; float r22915 = r22913 + r22914; float r22916 = log(r22915); float r22917 = log(r22913); float r22918 = r22916 - r22917; float r22919 = 9.045593907618368e-32; bool r22920 = r22918 <= r22919; float r22921 = r22913 * r22913; float r22922 = r22914 / r22921; float r22923 = 1/3; float r22924 = r22923 / r22913; float r22925 = 1/2; float r22926 = r22924 - r22925; float r22927 = r22914 / r22913; float r22928 = fma(r22922, r22926, r22927); float r22929 = log1p(r22913); float r22930 = exp(r22929); float r22931 = r22930 / r22913; float r22932 = log(r22931); float r22933 = r22920 ? r22928 : r22932; return r22933; } double f_od(double N) { double r22934 = N; double r22935 = 1; double r22936 = r22934 + r22935; double r22937 = log(r22936); double r22938 = log(r22934); double r22939 = r22937 - r22938; double r22940 = 9.045593907618368e-32; bool r22941 = r22939 <= r22940; double r22942 = r22934 * r22934; double r22943 = r22935 / r22942; double r22944 = 1/3; double r22945 = r22944 / r22934; double r22946 = 1/2; double r22947 = r22945 - r22946; double r22948 = r22935 / r22934; double r22949 = fma(r22943, r22947, r22948); double r22950 = log1p(r22934); double r22951 = exp(r22950); double r22952 = r22951 / r22934; double r22953 = log(r22952); double r22954 = r22941 ? r22949 : r22953; return r22954; } 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 r22955, r22956, r22957, r22958, r22959, r22960; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r22955); mpfr_init_set_str(r22956, "1", 10, MPFR_RNDN); mpfr_init(r22957); mpfr_init(r22958); mpfr_init(r22959); mpfr_init(r22960); } double f_im(double N) { mpfr_set_d(r22955, N, MPFR_RNDN); ; mpfr_add(r22957, r22955, r22956, MPFR_RNDN); mpfr_log(r22958, r22957, MPFR_RNDN); mpfr_log(r22959, r22955, MPFR_RNDN); mpfr_sub(r22960, r22958, r22959, MPFR_RNDN); return mpfr_get_d(r22960, MPFR_RNDN); } static mpfr_t r22961, r22962, r22963, r22964, r22965, r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980, r22981; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r22961); mpfr_init_set_str(r22962, "1", 10, MPFR_RNDN); mpfr_init(r22963); mpfr_init(r22964); mpfr_init(r22965); mpfr_init(r22966); mpfr_init_set_str(r22967, "9.045593907618368e-32", 10, MPFR_RNDN); mpfr_init(r22968); mpfr_init(r22969); mpfr_init(r22970); mpfr_init_set_str(r22971, "1/3", 10, MPFR_RNDN); mpfr_init(r22972); mpfr_init_set_str(r22973, "1/2", 10, MPFR_RNDN); mpfr_init(r22974); mpfr_init(r22975); mpfr_init(r22976); mpfr_init(r22977); mpfr_init(r22978); mpfr_init(r22979); mpfr_init(r22980); mpfr_init(r22981); } double f_fm(double N) { mpfr_set_d(r22961, N, MPFR_RNDN); ; mpfr_add(r22963, r22961, r22962, MPFR_RNDN); mpfr_log(r22964, r22963, MPFR_RNDN); mpfr_log(r22965, r22961, MPFR_RNDN); mpfr_sub(r22966, r22964, r22965, MPFR_RNDN); ; mpfr_set_si(r22968, mpfr_cmp(r22966, r22967) <= 0, MPFR_RNDN); mpfr_mul(r22969, r22961, r22961, MPFR_RNDN); mpfr_div(r22970, r22962, r22969, MPFR_RNDN); ; mpfr_div(r22972, r22971, r22961, MPFR_RNDN); ; mpfr_sub(r22974, r22972, r22973, MPFR_RNDN); mpfr_div(r22975, r22962, r22961, MPFR_RNDN); mpfr_fma(r22976, r22970, r22974, r22975, MPFR_RNDN); mpfr_log1p(r22977, r22961, MPFR_RNDN); mpfr_exp(r22978, r22977, MPFR_RNDN); mpfr_div(r22979, r22978, r22961, MPFR_RNDN); mpfr_log(r22980, r22979, MPFR_RNDN); if (mpfr_get_si(r22968, MPFR_RNDN)) { mpfr_set(r22981, r22976, MPFR_RNDN); } else { mpfr_set(r22981, r22980, MPFR_RNDN); }; return mpfr_get_d(r22981, MPFR_RNDN); } static mpfr_t r22982, r22983, r22984, r22985, r22986, r22987, r22988, r22989, r22990, r22991, r22992, r22993, r22994, r22995, r22996, r22997, r22998, r22999, r23000, r23001, r23002; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r22982); mpfr_init_set_str(r22983, "1", 10, MPFR_RNDN); mpfr_init(r22984); mpfr_init(r22985); mpfr_init(r22986); mpfr_init(r22987); mpfr_init_set_str(r22988, "9.045593907618368e-32", 10, MPFR_RNDN); mpfr_init(r22989); mpfr_init(r22990); mpfr_init(r22991); mpfr_init_set_str(r22992, "1/3", 10, MPFR_RNDN); mpfr_init(r22993); mpfr_init_set_str(r22994, "1/2", 10, MPFR_RNDN); mpfr_init(r22995); mpfr_init(r22996); mpfr_init(r22997); mpfr_init(r22998); mpfr_init(r22999); mpfr_init(r23000); mpfr_init(r23001); mpfr_init(r23002); } double f_dm(double N) { mpfr_set_d(r22982, N, MPFR_RNDN); ; mpfr_add(r22984, r22982, r22983, MPFR_RNDN); mpfr_log(r22985, r22984, MPFR_RNDN); mpfr_log(r22986, r22982, MPFR_RNDN); mpfr_sub(r22987, r22985, r22986, MPFR_RNDN); ; mpfr_set_si(r22989, mpfr_cmp(r22987, r22988) <= 0, MPFR_RNDN); mpfr_mul(r22990, r22982, r22982, MPFR_RNDN); mpfr_div(r22991, r22983, r22990, MPFR_RNDN); ; mpfr_div(r22993, r22992, r22982, MPFR_RNDN); ; mpfr_sub(r22995, r22993, r22994, MPFR_RNDN); mpfr_div(r22996, r22983, r22982, MPFR_RNDN); mpfr_fma(r22997, r22991, r22995, r22996, MPFR_RNDN); mpfr_log1p(r22998, r22982, MPFR_RNDN); mpfr_exp(r22999, r22998, MPFR_RNDN); mpfr_div(r23000, r22999, r22982, MPFR_RNDN); mpfr_log(r23001, r23000, MPFR_RNDN); if (mpfr_get_si(r22989, MPFR_RNDN)) { mpfr_set(r23002, r22997, MPFR_RNDN); } else { mpfr_set(r23002, r23001, MPFR_RNDN); }; return mpfr_get_d(r23002, MPFR_RNDN); }