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