#include #include #include #include #include char *name = "math.sin on complex, imaginary part"; double f_if(float re, float im) { float r19899 = 0.5f; float r19900 = re; float r19901 = cos(r19900); float r19902 = r19899 * r19901; float r19903 = 0.0f; float r19904 = im; float r19905 = r19903 - r19904; float r19906 = exp(r19905); float r19907 = exp(r19904); float r19908 = r19906 - r19907; float r19909 = r19902 * r19908; return r19909; } double f_id(double re, double im) { double r19910 = 0.5; double r19911 = re; double r19912 = cos(r19911); double r19913 = r19910 * r19912; double r19914 = 0.0; double r19915 = im; double r19916 = r19914 - r19915; double r19917 = exp(r19916); double r19918 = exp(r19915); double r19919 = r19917 - r19918; double r19920 = r19913 * r19919; return r19920; } double f_of(float re, float im) { float r19921 = 0.016666666666666666f; float r19922 = im; float r19923 = 5.0f; float r19924 = pow(r19922, r19923); float r19925 = r19921 * r19924; float r19926 = 2.0f; float r19927 = r19926 * r19922; float r19928 = 0.3333333333333333f; float r19929 = 3.0f; float r19930 = pow(r19922, r19929); float r19931 = r19928 * r19930; float r19932 = r19927 + r19931; float r19933 = r19925 + r19932; float r19934 = -r19933; float r19935 = re; float r19936 = cos(r19935); float r19937 = 0.5f; float r19938 = r19936 * r19937; float r19939 = r19934 * r19938; return r19939; } double f_od(double re, double im) { double r19940 = 0.016666666666666666; double r19941 = im; double r19942 = 5.0; double r19943 = pow(r19941, r19942); double r19944 = r19940 * r19943; double r19945 = 2.0; double r19946 = r19945 * r19941; double r19947 = 0.3333333333333333; double r19948 = 3.0; double r19949 = pow(r19941, r19948); double r19950 = r19947 * r19949; double r19951 = r19946 + r19950; double r19952 = r19944 + r19951; double r19953 = -r19952; double r19954 = re; double r19955 = cos(r19954); double r19956 = 0.5; double r19957 = r19955 * r19956; double r19958 = r19953 * r19957; return r19958; } 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 r19959, r19960, r19961, r19962, r19963, r19964, r19965, r19966, r19967, r19968, r19969; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r19959, "0.5", 10, MPFR_RNDN); mpfr_init(r19960); mpfr_init(r19961); mpfr_init(r19962); mpfr_init_set_str(r19963, "0", 10, MPFR_RNDN); mpfr_init(r19964); mpfr_init(r19965); mpfr_init(r19966); mpfr_init(r19967); mpfr_init(r19968); mpfr_init(r19969); } double f_im(double re, double im) { ; mpfr_set_d(r19960, re, MPFR_RNDN); mpfr_cos(r19961, r19960, MPFR_RNDN); mpfr_mul(r19962, r19959, r19961, MPFR_RNDN); ; mpfr_set_d(r19964, im, MPFR_RNDN); mpfr_sub(r19965, r19963, r19964, MPFR_RNDN); mpfr_exp(r19966, r19965, MPFR_RNDN); mpfr_exp(r19967, r19964, MPFR_RNDN); mpfr_sub(r19968, r19966, r19967, MPFR_RNDN); mpfr_mul(r19969, r19962, r19968, MPFR_RNDN); return mpfr_get_d(r19969, MPFR_RNDN); } static mpfr_t r19970, r19971, r19972, r19973, r19974, r19975, r19976, r19977, r19978, r19979, r19980, r19981, r19982, r19983, r19984, r19985, r19986, r19987, r19988; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19970, "1/60", 10, MPFR_RNDN); mpfr_init(r19971); mpfr_init_set_str(r19972, "5", 10, MPFR_RNDN); mpfr_init(r19973); mpfr_init(r19974); mpfr_init_set_str(r19975, "2", 10, MPFR_RNDN); mpfr_init(r19976); mpfr_init_set_str(r19977, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19978, "3", 10, MPFR_RNDN); mpfr_init(r19979); mpfr_init(r19980); mpfr_init(r19981); mpfr_init(r19982); mpfr_init(r19983); mpfr_init(r19984); mpfr_init(r19985); mpfr_init_set_str(r19986, "0.5", 10, MPFR_RNDN); mpfr_init(r19987); mpfr_init(r19988); } double f_fm(double re, double im) { ; mpfr_set_d(r19971, im, MPFR_RNDN); ; mpfr_pow(r19973, r19971, r19972, MPFR_RNDN); mpfr_mul(r19974, r19970, r19973, MPFR_RNDN); ; mpfr_mul(r19976, r19975, r19971, MPFR_RNDN); ; ; mpfr_pow(r19979, r19971, r19978, MPFR_RNDN); mpfr_mul(r19980, r19977, r19979, MPFR_RNDN); mpfr_add(r19981, r19976, r19980, MPFR_RNDN); mpfr_add(r19982, r19974, r19981, MPFR_RNDN); mpfr_neg(r19983, r19982, MPFR_RNDN); mpfr_set_d(r19984, re, MPFR_RNDN); mpfr_cos(r19985, r19984, MPFR_RNDN); ; mpfr_mul(r19987, r19985, r19986, MPFR_RNDN); mpfr_mul(r19988, r19983, r19987, MPFR_RNDN); return mpfr_get_d(r19988, MPFR_RNDN); } static mpfr_t r19989, r19990, r19991, r19992, r19993, r19994, r19995, r19996, r19997, r19998, r19999, r20000, r20001, r20002, r20003, r20004, r20005, r20006, r20007; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r19989, "1/60", 10, MPFR_RNDN); mpfr_init(r19990); mpfr_init_set_str(r19991, "5", 10, MPFR_RNDN); mpfr_init(r19992); mpfr_init(r19993); mpfr_init_set_str(r19994, "2", 10, MPFR_RNDN); mpfr_init(r19995); mpfr_init_set_str(r19996, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r19997, "3", 10, MPFR_RNDN); mpfr_init(r19998); mpfr_init(r19999); mpfr_init(r20000); mpfr_init(r20001); mpfr_init(r20002); mpfr_init(r20003); mpfr_init(r20004); mpfr_init_set_str(r20005, "0.5", 10, MPFR_RNDN); mpfr_init(r20006); mpfr_init(r20007); } double f_dm(double re, double im) { ; mpfr_set_d(r19990, im, MPFR_RNDN); ; mpfr_pow(r19992, r19990, r19991, MPFR_RNDN); mpfr_mul(r19993, r19989, r19992, MPFR_RNDN); ; mpfr_mul(r19995, r19994, r19990, MPFR_RNDN); ; ; mpfr_pow(r19998, r19990, r19997, MPFR_RNDN); mpfr_mul(r19999, r19996, r19998, MPFR_RNDN); mpfr_add(r20000, r19995, r19999, MPFR_RNDN); mpfr_add(r20001, r19993, r20000, MPFR_RNDN); mpfr_neg(r20002, r20001, MPFR_RNDN); mpfr_set_d(r20003, re, MPFR_RNDN); mpfr_cos(r20004, r20003, MPFR_RNDN); ; mpfr_mul(r20006, r20004, r20005, MPFR_RNDN); mpfr_mul(r20007, r20002, r20006, MPFR_RNDN); return mpfr_get_d(r20007, MPFR_RNDN); }