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