#include #include #include #include #include char *name = "2tan (problem 3.3.2)"; double f_if(float x, float eps) { float r24243 = x; float r24244 = eps; float r24245 = r24243 + r24244; float r24246 = tan(r24245); float r24247 = tan(r24243); float r24248 = r24246 - r24247; return r24248; } double f_id(double x, double eps) { double r24249 = x; double r24250 = eps; double r24251 = r24249 + r24250; double r24252 = tan(r24251); double r24253 = tan(r24249); double r24254 = r24252 - r24253; return r24254; } double f_of(float x, float eps) { float r24255 = eps; float r24256 = sin(r24255); float r24257 = cos(r24255); float r24258 = x; float r24259 = sin(r24258); float r24260 = r24259 * r24256; float r24261 = cos(r24258); float r24262 = r24257 * r24261; float r24263 = r24260 / r24262; float r24264 = 1; float r24265 = r24263 - r24264; float r24266 = r24257 * r24265; float r24267 = r24256 / r24266; float r24268 = r24261 * r24259; float r24269 = r24259 / r24261; float r24270 = r24256 / r24257; float r24271 = -r24264; float r24272 = fma(r24269, r24270, r24271); float r24273 = fma(r24268, r24272, r24268); float r24274 = r24261 * r24265; float r24275 = r24261 * r24274; float r24276 = r24273 / r24275; float r24277 = r24267 + r24276; float r24278 = -r24277; return r24278; } double f_od(double x, double eps) { double r24279 = eps; double r24280 = sin(r24279); double r24281 = cos(r24279); double r24282 = x; double r24283 = sin(r24282); double r24284 = r24283 * r24280; double r24285 = cos(r24282); double r24286 = r24281 * r24285; double r24287 = r24284 / r24286; double r24288 = 1; double r24289 = r24287 - r24288; double r24290 = r24281 * r24289; double r24291 = r24280 / r24290; double r24292 = r24285 * r24283; double r24293 = r24283 / r24285; double r24294 = r24280 / r24281; double r24295 = -r24288; double r24296 = fma(r24293, r24294, r24295); double r24297 = fma(r24292, r24296, r24292); double r24298 = r24285 * r24289; double r24299 = r24285 * r24298; double r24300 = r24297 / r24299; double r24301 = r24291 + r24300; double r24302 = -r24301; return r24302; } 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 r24303, r24304, r24305, r24306, r24307, r24308; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r24303); mpfr_init(r24304); mpfr_init(r24305); mpfr_init(r24306); mpfr_init(r24307); mpfr_init(r24308); } double f_im(double x, double eps) { mpfr_set_d(r24303, x, MPFR_RNDN); mpfr_set_d(r24304, eps, MPFR_RNDN); mpfr_add(r24305, r24303, r24304, MPFR_RNDN); mpfr_tan(r24306, r24305, MPFR_RNDN); mpfr_tan(r24307, r24303, MPFR_RNDN); mpfr_sub(r24308, r24306, r24307, MPFR_RNDN); return mpfr_get_d(r24308, MPFR_RNDN); } static mpfr_t r24309, r24310, r24311, r24312, r24313, r24314, r24315, r24316, r24317, r24318, r24319, r24320, r24321, r24322, r24323, r24324, r24325, r24326, r24327, r24328, r24329, r24330, r24331, r24332; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init(r24309); mpfr_init(r24310); mpfr_init(r24311); mpfr_init(r24312); mpfr_init(r24313); mpfr_init(r24314); mpfr_init(r24315); mpfr_init(r24316); mpfr_init(r24317); mpfr_init_set_str(r24318, "1", 10, MPFR_RNDN); mpfr_init(r24319); mpfr_init(r24320); mpfr_init(r24321); mpfr_init(r24322); mpfr_init(r24323); mpfr_init(r24324); mpfr_init(r24325); mpfr_init(r24326); mpfr_init(r24327); mpfr_init(r24328); mpfr_init(r24329); mpfr_init(r24330); mpfr_init(r24331); mpfr_init(r24332); } double f_fm(double x, double eps) { mpfr_set_d(r24309, eps, MPFR_RNDN); mpfr_sin(r24310, r24309, MPFR_RNDN); mpfr_cos(r24311, r24309, MPFR_RNDN); mpfr_set_d(r24312, x, MPFR_RNDN); mpfr_sin(r24313, r24312, MPFR_RNDN); mpfr_mul(r24314, r24313, r24310, MPFR_RNDN); mpfr_cos(r24315, r24312, MPFR_RNDN); mpfr_mul(r24316, r24311, r24315, MPFR_RNDN); mpfr_div(r24317, r24314, r24316, MPFR_RNDN); ; mpfr_sub(r24319, r24317, r24318, MPFR_RNDN); mpfr_mul(r24320, r24311, r24319, MPFR_RNDN); mpfr_div(r24321, r24310, r24320, MPFR_RNDN); mpfr_mul(r24322, r24315, r24313, MPFR_RNDN); mpfr_div(r24323, r24313, r24315, MPFR_RNDN); mpfr_div(r24324, r24310, r24311, MPFR_RNDN); mpfr_neg(r24325, r24318, MPFR_RNDN); mpfr_fma(r24326, r24323, r24324, r24325, MPFR_RNDN); mpfr_fma(r24327, r24322, r24326, r24322, MPFR_RNDN); mpfr_mul(r24328, r24315, r24319, MPFR_RNDN); mpfr_mul(r24329, r24315, r24328, MPFR_RNDN); mpfr_div(r24330, r24327, r24329, MPFR_RNDN); mpfr_add(r24331, r24321, r24330, MPFR_RNDN); mpfr_neg(r24332, r24331, MPFR_RNDN); return mpfr_get_d(r24332, MPFR_RNDN); } static mpfr_t r24333, r24334, r24335, r24336, r24337, r24338, r24339, r24340, r24341, r24342, r24343, r24344, r24345, r24346, r24347, r24348, r24349, r24350, r24351, r24352, r24353, r24354, r24355, r24356; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init(r24333); mpfr_init(r24334); mpfr_init(r24335); mpfr_init(r24336); mpfr_init(r24337); mpfr_init(r24338); mpfr_init(r24339); mpfr_init(r24340); mpfr_init(r24341); mpfr_init_set_str(r24342, "1", 10, MPFR_RNDN); mpfr_init(r24343); mpfr_init(r24344); mpfr_init(r24345); mpfr_init(r24346); mpfr_init(r24347); mpfr_init(r24348); mpfr_init(r24349); mpfr_init(r24350); mpfr_init(r24351); mpfr_init(r24352); mpfr_init(r24353); mpfr_init(r24354); mpfr_init(r24355); mpfr_init(r24356); } double f_dm(double x, double eps) { mpfr_set_d(r24333, eps, MPFR_RNDN); mpfr_sin(r24334, r24333, MPFR_RNDN); mpfr_cos(r24335, r24333, MPFR_RNDN); mpfr_set_d(r24336, x, MPFR_RNDN); mpfr_sin(r24337, r24336, MPFR_RNDN); mpfr_mul(r24338, r24337, r24334, MPFR_RNDN); mpfr_cos(r24339, r24336, MPFR_RNDN); mpfr_mul(r24340, r24335, r24339, MPFR_RNDN); mpfr_div(r24341, r24338, r24340, MPFR_RNDN); ; mpfr_sub(r24343, r24341, r24342, MPFR_RNDN); mpfr_mul(r24344, r24335, r24343, MPFR_RNDN); mpfr_div(r24345, r24334, r24344, MPFR_RNDN); mpfr_mul(r24346, r24339, r24337, MPFR_RNDN); mpfr_div(r24347, r24337, r24339, MPFR_RNDN); mpfr_div(r24348, r24334, r24335, MPFR_RNDN); mpfr_neg(r24349, r24342, MPFR_RNDN); mpfr_fma(r24350, r24347, r24348, r24349, MPFR_RNDN); mpfr_fma(r24351, r24346, r24350, r24346, MPFR_RNDN); mpfr_mul(r24352, r24339, r24343, MPFR_RNDN); mpfr_mul(r24353, r24339, r24352, MPFR_RNDN); mpfr_div(r24354, r24351, r24353, MPFR_RNDN); mpfr_add(r24355, r24345, r24354, MPFR_RNDN); mpfr_neg(r24356, r24355, MPFR_RNDN); return mpfr_get_d(r24356, MPFR_RNDN); }