#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r19869 = 0.5f; float r19870 = re; float r19871 = sin(r19870); float r19872 = r19869 * r19871; float r19873 = im; float r19874 = -r19873; float r19875 = exp(r19874); float r19876 = exp(r19873); float r19877 = r19875 - r19876; float r19878 = r19872 * r19877; return r19878; } double f_id(double re, double im) { double r19879 = 0.5; double r19880 = re; double r19881 = sin(r19880); double r19882 = r19879 * r19881; double r19883 = im; double r19884 = -r19883; double r19885 = exp(r19884); double r19886 = exp(r19883); double r19887 = r19885 - r19886; double r19888 = r19882 * r19887; return r19888; } double f_of(float re, float im) { float r19889 = 0.5f; float r19890 = re; float r19891 = sin(r19890); float r19892 = r19889 * r19891; float r19893 = 0.016666666666666666f; float r19894 = im; float r19895 = 5.0f; float r19896 = pow(r19894, r19895); float r19897 = r19893 * r19896; float r19898 = 2.0f; float r19899 = r19898 * r19894; float r19900 = 0.3333333333333333f; float r19901 = 3.0f; float r19902 = pow(r19894, r19901); float r19903 = r19900 * r19902; float r19904 = r19899 + r19903; float r19905 = r19897 + r19904; float r19906 = -r19905; float r19907 = r19892 * r19906; return r19907; } double f_od(double re, double im) { double r19908 = 0.5; double r19909 = re; double r19910 = sin(r19909); double r19911 = r19908 * r19910; double r19912 = 0.016666666666666666; double r19913 = im; double r19914 = 5.0; double r19915 = pow(r19913, r19914); double r19916 = r19912 * r19915; double r19917 = 2.0; double r19918 = r19917 * r19913; double r19919 = 0.3333333333333333; double r19920 = 3.0; double r19921 = pow(r19913, r19920); double r19922 = r19919 * r19921; double r19923 = r19918 + r19922; double r19924 = r19916 + r19923; double r19925 = -r19924; double r19926 = r19911 * r19925; return r19926; } 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 r19927, r19928, r19929, r19930, r19931, r19932, r19933, r19934, r19935, r19936; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19927, "0.5", 10, MPFR_RNDN); mpfr_init(r19928); mpfr_init(r19929); mpfr_init(r19930); mpfr_init(r19931); mpfr_init(r19932); mpfr_init(r19933); mpfr_init(r19934); mpfr_init(r19935); mpfr_init(r19936); } double f_im(double re, double im) { ; mpfr_set_d(r19928, re, MPFR_RNDN); mpfr_sin(r19929, r19928, MPFR_RNDN); mpfr_mul(r19930, r19927, r19929, MPFR_RNDN); mpfr_set_d(r19931, im, MPFR_RNDN); mpfr_neg(r19932, r19931, MPFR_RNDN); mpfr_exp(r19933, r19932, MPFR_RNDN); mpfr_exp(r19934, r19931, MPFR_RNDN); mpfr_sub(r19935, r19933, r19934, MPFR_RNDN); mpfr_mul(r19936, r19930, r19935, MPFR_RNDN); return mpfr_get_d(r19936, MPFR_RNDN); } static mpfr_t r19937, r19938, r19939, r19940, r19941, r19942, r19943, r19944, r19945, r19946, r19947, r19948, r19949, r19950, r19951, r19952, r19953, r19954, r19955; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19937, "0.5", 10, MPFR_RNDN); mpfr_init(r19938); mpfr_init(r19939); mpfr_init(r19940); mpfr_init_set_str(r19941, "1/60", 10, MPFR_RNDN); mpfr_init(r19942); mpfr_init_set_str(r19943, "5", 10, MPFR_RNDN); mpfr_init(r19944); mpfr_init(r19945); mpfr_init_set_str(r19946, "2", 10, MPFR_RNDN); mpfr_init(r19947); mpfr_init_set_str(r19948, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19949, "3", 10, MPFR_RNDN); mpfr_init(r19950); mpfr_init(r19951); mpfr_init(r19952); mpfr_init(r19953); mpfr_init(r19954); mpfr_init(r19955); } double f_fm(double re, double im) { ; mpfr_set_d(r19938, re, MPFR_RNDN); mpfr_sin(r19939, r19938, MPFR_RNDN); mpfr_mul(r19940, r19937, r19939, MPFR_RNDN); ; mpfr_set_d(r19942, im, MPFR_RNDN); ; mpfr_pow(r19944, r19942, r19943, MPFR_RNDN); mpfr_mul(r19945, r19941, r19944, MPFR_RNDN); ; mpfr_mul(r19947, r19946, r19942, MPFR_RNDN); ; ; mpfr_pow(r19950, r19942, r19949, MPFR_RNDN); mpfr_mul(r19951, r19948, r19950, MPFR_RNDN); mpfr_add(r19952, r19947, r19951, MPFR_RNDN); mpfr_add(r19953, r19945, r19952, MPFR_RNDN); mpfr_neg(r19954, r19953, MPFR_RNDN); mpfr_mul(r19955, r19940, r19954, MPFR_RNDN); return mpfr_get_d(r19955, MPFR_RNDN); } static mpfr_t r19956, r19957, r19958, r19959, r19960, r19961, r19962, r19963, r19964, r19965, r19966, r19967, r19968, r19969, r19970, r19971, r19972, r19973, r19974; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19956, "0.5", 10, MPFR_RNDN); mpfr_init(r19957); mpfr_init(r19958); mpfr_init(r19959); mpfr_init_set_str(r19960, "1/60", 10, MPFR_RNDN); mpfr_init(r19961); mpfr_init_set_str(r19962, "5", 10, MPFR_RNDN); mpfr_init(r19963); mpfr_init(r19964); mpfr_init_set_str(r19965, "2", 10, MPFR_RNDN); mpfr_init(r19966); mpfr_init_set_str(r19967, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19968, "3", 10, MPFR_RNDN); mpfr_init(r19969); mpfr_init(r19970); mpfr_init(r19971); mpfr_init(r19972); mpfr_init(r19973); mpfr_init(r19974); } double f_dm(double re, double im) { ; mpfr_set_d(r19957, re, MPFR_RNDN); mpfr_sin(r19958, r19957, MPFR_RNDN); mpfr_mul(r19959, r19956, r19958, MPFR_RNDN); ; mpfr_set_d(r19961, im, MPFR_RNDN); ; mpfr_pow(r19963, r19961, r19962, MPFR_RNDN); mpfr_mul(r19964, r19960, r19963, MPFR_RNDN); ; mpfr_mul(r19966, r19965, r19961, MPFR_RNDN); ; ; mpfr_pow(r19969, r19961, r19968, MPFR_RNDN); mpfr_mul(r19970, r19967, r19969, MPFR_RNDN); mpfr_add(r19971, r19966, r19970, MPFR_RNDN); mpfr_add(r19972, r19964, r19971, MPFR_RNDN); mpfr_neg(r19973, r19972, MPFR_RNDN); mpfr_mul(r19974, r19959, r19973, MPFR_RNDN); return mpfr_get_d(r19974, MPFR_RNDN); }