#include #include #include #include #include char *name = "math.cos on complex, imaginary part"; double f_if(float re, float im) { float r20045 = 0.5f; float r20046 = re; float r20047 = sin(r20046); float r20048 = r20045 * r20047; float r20049 = im; float r20050 = -r20049; float r20051 = exp(r20050); float r20052 = exp(r20049); float r20053 = r20051 - r20052; float r20054 = r20048 * r20053; return r20054; } double f_id(double re, double im) { double r20055 = 0.5; double r20056 = re; double r20057 = sin(r20056); double r20058 = r20055 * r20057; double r20059 = im; double r20060 = -r20059; double r20061 = exp(r20060); double r20062 = exp(r20059); double r20063 = r20061 - r20062; double r20064 = r20058 * r20063; return r20064; } double f_of(float re, float im) { float r20065 = 0.5f; float r20066 = re; float r20067 = sin(r20066); float r20068 = r20065 * r20067; float r20069 = 0.01666666753590107f; float r20070 = im; float r20071 = 5.0f; float r20072 = pow(r20070, r20071); float r20073 = r20069 * r20072; float r20074 = 2.0f; float r20075 = r20074 * r20070; float r20076 = 0.3333333432674408f; float r20077 = 3.0f; float r20078 = pow(r20070, r20077); float r20079 = r20076 * r20078; float r20080 = r20075 + r20079; float r20081 = r20073 + r20080; float r20082 = -r20081; float r20083 = r20068 * r20082; return r20083; } double f_od(double re, double im) { double r20084 = 0.5; double r20085 = re; double r20086 = sin(r20085); double r20087 = r20084 * r20086; double r20088 = 0.01666666753590107; double r20089 = im; double r20090 = 5.0; double r20091 = pow(r20089, r20090); double r20092 = r20088 * r20091; double r20093 = 2.0; double r20094 = r20093 * r20089; double r20095 = 0.3333333432674408; double r20096 = 3.0; double r20097 = pow(r20089, r20096); double r20098 = r20095 * r20097; double r20099 = r20094 + r20098; double r20100 = r20092 + r20099; double r20101 = -r20100; double r20102 = r20087 * r20101; return r20102; } 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 r20103, r20104, r20105, r20106, r20107, r20108, r20109, r20110, r20111, r20112; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r20103, "0.5", 10, MPFR_RNDN); mpfr_init(r20104); mpfr_init(r20105); mpfr_init(r20106); mpfr_init(r20107); mpfr_init(r20108); mpfr_init(r20109); mpfr_init(r20110); mpfr_init(r20111); mpfr_init(r20112); } double f_im(double re, double im) { ; mpfr_set_d(r20104, re, MPFR_RNDN); mpfr_sin(r20105, r20104, MPFR_RNDN); mpfr_mul(r20106, r20103, r20105, MPFR_RNDN); mpfr_set_d(r20107, im, MPFR_RNDN); mpfr_neg(r20108, r20107, MPFR_RNDN); mpfr_exp(r20109, r20108, MPFR_RNDN); mpfr_exp(r20110, r20107, MPFR_RNDN); mpfr_sub(r20111, r20109, r20110, MPFR_RNDN); mpfr_mul(r20112, r20106, r20111, MPFR_RNDN); return mpfr_get_d(r20112, MPFR_RNDN); } static mpfr_t r20113, r20114, r20115, r20116, r20117, r20118, r20119, r20120, r20121, r20122, r20123, r20124, r20125, r20126, r20127, r20128, r20129, r20130, r20131; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20113, "0.5", 10, MPFR_RNDN); mpfr_init(r20114); mpfr_init(r20115); mpfr_init(r20116); mpfr_init_set_str(r20117, "1/60", 10, MPFR_RNDN); mpfr_init(r20118); mpfr_init_set_str(r20119, "5", 10, MPFR_RNDN); mpfr_init(r20120); mpfr_init(r20121); mpfr_init_set_str(r20122, "2", 10, MPFR_RNDN); mpfr_init(r20123); mpfr_init_set_str(r20124, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20125, "3", 10, MPFR_RNDN); mpfr_init(r20126); mpfr_init(r20127); mpfr_init(r20128); mpfr_init(r20129); mpfr_init(r20130); mpfr_init(r20131); } double f_fm(double re, double im) { ; mpfr_set_d(r20114, re, MPFR_RNDN); mpfr_sin(r20115, r20114, MPFR_RNDN); mpfr_mul(r20116, r20113, r20115, MPFR_RNDN); ; mpfr_set_d(r20118, im, MPFR_RNDN); ; mpfr_pow(r20120, r20118, r20119, MPFR_RNDN); mpfr_mul(r20121, r20117, r20120, MPFR_RNDN); ; mpfr_mul(r20123, r20122, r20118, MPFR_RNDN); ; ; mpfr_pow(r20126, r20118, r20125, MPFR_RNDN); mpfr_mul(r20127, r20124, r20126, MPFR_RNDN); mpfr_add(r20128, r20123, r20127, MPFR_RNDN); mpfr_add(r20129, r20121, r20128, MPFR_RNDN); mpfr_neg(r20130, r20129, MPFR_RNDN); mpfr_mul(r20131, r20116, r20130, MPFR_RNDN); return mpfr_get_d(r20131, MPFR_RNDN); } static mpfr_t r20132, r20133, r20134, r20135, r20136, r20137, r20138, r20139, r20140, r20141, r20142, r20143, r20144, r20145, r20146, r20147, r20148, r20149, r20150; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r20132, "0.5", 10, MPFR_RNDN); mpfr_init(r20133); mpfr_init(r20134); mpfr_init(r20135); mpfr_init_set_str(r20136, "1/60", 10, MPFR_RNDN); mpfr_init(r20137); mpfr_init_set_str(r20138, "5", 10, MPFR_RNDN); mpfr_init(r20139); mpfr_init(r20140); mpfr_init_set_str(r20141, "2", 10, MPFR_RNDN); mpfr_init(r20142); mpfr_init_set_str(r20143, "1/3", 10, MPFR_RNDN); mpfr_init_set_str(r20144, "3", 10, MPFR_RNDN); mpfr_init(r20145); mpfr_init(r20146); mpfr_init(r20147); mpfr_init(r20148); mpfr_init(r20149); mpfr_init(r20150); } double f_dm(double re, double im) { ; mpfr_set_d(r20133, re, MPFR_RNDN); mpfr_sin(r20134, r20133, MPFR_RNDN); mpfr_mul(r20135, r20132, r20134, MPFR_RNDN); ; mpfr_set_d(r20137, im, MPFR_RNDN); ; mpfr_pow(r20139, r20137, r20138, MPFR_RNDN); mpfr_mul(r20140, r20136, r20139, MPFR_RNDN); ; mpfr_mul(r20142, r20141, r20137, MPFR_RNDN); ; ; mpfr_pow(r20145, r20137, r20144, MPFR_RNDN); mpfr_mul(r20146, r20143, r20145, MPFR_RNDN); mpfr_add(r20147, r20142, r20146, MPFR_RNDN); mpfr_add(r20148, r20140, r20147, MPFR_RNDN); mpfr_neg(r20149, r20148, MPFR_RNDN); mpfr_mul(r20150, r20135, r20149, MPFR_RNDN); return mpfr_get_d(r20150, MPFR_RNDN); }