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