#include #include #include #include #include char *name = "2log (problem 3.3.6)"; double f_if(float N) { float r22982 = N; float r22983 = 1; float r22984 = r22982 + r22983; float r22985 = log(r22984); float r22986 = log(r22982); float r22987 = r22985 - r22986; return r22987; } double f_id(double N) { double r22988 = N; double r22989 = 1; double r22990 = r22988 + r22989; double r22991 = log(r22990); double r22992 = log(r22988); double r22993 = r22991 - r22992; return r22993; } double f_of(float N) { float r22994 = N; float r22995 = 1; float r22996 = r22994 + r22995; float r22997 = log(r22996); float r22998 = log(r22994); float r22999 = r22997 - r22998; float r23000 = 4.4676756717265675e-05; bool r23001 = r22999 <= r23000; float r23002 = r22994 * r22994; float r23003 = r22995 / r23002; float r23004 = 1/3; float r23005 = r23004 / r22994; float r23006 = 1/2; float r23007 = r23005 - r23006; float r23008 = r22995 / r22994; float r23009 = fma(r23003, r23007, r23008); float r23010 = r22996 / r22994; float r23011 = log(r23010); float r23012 = r23001 ? r23009 : r23011; return r23012; } double f_od(double N) { double r23013 = N; double r23014 = 1; double r23015 = r23013 + r23014; double r23016 = log(r23015); double r23017 = log(r23013); double r23018 = r23016 - r23017; double r23019 = 4.4676756717265675e-05; bool r23020 = r23018 <= r23019; double r23021 = r23013 * r23013; double r23022 = r23014 / r23021; double r23023 = 1/3; double r23024 = r23023 / r23013; double r23025 = 1/2; double r23026 = r23024 - r23025; double r23027 = r23014 / r23013; double r23028 = fma(r23022, r23026, r23027); double r23029 = r23015 / r23013; double r23030 = log(r23029); double r23031 = r23020 ? r23028 : r23030; return r23031; } 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 r23032, r23033, r23034, r23035, r23036, r23037; void setup_mpfr_f_im() { mpfr_set_default_prec(1616); mpfr_init(r23032); mpfr_init_set_str(r23033, "1", 10, MPFR_RNDN); mpfr_init(r23034); mpfr_init(r23035); mpfr_init(r23036); mpfr_init(r23037); } double f_im(double N) { mpfr_set_d(r23032, N, MPFR_RNDN); ; mpfr_add(r23034, r23032, r23033, MPFR_RNDN); mpfr_log(r23035, r23034, MPFR_RNDN); mpfr_log(r23036, r23032, MPFR_RNDN); mpfr_sub(r23037, r23035, r23036, MPFR_RNDN); return mpfr_get_d(r23037, MPFR_RNDN); } static mpfr_t r23038, r23039, r23040, r23041, r23042, r23043, r23044, r23045, r23046, r23047, r23048, r23049, r23050, r23051, r23052, r23053, r23054, r23055, r23056; void setup_mpfr_f_fm() { mpfr_set_default_prec(1616); mpfr_init(r23038); mpfr_init_set_str(r23039, "1", 10, MPFR_RNDN); mpfr_init(r23040); mpfr_init(r23041); mpfr_init(r23042); mpfr_init(r23043); mpfr_init_set_str(r23044, "4.4676756717265675e-05", 10, MPFR_RNDN); mpfr_init(r23045); mpfr_init(r23046); mpfr_init(r23047); mpfr_init_set_str(r23048, "1/3", 10, MPFR_RNDN); mpfr_init(r23049); mpfr_init_set_str(r23050, "1/2", 10, MPFR_RNDN); mpfr_init(r23051); mpfr_init(r23052); mpfr_init(r23053); mpfr_init(r23054); mpfr_init(r23055); mpfr_init(r23056); } double f_fm(double N) { mpfr_set_d(r23038, N, MPFR_RNDN); ; mpfr_add(r23040, r23038, r23039, MPFR_RNDN); mpfr_log(r23041, r23040, MPFR_RNDN); mpfr_log(r23042, r23038, MPFR_RNDN); mpfr_sub(r23043, r23041, r23042, MPFR_RNDN); ; mpfr_set_si(r23045, mpfr_cmp(r23043, r23044) <= 0, MPFR_RNDN); mpfr_mul(r23046, r23038, r23038, MPFR_RNDN); mpfr_div(r23047, r23039, r23046, MPFR_RNDN); ; mpfr_div(r23049, r23048, r23038, MPFR_RNDN); ; mpfr_sub(r23051, r23049, r23050, MPFR_RNDN); mpfr_div(r23052, r23039, r23038, MPFR_RNDN); mpfr_fma(r23053, r23047, r23051, r23052, MPFR_RNDN); mpfr_div(r23054, r23040, r23038, MPFR_RNDN); mpfr_log(r23055, r23054, MPFR_RNDN); if (mpfr_get_si(r23045, MPFR_RNDN)) { mpfr_set(r23056, r23053, MPFR_RNDN); } else { mpfr_set(r23056, r23055, MPFR_RNDN); }; return mpfr_get_d(r23056, MPFR_RNDN); } static mpfr_t r23057, r23058, r23059, r23060, r23061, r23062, r23063, r23064, r23065, r23066, r23067, r23068, r23069, r23070, r23071, r23072, r23073, r23074, r23075; void setup_mpfr_f_dm() { mpfr_set_default_prec(1616); mpfr_init(r23057); mpfr_init_set_str(r23058, "1", 10, MPFR_RNDN); mpfr_init(r23059); mpfr_init(r23060); mpfr_init(r23061); mpfr_init(r23062); mpfr_init_set_str(r23063, "4.4676756717265675e-05", 10, MPFR_RNDN); mpfr_init(r23064); mpfr_init(r23065); mpfr_init(r23066); mpfr_init_set_str(r23067, "1/3", 10, MPFR_RNDN); mpfr_init(r23068); mpfr_init_set_str(r23069, "1/2", 10, MPFR_RNDN); mpfr_init(r23070); mpfr_init(r23071); mpfr_init(r23072); mpfr_init(r23073); mpfr_init(r23074); mpfr_init(r23075); } double f_dm(double N) { mpfr_set_d(r23057, N, MPFR_RNDN); ; mpfr_add(r23059, r23057, r23058, MPFR_RNDN); mpfr_log(r23060, r23059, MPFR_RNDN); mpfr_log(r23061, r23057, MPFR_RNDN); mpfr_sub(r23062, r23060, r23061, MPFR_RNDN); ; mpfr_set_si(r23064, mpfr_cmp(r23062, r23063) <= 0, MPFR_RNDN); mpfr_mul(r23065, r23057, r23057, MPFR_RNDN); mpfr_div(r23066, r23058, r23065, MPFR_RNDN); ; mpfr_div(r23068, r23067, r23057, MPFR_RNDN); ; mpfr_sub(r23070, r23068, r23069, MPFR_RNDN); mpfr_div(r23071, r23058, r23057, MPFR_RNDN); mpfr_fma(r23072, r23066, r23070, r23071, MPFR_RNDN); mpfr_div(r23073, r23059, r23057, MPFR_RNDN); mpfr_log(r23074, r23073, MPFR_RNDN); if (mpfr_get_si(r23064, MPFR_RNDN)) { mpfr_set(r23075, r23072, MPFR_RNDN); } else { mpfr_set(r23075, r23074, MPFR_RNDN); }; return mpfr_get_d(r23075, MPFR_RNDN); }