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