#include #include #include #include #include char *name = "2tan (problem 3.3.2)"; double f_if(float x, float eps) { float r24370 = x; float r24371 = eps; float r24372 = r24370 + r24371; float r24373 = tan(r24372); float r24374 = tan(r24370); float r24375 = r24373 - r24374; return r24375; } double f_id(double x, double eps) { double r24376 = x; double r24377 = eps; double r24378 = r24376 + r24377; double r24379 = tan(r24378); double r24380 = tan(r24376); double r24381 = r24379 - r24380; return r24381; } double f_of(float x, float eps) { float r24382 = eps; float r24383 = sin(r24382); float r24384 = cos(r24382); float r24385 = x; float r24386 = sin(r24385); float r24387 = 3; float r24388 = pow(r24386, r24387); float r24389 = pow(r24383, r24387); float r24390 = r24388 * r24389; float r24391 = cos(r24385); float r24392 = r24391 * r24384; float r24393 = pow(r24392, r24387); float r24394 = r24390 / r24393; float r24395 = cbrt(r24394); float r24396 = 1; float r24397 = r24395 - r24396; float r24398 = r24384 * r24397; float r24399 = r24383 / r24398; float r24400 = r24386 / r24391; float r24401 = r24386 * r24383; float r24402 = r24384 * r24391; float r24403 = r24401 / r24402; float r24404 = r24403 - r24396; float r24405 = r24391 * r24404; float r24406 = r24386 / r24405; float r24407 = r24400 + r24406; float r24408 = r24399 + r24407; float r24409 = -r24408; return r24409; } double f_od(double x, double eps) { double r24410 = eps; double r24411 = sin(r24410); double r24412 = cos(r24410); double r24413 = x; double r24414 = sin(r24413); double r24415 = 3; double r24416 = pow(r24414, r24415); double r24417 = pow(r24411, r24415); double r24418 = r24416 * r24417; double r24419 = cos(r24413); double r24420 = r24419 * r24412; double r24421 = pow(r24420, r24415); double r24422 = r24418 / r24421; double r24423 = cbrt(r24422); double r24424 = 1; double r24425 = r24423 - r24424; double r24426 = r24412 * r24425; double r24427 = r24411 / r24426; double r24428 = r24414 / r24419; double r24429 = r24414 * r24411; double r24430 = r24412 * r24419; double r24431 = r24429 / r24430; double r24432 = r24431 - r24424; double r24433 = r24419 * r24432; double r24434 = r24414 / r24433; double r24435 = r24428 + r24434; double r24436 = r24427 + r24435; double r24437 = -r24436; return r24437; } 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 r24438, r24439, r24440, r24441, r24442, r24443; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r24438); mpfr_init(r24439); mpfr_init(r24440); mpfr_init(r24441); mpfr_init(r24442); mpfr_init(r24443); } double f_im(double x, double eps) { mpfr_set_d(r24438, x, MPFR_RNDN); mpfr_set_d(r24439, eps, MPFR_RNDN); mpfr_add(r24440, r24438, r24439, MPFR_RNDN); mpfr_tan(r24441, r24440, MPFR_RNDN); mpfr_tan(r24442, r24438, MPFR_RNDN); mpfr_sub(r24443, r24441, r24442, MPFR_RNDN); return mpfr_get_d(r24443, MPFR_RNDN); } static mpfr_t r24444, r24445, r24446, r24447, r24448, r24449, r24450, r24451, r24452, r24453, r24454, r24455, r24456, r24457, r24458, r24459, r24460, r24461, r24462, r24463, r24464, r24465, r24466, r24467, r24468, r24469, r24470, r24471; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r24444); mpfr_init(r24445); mpfr_init(r24446); mpfr_init(r24447); mpfr_init(r24448); mpfr_init_set_str(r24449, "3", 10, MPFR_RNDN); mpfr_init(r24450); mpfr_init(r24451); mpfr_init(r24452); mpfr_init(r24453); mpfr_init(r24454); mpfr_init(r24455); mpfr_init(r24456); mpfr_init(r24457); mpfr_init_set_str(r24458, "1", 10, MPFR_RNDN); mpfr_init(r24459); mpfr_init(r24460); mpfr_init(r24461); mpfr_init(r24462); mpfr_init(r24463); mpfr_init(r24464); mpfr_init(r24465); mpfr_init(r24466); mpfr_init(r24467); mpfr_init(r24468); mpfr_init(r24469); mpfr_init(r24470); mpfr_init(r24471); } double f_fm(double x, double eps) { mpfr_set_d(r24444, eps, MPFR_RNDN); mpfr_sin(r24445, r24444, MPFR_RNDN); mpfr_cos(r24446, r24444, MPFR_RNDN); mpfr_set_d(r24447, x, MPFR_RNDN); mpfr_sin(r24448, r24447, MPFR_RNDN); ; mpfr_pow(r24450, r24448, r24449, MPFR_RNDN); mpfr_pow(r24451, r24445, r24449, MPFR_RNDN); mpfr_mul(r24452, r24450, r24451, MPFR_RNDN); mpfr_cos(r24453, r24447, MPFR_RNDN); mpfr_mul(r24454, r24453, r24446, MPFR_RNDN); mpfr_pow(r24455, r24454, r24449, MPFR_RNDN); mpfr_div(r24456, r24452, r24455, MPFR_RNDN); mpfr_cbrt(r24457, r24456, MPFR_RNDN); ; mpfr_sub(r24459, r24457, r24458, MPFR_RNDN); mpfr_mul(r24460, r24446, r24459, MPFR_RNDN); mpfr_div(r24461, r24445, r24460, MPFR_RNDN); mpfr_div(r24462, r24448, r24453, MPFR_RNDN); mpfr_mul(r24463, r24448, r24445, MPFR_RNDN); mpfr_mul(r24464, r24446, r24453, MPFR_RNDN); mpfr_div(r24465, r24463, r24464, MPFR_RNDN); mpfr_sub(r24466, r24465, r24458, MPFR_RNDN); mpfr_mul(r24467, r24453, r24466, MPFR_RNDN); mpfr_div(r24468, r24448, r24467, MPFR_RNDN); mpfr_add(r24469, r24462, r24468, MPFR_RNDN); mpfr_add(r24470, r24461, r24469, MPFR_RNDN); mpfr_neg(r24471, r24470, MPFR_RNDN); return mpfr_get_d(r24471, MPFR_RNDN); } static mpfr_t r24472, r24473, r24474, r24475, r24476, r24477, r24478, r24479, r24480, r24481, r24482, r24483, r24484, r24485, r24486, r24487, r24488, r24489, r24490, r24491, r24492, r24493, r24494, r24495, r24496, r24497, r24498, r24499; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r24472); mpfr_init(r24473); mpfr_init(r24474); mpfr_init(r24475); mpfr_init(r24476); mpfr_init_set_str(r24477, "3", 10, MPFR_RNDN); mpfr_init(r24478); mpfr_init(r24479); mpfr_init(r24480); mpfr_init(r24481); mpfr_init(r24482); mpfr_init(r24483); mpfr_init(r24484); mpfr_init(r24485); mpfr_init_set_str(r24486, "1", 10, MPFR_RNDN); mpfr_init(r24487); mpfr_init(r24488); mpfr_init(r24489); mpfr_init(r24490); mpfr_init(r24491); mpfr_init(r24492); mpfr_init(r24493); mpfr_init(r24494); mpfr_init(r24495); mpfr_init(r24496); mpfr_init(r24497); mpfr_init(r24498); mpfr_init(r24499); } double f_dm(double x, double eps) { mpfr_set_d(r24472, eps, MPFR_RNDN); mpfr_sin(r24473, r24472, MPFR_RNDN); mpfr_cos(r24474, r24472, MPFR_RNDN); mpfr_set_d(r24475, x, MPFR_RNDN); mpfr_sin(r24476, r24475, MPFR_RNDN); ; mpfr_pow(r24478, r24476, r24477, MPFR_RNDN); mpfr_pow(r24479, r24473, r24477, MPFR_RNDN); mpfr_mul(r24480, r24478, r24479, MPFR_RNDN); mpfr_cos(r24481, r24475, MPFR_RNDN); mpfr_mul(r24482, r24481, r24474, MPFR_RNDN); mpfr_pow(r24483, r24482, r24477, MPFR_RNDN); mpfr_div(r24484, r24480, r24483, MPFR_RNDN); mpfr_cbrt(r24485, r24484, MPFR_RNDN); ; mpfr_sub(r24487, r24485, r24486, MPFR_RNDN); mpfr_mul(r24488, r24474, r24487, MPFR_RNDN); mpfr_div(r24489, r24473, r24488, MPFR_RNDN); mpfr_div(r24490, r24476, r24481, MPFR_RNDN); mpfr_mul(r24491, r24476, r24473, MPFR_RNDN); mpfr_mul(r24492, r24474, r24481, MPFR_RNDN); mpfr_div(r24493, r24491, r24492, MPFR_RNDN); mpfr_sub(r24494, r24493, r24486, MPFR_RNDN); mpfr_mul(r24495, r24481, r24494, MPFR_RNDN); mpfr_div(r24496, r24476, r24495, MPFR_RNDN); mpfr_add(r24497, r24490, r24496, MPFR_RNDN); mpfr_add(r24498, r24489, r24497, MPFR_RNDN); mpfr_neg(r24499, r24498, MPFR_RNDN); return mpfr_get_d(r24499, MPFR_RNDN); }