#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19871 = 0.5f; float r19872 = re; float r19873 = cos(r19872); float r19874 = r19871 * r19873; float r19875 = 0.0f; float r19876 = im; float r19877 = r19875 - r19876; float r19878 = exp(r19877); float r19879 = exp(r19876); float r19880 = r19878 - r19879; float r19881 = r19874 * r19880; return r19881; } double f_id(double re, double im) { double r19882 = 0.5; double r19883 = re; double r19884 = cos(r19883); double r19885 = r19882 * r19884; double r19886 = 0.0; double r19887 = im; double r19888 = r19886 - r19887; double r19889 = exp(r19888); double r19890 = exp(r19887); double r19891 = r19889 - r19890; double r19892 = r19885 * r19891; return r19892; } double f_of(float re, float im) { 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 = re; float r19908 = cos(r19907); float r19909 = 0.5f; float r19910 = r19908 * r19909; float r19911 = r19906 * r19910; return r19911; } double f_od(double re, double im) { 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 = re; double r19927 = cos(r19926); double r19928 = 0.5; double r19929 = r19927 * r19928; double r19930 = r19925 * r19929; return r19930; } 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 r19931, r19932, r19933, r19934, r19935, r19936, r19937, r19938, r19939, r19940, r19941; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19931, "0.5", 10, MPFR_RNDN); mpfr_init(r19932); mpfr_init(r19933); mpfr_init(r19934); mpfr_init_set_str(r19935, "0", 10, MPFR_RNDN); mpfr_init(r19936); mpfr_init(r19937); mpfr_init(r19938); mpfr_init(r19939); mpfr_init(r19940); mpfr_init(r19941); } double f_im(double re, double im) { ; mpfr_set_d(r19932, re, MPFR_RNDN); mpfr_cos(r19933, r19932, MPFR_RNDN); mpfr_mul(r19934, r19931, r19933, MPFR_RNDN); ; mpfr_set_d(r19936, im, MPFR_RNDN); mpfr_sub(r19937, r19935, r19936, MPFR_RNDN); mpfr_exp(r19938, r19937, MPFR_RNDN); mpfr_exp(r19939, r19936, MPFR_RNDN); mpfr_sub(r19940, r19938, r19939, MPFR_RNDN); mpfr_mul(r19941, r19934, r19940, MPFR_RNDN); return mpfr_get_d(r19941, MPFR_RNDN); } static mpfr_t r19942, r19943, r19944, r19945, r19946, r19947, r19948, r19949, r19950, r19951, r19952, r19953, r19954, r19955, r19956, r19957, r19958, r19959, r19960; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19942, "1/60", 10, MPFR_RNDN); mpfr_init(r19943); mpfr_init_set_str(r19944, "5", 10, MPFR_RNDN); mpfr_init(r19945); mpfr_init(r19946); mpfr_init_set_str(r19947, "2", 10, MPFR_RNDN); mpfr_init(r19948); mpfr_init_set_str(r19949, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19950, "3", 10, MPFR_RNDN); mpfr_init(r19951); mpfr_init(r19952); mpfr_init(r19953); mpfr_init(r19954); mpfr_init(r19955); mpfr_init(r19956); mpfr_init(r19957); mpfr_init_set_str(r19958, "0.5", 10, MPFR_RNDN); mpfr_init(r19959); mpfr_init(r19960); } double f_fm(double re, double im) { ; mpfr_set_d(r19943, im, MPFR_RNDN); ; mpfr_pow(r19945, r19943, r19944, MPFR_RNDN); mpfr_mul(r19946, r19942, r19945, MPFR_RNDN); ; mpfr_mul(r19948, r19947, r19943, MPFR_RNDN); ; ; mpfr_pow(r19951, r19943, r19950, MPFR_RNDN); mpfr_mul(r19952, r19949, r19951, MPFR_RNDN); mpfr_add(r19953, r19948, r19952, MPFR_RNDN); mpfr_add(r19954, r19946, r19953, MPFR_RNDN); mpfr_neg(r19955, r19954, MPFR_RNDN); mpfr_set_d(r19956, re, MPFR_RNDN); mpfr_cos(r19957, r19956, MPFR_RNDN); ; mpfr_mul(r19959, r19957, r19958, MPFR_RNDN); mpfr_mul(r19960, r19955, r19959, MPFR_RNDN); return mpfr_get_d(r19960, MPFR_RNDN); } static mpfr_t r19961, r19962, r19963, r19964, r19965, r19966, r19967, r19968, r19969, r19970, r19971, r19972, r19973, r19974, r19975, r19976, r19977, r19978, r19979; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19961, "1/60", 10, MPFR_RNDN); mpfr_init(r19962); mpfr_init_set_str(r19963, "5", 10, MPFR_RNDN); mpfr_init(r19964); mpfr_init(r19965); mpfr_init_set_str(r19966, "2", 10, MPFR_RNDN); mpfr_init(r19967); mpfr_init_set_str(r19968, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19969, "3", 10, MPFR_RNDN); mpfr_init(r19970); mpfr_init(r19971); mpfr_init(r19972); mpfr_init(r19973); mpfr_init(r19974); mpfr_init(r19975); mpfr_init(r19976); mpfr_init_set_str(r19977, "0.5", 10, MPFR_RNDN); mpfr_init(r19978); mpfr_init(r19979); } double f_dm(double re, double im) { ; mpfr_set_d(r19962, im, MPFR_RNDN); ; mpfr_pow(r19964, r19962, r19963, MPFR_RNDN); mpfr_mul(r19965, r19961, r19964, MPFR_RNDN); ; mpfr_mul(r19967, r19966, r19962, MPFR_RNDN); ; ; mpfr_pow(r19970, r19962, r19969, MPFR_RNDN); mpfr_mul(r19971, r19968, r19970, MPFR_RNDN); mpfr_add(r19972, r19967, r19971, MPFR_RNDN); mpfr_add(r19973, r19965, r19972, MPFR_RNDN); mpfr_neg(r19974, r19973, MPFR_RNDN); mpfr_set_d(r19975, re, MPFR_RNDN); mpfr_cos(r19976, r19975, MPFR_RNDN); ; mpfr_mul(r19978, r19976, r19977, MPFR_RNDN); mpfr_mul(r19979, r19974, r19978, MPFR_RNDN); return mpfr_get_d(r19979, MPFR_RNDN); }