#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r18879 = 0.5f; float r18880 = re; float r18881 = cos(r18880); float r18882 = r18879 * r18881; float r18883 = 0.0f; float r18884 = im; float r18885 = r18883 - r18884; float r18886 = exp(r18885); float r18887 = exp(r18884); float r18888 = r18886 - r18887; float r18889 = r18882 * r18888; return r18889; } double f_id(double re, double im) { double r18890 = 0.5; double r18891 = re; double r18892 = cos(r18891); double r18893 = r18890 * r18892; double r18894 = 0.0; double r18895 = im; double r18896 = r18894 - r18895; double r18897 = exp(r18896); double r18898 = exp(r18895); double r18899 = r18897 - r18898; double r18900 = r18893 * r18899; return r18900; } double f_of(float re, float im) { float r18901 = im; float r18902 = r18901 * (r18901 * r18901); float r18903 = 0.3333333333333333f; float r18904 = r18902 * r18903; float r18905 = 0.016666666666666666f; float r18906 = 5.0f; float r18907 = pow(r18901, r18906); float r18908 = r18905 * r18907; float r18909 = 2.0f; float r18910 = r18901 * r18909; float r18911 = r18908 + r18910; float r18912 = r18904 + r18911; float r18913 = 1.0f; float r18914 = r18912 / r18913; float r18915 = 0.5f; float r18916 = -r18915; float r18917 = re; float r18918 = cos(r18917); float r18919 = r18916 * r18918; float r18920 = r18914 * r18919; return r18920; } double f_od(double re, double im) { double r18921 = im; double r18922 = r18921 * (r18921 * r18921); double r18923 = 0.3333333333333333; double r18924 = r18922 * r18923; double r18925 = 0.016666666666666666; double r18926 = 5.0; double r18927 = pow(r18921, r18926); double r18928 = r18925 * r18927; double r18929 = 2.0; double r18930 = r18921 * r18929; double r18931 = r18928 + r18930; double r18932 = r18924 + r18931; double r18933 = 1.0; double r18934 = r18932 / r18933; double r18935 = 0.5; double r18936 = -r18935; double r18937 = re; double r18938 = cos(r18937); double r18939 = r18936 * r18938; double r18940 = r18934 * r18939; return r18940; } 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 r18941, r18942, r18943, r18944, r18945, r18946, r18947, r18948, r18949, r18950, r18951; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r18941, "0.5", 10, MPFR_RNDN); mpfr_init(r18942); mpfr_init(r18943); mpfr_init(r18944); mpfr_init_set_str(r18945, "0", 10, MPFR_RNDN); mpfr_init(r18946); mpfr_init(r18947); mpfr_init(r18948); mpfr_init(r18949); mpfr_init(r18950); mpfr_init(r18951); } double f_im(double re, double im) { ; mpfr_set_d(r18942, re, MPFR_RNDN); mpfr_cos(r18943, r18942, MPFR_RNDN); mpfr_mul(r18944, r18941, r18943, MPFR_RNDN); ; mpfr_set_d(r18946, im, MPFR_RNDN); mpfr_sub(r18947, r18945, r18946, MPFR_RNDN); mpfr_exp(r18948, r18947, MPFR_RNDN); mpfr_exp(r18949, r18946, MPFR_RNDN); mpfr_sub(r18950, r18948, r18949, MPFR_RNDN); mpfr_mul(r18951, r18944, r18950, MPFR_RNDN); return mpfr_get_d(r18951, MPFR_RNDN); } static mpfr_t r18952, r18953, r18954, r18955, r18956, r18957, r18958, r18959, r18960, r18961, r18962, r18963, r18964, r18965, r18966, r18967, r18968, r18969, r18970, r18971; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18952); mpfr_init(r18953); mpfr_init_set_str(r18954, "1/3", 10, MPFR_RNDN); mpfr_init(r18955); mpfr_init_set_str(r18956, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18957, "5", 10, MPFR_RNDN); mpfr_init(r18958); mpfr_init(r18959); mpfr_init_set_str(r18960, "2", 10, MPFR_RNDN); mpfr_init(r18961); mpfr_init(r18962); mpfr_init(r18963); mpfr_init_set_str(r18964, "1", 10, MPFR_RNDN); mpfr_init(r18965); mpfr_init_set_str(r18966, "0.5", 10, MPFR_RNDN); mpfr_init(r18967); mpfr_init(r18968); mpfr_init(r18969); mpfr_init(r18970); mpfr_init(r18971); } double f_fm(double re, double im) { mpfr_set_d(r18952, im, MPFR_RNDN); mpfr_mul(r18953, r18952, r18952, MPFR_RNDN); mpfr_mul(r18953, r18953, r18952, MPFR_RNDN); ; mpfr_mul(r18955, r18953, r18954, MPFR_RNDN); ; ; mpfr_pow(r18958, r18952, r18957, MPFR_RNDN); mpfr_mul(r18959, r18956, r18958, MPFR_RNDN); ; mpfr_mul(r18961, r18952, r18960, MPFR_RNDN); mpfr_add(r18962, r18959, r18961, MPFR_RNDN); mpfr_add(r18963, r18955, r18962, MPFR_RNDN); ; mpfr_div(r18965, r18963, r18964, MPFR_RNDN); ; mpfr_neg(r18967, r18966, MPFR_RNDN); mpfr_set_d(r18968, re, MPFR_RNDN); mpfr_cos(r18969, r18968, MPFR_RNDN); mpfr_mul(r18970, r18967, r18969, MPFR_RNDN); mpfr_mul(r18971, r18965, r18970, MPFR_RNDN); return mpfr_get_d(r18971, MPFR_RNDN); } static mpfr_t r18972, r18973, r18974, r18975, r18976, r18977, r18978, r18979, r18980, r18981, r18982, r18983, r18984, r18985, r18986, r18987, r18988, r18989, r18990, r18991; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18972); mpfr_init(r18973); mpfr_init_set_str(r18974, "1/3", 10, MPFR_RNDN); mpfr_init(r18975); mpfr_init_set_str(r18976, "1/60", 10, MPFR_RNDN); mpfr_init_set_str(r18977, "5", 10, MPFR_RNDN); mpfr_init(r18978); mpfr_init(r18979); mpfr_init_set_str(r18980, "2", 10, MPFR_RNDN); mpfr_init(r18981); mpfr_init(r18982); mpfr_init(r18983); mpfr_init_set_str(r18984, "1", 10, MPFR_RNDN); mpfr_init(r18985); mpfr_init_set_str(r18986, "0.5", 10, MPFR_RNDN); mpfr_init(r18987); mpfr_init(r18988); mpfr_init(r18989); mpfr_init(r18990); mpfr_init(r18991); } double f_dm(double re, double im) { mpfr_set_d(r18972, im, MPFR_RNDN); mpfr_mul(r18973, r18972, r18972, MPFR_RNDN); mpfr_mul(r18973, r18973, r18972, MPFR_RNDN); ; mpfr_mul(r18975, r18973, r18974, MPFR_RNDN); ; ; mpfr_pow(r18978, r18972, r18977, MPFR_RNDN); mpfr_mul(r18979, r18976, r18978, MPFR_RNDN); ; mpfr_mul(r18981, r18972, r18980, MPFR_RNDN); mpfr_add(r18982, r18979, r18981, MPFR_RNDN); mpfr_add(r18983, r18975, r18982, MPFR_RNDN); ; mpfr_div(r18985, r18983, r18984, MPFR_RNDN); ; mpfr_neg(r18987, r18986, MPFR_RNDN); mpfr_set_d(r18988, re, MPFR_RNDN); mpfr_cos(r18989, r18988, MPFR_RNDN); mpfr_mul(r18990, r18987, r18989, MPFR_RNDN); mpfr_mul(r18991, r18985, r18990, MPFR_RNDN); return mpfr_get_d(r18991, MPFR_RNDN); }