#include #include #include #include #include char *name = "2cos (problem 3.3.5)"; double f_if(float x, float eps) { float r25865 = x; float r25866 = eps; float r25867 = r25865 + r25866; float r25868 = cos(r25867); float r25869 = cos(r25865); float r25870 = r25868 - r25869; return r25870; } double f_id(double x, double eps) { double r25871 = x; double r25872 = eps; double r25873 = r25871 + r25872; double r25874 = cos(r25873); double r25875 = cos(r25871); double r25876 = r25874 - r25875; return r25876; } double f_of(float x, float eps) { float r25877 = -2; float r25878 = eps; float r25879 = 2; float r25880 = r25878 / r25879; float r25881 = sin(r25880); float r25882 = 1/2; float r25883 = r25882 * r25878; float r25884 = sin(r25883); float r25885 = x; float r25886 = cos(r25885); float r25887 = r25884 * r25886; float r25888 = r25887 * r25887; float r25889 = cos(r25883); float r25890 = sin(r25885); float r25891 = r25889 * r25890; float r25892 = r25891 * r25891; float r25893 = r25888 - r25892; float r25894 = r25887 - r25891; float r25895 = r25893 / r25894; float r25896 = r25881 * r25895; float r25897 = r25877 * r25896; return r25897; } double f_od(double x, double eps) { double r25898 = -2; double r25899 = eps; double r25900 = 2; double r25901 = r25899 / r25900; double r25902 = sin(r25901); double r25903 = 1/2; double r25904 = r25903 * r25899; double r25905 = sin(r25904); double r25906 = x; double r25907 = cos(r25906); double r25908 = r25905 * r25907; double r25909 = r25908 * r25908; double r25910 = cos(r25904); double r25911 = sin(r25906); double r25912 = r25910 * r25911; double r25913 = r25912 * r25912; double r25914 = r25909 - r25913; double r25915 = r25908 - r25912; double r25916 = r25914 / r25915; double r25917 = r25902 * r25916; double r25918 = r25898 * r25917; return r25918; } 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 r25919, r25920, r25921, r25922, r25923, r25924; void setup_mpfr_f_im() { mpfr_set_default_prec(2384); mpfr_init(r25919); mpfr_init(r25920); mpfr_init(r25921); mpfr_init(r25922); mpfr_init(r25923); mpfr_init(r25924); } double f_im(double x, double eps) { mpfr_set_d(r25919, x, MPFR_RNDN); mpfr_set_d(r25920, eps, MPFR_RNDN); mpfr_add(r25921, r25919, r25920, MPFR_RNDN); mpfr_cos(r25922, r25921, MPFR_RNDN); mpfr_cos(r25923, r25919, MPFR_RNDN); mpfr_sub(r25924, r25922, r25923, MPFR_RNDN); return mpfr_get_d(r25924, MPFR_RNDN); } static mpfr_t r25925, r25926, r25927, r25928, r25929, r25930, r25931, r25932, r25933, r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945; void setup_mpfr_f_fm() { mpfr_set_default_prec(2384); mpfr_init_set_str(r25925, "-2", 10, MPFR_RNDN); mpfr_init(r25926); mpfr_init_set_str(r25927, "2", 10, MPFR_RNDN); mpfr_init(r25928); mpfr_init(r25929); mpfr_init_set_str(r25930, "1/2", 10, MPFR_RNDN); mpfr_init(r25931); mpfr_init(r25932); mpfr_init(r25933); mpfr_init(r25934); mpfr_init(r25935); mpfr_init(r25936); mpfr_init(r25937); mpfr_init(r25938); mpfr_init(r25939); mpfr_init(r25940); mpfr_init(r25941); mpfr_init(r25942); mpfr_init(r25943); mpfr_init(r25944); mpfr_init(r25945); } double f_fm(double x, double eps) { ; mpfr_set_d(r25926, eps, MPFR_RNDN); ; mpfr_div(r25928, r25926, r25927, MPFR_RNDN); mpfr_sin(r25929, r25928, MPFR_RNDN); ; mpfr_mul(r25931, r25930, r25926, MPFR_RNDN); mpfr_sin(r25932, r25931, MPFR_RNDN); mpfr_set_d(r25933, x, MPFR_RNDN); mpfr_cos(r25934, r25933, MPFR_RNDN); mpfr_mul(r25935, r25932, r25934, MPFR_RNDN); mpfr_mul(r25936, r25935, r25935, MPFR_RNDN); mpfr_cos(r25937, r25931, MPFR_RNDN); mpfr_sin(r25938, r25933, MPFR_RNDN); mpfr_mul(r25939, r25937, r25938, MPFR_RNDN); mpfr_mul(r25940, r25939, r25939, MPFR_RNDN); mpfr_sub(r25941, r25936, r25940, MPFR_RNDN); mpfr_sub(r25942, r25935, r25939, MPFR_RNDN); mpfr_div(r25943, r25941, r25942, MPFR_RNDN); mpfr_mul(r25944, r25929, r25943, MPFR_RNDN); mpfr_mul(r25945, r25925, r25944, MPFR_RNDN); return mpfr_get_d(r25945, MPFR_RNDN); } static mpfr_t r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955, r25956, r25957, r25958, r25959, r25960, r25961, r25962, r25963, r25964, r25965, r25966; void setup_mpfr_f_dm() { mpfr_set_default_prec(2384); mpfr_init_set_str(r25946, "-2", 10, MPFR_RNDN); mpfr_init(r25947); mpfr_init_set_str(r25948, "2", 10, MPFR_RNDN); mpfr_init(r25949); mpfr_init(r25950); mpfr_init_set_str(r25951, "1/2", 10, MPFR_RNDN); mpfr_init(r25952); mpfr_init(r25953); mpfr_init(r25954); mpfr_init(r25955); mpfr_init(r25956); mpfr_init(r25957); mpfr_init(r25958); mpfr_init(r25959); mpfr_init(r25960); mpfr_init(r25961); mpfr_init(r25962); mpfr_init(r25963); mpfr_init(r25964); mpfr_init(r25965); mpfr_init(r25966); } double f_dm(double x, double eps) { ; mpfr_set_d(r25947, eps, MPFR_RNDN); ; mpfr_div(r25949, r25947, r25948, MPFR_RNDN); mpfr_sin(r25950, r25949, MPFR_RNDN); ; mpfr_mul(r25952, r25951, r25947, MPFR_RNDN); mpfr_sin(r25953, r25952, MPFR_RNDN); mpfr_set_d(r25954, x, MPFR_RNDN); mpfr_cos(r25955, r25954, MPFR_RNDN); mpfr_mul(r25956, r25953, r25955, MPFR_RNDN); mpfr_mul(r25957, r25956, r25956, MPFR_RNDN); mpfr_cos(r25958, r25952, MPFR_RNDN); mpfr_sin(r25959, r25954, MPFR_RNDN); mpfr_mul(r25960, r25958, r25959, MPFR_RNDN); mpfr_mul(r25961, r25960, r25960, MPFR_RNDN); mpfr_sub(r25962, r25957, r25961, MPFR_RNDN); mpfr_sub(r25963, r25956, r25960, MPFR_RNDN); mpfr_div(r25964, r25962, r25963, MPFR_RNDN); mpfr_mul(r25965, r25950, r25964, MPFR_RNDN); mpfr_mul(r25966, r25946, r25965, MPFR_RNDN); return mpfr_get_d(r25966, MPFR_RNDN); }