#include #include #include #include #include char *name = "2tan (problem 3.3.2)"; double f_if(float x, float eps) { float r24495 = x; float r24496 = eps; float r24497 = r24495 + r24496; float r24498 = tan(r24497); float r24499 = tan(r24495); float r24500 = r24498 - r24499; return r24500; } double f_id(double x, double eps) { double r24501 = x; double r24502 = eps; double r24503 = r24501 + r24502; double r24504 = tan(r24503); double r24505 = tan(r24501); double r24506 = r24504 - r24505; return r24506; } double f_of(float x, float eps) { float r24507 = eps; float r24508 = sin(r24507); float r24509 = 1; float r24510 = x; float r24511 = sin(r24510); float r24512 = cos(r24510); float r24513 = r24511 / r24512; float r24514 = cos(r24507); float r24515 = r24508 / r24514; float r24516 = r24513 * r24515; float r24517 = r24509 - r24516; float r24518 = r24517 * r24514; float r24519 = r24508 / r24518; float r24520 = r24508 * r24511; float r24521 = r24512 * r24514; float r24522 = r24520 / r24521; float r24523 = r24509 - r24522; float r24524 = r24513 / r24523; float r24525 = r24513 - r24524; float r24526 = exp(r24525); float r24527 = log(r24526); float r24528 = r24519 - r24527; return r24528; } double f_od(double x, double eps) { double r24529 = eps; double r24530 = sin(r24529); double r24531 = 1; double r24532 = x; double r24533 = sin(r24532); double r24534 = cos(r24532); double r24535 = r24533 / r24534; double r24536 = cos(r24529); double r24537 = r24530 / r24536; double r24538 = r24535 * r24537; double r24539 = r24531 - r24538; double r24540 = r24539 * r24536; double r24541 = r24530 / r24540; double r24542 = r24530 * r24533; double r24543 = r24534 * r24536; double r24544 = r24542 / r24543; double r24545 = r24531 - r24544; double r24546 = r24535 / r24545; double r24547 = r24535 - r24546; double r24548 = exp(r24547); double r24549 = log(r24548); double r24550 = r24541 - r24549; return r24550; } 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 r24551, r24552, r24553, r24554, r24555, r24556; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r24551); mpfr_init(r24552); mpfr_init(r24553); mpfr_init(r24554); mpfr_init(r24555); mpfr_init(r24556); } double f_im(double x, double eps) { mpfr_set_d(r24551, x, MPFR_RNDN); mpfr_set_d(r24552, eps, MPFR_RNDN); mpfr_add(r24553, r24551, r24552, MPFR_RNDN); mpfr_tan(r24554, r24553, MPFR_RNDN); mpfr_tan(r24555, r24551, MPFR_RNDN); mpfr_sub(r24556, r24554, r24555, MPFR_RNDN); return mpfr_get_d(r24556, MPFR_RNDN); } static mpfr_t r24557, r24558, r24559, r24560, r24561, r24562, r24563, r24564, r24565, r24566, r24567, r24568, r24569, r24570, r24571, r24572, r24573, r24574, r24575, r24576, r24577, r24578; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r24557); mpfr_init(r24558); mpfr_init_set_str(r24559, "1", 10, MPFR_RNDN); mpfr_init(r24560); mpfr_init(r24561); mpfr_init(r24562); mpfr_init(r24563); mpfr_init(r24564); mpfr_init(r24565); mpfr_init(r24566); mpfr_init(r24567); mpfr_init(r24568); mpfr_init(r24569); mpfr_init(r24570); mpfr_init(r24571); mpfr_init(r24572); mpfr_init(r24573); mpfr_init(r24574); mpfr_init(r24575); mpfr_init(r24576); mpfr_init(r24577); mpfr_init(r24578); } double f_fm(double x, double eps) { mpfr_set_d(r24557, eps, MPFR_RNDN); mpfr_sin(r24558, r24557, MPFR_RNDN); ; mpfr_set_d(r24560, x, MPFR_RNDN); mpfr_sin(r24561, r24560, MPFR_RNDN); mpfr_cos(r24562, r24560, MPFR_RNDN); mpfr_div(r24563, r24561, r24562, MPFR_RNDN); mpfr_cos(r24564, r24557, MPFR_RNDN); mpfr_div(r24565, r24558, r24564, MPFR_RNDN); mpfr_mul(r24566, r24563, r24565, MPFR_RNDN); mpfr_sub(r24567, r24559, r24566, MPFR_RNDN); mpfr_mul(r24568, r24567, r24564, MPFR_RNDN); mpfr_div(r24569, r24558, r24568, MPFR_RNDN); mpfr_mul(r24570, r24558, r24561, MPFR_RNDN); mpfr_mul(r24571, r24562, r24564, MPFR_RNDN); mpfr_div(r24572, r24570, r24571, MPFR_RNDN); mpfr_sub(r24573, r24559, r24572, MPFR_RNDN); mpfr_div(r24574, r24563, r24573, MPFR_RNDN); mpfr_sub(r24575, r24563, r24574, MPFR_RNDN); mpfr_exp(r24576, r24575, MPFR_RNDN); mpfr_log(r24577, r24576, MPFR_RNDN); mpfr_sub(r24578, r24569, r24577, MPFR_RNDN); return mpfr_get_d(r24578, MPFR_RNDN); } static mpfr_t r24579, r24580, r24581, r24582, r24583, r24584, r24585, r24586, r24587, r24588, r24589, r24590, r24591, r24592, r24593, r24594, r24595, r24596, r24597, r24598, r24599, r24600; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r24579); mpfr_init(r24580); mpfr_init_set_str(r24581, "1", 10, MPFR_RNDN); mpfr_init(r24582); mpfr_init(r24583); mpfr_init(r24584); mpfr_init(r24585); mpfr_init(r24586); mpfr_init(r24587); mpfr_init(r24588); mpfr_init(r24589); mpfr_init(r24590); mpfr_init(r24591); mpfr_init(r24592); mpfr_init(r24593); mpfr_init(r24594); mpfr_init(r24595); mpfr_init(r24596); mpfr_init(r24597); mpfr_init(r24598); mpfr_init(r24599); mpfr_init(r24600); } double f_dm(double x, double eps) { mpfr_set_d(r24579, eps, MPFR_RNDN); mpfr_sin(r24580, r24579, MPFR_RNDN); ; mpfr_set_d(r24582, x, MPFR_RNDN); mpfr_sin(r24583, r24582, MPFR_RNDN); mpfr_cos(r24584, r24582, MPFR_RNDN); mpfr_div(r24585, r24583, r24584, MPFR_RNDN); mpfr_cos(r24586, r24579, MPFR_RNDN); mpfr_div(r24587, r24580, r24586, MPFR_RNDN); mpfr_mul(r24588, r24585, r24587, MPFR_RNDN); mpfr_sub(r24589, r24581, r24588, MPFR_RNDN); mpfr_mul(r24590, r24589, r24586, MPFR_RNDN); mpfr_div(r24591, r24580, r24590, MPFR_RNDN); mpfr_mul(r24592, r24580, r24583, MPFR_RNDN); mpfr_mul(r24593, r24584, r24586, MPFR_RNDN); mpfr_div(r24594, r24592, r24593, MPFR_RNDN); mpfr_sub(r24595, r24581, r24594, MPFR_RNDN); mpfr_div(r24596, r24585, r24595, MPFR_RNDN); mpfr_sub(r24597, r24585, r24596, MPFR_RNDN); mpfr_exp(r24598, r24597, MPFR_RNDN); mpfr_log(r24599, r24598, MPFR_RNDN); mpfr_sub(r24600, r24591, r24599, MPFR_RNDN); return mpfr_get_d(r24600, MPFR_RNDN); }