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