#include #include #include #include #include char *name = "2cos (problem 3.3.5)"; double f_if(float x, float eps) { float r23242 = x; float r23243 = eps; float r23244 = r23242 + r23243; float r23245 = cos(r23244); float r23246 = cos(r23242); float r23247 = r23245 - r23246; return r23247; } double f_id(double x, double eps) { double r23248 = x; double r23249 = eps; double r23250 = r23248 + r23249; double r23251 = cos(r23250); double r23252 = cos(r23248); double r23253 = r23251 - r23252; return r23253; } double f_of(float x, float eps) { float r23254 = eps; float r23255 = -0.15281236229617784; bool r23256 = r23254 <= r23255; float r23257 = 0.00023034535512454643; bool r23258 = r23254 <= r23257; float r23259 = !r23258; bool r23260 = r23256 || r23259; float r23261 = cos(r23254); float r23262 = x; float r23263 = cos(r23262); float r23264 = r23261 * r23263; float r23265 = sin(r23254); float r23266 = sin(r23262); float r23267 = r23265 * r23266; float r23268 = r23264 - r23267; float r23269 = r23268 - r23263; float r23270 = -2; float r23271 = r23254 + r23262; float r23272 = r23271 + r23262; float r23273 = 2; float r23274 = r23272 / r23273; float r23275 = sin(r23274); float r23276 = r23254 / r23273; float r23277 = sin(r23276); float r23278 = r23275 * r23277; float r23279 = r23270 * r23278; float r23280 = r23260 ? r23269 : r23279; return r23280; } double f_od(double x, double eps) { double r23281 = eps; double r23282 = -0.15281236229617784; bool r23283 = r23281 <= r23282; double r23284 = 0.00023034535512454643; bool r23285 = r23281 <= r23284; double r23286 = !r23285; bool r23287 = r23283 || r23286; double r23288 = cos(r23281); double r23289 = x; double r23290 = cos(r23289); double r23291 = r23288 * r23290; double r23292 = sin(r23281); double r23293 = sin(r23289); double r23294 = r23292 * r23293; double r23295 = r23291 - r23294; double r23296 = r23295 - r23290; double r23297 = -2; double r23298 = r23281 + r23289; double r23299 = r23298 + r23289; double r23300 = 2; double r23301 = r23299 / r23300; double r23302 = sin(r23301); double r23303 = r23281 / r23300; double r23304 = sin(r23303); double r23305 = r23302 * r23304; double r23306 = r23297 * r23305; double r23307 = r23287 ? r23296 : r23306; return r23307; } 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 r23308, r23309, r23310, r23311, r23312, r23313; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r23308); mpfr_init(r23309); mpfr_init(r23310); mpfr_init(r23311); mpfr_init(r23312); mpfr_init(r23313); } double f_im(double x, double eps) { mpfr_set_d(r23308, x, MPFR_RNDN); mpfr_set_d(r23309, eps, MPFR_RNDN); mpfr_add(r23310, r23308, r23309, MPFR_RNDN); mpfr_cos(r23311, r23310, MPFR_RNDN); mpfr_cos(r23312, r23308, MPFR_RNDN); mpfr_sub(r23313, r23311, r23312, MPFR_RNDN); return mpfr_get_d(r23313, MPFR_RNDN); } static mpfr_t r23314, r23315, r23316, r23317, r23318, r23319, r23320, r23321, r23322, r23323, r23324, r23325, r23326, r23327, r23328, r23329, r23330, r23331, r23332, r23333, r23334, r23335, r23336, r23337, r23338, r23339, r23340; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r23314); mpfr_init_set_str(r23315, "-0.15281236229617784", 10, MPFR_RNDN); mpfr_init(r23316); mpfr_init_set_str(r23317, "0.00023034535512454643", 10, MPFR_RNDN); mpfr_init(r23318); mpfr_init(r23319); mpfr_init(r23320); mpfr_init(r23321); mpfr_init(r23322); mpfr_init(r23323); mpfr_init(r23324); mpfr_init(r23325); mpfr_init(r23326); mpfr_init(r23327); mpfr_init(r23328); mpfr_init(r23329); mpfr_init_set_str(r23330, "-2", 10, MPFR_RNDN); mpfr_init(r23331); mpfr_init(r23332); mpfr_init_set_str(r23333, "2", 10, MPFR_RNDN); mpfr_init(r23334); mpfr_init(r23335); mpfr_init(r23336); mpfr_init(r23337); mpfr_init(r23338); mpfr_init(r23339); mpfr_init(r23340); } double f_fm(double x, double eps) { mpfr_set_d(r23314, eps, MPFR_RNDN); ; mpfr_set_si(r23316, mpfr_cmp(r23314, r23315) <= 0, MPFR_RNDN); ; mpfr_set_si(r23318, mpfr_cmp(r23314, r23317) <= 0, MPFR_RNDN); mpfr_set_si(r23319, !mpfr_get_si(r23318, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r23320, mpfr_get_si(r23316, MPFR_RNDN) || mpfr_get_si(r23319, MPFR_RNDN), MPFR_RNDN); mpfr_cos(r23321, r23314, MPFR_RNDN); mpfr_set_d(r23322, x, MPFR_RNDN); mpfr_cos(r23323, r23322, MPFR_RNDN); mpfr_mul(r23324, r23321, r23323, MPFR_RNDN); mpfr_sin(r23325, r23314, MPFR_RNDN); mpfr_sin(r23326, r23322, MPFR_RNDN); mpfr_mul(r23327, r23325, r23326, MPFR_RNDN); mpfr_sub(r23328, r23324, r23327, MPFR_RNDN); mpfr_sub(r23329, r23328, r23323, MPFR_RNDN); ; mpfr_add(r23331, r23314, r23322, MPFR_RNDN); mpfr_add(r23332, r23331, r23322, MPFR_RNDN); ; mpfr_div(r23334, r23332, r23333, MPFR_RNDN); mpfr_sin(r23335, r23334, MPFR_RNDN); mpfr_div(r23336, r23314, r23333, MPFR_RNDN); mpfr_sin(r23337, r23336, MPFR_RNDN); mpfr_mul(r23338, r23335, r23337, MPFR_RNDN); mpfr_mul(r23339, r23330, r23338, MPFR_RNDN); if (mpfr_get_si(r23320, MPFR_RNDN)) { mpfr_set(r23340, r23329, MPFR_RNDN); } else { mpfr_set(r23340, r23339, MPFR_RNDN); }; return mpfr_get_d(r23340, MPFR_RNDN); } static mpfr_t r23341, r23342, r23343, r23344, r23345, r23346, r23347, r23348, r23349, r23350, r23351, r23352, r23353, r23354, r23355, r23356, r23357, r23358, r23359, r23360, r23361, r23362, r23363, r23364, r23365, r23366, r23367; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r23341); mpfr_init_set_str(r23342, "-0.15281236229617784", 10, MPFR_RNDN); mpfr_init(r23343); mpfr_init_set_str(r23344, "0.00023034535512454643", 10, MPFR_RNDN); mpfr_init(r23345); mpfr_init(r23346); mpfr_init(r23347); mpfr_init(r23348); mpfr_init(r23349); mpfr_init(r23350); mpfr_init(r23351); mpfr_init(r23352); mpfr_init(r23353); mpfr_init(r23354); mpfr_init(r23355); mpfr_init(r23356); mpfr_init_set_str(r23357, "-2", 10, MPFR_RNDN); mpfr_init(r23358); mpfr_init(r23359); mpfr_init_set_str(r23360, "2", 10, MPFR_RNDN); mpfr_init(r23361); mpfr_init(r23362); mpfr_init(r23363); mpfr_init(r23364); mpfr_init(r23365); mpfr_init(r23366); mpfr_init(r23367); } double f_dm(double x, double eps) { mpfr_set_d(r23341, eps, MPFR_RNDN); ; mpfr_set_si(r23343, mpfr_cmp(r23341, r23342) <= 0, MPFR_RNDN); ; mpfr_set_si(r23345, mpfr_cmp(r23341, r23344) <= 0, MPFR_RNDN); mpfr_set_si(r23346, !mpfr_get_si(r23345, MPFR_RNDN), MPFR_RNDN); mpfr_set_si(r23347, mpfr_get_si(r23343, MPFR_RNDN) || mpfr_get_si(r23346, MPFR_RNDN), MPFR_RNDN); mpfr_cos(r23348, r23341, MPFR_RNDN); mpfr_set_d(r23349, x, MPFR_RNDN); mpfr_cos(r23350, r23349, MPFR_RNDN); mpfr_mul(r23351, r23348, r23350, MPFR_RNDN); mpfr_sin(r23352, r23341, MPFR_RNDN); mpfr_sin(r23353, r23349, MPFR_RNDN); mpfr_mul(r23354, r23352, r23353, MPFR_RNDN); mpfr_sub(r23355, r23351, r23354, MPFR_RNDN); mpfr_sub(r23356, r23355, r23350, MPFR_RNDN); ; mpfr_add(r23358, r23341, r23349, MPFR_RNDN); mpfr_add(r23359, r23358, r23349, MPFR_RNDN); ; mpfr_div(r23361, r23359, r23360, MPFR_RNDN); mpfr_sin(r23362, r23361, MPFR_RNDN); mpfr_div(r23363, r23341, r23360, MPFR_RNDN); mpfr_sin(r23364, r23363, MPFR_RNDN); mpfr_mul(r23365, r23362, r23364, MPFR_RNDN); mpfr_mul(r23366, r23357, r23365, MPFR_RNDN); if (mpfr_get_si(r23347, MPFR_RNDN)) { mpfr_set(r23367, r23356, MPFR_RNDN); } else { mpfr_set(r23367, r23366, MPFR_RNDN); }; return mpfr_get_d(r23367, MPFR_RNDN); }