#include #include #include #include #include char *name = "x / (x^2 + 1)"; double f_if(float x) { float r22950 = x; float r22951 = r22950 * r22950; float r22952 = 1; float r22953 = r22951 + r22952; float r22954 = r22950 / r22953; return r22954; } double f_id(double x) { double r22955 = x; double r22956 = r22955 * r22955; double r22957 = 1; double r22958 = r22956 + r22957; double r22959 = r22955 / r22958; return r22959; } double f_of(float x) { float r22960 = x; float r22961 = -23253.23141424013; bool r22962 = r22960 <= r22961; float r22963 = 5842628.826673498; bool r22964 = r22960 <= r22963; float r22965 = !r22964; bool r22966 = r22962 || r22965; float r22967 = 1; float r22968 = 5; float r22969 = pow(r22960, r22968); float r22970 = r22967 / r22969; float r22971 = r22967 / r22960; float r22972 = r22970 + r22971; float r22973 = 3; float r22974 = pow(r22960, r22973); float r22975 = r22967 / r22974; float r22976 = r22972 - r22975; float r22977 = r22960 * r22960; float r22978 = r22977 + r22967; float r22979 = r22967 / r22978; float r22980 = r22979 * r22960; float r22981 = r22966 ? r22976 : r22980; return r22981; } double f_od(double x) { double r22982 = x; double r22983 = -23253.23141424013; bool r22984 = r22982 <= r22983; double r22985 = 5842628.826673498; bool r22986 = r22982 <= r22985; double r22987 = !r22986; bool r22988 = r22984 || r22987; double r22989 = 1; double r22990 = 5; double r22991 = pow(r22982, r22990); double r22992 = r22989 / r22991; double r22993 = r22989 / r22982; double r22994 = r22992 + r22993; double r22995 = 3; double r22996 = pow(r22982, r22995); double r22997 = r22989 / r22996; double r22998 = r22994 - r22997; double r22999 = r22982 * r22982; double r23000 = r22999 + r22989; double r23001 = r22989 / r23000; double r23002 = r23001 * r22982; double r23003 = r22988 ? r22998 : r23002; return r23003; } 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 r23004, r23005, r23006, r23007, r23008; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23004); mpfr_init(r23005); mpfr_init_set_str(r23006, "1", 10, MPFR_RNDN); mpfr_init(r23007); mpfr_init(r23008); } double f_im(double x) { mpfr_set_d(r23004, x, MPFR_RNDN); mpfr_mul(r23005, r23004, r23004, MPFR_RNDN); ; mpfr_add(r23007, r23005, r23006, MPFR_RNDN); mpfr_div(r23008, r23004, r23007, MPFR_RNDN); return mpfr_get_d(r23008, MPFR_RNDN); } static mpfr_t r23009, r23010, r23011, r23012, r23013, r23014, r23015, r23016, r23017, r23018, r23019, r23020, r23021, r23022, r23023, r23024, r23025, r23026, r23027, r23028, r23029, r23030; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23009); mpfr_init_set_str(r23010, "-23253.23141424013", 10, MPFR_RNDN); mpfr_init(r23011); mpfr_init_set_str(r23012, "5842628.826673498", 10, MPFR_RNDN); mpfr_init(r23013); mpfr_init(r23014); mpfr_init(r23015); mpfr_init_set_str(r23016, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23017, "5", 10, MPFR_RNDN); mpfr_init(r23018); mpfr_init(r23019); mpfr_init(r23020); mpfr_init(r23021); mpfr_init_set_str(r23022, "3", 10, MPFR_RNDN); mpfr_init(r23023); mpfr_init(r23024); mpfr_init(r23025); mpfr_init(r23026); mpfr_init(r23027); mpfr_init(r23028); mpfr_init(r23029); mpfr_init(r23030); } double f_fm(double x) { mpfr_set_d(r23009, x, MPFR_RNDN); ; mpfr_set_si(r23011, mpfr_cmp(r23009, r23010) <= 0, MPFR_RNDN); ; mpfr_set_si(r23013, mpfr_cmp(r23009, r23012) <= 0, MPFR_RNDN); mpfr_set_si(r23014, !mpfr_get_si(r23013, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r23015, mpfr_get_si(r23011, MPFR_RNDN) || mpfr_get_si(r23014, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r23018, r23009, r23017, MPFR_RNDN); mpfr_div(r23019, r23016, r23018, MPFR_RNDN); mpfr_div(r23020, r23016, r23009, MPFR_RNDN); mpfr_add(r23021, r23019, r23020, MPFR_RNDN); ; mpfr_pow(r23023, r23009, r23022, MPFR_RNDN); mpfr_div(r23024, r23016, r23023, MPFR_RNDN); mpfr_sub(r23025, r23021, r23024, MPFR_RNDN); mpfr_mul(r23026, r23009, r23009, MPFR_RNDN); mpfr_add(r23027, r23026, r23016, MPFR_RNDN); mpfr_div(r23028, r23016, r23027, MPFR_RNDN); mpfr_mul(r23029, r23028, r23009, MPFR_RNDN); if (mpfr_get_si(r23015, MPFR_RNDN)) { mpfr_set(r23030, r23025, MPFR_RNDN); } else { mpfr_set(r23030, r23029, MPFR_RNDN); }; return mpfr_get_d(r23030, MPFR_RNDN); } static mpfr_t r23031, r23032, r23033, r23034, r23035, r23036, r23037, r23038, r23039, r23040, r23041, r23042, r23043, r23044, r23045, r23046, r23047, r23048, r23049, r23050, r23051, r23052; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23031); mpfr_init_set_str(r23032, "-23253.23141424013", 10, MPFR_RNDN); mpfr_init(r23033); mpfr_init_set_str(r23034, "5842628.826673498", 10, MPFR_RNDN); mpfr_init(r23035); mpfr_init(r23036); mpfr_init(r23037); mpfr_init_set_str(r23038, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23039, "5", 10, MPFR_RNDN); mpfr_init(r23040); mpfr_init(r23041); mpfr_init(r23042); mpfr_init(r23043); mpfr_init_set_str(r23044, "3", 10, MPFR_RNDN); mpfr_init(r23045); mpfr_init(r23046); mpfr_init(r23047); mpfr_init(r23048); mpfr_init(r23049); mpfr_init(r23050); mpfr_init(r23051); mpfr_init(r23052); } double f_dm(double x) { mpfr_set_d(r23031, x, MPFR_RNDN); ; mpfr_set_si(r23033, mpfr_cmp(r23031, r23032) <= 0, MPFR_RNDN); ; mpfr_set_si(r23035, mpfr_cmp(r23031, r23034) <= 0, MPFR_RNDN); mpfr_set_si(r23036, !mpfr_get_si(r23035, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r23037, mpfr_get_si(r23033, MPFR_RNDN) || mpfr_get_si(r23036, MPFR_RNDN), MPFR_RNDN); ; ; mpfr_pow(r23040, r23031, r23039, MPFR_RNDN); mpfr_div(r23041, r23038, r23040, MPFR_RNDN); mpfr_div(r23042, r23038, r23031, MPFR_RNDN); mpfr_add(r23043, r23041, r23042, MPFR_RNDN); ; mpfr_pow(r23045, r23031, r23044, MPFR_RNDN); mpfr_div(r23046, r23038, r23045, MPFR_RNDN); mpfr_sub(r23047, r23043, r23046, MPFR_RNDN); mpfr_mul(r23048, r23031, r23031, MPFR_RNDN); mpfr_add(r23049, r23048, r23038, MPFR_RNDN); mpfr_div(r23050, r23038, r23049, MPFR_RNDN); mpfr_mul(r23051, r23050, r23031, MPFR_RNDN); if (mpfr_get_si(r23037, MPFR_RNDN)) { mpfr_set(r23052, r23047, MPFR_RNDN); } else { mpfr_set(r23052, r23051, MPFR_RNDN); }; return mpfr_get_d(r23052, MPFR_RNDN); }