#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r25959 = N; float r25960 = 1; float r25961 = r25959 + r25960; float r25962 = log(r25961); float r25963 = log(r25959); float r25964 = r25962 - r25963; return r25964; } double f_id(double N) { double r25965 = N; double r25966 = 1; double r25967 = r25965 + r25966; double r25968 = log(r25967); double r25969 = log(r25965); double r25970 = r25968 - r25969; return r25970; } double f_of(float N) { float r25971 = N; float r25972 = 1; float r25973 = r25971 + r25972; float r25974 = log(r25973); float r25975 = log(r25971); float r25976 = r25974 - r25975; float r25977 = 1.0007174098329816e-09; bool r25978 = r25976 <= r25977; float r25979 = r25971 * r25971; float r25980 = r25972 / r25979; float r25981 = 1/3; float r25982 = r25981 / r25971; float r25983 = 1/2; float r25984 = r25982 - r25983; float r25985 = r25972 / r25971; float r25986 = fma(r25980, r25984, r25985); float r25987 = log1p(r25971); float r25988 = cbrt(r25987); float r25989 = r25988 * r25988; float r25990 = r25989 * r25988; float r25991 = r25990 - r25975; float r25992 = r25978 ? r25986 : r25991; return r25992; } double f_od(double N) { double r25993 = N; double r25994 = 1; double r25995 = r25993 + r25994; double r25996 = log(r25995); double r25997 = log(r25993); double r25998 = r25996 - r25997; double r25999 = 1.0007174098329816e-09; bool r26000 = r25998 <= r25999; double r26001 = r25993 * r25993; double r26002 = r25994 / r26001; double r26003 = 1/3; double r26004 = r26003 / r25993; double r26005 = 1/2; double r26006 = r26004 - r26005; double r26007 = r25994 / r25993; double r26008 = fma(r26002, r26006, r26007); double r26009 = log1p(r25993); double r26010 = cbrt(r26009); double r26011 = r26010 * r26010; double r26012 = r26011 * r26010; double r26013 = r26012 - r25997; double r26014 = r26000 ? r26008 : r26013; return r26014; } 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 r26015, r26016, r26017, r26018, r26019, r26020; void setup_mpfr_f_im() { mpfr_set_default_prec(1360); mpfr_init(r26015); mpfr_init_set_str(r26016, "1", 10, MPFR_RNDN); mpfr_init(r26017); mpfr_init(r26018); mpfr_init(r26019); mpfr_init(r26020); } double f_im(double N) { mpfr_set_d(r26015, N, MPFR_RNDN); ; mpfr_add(r26017, r26015, r26016, MPFR_RNDN); mpfr_log(r26018, r26017, MPFR_RNDN); mpfr_log(r26019, r26015, MPFR_RNDN); mpfr_sub(r26020, r26018, r26019, MPFR_RNDN); return mpfr_get_d(r26020, MPFR_RNDN); } static mpfr_t r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042; void setup_mpfr_f_fm() { mpfr_set_default_prec(1360); mpfr_init(r26021); mpfr_init_set_str(r26022, "1", 10, MPFR_RNDN); mpfr_init(r26023); mpfr_init(r26024); mpfr_init(r26025); mpfr_init(r26026); mpfr_init_set_str(r26027, "1.0007174098329816e-09", 10, MPFR_RNDN); mpfr_init(r26028); mpfr_init(r26029); mpfr_init(r26030); mpfr_init_set_str(r26031, "1/3", 10, MPFR_RNDN); mpfr_init(r26032); mpfr_init_set_str(r26033, "1/2", 10, MPFR_RNDN); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init(r26038); mpfr_init(r26039); mpfr_init(r26040); mpfr_init(r26041); mpfr_init(r26042); } double f_fm(double N) { mpfr_set_d(r26021, N, MPFR_RNDN); ; mpfr_add(r26023, r26021, r26022, MPFR_RNDN); mpfr_log(r26024, r26023, MPFR_RNDN); mpfr_log(r26025, r26021, MPFR_RNDN); mpfr_sub(r26026, r26024, r26025, MPFR_RNDN); ; mpfr_set_si(r26028, mpfr_cmp(r26026, r26027) <= 0, MPFR_RNDN); mpfr_mul(r26029, r26021, r26021, MPFR_RNDN); mpfr_div(r26030, r26022, r26029, MPFR_RNDN); ; mpfr_div(r26032, r26031, r26021, MPFR_RNDN); ; mpfr_sub(r26034, r26032, r26033, MPFR_RNDN); mpfr_div(r26035, r26022, r26021, MPFR_RNDN); mpfr_fma(r26036, r26030, r26034, r26035, MPFR_RNDN); mpfr_log1p(r26037, r26021, MPFR_RNDN); mpfr_cbrt(r26038, r26037, MPFR_RNDN); mpfr_mul(r26039, r26038, r26038, MPFR_RNDN); mpfr_mul(r26040, r26039, r26038, MPFR_RNDN); mpfr_sub(r26041, r26040, r26025, MPFR_RNDN); if (mpfr_get_si(r26028, MPFR_RNDN)) { mpfr_set(r26042, r26036, MPFR_RNDN); } else { mpfr_set(r26042, r26041, MPFR_RNDN); }; return mpfr_get_d(r26042, MPFR_RNDN); } static mpfr_t r26043, r26044, r26045, r26046, r26047, r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064; void setup_mpfr_f_dm() { mpfr_set_default_prec(1360); mpfr_init(r26043); mpfr_init_set_str(r26044, "1", 10, MPFR_RNDN); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); mpfr_init(r26048); mpfr_init_set_str(r26049, "1.0007174098329816e-09", 10, MPFR_RNDN); mpfr_init(r26050); mpfr_init(r26051); mpfr_init(r26052); mpfr_init_set_str(r26053, "1/3", 10, MPFR_RNDN); mpfr_init(r26054); mpfr_init_set_str(r26055, "1/2", 10, MPFR_RNDN); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); } double f_dm(double N) { mpfr_set_d(r26043, N, MPFR_RNDN); ; mpfr_add(r26045, r26043, r26044, MPFR_RNDN); mpfr_log(r26046, r26045, MPFR_RNDN); mpfr_log(r26047, r26043, MPFR_RNDN); mpfr_sub(r26048, r26046, r26047, MPFR_RNDN); ; mpfr_set_si(r26050, mpfr_cmp(r26048, r26049) <= 0, MPFR_RNDN); mpfr_mul(r26051, r26043, r26043, MPFR_RNDN); mpfr_div(r26052, r26044, r26051, MPFR_RNDN); ; mpfr_div(r26054, r26053, r26043, MPFR_RNDN); ; mpfr_sub(r26056, r26054, r26055, MPFR_RNDN); mpfr_div(r26057, r26044, r26043, MPFR_RNDN); mpfr_fma(r26058, r26052, r26056, r26057, MPFR_RNDN); mpfr_log1p(r26059, r26043, MPFR_RNDN); mpfr_cbrt(r26060, r26059, MPFR_RNDN); mpfr_mul(r26061, r26060, r26060, MPFR_RNDN); mpfr_mul(r26062, r26061, r26060, MPFR_RNDN); mpfr_sub(r26063, r26062, r26047, MPFR_RNDN); if (mpfr_get_si(r26050, MPFR_RNDN)) { mpfr_set(r26064, r26058, MPFR_RNDN); } else { mpfr_set(r26064, r26063, MPFR_RNDN); }; return mpfr_get_d(r26064, MPFR_RNDN); }