#include #include #include #include #include char *name = "2tan (problem 3.3.2)"; double f_if(float x, float eps) { float r22644 = x; float r22645 = eps; float r22646 = r22644 + r22645; float r22647 = tan(r22646); float r22648 = tan(r22644); float r22649 = r22647 - r22648; return r22649; } double f_id(double x, double eps) { double r22650 = x; double r22651 = eps; double r22652 = r22650 + r22651; double r22653 = tan(r22652); double r22654 = tan(r22650); double r22655 = r22653 - r22654; return r22655; } double f_of(float x, float eps) { float r22656 = eps; float r22657 = sin(r22656); float r22658 = 1; float r22659 = cos(r22656); float r22660 = r22657 / r22659; float r22661 = x; float r22662 = sin(r22661); float r22663 = cos(r22661); float r22664 = r22662 / r22663; float r22665 = r22660 * r22664; float r22666 = r22658 - r22665; float r22667 = r22666 * r22659; float r22668 = r22657 / r22667; float r22669 = r22664 * r22664; float r22670 = r22662 / r22666; float r22671 = r22670 / r22663; float r22672 = r22671 * r22671; float r22673 = r22669 - r22672; float r22674 = r22664 + r22671; float r22675 = r22673 / r22674; float r22676 = r22668 - r22675; return r22676; } double f_od(double x, double eps) { double r22677 = eps; double r22678 = sin(r22677); double r22679 = 1; double r22680 = cos(r22677); double r22681 = r22678 / r22680; double r22682 = x; double r22683 = sin(r22682); double r22684 = cos(r22682); double r22685 = r22683 / r22684; double r22686 = r22681 * r22685; double r22687 = r22679 - r22686; double r22688 = r22687 * r22680; double r22689 = r22678 / r22688; double r22690 = r22685 * r22685; double r22691 = r22683 / r22687; double r22692 = r22691 / r22684; double r22693 = r22692 * r22692; double r22694 = r22690 - r22693; double r22695 = r22685 + r22692; double r22696 = r22694 / r22695; double r22697 = r22689 - r22696; return r22697; } 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 r22698, r22699, r22700, r22701, r22702, r22703; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init(r22698); mpfr_init(r22699); mpfr_init(r22700); mpfr_init(r22701); mpfr_init(r22702); mpfr_init(r22703); } double f_im(double x, double eps) { mpfr_set_d(r22698, x, MPFR_RNDN); mpfr_set_d(r22699, eps, MPFR_RNDN); mpfr_add(r22700, r22698, r22699, MPFR_RNDN); mpfr_tan(r22701, r22700, MPFR_RNDN); mpfr_tan(r22702, r22698, MPFR_RNDN); mpfr_sub(r22703, r22701, r22702, MPFR_RNDN); return mpfr_get_d(r22703, MPFR_RNDN); } static mpfr_t r22704, r22705, r22706, r22707, r22708, r22709, r22710, r22711, r22712, r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r22704); mpfr_init(r22705); mpfr_init_set_str(r22706, "1", 10, MPFR_RNDN); mpfr_init(r22707); mpfr_init(r22708); mpfr_init(r22709); mpfr_init(r22710); mpfr_init(r22711); mpfr_init(r22712); mpfr_init(r22713); mpfr_init(r22714); mpfr_init(r22715); mpfr_init(r22716); mpfr_init(r22717); mpfr_init(r22718); mpfr_init(r22719); mpfr_init(r22720); mpfr_init(r22721); mpfr_init(r22722); mpfr_init(r22723); mpfr_init(r22724); } double f_fm(double x, double eps) { mpfr_set_d(r22704, eps, MPFR_RNDN); mpfr_sin(r22705, r22704, MPFR_RNDN); ; mpfr_cos(r22707, r22704, MPFR_RNDN); mpfr_div(r22708, r22705, r22707, MPFR_RNDN); mpfr_set_d(r22709, x, MPFR_RNDN); mpfr_sin(r22710, r22709, MPFR_RNDN); mpfr_cos(r22711, r22709, MPFR_RNDN); mpfr_div(r22712, r22710, r22711, MPFR_RNDN); mpfr_mul(r22713, r22708, r22712, MPFR_RNDN); mpfr_sub(r22714, r22706, r22713, MPFR_RNDN); mpfr_mul(r22715, r22714, r22707, MPFR_RNDN); mpfr_div(r22716, r22705, r22715, MPFR_RNDN); mpfr_mul(r22717, r22712, r22712, MPFR_RNDN); mpfr_div(r22718, r22710, r22714, MPFR_RNDN); mpfr_div(r22719, r22718, r22711, MPFR_RNDN); mpfr_mul(r22720, r22719, r22719, MPFR_RNDN); mpfr_sub(r22721, r22717, r22720, MPFR_RNDN); mpfr_add(r22722, r22712, r22719, MPFR_RNDN); mpfr_div(r22723, r22721, r22722, MPFR_RNDN); mpfr_sub(r22724, r22716, r22723, MPFR_RNDN); return mpfr_get_d(r22724, MPFR_RNDN); } static mpfr_t r22725, r22726, r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r22725); mpfr_init(r22726); mpfr_init_set_str(r22727, "1", 10, MPFR_RNDN); mpfr_init(r22728); mpfr_init(r22729); mpfr_init(r22730); mpfr_init(r22731); mpfr_init(r22732); mpfr_init(r22733); mpfr_init(r22734); mpfr_init(r22735); mpfr_init(r22736); mpfr_init(r22737); mpfr_init(r22738); mpfr_init(r22739); mpfr_init(r22740); mpfr_init(r22741); mpfr_init(r22742); mpfr_init(r22743); mpfr_init(r22744); mpfr_init(r22745); } double f_dm(double x, double eps) { mpfr_set_d(r22725, eps, MPFR_RNDN); mpfr_sin(r22726, r22725, MPFR_RNDN); ; mpfr_cos(r22728, r22725, MPFR_RNDN); mpfr_div(r22729, r22726, r22728, MPFR_RNDN); mpfr_set_d(r22730, x, MPFR_RNDN); mpfr_sin(r22731, r22730, MPFR_RNDN); mpfr_cos(r22732, r22730, MPFR_RNDN); mpfr_div(r22733, r22731, r22732, MPFR_RNDN); mpfr_mul(r22734, r22729, r22733, MPFR_RNDN); mpfr_sub(r22735, r22727, r22734, MPFR_RNDN); mpfr_mul(r22736, r22735, r22728, MPFR_RNDN); mpfr_div(r22737, r22726, r22736, MPFR_RNDN); mpfr_mul(r22738, r22733, r22733, MPFR_RNDN); mpfr_div(r22739, r22731, r22735, MPFR_RNDN); mpfr_div(r22740, r22739, r22732, MPFR_RNDN); mpfr_mul(r22741, r22740, r22740, MPFR_RNDN); mpfr_sub(r22742, r22738, r22741, MPFR_RNDN); mpfr_add(r22743, r22733, r22740, MPFR_RNDN); mpfr_div(r22744, r22742, r22743, MPFR_RNDN); mpfr_sub(r22745, r22737, r22744, MPFR_RNDN); return mpfr_get_d(r22745, MPFR_RNDN); }